PSK系统设计与仿真

课程设计报告

课程设计题目：PSK系统设计与仿真

学号：123

学生姓名：陈

专业：通信工程

班级：1221

指导教师：涂老师

2015年1 月12 日

目录

一引言

1.1 课程设计目的 (3)

1.2 课程设计要求 (3)

1.3 课程设计注意事项 (3)

二 PSK信号调制解调模型的建立及分析

2.1 PSK信号解调模型的建立 (4)

2.2 PSK信号解调模型的建立 (4)

2.3 PSK调制过程分析 (5)

2.4 PSK解调过程分析 (7)

三高斯白噪声对系统影响分析............................. .......... .. (10)

四PSK调制系统的抗噪声性能分析 (12)

五PSK系统的仿真 (13)

六实验总结 (17)

.

一引言

本课程设计用于实现PSK系统设计与仿真,移动通信的迅速发展，离不开很多关键技术的支持与应用，数字调制在通信领域中就发挥着重大的作用,为了使数字信号在带通信道中传输，必须使用数字基带信号对载波进行调制，以使信号与信道的特性匹配.由于PSK在生活中有着广泛的应用，本课题主要介绍了PSK波形的产生和仿真过程.

1.1、课程设计目的

本课程设计是实现PSK系统仿真。在此次课程设计中，我将通过多方搜集资料与分析，来理解PSK系统仿真的具体过程和它在MATLAB中的实现方法。预期通过这个阶段的研习，更清晰地认识PSK系统仿真原理，同时加深对MATLAB这款通信仿真软件操作的熟练度，并在使用中去感受MATLAB的应用方式与特色。利用自主的设计过程来锻炼自己独立思考，分析和解决问题的能力，为我今后的自主学习研究提供具有实用性的经验

1.2、课程设计要求

（1）熟悉MATLAB中M文件的使用方法，掌握PSK系统仿真原理，以此为基础用M文件编程实现PSK系统仿真。

（2）绘制出PSK系统仿真在时域和频域中的波形，观察两者在解调前后的变化，通过对分析结果来加强对PSK系统仿真原理的理解。

（3）对信号分别叠加大小不同的噪声后再进行解调，绘制出解调前后信号的时域和频域波形，比较未叠加噪声时和分别叠加大小噪声时解调信号的波形有何区别，由所得结果来分析噪声对信号解调造成的影响。

（4）在老师的指导下，独立完成课程设计的全部内容，并按要求编写课程设计论文，文中能正确阐述和分析设计和实验结果。

1.3、课程设计注意事项

(1)所有的仿真用MATLAB或VC程序实现（如用MATLAB则只能用代码的形式，不能用SIMULINK实现）

(2)系统经过的信道都假设为高斯白噪声信道。

（3）模拟调制要求用程序画出调制信号，载波，已调信号、解调信号的波形，数字调制要求画出误码率随信噪比的变化曲线。

二 PSK 信号调制解调模型的建立

2.1、PSK 信号调制模型的建立

相移键控（PSK ）是利用载波的相位变化来传递数字信息，而振幅和频率保持不变的一种数字信号传递方法。PSK 的调制原理框图如下图所示，与ASK 信号的产生方法比较，只是对s 的要求不同，在ASK 中s 是单极性的，而在PSK 中S 是双极性的基带信号。

)(t s 双极性不归零 )(2t e P S K

t c ωc o s

图1 模拟调制方法

开关电路

0 )(2t e PSK

π

)(t s

图2 键控法

2.2、PSK 信号解调模型的建立

PSK 信号的解调通常采用相干解调法，解调器原理框图如下图。在相干解调中，怎样得到与接收的PSK 信号同频同相的相干载波是一个关键的问题。

)(2t e P S K

a c d e

t c ωc o s b 定时脉冲

图3 PSK 信号的解调原理框图

码型变换

乘法器 t c ωcos o

180相移 带通 滤波器 相乘器 低通 滤波器 抽样 判决器

2.3、PSK 调制过程分析

根据PSK 调制的定义，设初始相位0和π分别表示二进制“1”和“0”。因此，PSK 信号的时域表达式为

2()cos()PSK c n e t A t ωϕ=+，其中，n ϕ表示第n 个符号的绝对相位：

因此，可得到下式

典型波形如下图所示

图4 PSK 信号的时间波形

由于两种码元的波形相同，极性相反，故2PSK 信号可以表述为一个双极性全占空矩形脉冲序列与一个正弦载波的相乘，即

其中，

这里，g(t)是脉宽为TS 的单个矩形脉冲，而an 的统计特性为

即发送二进制符号“0”时（an 取+1），)(2t e P S K 取0相位；

发送二进制符号“1”时（ an 取 -1），)(2t e P S K 取π相位。这种以载波的不同相位直接去表示相应二进制数字信号的调制方式，称为二进制绝对相移方式，且其带宽为基带信号的两倍。

调制过程产生的代码和波形如下

2PSK Acos ,()Acos ,1c c t P

e t t P ωω⎧=⎨--⎩概率为概率为00,1n ϕπ⎧=⎨⎩，发送“”时发送“”时1101t 0s T ()2PSK

()cos c e t s t t ω=()()n s n s t a g t nT =-∑1,1,1n P a P ⎧=⎨--⎩概率为概率为

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);

end

cp=[cp A]; %s(t),码元宽度200

c=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); %载波信号

end

cm=[cm B]; %s(t),码元宽度200

mod=[mod c]; %与s(t)等长的载波信号

end

tiaoz=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信号频谱');