基带传输实验图形.doc

正交信号在AWGN信道下的传输性能 代码： clear all nsamp=10; s0=ones(1,nsamp); s1=[ones(1,nsamp/2) -ones(1,nsamp/2)]; nsymbol=100000; %每种信噪比下的发送符号数 EbN0=0:12; %信噪比，E/N0 msg=randint(1,nsymbol); %消息数据 s00=zeros(nsymbol,1); s11=zeros(nsymbol,1); indx=find(msg==0); %比特0在发送消息中的位置 s00(indx)=1; s00=s00*s0; %比特0影射为发送波形s0 indx1=find(msg==1); %比特1在发送消息中的位置 s11(indx1)=1; s11=s11*s1; %比特1映射为发送波形s1 s=s00+s11; %总的发送波形 s=s.; %数据转置，方便接收端处理 for indx=1:length(EbN0) decmsg=zeros(1,nsymbol); r=awgn(s,EbN0(indx)-7); %通过AWGN信道 r00=s0*r; %与s0相关 r11=s1*r; %与s1相关 indx1=find(r11=r00); decmsg(indx1)=1; %判决 [err,ber(indx)]=biterr(msg,decmsg); end semilogy(EbN0,ber,-ko,EbN0,qfunc(sqrt(10.^(EbN0/10)))); title(二进制正交信号在AWGN信道下的误比特率性能) xlabel(Eb/N0);ylabel(误比特率Pe) legend(仿真结果,理论结果) 图形 改善代码可以不考虑发送信号波形的影响： clear all nsymbol=100000; %发送符号数 EbN0=0:12; %信噪比 msg=randint(1,nsymbol); %消息数据 E=1; r0=zeros(1,nsymbol); r1=zeros(1,nsymbol); indx=find(msg==0); r0(indx)=E; indx1=find(msg==1); r1(indx1)=E; for indx=1:length(EbN0) dec=zeros(1,length(msg)); snr=10.^(EbN0(indx)/10); %dB转换为线性值 sigma=1/(2*snr); %噪声方差 r00=r0+sqrt(sigma)*randn(1,length(msg)); %相关器的输出 r11=r1+sqrt(sigma)*randn(1,length(msg)); indx1=find(r11=r00); %判决 dec(indx1)=1; [err,ber(indx)]=biterr(msg,dec); end figure semilogy(EbN0,ber,-ko,EbN0,qfunc(sqrt(10.^(EbN0/10)))); title(二进制正交信号在AWGN信道下的误比特率性能) xlabel(Eb/N0);ylabel(误比特率Pe) legend