通信原理课程设计代码

1、PSK调制代码clearall;closeall;clf;%清除窗口中的图形max=20%定义max长度g=zeros(1,max);g=randint(1,max);%长度为max的随机二进制序列cp=;mod1=;f=2*2*pi;t=0:2*pi/199:2*pi;forn=1:length(g);ifg(n)=0;A=zeros(1,200);%每个值200个点elseg(n)=1A=ones(1,200);endcp=cpA;%s(t),码元宽度200c=cos(f*t);%载波信号mod1=mod1c;%与s(t)等长的载波信号,变为矩阵形式endfigure(1);subplot

2、(3,2,1);plot(cp);gridon;axis(0200*length(g)-22);title(随机二进制信号序列);cm=;mod=;forn=1:length(g);ifg(n)=0;B=ones(1,200);%每个值200个点c=cos(f*t);%载波信号elseg(n)=1;B=ones(1,200);c=cos(f*t+pi);%载波信号endcm=cmB;%s(t),码元宽度200mod=modc;%与s(t)等长的载波信号endtiaoz=cm.*mod;%e(t)调制figure(1);subplot(3,2,2);plot(tiaoz);gridon;axis

3、(0200*length(g)-22);title(2PSK调制信号);figure(2);subplot(3,2,1);plot(abs(fft(cp);axis(0200*length(g)0400);title(原始信号频谱);figure(2);subplot(3,2,2);plot(abs(fft(tiaoz);axis(0200*length(g)0400);title(2PSK信号频谱);运行结果PSK解调代码clearall;closeall;clf;%清除窗口中的图形max=20%定义max长度g=zeros(1,max);g=randint(1,max);%长度为max的随

4、机二进制序列cp=;mod1=;f=2*2*pi;t=0:2*pi/199:2*pi;forn=1:length(g);ifg(n)=0;A=zeros(1,200);%每个值200个点elseg(n)=1;A=ones(1,200);endcp=cpA;%s(t),码元宽度200c=cos(f*t);%载波信号mod1=mod1c;%与s(t)等长的载波信号,变为矩阵形式end%figure(1);subplot(3,2,1);plot(cp);gridon;%axis(0200*length(g)-22);title(随机二进制信号序列);cm=;mod=;forn=1:length(g)

5、;ifg(n)=0;B=ones(1,200);%每个值200个点c=cos(f*t);%载波信号elseg(n)=1;B=ones(1,200);c=cos(f*t+pi);%载波信号endcm=cmB;%s(t),码元宽度200mod=modc;%与s(t)等长的载波信号endtiaoz=cm.*mod;%e(t)调制tz=awgn(tiaoz,10);%信号调制中加入白噪声，信噪比为10jiet=2*mod1.*tz;%同步解调figure(1);subplot(3,2,1);plot(jiet);gridonaxis(0200*length(g)-22);title(相乘后信号波形)f

6、igure(1);subplot(3,2,2);plot(abs(fft(jiet);axis(0200*length(g)0400);title(相乘后信号频谱);%低通滤波器fp=500;fs=700;rp=3;rs=20;fn=11025;ws=fs/(fn/2);wp=fp/(fn/2);%计算归一化角频率n,wn=buttord(wp,ws,rp,rs);%计算阶数和截止频率b,a=butter(n,wn);%计算H(z）jt=filter(b,a,jiet);figure(1);subplot(3,2,3);plot(jt);gridonaxis(0200*length(g)-22

7、);title(经低通滤波器后信号波形)figure(1);subplot(3,2,4);plot(abs(fft(jt);axis(0200*length(g)0400);title(经低通滤波器后信号频谱);%抽样判决form=1:200*length(g);ifjt(m)=0;jt(m)=0;endendfigure(1);subplot(3,2,5);plot(jt);gridonaxis(0200*length(g)-22);title(经抽样判决后信号波形)figure(1);subplot(3,2,6);plot(abs(fft(jt);axis(0200*length(g)04

8、00);title(经抽样判决后信号频谱);运行结果PSK仿真代码clearall;closeall;clf;%清除窗口中的图形max=50%定义max长度g=zeros(1,max);g=randint(1,max);%长度为max的随机二进制序列cp=;mod1=;f=2*2*pi;t=0:2*pi/199:2*pi;forn=1:length(g);ifg(n)=0;A=zeros(1,200);%每个值200个点elseg(n)=1;A=ones(1,200);endcp=cpA;%s(t),码元宽度200c=cos(f*t);%载波信号mod1=mod1c;%与s(t)等长的载波信号

9、,变为矩阵形式endfigure(1);subplot(3,2,1);plot(cp);gridon;axis(020*length(g)-22);title(随机二进制信号序列);cm=;mod=;forn=1:length(g);ifg(n)=0;B=ones(1,200);%每个值200个点c=cos(f*t);%载波信号elseg(n)=1;B=ones(1,200);c=cos(f*t+pi);%载波信号endcm=cmB;%s(t),码元宽度200mod=modc;%与s(t)等长的载波信号endtiaoz=cm.*mod;%e(t)调制figure(1);subplot(3,2,

10、2);plot(tiaoz);gridon;axis(020*length(g)-22);title(2PSK调制信号);figure(2);subplot(3,2,1);plot(abs(fft(cp);axis(020*length(g)0400);title(原始信号频谱);figure(2);subplot(3,2,2);plot(abs(fft(tiaoz);axis(020*length(g)0400);title(2PSK信号频谱);%带有高斯白噪声的信道tz=awgn(tiaoz,10);%信号调制中加入白噪声，信噪比为10figure(1);subplot(3,2,3);pl

11、ot(tz);gridonaxis(020*length(g)-22);title(通过高斯白噪声信道后的2PSK信号);figure(2);subplot(3,2,3);plot(abs(fft(tz);axis(020*length(g)0400);title(加入白噪声的2PSK信号频谱);jiet=2*mod1.*tz;%同步解调figure(1);subplot(3,2,4);plot(jiet);gridonaxis(020*length(g)-22);title(相乘后信号波形)figure(2);subplot(3,2,4);plot(abs(fft(jiet);axis(02

12、0*length(g)0400);title(相乘后信号频谱);%低通滤波器fp=500;fs=700;rp=3;rs=20;fn=11025;ws=fs/(fn/2);wp=fp/(fn/2);%计算归一化角频率n,wn=buttord(wp,ws,rp,rs);%计算阶数和截止频率b,a=butter(n,wn);%计算H(z）jt=filter(b,a,jiet);figure(1);subplot(3,2,5);plot(jt);gridonaxis(020*length(g)-22);title(经低通滤波器后信号波形)figure(2);subplot(3,2,5);plot(ab

13、s(fft(jt);axis(020*length(g)0400);title(经低通滤波器后信号频谱);%抽样判决form=1:200*length(g);ifjt(m)=0;jt(m)=0;endendfigure(1);subplot(3,2,6);plot(jt);gridonaxis(020*length(g)-22);title(经抽样判决后信号波形)figure(2);subplot(3,2,6);plot(abs(fft(jt);axis(020*length(g)0400);title(经抽样判决后信号频谱);%误码率随信噪比的变化曲线snr=0:2:20;%信噪比范围len

14、_snr=length(snr);fori=1:len_snrSNR_eb=exp(snr(i)*log(10)/10);theo_err_prb(i)=(1/2)*erfc(sqrt(SNR_eb);%PSK系统的误码率与信噪比的关系endfigure(3);semilogy(snr,theo_err_prb,ko-);gridon;title(误码率随信噪比的变化曲线);仿真结果心得体会经过这一周的通信原理实验课的学习，让我收获多多。但在这中间，我也发现了我存在的很多不足。我的动手能力还不够强，当有些实验需要很强的动手能力时我还不能从容应对；我的探索方式还有待改善，当面对一些复杂的实验时我

15、还不能很快很好的完成；我的数据处理能力还得提高，当眼前摆着一大堆复杂数据时我处理的方式及能力还不足，不能用最佳的处理手段使实验误差减小到最小程度在整整一个星期的日子里，可以说得是苦多于甜，但是可以学到很多很多的的东西，同时不仅可以巩固了以前所学过的知识，而且学到了很多在书本上所没有学到过的知识。通过这次课程设计使我懂得了理论与实际相结合是很重要的，只有理论知识是远远不够的，只有把所学的理论知识与实践相结合起来，从理论中得出结论，才能真正为社会服务，从而提高自己的实际动手能力和独立思考的能力总之，通信原理课程设计让我收获颇丰，同时也让我发现了自身的不足。在实验课上学得的，我将发挥到其它中去，也将

16、在今后的学习和工作中不断提高、完善；在此间发现的不足，我将努力改善，通过学习、实践等方式不断提高，克服那些不应成为学习、获得知识的障碍。在今后的学习、工作中有更大的收获，在不断地探索中、在无私的学习、奉献中实现自己的人身价值！2方案分析与方案确定2.1数字基带信号的传输本次课设中涉及到的通信技术为信息传输技术中“数字调制与解调。通过数字调制与解调，可以把信号变成适合在信道中传输的信号。调制与解调技术可以提高信号传输的可靠性和有效性，增强信号的抗干扰能力，减小发射天线的尺寸，削减传输成本。基于以上特点，数字信号的调制与借条在通信领域有广泛的应用。调制一般可分为两种类型：(1)利用模拟方法去实现数

17、字调制，即把数字基带信号当作模拟信号的特殊情况来处理；(2)利用数字信号的离散取值特点键控载波，从而实现数字调制键控法。如：2ASK、2FSK、2PSK等。键控法的特点：数字电路实现，调制变换速率快，调整测试方便，体积小和设备可靠性高。移相键控PSK在数字通信系统中是一种重要的调制方式，其抗噪性能和信道频带利用率均优于ASK和FSK，即使在有衰落的信道中也能获得很好的效果。因而在实际的高速数据传输系统中得到广泛的应用。基于以上分析，于上述分析，本次课设的内容大致可以分为调制和解调两部分。2.2数字基带信号的调制调制是用调制信号去控制载波一些参数的过程。在数字信号的调制类型里，大致有FSK,AS

18、K,PSK三种。其中移相键控PSK在数字通信系统中是一种重要的调制方式，其抗噪性能和信道频带利用率均优于ASK和FSK，即使在有衰落的信道中也能获得很好的效果。因而在实际的高速数据传输系统中得到广泛的应用。数字相位调制又称频相键控，简记PSK,二进制移相键控记作2PSK。它是利用载相位的变化来传送数字信息的。通常有两种类型：(1)绝对相移(2PSK或BPSK)(2)相对相移(差分相移/2DPSK或DBPSK)此次课设中，我们用到的是绝对相移法。所谓绝对移相是以载波的不同相位的绝对值来直接传送相应二进制数字信号的一种调制方式，简称2PSK。通常用已调载波的“0”相和“”相，分别表示二进制数字的“1”和“0”绝对相移2PSK的电路主要由五部分组成。载波振荡器产生高频载波，可知直接输入到开关电路。高频载波有一路需要输入到反相器，使其产生180度相移，然后输入到开关电路。这样便有两路信号输入到开关电路。实际使用调制器时，用乘法器和加法器代替，将输入的M序列信号一分为二，一路为原信号，另一路为取反后的M序列。两路信号分别输入到乘法器，经过加法器相加后，便可输出2PSK信号。绝对相移2PSK的工作原理：条件：以未调载波的“0”相表示“1