基于AF功率优化和QPSK调制解调的通信系统误码率仿真,含仿真操作录像

clc; clear; close all; warning off; addpath(genpath(pwd)); fprintf( 'AF仿真\n') ; %产生随机序列 N=120000; maxsnr=40; BerSnrTable1 = zeros(maxsnr+1,3); for SNR=0:maxsnr SNR BerSnrTable1(SNR+1,1) = SNR; BitTx=floor(rand(1,N)*2); %星座图映射，QPSK SymQpsk=QpskMapping(BitTx); SymQpskd=[]; M=length(SymQpsk); %噪声能量（已归一化） zsnl=sqrt(1/(10^((SNR)/10))); %========================================== for i=1:M %生成服从瑞利分布的信道系数 H0=(randn+j*randn)/sqrt(2); H1=(randn+j*randn)/sqrt(2); H2=(randn+j*randn)/sqrt(2); %生成加性高斯白噪声 N0=zsnl*(randn+j*randn); N2=zsnl*(randn+j*randn); N1=zsnl*(randn+j*randn); a=abs(H1); G=sqrt(2/(4*a*a+3*zsnl*zsnl));%E0=4/3，E1=2/3,根据<无线合作分集网络中的资源分配研究> %接收信号 R1=H0*sqrt(4/3)*SymQpsk(i)+N0; R2=H2*G*(H1*sqrt(4/3)*SymQpsk(i)+N1)+N2; %信号合成 S1=conj(H0)*sqrt(4/3)*R1/(zsnl*zsnl); S2=conj(H2)*conj(H1)*sqrt(4/3)*G*R2/(zsnl*zsnl); SS=S1+S2; %解调 dh = [1+j -1+j -1-j 1-j]/sqrt(2); D1=abs(SS*[1 1 1 1]-dh).^2; [minScale1 positionmin1]=min(D1); SymQpskd=[SymQpskd dh(positionmin1)]; end BitRx=QpskInverseMapping(SymQpskd); [Num,Ber]=symerr(BitTx,BitRx); BerSnrTable1(SNR+1,2) = Num ; BerSnrTable1(SNR+1,3) = Ber ; end figure(1); semilogy(BerSnrTable1(1:4:maxsnr+1,1),BerSnrTable1(1:4:maxsnr+1,3),'b-o');