青岛农业大学理学与信息科学学院通信系统仿真课程设计报告设计题目 PSK调制和解调系统学生专业班级学生姓名（学号）指导教师完成时间实习（设计）地点理信学院机房2015年 09月01日一、课程设计目的和任务学生通过本课程的实践，能进一步掌握高级语言程序设计基本概念，掌握基本的程序设计方法；通过设计一个完整的小型程序，初步掌握开发软件所需的需求定义能力、功能分解能力和程序设计能力、代码调试技能；学习编写软件设计文档；为未来的软件设计打下良好的基础。
要求学生掌握所学的程序设计方法的基本知识。
结合实际深入理解高级语言程序设计的基本概念、原理及方法。
运用所学的基础知识开发一个小型的程序，能根据问题的需要构造所需的数据结构，设计适合的算法，解决问题。
掌握设计任务的具体要求，进行设计、调试软件的具体方法、步骤和技巧。
对一个实际课题的软件设计有基本了解，拓展知识面，激发在此领域中继续学习和研究的兴趣，为学习后续课程做准备。
二 PSK 信号调制解调模型的建立2.1、PSK 信号调制模型的建立相移键控（PSK ）是利用载波的相位变化来传递数字信息，而振幅和频率保持不变的一种数字信号传递方法。
PSK 的调制原理框图如下图所示，与ASK 信号的产生方法比较，只是对s 的要求不同，在ASK 中s 是单极性的，而在PSK 中S 是双极性的基带信号。
)(t st c ωcos图1 模拟调制方法开关电路)(2t e PSK图2 键控法2.2、PSK 信号解调模型的建立PSK 信号的解调通常采用相干解调法，解调器原理框图如下图。
在相干解调中，怎样得到与接收的PSK 信号同频同相的相干载波是一个关键的问题。
)(2te PSK e图3 PSK 信号的解调原理框图2.3、PSK 调制过程分析根据PSK 调制的定义，设初始相位0和π分别表示二进制“1”和“0”。
因此，PSK 信号的时域表达式为2()cos()PSK c n e t A t ωϕ=+，其中，n ϕ表示第n 个符号的绝对相位：因此，可得到下式典型波形如下图所示图4 PSK 信号的时间波形由于两种码元的波形相同，极性相反，故2PSK 信号可以表述为一个双极性全占空矩形脉冲序列与一个正弦载波的相乘，即其中，2PSK Acos ,()Acos ,1c c t Pe t t P ωω⎧=⎨--⎩概率为概率为00,1n ϕπ⎧=⎨⎩，发送“”时发送“”时()2PSK()cos c e t s t t ω=()()n s n s t a g t nT =-∑这里，g(t)是脉宽为TS 的单个矩形脉冲，而an 的统计特性为即发送二进制符号“0”时（an 取+1），)(2t e PSK 取0相位；发送二进制符号“1”时（ an 取 -1），)(2t e PSK 取π相位。
这种以载波的不同相位直接去表示相应二进制数字信号的调制方式，称为二进制绝对相移方式，且其带宽为基带信号的两倍。
调制过程产生的代码和波形如下：clear all;close all;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;for n=1:length(g);if g(n)==0;A=zeros(1,200);%每个值200个点else g(n)==1;A=ones(1,200);endcp=[cp A]; %s(t),码元宽度200c=cos(f*t);%载波信号mod1=[mod1 c];%与s(t)等长的载波信号,变为矩阵形式endfigure(1);subplot(3,2,1);plot(cp);grid on;axis([0 200*length(g) -2 2]);title('随机二进制信号序列');cm=[];mod=[];for n=1:length(g);if g(n)==0;B=ones(1,200);%每个值200个点c=cos(f*t); %载波信号else g(n)==1;B=ones(1,200);c=cos(f*t+pi); %载波信号endcm=[cm B]; %s(t),码元宽度200mod=[mod c]; %与s(t)等长的载波信号1,1,1n P a P ⎧=⎨--⎩概率为概率为endtiaoz=cm.*mod;%e(t)调制figure(1);subplot(3,2,2);plot(tiaoz);grid on;axis([0 200*length(g) -2 2]);title('2PSK 调制信号');figure(2);subplot(3,2,1);plot(abs(fft(cp)));axis([0 200*length(g) 0 400]);title('原始信号频谱');figure(2);subplot(3,2,2);plot(abs(fft(tiaoz)));axis([0 200*length(g) 0 400]);title('2PSK 信号频谱');运行结果：图5 二进制、PSK 信号波形和频谱图2.4、PSK 解调过程分析根据PSK 调制的原理框图图6 PSK 信号的解调原理框图带通滤波器的意义是让有用信号（已调信号）通过，滤除一部分噪声，所以有用信号在a 处得到信号为()()cos c a t s t tω=假设相干载波的基准相位与2PSK 信号的调制载波的基准相位一致（通常默认为0相位）。
所以得到下式 带通滤波器相乘器低通滤波器抽样判决器定时脉冲输出)(2t e PSK t c ωcos a b c d e222211()()cos ()cos ()()cos 222PKS c c c c t e t t s t t s t s t t ωωω===+通过低通滤波器后21()()2d t s t =最后通过抽样判决器恢复出数字信号。
但是，由于在2PSK 信号的载波恢复过程中存在着的相位模糊，即恢复的本地载波与所需的相干载波可能同相，也可能反相，这种相位关系的不确定性将会造成解调出的数字基带信号与发送的数字基带信号正好相反，即“1”变为“0”，“0”变为“1”，判决器输出数字信号全部出错。
这种现象称为2PSK 方式的“倒π”现象或“反相工作”。
这也是2PSK 方式在实际中很少采用的主要原因。
另外，在随机信号码元序列中，信号波形有可能出现长时间连续的正弦波形，致使在接收端无法辨认信号码元的起止时刻。
2PSK 信号相干解调各点时间波形如下图所示：图7 2PSK 信号相干解调时各点时间波形为了更直观的了解解调过程，我用MATLAB 绘出解调过程的相关波形，代码和波形如下：clear all;close all;clf; %清除窗口中的图形max=20 %定义max 长度g=zeros(1,max);g=randint(1,max); %长度为max 的随机二进制序列cp=[];mod1=[];f=2*2*pi;a b c d et=0:2*pi/199:2*pi;for n=1:length(g);if g(n)==0;A=zeros(1,200); %每个值200个点else g(n)==1;A=ones(1,200);endcp=[cp A]; %s(t),码元宽度200c=cos(f*t); %载波信号mod1=[mod1 c]; %与s(t)等长的载波信号,变为矩阵形式end%figure(1);subplot(3,2,1);plot(cp);grid on;%axis([0 200*length(g) -2 2]);title('随机二进制信号序列');cm=[];mod=[];for n=1:length(g);if g(n)==0;B=ones(1,200); %每个值200个点c=cos(f*t); %载波信号else g(n)==1;B=ones(1,200);c=cos(f*t+pi); %载波信号endcm=[cm B]; %s(t),码元宽度200mod=[mod c]; %与s(t)等长的载波信号endtiaoz=cm.*mod; %e(t)调制tz=awgn(tiaoz,10); %信号调制中加入白噪声，信噪比为10 jiet=2*mod1.*tz; %同步解调figure(1);subplot(3,2,1);plot(jiet);grid onaxis([0 200*length(g) -2 2]);title('相乘后信号波形')figure(1);subplot(3,2,2);plot(abs(fft(jiet)));axis([0 200*length(g) 0 400]);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);grid onaxis([0 200*length(g) -2 2]);title('经低通滤波器后信号波形') figure(1);subplot(3,2,4);plot(abs(fft(jt)));axis([0 200*length(g) 0 400]);title('经低通滤波器后信号频谱'); %抽样判决for m=1:200*length(g);if jt(m)<0;jt(m)=1;else jt(m)>=0;jt(m)=0;endendfigure(1);subplot(3,2,5);plot(jt);grid onaxis([0 200*length(g) -2 2]);title('经抽样判决后信号波形')figure(1);subplot(3,2,6);plot(abs(fft(jt)));axis([0 200*length(g) 0 400]);title('经抽样判决后信号频谱');运行结果：图8 解调相关波形三、高斯白噪声对系统影响分析当信号经过信道传输时会受到噪声的影响，这是不可避免的。
