1.算法描述 16QAM全称正交幅度调制是英文Quadrature Amplitude Modulation的缩略语简称,意思是正交幅度调制,是一种数字调制方式。产生的方法有正交调幅法和复合相移法。
波束形成是阵列信号处理中非常重要的任务,基本处理过程为:
1.采用空间分布的传感器阵列采集场数据
2.对所采集的阵列数据进行线性加权组合处理得到一个标量波束输出。
自适应阵列天线具有抗干扰的作用,他的基本原理是利用算法,在保证信号大增益接收的前提下,自适应的使天线的方向图零陷对准干扰方向,从而抑制干扰或者降低干扰信号的强度,使有用接收信号的增益最大,干扰的增益最小。自适应波束形成通过不同的准则确定自适应加权,主要准则有:1.MSE(最小均方误差)2.SNR(最大信噪比)3.LH(最大似然比)
这些准则都是采用一定的算法,调整波束方向图,从而实现自适应控制。
QAM调制中,数据信号由相互正交的两个载波的幅度变化表示。QAM是一种矢量调制,将输入比特先映射(一般采用格雷码)到一个复平面(星座)上,形成复数调制符号,然后将该符号的I、Q分量(即该复数的实部和虚部)采用幅度调制,分别对应调制在相互正交(时域正交)的两个载波coswt和sinwt上。
正交振幅调制可以表示为:
式中:两个相互正交的载波分量中,每个载波被一组离散的振幅{Am}、{Bm}所调制,故称这种方式为正交振幅调制;Ts是码元宽度,m=1、2、...、M,M为Am和Bm的电平数。振幅Am和Bm可以表示成:
式中:A是固定振幅,dm、em由输入信号确定。dm、em决定已调QAM信号在星座图中的坐标点。QAM是幅度、相位联合调制的技术,它同时利用载波的幅度和相位来传递信息比特,因此在最小距离相同的条件下可实现更高的频带利用率。 根据802.11a规定,16QAM编码表如表4-14所示。经过表4-14映射后得到的I/Q数据再乘上1/sqrt(10)进行归一化,即得到调制后的I/Q值。I/Q数据分别进行A/D变换,得到两路模拟电平信号,用于和coswt、-sinwt相乘,从而实现调制。
2.仿真效果预览 matlab2022a仿真如下:
3.MATLAB核心程序
N0 = 1/(10^(SNR/10));
delta2 = N0;
deltab = 0.5;
error_count = 0;
bit_count = 0;
index = 0;
ERR_NUM = [];
tmps = 0;
while error_count < 1000
index = index+1;
for kk=1:nTx
bits(kk,:) = round(rand(1,frame_length));
symbols(kk,:) = qammod(bits(kk,:),16);
end
%transmit signal
s = symbols;
u = reshape(s,nTx,nRx,length(s)/nRx);
%Channel
h = 1/sqrt(2)*[randn(nRx,nTx,length(s)/nRx) + j*randn(nRx,nTx,length(s)/nRx)];
for ij = 1:nTx
p(:,ij,:) = h(:,ij,:).*exp(j*alpha(ij)*pi/180);
end
%mmse beamforming
LMS = zeros(1,length(s)/nRx);
for i=1:length(s)/nRx
for ii = 1:length(SIR)
u2(:,ii) = u(:,ii,i)*10^(SIR(ii)/20);
end
XN(:,:,i) = awgn(u2,SNR,'measured');
%定义接收信号
w(:,:,i) = inv(p(:,:,i)*p(:,:,i)'+2*delta2/deltab*eye(nRx))*p(:,1,i);
WK = w(:,:,i)'*p(:,:,i);
WK2 = WK/max(abs(WK));
yhat(:,:,i) = WK2*XN(:,:,i);
R1 = real(yhat(:,:,i));
I1 = imag(yhat(:,:,i));
s_hat(:,:,i)= qamdemod(R1+sqrt(-1)*I1,16);
end
s_hat1 = squeeze(s_hat);
recovered_bits = reshape(s_hat1,1,length(s));
ERR_NUM = sum(recovered_bits ~= bits(1,:));
if index <= 50
tmps = tmps + ERR_NUM;
else
if ERR_NUM/(tmps/50) < 10
SNR
error_count
error_count = error_count + ERR_NUM;
bit_count = bit_count + frame_length;
end
end
end
%Calculate the BER
BER = error_count/bit_count;
Bers = [Bers,BER];
end
01_104m
标签:nTx,正交,length,nRx,matlab,16QAM,波束,bits,调制
From: https://blog.51cto.com/matworld/6033873