基带传输实验图形.doc

正交信号在AWGN信道下的传输性能代码：clear allnsamp=10; s0=ones(1,nsamp); s1=ones(1,nsamp/2) -ones(1,nsamp/2); nsymbol=100000; %每种信噪比下的发送符号数EbN0=0:12; %信噪比，E/N0msg=randint(1,nsymbol); %消息数据s00=zeros(nsymbol,1); s11=zeros(nsymbol,1);indx=find(msg=0); %比特0在发送消息中的位置s00(indx)=1;s00=s00*s0; %比特0影射为发送波形s0indx1=find(msg=1); %比特1在发送消息中的位置s11(indx1)=1; s11=s11*s1; %比特1映射为发送波形s1s=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);endsemilogy(EbN0,ber,-ko,EbN0,qfunc(sqrt(10.(EbN0/10);title(二进制正交信号在AWGN信道下的误比特率性能)xlabel(Eb/N0);ylabel(误比特率Pe)legend(仿真结果,理论结果)图形改善代码可以不考虑发送信号波形的影响：clear allnsymbol=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);endfiguresemilogy(EbN0,ber,-ko,EbN0,qfunc(sqrt(10.(EbN0/10);title(二进制正交信号在AWGN信道下的误比特率性能)xlabel(Eb/N0);ylabel(误比特率Pe)legend(仿真结果,理论结果)图形：这种修改可以使仿真的运行时间变短，效率高用simulink仿真上边的代码：仿真时间设为100000双极性信号在AWGN信道下的传输性能程序代码如下clear allnsamp=10; %每个脉冲信号的抽样点数s0=ones(1,nsamp); %基带脉冲信号s1=-s0; nsymbol=100000; %每种信噪比下的发送符号数EbN0=0:10; %信噪比，E/N0msg=randint(1,nsymbol); %消息数据s00=zeros(nsymbol,1); s11=zeros(nsymbol,1);indx=find(msg=0); %比特0在发送消息中的位置s00(indx)=1;s00=s00*s0; %比特0影射为发送波形s0indx1=find(msg=1); %比特1在发送消息中的位置s11(indx1)=1; s11=s11*s1; %比特1映射为发送波形s1s=s00+s11; %总的发送波形s=s.; %数据转置，方便接收端处理for indx=1:length(EbN0) decmsg=zeros(1,nsymbol); r=awgn(s,EbN0(indx)-7); %通过AWGN信道 r00=s0*r; %与s0相关 indx1=find(r005); decmsg(indx1)=1; %判决 err,ber(indx)=biterr(msg,decmsg);endsemilogy(EbN0,ber,-ko,EbN0,qfunc(sqrt(10.(EbN0/10)/2),EbN0,qfunc(sqrt(10.(EbN0/10),-k*,EbN0,qfunc(sqrt(2*10.(EbN0/10),-kv);title(单极性信号在AWGN信道下的误比特率性能)xlabel(Eb/N0);ylabel(误比特率Pe)legend(单极性信号仿真结果,单极性信号理论误比特率,正交信号理论误比特率,双极性信号理论误比特率)程序运行结果3基带PAM信号传输3.1基带4PAM的信号波形3.2基带4PAM信号在AWGN信道下的最佳接收3.3基带4PAM信号在AWGN信道下的传输性能程序代码：clear allnsymbol=100000; %每种信噪比下的发送符号数nsamp=10; %每个脉冲信号的抽样点数M=4; %4-PAMgraycode=0 1 3 2; %Gray编码规则 EsN0=0:15; %信噪比，E/N0msg=randint(1,nsymbol,4); %消息数据msg1=graycode(msg+1); %Gray映射msg2=pammod(msg1,M); %4-PAM调制s=rectpulse(msg2,nsamp); %矩形脉冲成形for indx=1:length(EsN0) decmsg=zeros(1,nsymbol); r=awgn(real(s),EsN0(indx)-7,measured); %通过AWGN信道 r1=intdump(r,nsamp); %相关器输出 msg_demod=pamdemod(r1,M); %判决 decmsg=graycode(msg_demod+1); %Gray逆映射 err,ber(indx)=biterr(msg,decmsg,log2(M); %求误比特率 err,ser(indx)=symerr(msg,decmsg);endsemilogy(EsN0,ber,-ko,EsN0,ser,-k*,EsN0,1.5*qfunc(sqrt(0.4*10.(EsN0/10);title(4-PAM信号在AWGN信道下的性能)xlabel(Es/N0);ylabel(误比特率和误符号率)legend(误比特率,误符号率,理论误符号率)程序运行结果： 4带限信道的信号传输4.1带限信道4.2带限信道信号无ISI的条件4.3带线信道信号传输的仿真程序代码： 原始数据与脉冲成型后的数据 滤波器的冲击响应程序代码clear allnsymbol=100000; %每种信噪比下的发送符号数Fd=1; %符号采样频率Fs=10; %滤波器采样频率rolloff=0.25; %滤波器滚降系数delay=5; %滤波器时延M=4; %4-PAMgraycode=0 1 3 2; %Gray编码规则 EsN0=0:15; %信噪比，E/N0msg=randint(1,nsymbol,4); %消息数据msg1=graycode(msg+1); %Gray映射msgmod=pammod(msg1,M); %4-PAM调制rrcfilter = rcosine(Fd,Fs,fir/sqrt,rolloff,delay); %设计根升余弦滤波器s=rcosflt(msgmod,Fd,Fs,filter,rrcfilter);for indx=1:length(EsN0) decmsg=zeros(1,nsymbol); r=awgn(real(s),EsN0(indx)-7,measured); rx=rcosflt(r,Fd,Fs,Fs/filter,rrcfilter); rx1=downsample(rx,Fs); rx2=rx1(2*delay+1:end-2*delay); msg_demod=pamdemod(rx2,M); %判决 decmsg=graycode(msg_demod+1); %Gray逆映射 err,ber(indx)=biterr(msg,decmsg,log2(M); %求误比特率 err,