2DPSK调制与解调的数字通信系统软件仿真.docx
2DPSK调制与解调的数字通信系统软件仿真一、试验内容1、设计系统整体框图及数学模型；2、产生离散二进制信源，进行信道编码（汉明码），产生2DPSK信号；3、加入信道噪声（高斯白噪声）；4、2DPSK信号相干解调，信道解码；5、系统性能分析（信号波形、频谱，白噪声的波形、频谱，信道编解码带来的系统误码率性能增益，误码率与理论曲线对比等）。二、试验原理2DPSK的调制原理在传输信号中，2PSK信号具有较好的误码率性能。但是，在2PSK信号的载波复原过程中存在着180度的相位模糊，即复原的本地载波与所需的相干载波可能同相也可能反相，这种相位关系的不确定性将会造成解调出的数字基带信号与发送出的数字基带信号正好相反，即“1”变为“0”，“0”变为“1”，判决器输出输出数字信号全部出错。为了克服此缺点提出二进制差分相移键控（2DPSK）方式。2DPSK信号的产生方法是先对二进制数字基带信号进行差分编码，即把表示数字信息序列的肯定码变换成相对码，然后再依据相对码进行肯定调相，从而变成二进制差分相移键控信号。2DPSK信号调制器原理框图如图2-1所示。开关舱光m图2-12DPSK信号调制至原理框图2DPSK的解调原理2DPSK的解调方法分为两种：一种是相干解调加码反变换即极性比较法，另一种是差分相干解调即相位比较法。这里采纳的是第一种其解调原理是：对2DPSK信号进行相干解调，复原出相对码，再经码反变换器变换为肯定码从而复原动身送的二进制数字信息。在解调过程中，由于载波相位模糊性的影响，使得解调出的相对码也可能是“1”和“0”倒置，但经差分译码得到的肯定码不会发生任何倒置的现象，从而解决了载波相位模糊带来的问题。2DPSK的相干解调原理框图如图2-3所示定时妹冲图2-32DPSK的相干IK语原理框图三、仿真设计S(t)2DPSK调制与解调原理框图四、试验结果&分析探讨试验仿真结果:基带信号Xx10,x10递新遍后的频谐同频同相载波波形024681012141618低通滤'波输出波形协羹点捶或阱I近拽分析：通过仿真得到的误符号率和误比特率曲线与理论曲线相比差别挺大，比如在理论曲线中在信噪比为12dB时的误符号率和误比特率可达到10一7数量级而在实际仿真中这一数据在IOTio-5之间。可能的缘由如下：仿真的码元数有点少，本次仿真中仿真的码元数为10万个，从曲线中可以看出在信噪比为9dB的时候，曲线起先变差，从前面的仿真来看跟仿真的码元数少有关系；另一方面，我们产生的是伪随机基带信号，可能带有肯定的周期性，在肯定程度上影响了误符号率和误比特率与信噪比的关系。程序源代码产生基带信号及码变换：functionCoded_Signal=CodeSignal(n/Num_Unit/Signal,bitRate)M=2;%二施制N三7j%(7,4)Hamming码m=3;DeCodeSignaI=;DeCodeSignaI(I)=0;CodeSignal=encode(SignalfNzN-m);%Hamming编码CodeSignal=reshape(CodeSignal.'zlog2(M),n*Nlog2(M).,;Signal=reshape(Signal.',log2(M),n*(N-m)log2(M).,;bit2=;fori=l:length(Signal)ifSignal(i)=Obit3=zeros(l,Num-Unit);elsebit3=ones(l,Num_Unit);endbit2=bit2,bit3;endfigure(l)SUbPlOt(3,1,1)plot(bit2);title('基带信号'),gridon;a×is(0,Num-Unit*(length(Signal)-4),-1.5,l5);1.enth_Of_CoSign=length(CodeSignaI);%>>>>>>>>>>>>>肯定码->相对码>>>>>>>>>>>>>>>DeCodeSignal(l)=0+CodeSignal(l);ifDeCodeSignaI(I)=ZDeCodeSignaI(I)=O;elseDecodeSignaI(I)=DeCodeSignaI(I);endforww=2:1.enth_Of_CoSignDeCodeSignal(ww)=DeCodeSignal(ww-l)+CodeSignal(ww);ifDeCodeSignal(ww)=2DeCodeSignaI(Ww)=O;elseDeCodeSignal(ww)=DeCodeSignal(ww);endendCoded_Signal=DeCodeSignaI;%>>>>>>>>>>>产生基带信号>>>>>>>>>>>>Signalrign=;fori=l:1.enth_Of_CoSignifCodeSignal(i)=0SignalOrignl=zeros(l,Num_Unit);elseSignaIOrignI=OneS(1,NUm_Unit);endSignalOrign=SignalOrign,SignalOrignl;endSUbPIOt(3,1,2)PloUSignaIOrign);tit®，编码后基带信号Igridon;axis(O,Num_Unit*1.enth_Of_CoSign厂1.5,1.5);SUbPIOt(3,1,3)%plot(abs(fft(SignalOrign);P,F=t2f(Signal0rign,Num_Unit*bitRate);Pot(F,P)tit®，基带信号频谱，)；end2DPSK调制:functionModuled_coSignal=Modu2DPSK(coSignal,bitRate,Freq_Carrier,N)t=linspace(O,lbitRate,N);cl=sin(2*pi*t*Freq-Carrier);c2=sin(2*pi*t*Freq_Carrier+pi);Moduled_coSignal=;forii=l:length(coSignal)ifcoSignal(ii)=lModuled_coSignal=Moduled-coSignalfcl;elseModuIed_coSignaI=Moduled-coSignabc2;endendfigure(9)P2,F2=t2f(Moduled_coSignal,bitRate*N);plot(F2zP2)tit同，调制后频谱)%figure(2)%画调制图%plot(l:length(transmittedSignal),transmittedSignal);title('ModulationofBPSK');gridon;%figure(3)%画频谱图%m=0:length(Moduled_Signal)-l;%F=fft(Moduled_Signal);%plot(m,abs(real(F),title('BPSK-frequency-domainanalysisreal');%gridon;EndfunctionP,f=t2f(s,fs)%s代表输入信号，S代表S的频谱，fs是采样频率N=length(三);%样点总数T=lfs*N;%视察时间f=-N/2:(N/2-l)/T;%频率采样点tmpl=fft(三)fs;tmp2=N*ifft(三)fs;S(l:N/2)=tmp2(N2+l:-1:2);S(N/2+l:N)=tmpl(l:N/2);P=abs(三)2T;EndfunctionSignal=SignalSet(nzNum-UnitzbitRate)Signal=randint(n,4,2);end编码：functioncoSignalAndNoise=lnforChannel(Moduled-coSignal,ebno,Num-UnitzbitRate,Freq-Carrier,n,conumbits,numbits)Eb=(l/conumbits)*sum(Moduled_coSignal.A2);EbNo=10.(ebno.10);No=Eb./EbNo;nstdv=sqrt(No.2);k=1;Noise(k,:)=nstdv.*randn(l,length(Moduled-coSignal);coSignalAndNoise(kz:)=NOiSe(k,:)+MOdUIed_coSignal;figure(2)SUbPIOt(3,1,1)plot(Moduled_coSignal);title(巳调信号,),gridon;%axis(0,Num-Unit*1.enth-Of_CoSign,-2.5,2.5);SUbPlOt(3,1,2)PlOt(INength(Noise),Noise),gridon;tit®，噪声叫%a×is(ONum_Unit*1.enth_Of_CoSignO400);SUbPlOt(3,1,3)Plot(Idength(CosignalAndNoise)zCoSignaIAndNoise),gridon;title。加噪声的已调信号%a×is(ONum_Unit*1.enth_Of_CoSignO400);figure(ll)subplot(311)p3,f3=t2f(coSignalAndNoise,Num_Unit*bitRate);plot(f3,p3)t=linspace(0,(7*n)bitRate,7*n*Nun-Unit)%7bits700份wp=2*pi*Freq_Carrier*0.92*pi*Freq_Carrier*l.l;ws=2*pi*Freq_Carrier*0.82*pi*Freq_Carrier*1.2;Rp=l;Rs=2;N,wc=buttord(wp,ws,Rp,Rss,);B,A=butter(N,wcs,);H,W=freqz(B,A);H=tf(B,A);coSignalAndNoise,tl=sim(H,CoSignaIAndNoise,t)；p4J4=t2f(coSignalAndNoise,Num_Unit*bitRate);subplot(313)plot(f4,p4);title('通过带通后的频谱)end相干解调：functionDemoduled-coSignal=demo2DPSK(coSignalAndNoise,bitRate,Freq_Carrier,n,Num-Unit)t=linspace(0,7bitRate,7*NumUnit)j%7bi
