MATLAB中V-Blast的ZF和MMSE检测算法仿真代码

MATLAB中V-Blast的ZF和MMSE检测算法仿真代码

发送天线m1，接受天线m2，高斯白噪声，ZF时的仿真程序

% 发射天线数tx,接收天线数rx,发射矩阵长度L(帧长) tx=m1;rx=m2;L=10000; Modulation='BPSK';

EbN0=[0:5:20];

B=30000;Ts=1/24300;

% 建立EbN0与SNR之间的换算关系 SNR=EbN0-10*log10(Ts*B); % 信源A

A=randint(tx*L,1);

% 经过BPSK调制的V-Blast发射矩阵X X=zeros(tx,L); for k=1:tx

X(k,:)=(-1).^(A(k:tx:end)+1); end

% 信道传输============================================================ % 快衰落高斯信道H

H=sqrt(1/2)*(randn(rx,tx,L));

% 均值为0方差为1的高斯白噪声n n=sqrt(1/2)*(randn(rx,L)); % 未叠加噪声的接收信号R R=zeros(rx,L); for k=1:L

R(:,k)=sqrt(1/tx)*H(:,:,k)*X(:,k); end

% 检测

%ZF================================================================== disp('berz'); berz=[];

% 在不同的信噪比下计算ZF接收机误比特率berz for m=SNR

m

% 每个子信道的平均信噪比为snr的接受信号R_noised snr=10^(m/10); R_noised=awgn(R,m,1); x=[];

a=zeros(tx*L,1);

% 逐时隙对接收符号矢量进行检测，合并得到一帧发射矩阵X的估计x for t=1:L

r=R_noised(:,t); % 迫零矩阵G G=pinv(H(:,:,t)); y=G*r; xtemp=(y>=0)-(y<0)+0; x=[x,xtemp]; end

% 从x求A的估计a for k=1:tx

a(k:tx:end)=(x(k:tx:end)+1)/2; end

% 比较A和a计算错值率temp_ber [errbit,temp_ber]=biterr(A,a); berz=[berz,temp_ber]; end

semilogy(EbN0,berz,'o- r'),grid on xlabel('Eb/N0(dB)'); ylabel('误比特率');

title('不同天线配置ZF系统仿真传输特性')

text(9,0.15,'\\fontsize{14}\\color{red}\\fontname{隶书}2*2') hold on

发送天线m1，接受天线m2，高斯白噪声，MMSE时的仿真程序

% 发射天线数tx,接收天线数rx,发射矩阵长度L(帧长) tx=m1;rx=m2;L=10000; Modulation='BPSK';

EbN0=[0:5:20];

B=30000;Ts=1/24300;

% 建立EbN0与SNR之间的换算关系 SNR=EbN0-10*log10(Ts*B); % 信源A

A=randint(tx*L,1);

% 经过BPSK调制的V-Blast发射矩阵X X=zeros(tx,L); for k=1:tx

X(k,:)=(-1).^(A(k:tx:end)+1);

end

% 信道传输============================================================ % 快衰落高斯信道H

H=sqrt(1/2)*(randn(rx,tx,L));

% 均值为0方差为1的高斯白噪声n n=sqrt(1/2)*(randn(rx,L));

% 未叠加噪声的接收信号R R=zeros(rx,L); for k=1:L

R(:,k)=sqrt(1/tx)*H(:,:,k)*X(:,k); end

% 检测

%MMSE=============================================================== disp('berz');

berz=[];

% 在不同的信噪比下计算ZF接收机误比特率berz for m=SNR

m

% 每个子信道的平均信噪比为snr的接受信号R_noised snr=10^(m/10);

R_noised=awgn(R,m,3);

x=[];

a=zeros(tx*L,1);

% 逐时隙对接收符号矢量进行检测，合并得到一帧发射矩阵X的估计x for t=1:L

r=R_noised(:,t);

HH=H(:,:,t); xtemp=zeros(tx,1);

w=inv(HH'*HH+(1/snr)*eye(tx))*HH'; y=w*r; xtemp=(y>=0)-(y<0)+0; x=[x,xtemp]; end

% 从x求A的估计a

for k=1:tx

a(k:tx:end)=(x(k:tx:end)+1)/2; end

% 比较A和a计算错值率temp_ber [errbit,temp_ber]=biterr(A,a); berz=[berz,temp_ber]; end

semilogy(EbN0,berz,'o- b'),grid on xlabel('Eb/N0(dB)'); ylabel('误比特率');

title('不同天线配置MMSE系统仿真传输性能') text(8,0.08,'\\fontsize{14}\\color{blue}\\fontname{隶书}1*2') hold on

本文来源：https://www.bwwdw.com/article/dbg8.html