1、基于MATLAB的DPSK通信系统仿真设计通信原理课程实验报告基于MATLAB的DPSK通信系统仿真院系: 班级: 学号: 姓名: 教师:2012年12月7日1设计题目 3 2设计要求 3 3设计原理 3 3.1DPSK 调希 9 3 3.2DPSK 解调 4 4设计流程 5 5代码实现 5 5.1基带信号的获取 5 5.2差分编码 6 5.3DPSK 调希 9 7 5.4AWGN 信道 8 5.5限噪处理 8 5.6差分相干解调 9 5.4 AWGN 信道 11 5.7功率谱密度图的绘制6测试结果 6.1调制解调正确性的测试 126.2系统抗噪性能的测试 14 -7总结 17 8心得体会 1
2、7 附：完整代码 18 -1设计题目用MATLAB仿真DPSK通信系统2设计要求基本参数如下：1)Fc=1800Hz2)Rb=1200bps3)考虑加入AWGN信道,Eb/N0=20dB4）自定义一小段信息序列结果要求：1） 绘制各阶段信号波形2） 绘制信号的PSD3设计原理二进制差分相移键控常简称为二相相对调相，记为2DPSK。它不是利用载波 相位的绝对数值传送数字信息，而是用前后码元的相对载波相位值传送数字信 息。所谓相对载波相位是只本码元初相与前一码元初相之差 。3.1DPSK 调制9DPSK方式即是利用前后相邻码元的相对载波相位值去表示数字信息的一种 方式。例如，假设相位值用相位偏移表
3、示（定义为本码兀初相与前一码 元初相只差），并设二n数字信息1=0 数字信息0则数字信息序列与2DPSK信号的码元相位关系可举例表示如如下：数字信息： 00111001012DPSK 信号相位：000 n 0 nnn 0 0 n或 nnn 0 n 0 0 0 n n 0画出的2PSK及DPSK信号的波形如图3-1所示。图3-1 2PSK及2DPSK信号的波形2DPSK的产生基本类似于2PSK,只是调制信号需要经过码型变换，将绝对码变为相对码。2DPSK产生的原理框图如图2-2所示1，图（a）为模拟调制图3-2 2DPSK信号的调制原理图从上面分析可见，无论接收信号是2DPSK还是2PSK信号，
4、单从接收端看 是区分不开的。因此2DPSK信号的功率谱密度和2PSK信号的功率谱密度是完 全一样的。3.2DPSK 解调DPSK信号可以采用相干解调法（极性比较法）和差分相干解调法（相位比 较法）。本次设计采用差分相干解调法，图3-3为差分相干解调法原理图。差分相干解调法的思想是直接比较前后码元的相位差 。由于此时的解调已同时完成了码变换作业，故无需另加码变换器。由于这种解调方法又无需专门4设计流程设计流程如下框图所示，箭头上的信号表示实际传递的信号名称bn dn mod modz( % L 、 差分编码 DPSK调制 AWGN信道 BPF对比Bnmtdemod 1抽样判决LPF延时相乘差分相
5、干解调5代码实现5.1基带信号的获取利用MATLAB自带函数ran di nt产生定长度的二进制随机序列sign ,为了将 sign转换为单极性NRZ信号，要对sign的每个随机码进行T点增样，最后得到 基带数字信号bn代码实现如下：%初始参量定义 %fc=1800;%载波频率Rb=1200;%基带数据传输速率Tb=1/Rb;T=50;%每个Tb内进行的增样点数M=8;%原始码字序列的长度N=9;%差分码字序列的长度sig n=randin t(1,M);%产生长度为M的随机二进制序列SNR=20;%信噪比%对sign进行T点增样后的基带数字信号 %bn=;for n=1:Mif sig n(
6、n)=0;F=zeros(1,T);elseF=on es(1,T);endbn=b n,F;end基带数字信号bn100 200 300 400210-15.2差分编码设差分码cfm （二进制序列），差分编码原理为cfmn = bn xor cfmn-1例如,女口果 sign= 1 0 0 1 0 0 1 1 , 那么 cfm=1 0 0 0 1 1 1 0 1 ;同时，还要对得到的差分码cfm进行T点增样，得到差分信号dn （单极性NRZ信号）。代码实现如下：% 差分编码得到差分码cfm %cfm=;cfm(1)=1; %设cfm的初值为1for i=1:Mcfm(i+1)=xor(sig
7、 n(i),cfm(i);end% 对cfm 进行 T点增样 %dn=; %差分信号for n=1:N;if cfm(n)=0;A=zeros(1,T); %每个symbol进行T点增样elseA=o nes(1,T);enddn=dn A;end差分信号血0 100 200 300 4005.3DPSK 调希9DPSK调制由于PSK调制的主要原因是DPSK调制解决了 不确定性反相”问 题，即DPSK信号通过相邻时隙载波相位的变化与否来 携带”信息。delta，易就能反映前在代码实现这一步的过程中，我假设了一个相邻时隙相位差变量知，当sign的一个码字为时，delta=pi，否则delta=O
8、.通过delta 后时隙的载波相位变化，进而实现DPSK调制。代码实现如下：% DPSK调制 %t=Tb/T:Tb/T:Tb; %对Tb=1/Rb间隔进行T点划分mod=cos(2*fc*pi*t); % 调制后的信号delta=0; %相位差for n=1:Mif sig n(n)=1;delta=delta+pi;elsedelta=delta;endc2=cos(2*fc*pi*t+delta);mod=mod c2;end绘图5.4 AWGN 信道利用MATLAB自带函数awgn实现对已调信号mod进行加噪处理代码实现如下：% 调制信号过高斯白噪声的信道 %modz=awgn(mod,
9、SNR); %过AWGN 信道5.5限噪处理利用MATLAB自带函数fir1构造FIR I型带通滤波器，使其对准信号频带，让信号几乎无失真通过，同时抑制带外噪声，此处取BPF的通带宽度为2*Rb， 中心频率为fc。(由于设计的是数字滤波器，所以不要忘记对频率进行归一化 处理)值得注意的是，由于设计的滤波器是100阶，所以通过BPF后的信号会有出 现接近于一个symbol宽度的时延。所以在信号modz接入BPF前，要对其进行 一个拓宽处理，对modz最后一个symbol进行复制，得到modz_broad ，再接 入至U BPF上。同时，不要忘记对从BPF出来的信号r0进行截取，因为r0的第一个s
10、ymbol 是滤波器时延带来的。代码实现如下：%BPF %modz_broad=modz modz(N-1)*T+1:N*T);fs=Rb*T;hBPF=fir1(100,2*pi*fc-Rb fc+Rb/fs);r0=filter(hBPF,1,modz_broad);r=r0(T+1:(N+1)*T); %对 r0 进行截取过日PF后的倍号210-120 100 200 300 4005.6差分相干解调差分相干解调的核心思想是比较两个相邻时隙上信号的相位 ，从而直接还原出信号Bn,分为三步：(1)延时相乘代码实现如下：% 延时相乘 %r_shift=r(T+1:N*T);demod=;fo
11、r i=1:M*T;a=r_shift(i) .* r(i);demod=demod a;end绘图LPF利用MATLAB自带函数firl构造FIR I型 低通滤波器,此处取LPF的通带宽度为 Rb。同BPF一样，由于设计的滤波器是100阶，所以通过LPF后的信号会有出现接 近于一个symbol宽度的时延。所以要对信号demod接入LPF前，进行拓宽处 理，得到demod_broad ，再接入到LPF上。同样，不要忘记对从 LPF出来的信号 mtO进行截取，因为mtO的第一个 symbol是滤波器时延带来的。代码实现如下：% LPF % demod_broad=demod demod(M-1)
12、*T+1:M*T) hLPF=fir1(100,2*pi*Rb/fs);mt0=filter(hLPF,1,demod_broad);%对mt0进行截取mt=mt0(T+1:(M+1)*T);抽样判决，如果sum大于0,则Bn实现思想：对mt 个symbol的T个采样点进行求和对应这个symbol内的所有点的值为1，反之为0.代码实现如下：% 抽样判决 %Bn=zeros(1,M*T);for i=1:M;if sum( mt( (i-1)*T+1 :i*T ) ) 0;Bn (i-1)*T+1:i*T)=0;elseBn (i-1)*T+1:i*T)=1;end end绘图经抽样判决后信号B
13、n210-1*2100 200300 400从上图可知，最终经DPSK通信系统调制解调后得到的信号 Bn与一开始的基带 数字信号bn完全一致。5.7功率谱密度图的绘制采用 Welch法绘制DPSK调制信号mod的功率谱密度图 (Power spectrummagnitude)。选择适当的窗函数 w(n)，并再周期图计算前直接加进去，力卩窗的 优点是无论什么样的窗函数均可使谱估计非负 。二是在分段时，可使各段之间 有重叠，这样会使方差减小。本设计中尝试了三种窗，分别是矩形窗、海明窗 和 Blacman 窗。代码实现如下：% PSD %win dow=boxcar(100);%矩形窗win dow
14、1=ham min g(100);%海明窗wi ndow2=blackma n(100);%Blackman 窗Nfft=N*T*Rb;noverlap=0; %数据无重叠range= onesided ; %计算单边带功率谱密度Pxx,f=pwelch(mod,wi ndow, no verlap,Nfft,fs,ra nge); Pxx1,f=pwelch(mod,wi ndow1, no verlap,Nfft,fs,ra nge); Pxx2,f=pwelch(mod,wi ndow2, no verlap,Nfft,fs,ra nge);figure(3)subplot(311)pl
15、ot(f,Pxx);axis(0 5000 0 0.0008);title( rectangular window PSD );xlabel( f/Hz);subplot(312)plot(f,Pxx1);axis(0 5000 0 0.0006); title( Hamming PSD );xlabel( f/Hz);subplot(313)plot(f,Pxx2);axis(0 5000 0 0.0006);title( Blackman PSD );xlabel( f/Hz);-4X 10recta ngular wi ndow PSD0 500 1000 1500 2000 2500
16、3000 3500 4000 4500 5000f/Hz-4x 10Hammi ng PSD0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000f/Hz-4x 10Blackma n PSD0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000f/Hz从第一幅图即rectangular PSD中可以看出,单边带功率谱密度图的中心频率大致在1800Hz即载波频率处，第一零点大致在 600Hz和3000Hz，即fc-Rb和 fc+Rb 处。6测试结果6.1调制解调正确性的测试在前面我取用了一个特例sign
17、=1 0 0 1 0 0 1 1 来说明每一步的正确性。下 面利用ran di nt产生定长度的二进制随机序列sign，多次运行结果都显示最终差 分相干解调得出的信号Bn与基带数字信号bn 一致，说明此次基于MATLAB设计 的DPSK信系统能够成功实现DPSK信号的调制与解调。某几次运行结果的截图如下：1 基带数字信号bn差分信号dnDPSK调制信号mod20-20 200 4000 200 400基带数字信号bn2 T0 -2 L0 200 400差分信号dn-2 0 200 400DPSK调制信号mod2 m0 200 4002 0 -2 匚0 200 400-2 0 200 400差分
18、信号dn DPSK调制信号mod基带数字信号bn2 p0 110-10 200 4002 0 1 -2 L0 200 4006.2系统抗噪性能的测试在前面调制解调正确性的测试中，设定信噪比SNR=20.下面通过检验在不同 信噪比值下，系统调制解调的正确性，来简单评定系统的抗噪性能。SNR=10时，正确解调基带数字信号bn差分信号dn|hj U 11 JDPSK调制信号mod20-20 200 400SNR=2时，正确解调SNR=1.5时，有误码出现基带数字信号bn0差分信号dn0DPSK调制信号mod2 V JI .ijl0-2-2-202004000 200 400过LPF后的信号mt0 2
19、00 400经抽样判决后信号Bn2-200200400所以，可以定性判断该系统的抗噪性能良好7总结经过两周的MATLAB课程设计，我收获多多。总结后主要有一下几点可谈1） 良好编程习惯的培养我们编程时要有良好的风格，源代码的逻辑简明清晰，易读易懂是好程序 的重要标准。我的几点小收获是：标识符命名尽量含义鲜明，适当的注解，注 意程序的视觉组织。2） 信息筛查能力的训练拿到一个课程设计题目后，应该从何下手，如何获取有效资源、筛查信 息，是对个人能力的一种考验。经过这次试验后，我认为首先需要把理论知识 弄透彻（尤其是课本上的），再阅读别人的设计报告，重复别人所做的工作是 一种不错学习方法，但是需要懂
20、得去粗取精，提取有效信息，否则会浪费大量 时间。3） 调试思想的建立这次试验中我被很多问题难住，却常常不知如何下手找出错误，经老师指 点后才懂得要按一定方法进行调试纠错，结果这些错误常常是一些拼写错误或 符号不对称等小错。明确整体的架构，有针对地进行测试。调试方法多种多 样，但调试思想的建立却是最为基础、根本。8心得体会课程实验是培养学生综合运用所学知识，发现，提出，分析和解决实际问题，锻炼 实践能力的重要环节，是对学生实际工作能力的具体训练和考察过程 。在这两周的时间里我切切实实的学到了许多知识，尤其是提高了在MATLAB中的编程能力。回顾起此次通信原理课程实验，至今我仍感慨颇多，在整整两周
21、 的时间里，我的心情和课程设计的进展起起伏伏 ，一度怀疑自己的能力开始否 定自己，直到最后把每处都弄通了、达到完美了，心情才真正舒畅开来。但是，我确实学到很多很多的的东西，同时巩固了所学知识。整个设计的过 程，可以说得是困难重重，这毕竟算是自己第一次独立完成，在设计的过程中 发现了自己的不足之处，对软件不熟悉，动手能力比较欠缺。设计的过程很纠 结，调试的过程也很痛苦，但是知识就是在这个过程中学习的，经验也是在这 个过程中积累的。经历过困难才会有更多的收获。掌握正确的方法，保持认真 的作风和平和的心态，这样才能真正的学到知识。这次课程实验终于完成了 ，在设计中遇到了很多问题，多谢李晓峰老师的耐
22、心讲解，从老师的身上我学到了很多实用的知识 ，也由衷地庆幸自己能遇到这 样一位好老师，不仅专业知识扎实，而且品格高尚、谦和有礼，这就是我理想 中的好老师。再次感谢李老师。同时，对给予我帮助的同学们也表示忠心的感 谢。附：完整代码% DPSK 调制解调 %初始参量定义 %fc=1800;%载波频率Rb=1200;%基带数据传输速率Tb=1/Rb;T=50;%每个Tb内进行的增样点数M=8;%原始码字序列的长度N=9;%差分码字序列的长度sig n=zeros(1,M);sig n=randin t(1,N);%产生长度为M的随机二进制序列SNR=1.5;%信噪比%对sign进行T点增样后的基带数
23、字信号 %bn=;for n=1:Mif sig n(n)=0;F=zeros(1,T);elseF=on es(1,T);endbn=b n,F;end% 差分编码得到差分码cfm %cfm=;cfm(1)=1; %设cfm的初值为1for i=1:Mcfm(i+1)=xor(sig n(i),cfm(i);end%对cfm进行T点增样 %dn=;%差分信号for n=1:N;if cfm(n)=0;A=zeros(1,T);%每个symbol进行T点增样elseA=o nes(1,T);enddn=dn A;end%DPSK调制 %t=Tb/T:Tb/T:Tb;%对Tb=1/Rb间隔进行T
24、点划分mod=cos(2*fc*pi*t);%调制后的信号delta=O;%相位差for n=1:Mif sig n(n)=1; delta=delta+pi;elsedelta=delta;endc2=cos(2*fc*pi*t+delta);mod=mod c2;end% 调制信号过高斯白噪声的信道 %modz=awg n(m od,SNR);%过AWGN信道% BPF %modz_broad=modz modz(N-1)*T+1:N*T);fs=Rb*T;hBPF=fir1(100,2*pi*fc-Rb fc+Rb/fs);rO=filter(hBPF,1,modz_broad);r=r
25、0(T+1:(N+1)*T); %对 rO 进行截取% 延时相乘 %r_shift=r(T+1:N*T);demod=;for i=1:M*T;a=r_shift(i) .* r(i);demod=demod a;end% LPF %demod_broad=demod demod(M-1)*T+1:M*T)hLPF=fir1(100,2*pi*Rb/fs);mtO=filter(hLPF,1,demod_broad);mt=mt0(T+1:(M+1)*T); %对mtO 进行截取% 抽样判决 %Bn=zeros(1,M*T);for i=1:M;if sum( mt( (i-1)*T+1 :i
26、*T ) ) 0;Bn (i-1)*T+1:i*T)=0;elseBn (i-1)*T+1:i*T)=1;endend% 画图 % figure(1) % figure(1);subplot(3,3,1);plot(bn);grid on ;axis(0 T*N -2 2);title(基带数字信号 bn);figure(1);subplot(3,3,2);plot(d n);grid on;axis(0 T*N -2 2);title(差分信号 dn);figure(1);subplot(3,3,3);plot(mod);grid on;axis(0 T*N -2 2);title( DPSK调制信号 mod);figure(1);subplot(3,3,4);plot(modz);grid onaxis(0 T*N -2 2);title(过AWGN信道后的modz信号);figure(1);subplot(3,3,5);plot(r);grid onaxis(0 T*N -2 2);title(过BPF后的信号r)figure(1);subplot(3,3,6);plot(demod);grid onaxis(0 T*N -1 1);title(
