PSK系统设计课程设计报告.docx
《PSK系统设计课程设计报告.docx》由会员分享，可在线阅读，更多相关《PSK系统设计课程设计报告.docx（18页珍藏版）》请在冰点文库上搜索。
PSK系统设计课程设计报告
华南理工大学
通信原理课程设计报告
题目：
2PSK系统仿真
专业：
班级：
姓名：
学号:
日期：
20XX年XX月
1、实验需要材料
MATLAB软件
2、实验要求
完成规定系统的MATLAB编程以及simulink的仿真，基本内容包括：
输入信号，系统中各个关键模块的输出情况。
并调整仿真的参数得到不同的仿真结果。
三、设计原理
2PSK汉语全称：
二进制相移键控。
2PSK是相移键控的最简单的一种形式，它用两个初相相隔为180的载波来传递二进制信息。
所以也被称为BPSK。
Simulink简介：
Simulink是Mathworks公司推出的基于Matlab平台的着名仿真环境Simulin作为一种专业和功能强大且操作简单的仿真工具，目前已被越来越多的工程技术人员所青睐，它搭建积木式的建模仿真方式既简单又直观，而且已经在各个领域得到了广泛的应用。
数字信号的传输方式分为基带传输和带通传输，在实际应用中，大多数信道具有带通特性而不能直接传输基带信号。
为了使数字信号在带通信道中传输，必须使用数字基带信号对载波进行调制，以使信号与信道的特性相匹配。
这种用数字基带信号控制载波，把数字基带信号变换为数字带通信号的过程称为数字调制。
数字调制技术的两种方法：
①利用模拟调制的方法去实现数字式调制，即把数字调制看成是模拟调制的一个特例，把数字基带信号当做模拟信号的特殊情况处理。
②利用数字信号的离散取值特点通过开关键控载波，从而实现数字调制。
这种方法通常称为键控法，比如对载波的相位进行键控，便可获得相移键控（PSK）基本的调制方式。
图1相应的信号波形的示例
101
调制原理：
在二进制数字调制中，当正弦载波的相位随二进制数字基带信号离散变化时，则产生二进制移相键控（2PSK）信号。
2PSK信号调制有两种方法，即模拟调制法和键控法。
通常用已调信号载波的0°和180°分别表示二进制数字基带信号的1和0，模拟调制法用两个反相的载波信号进行调制。
2PSK以载波的相位变化作为参考基准的，当基带信号为0时相位相对于初始相位为0°，当基带信号为1时相对于初始相位为180°。
键控法，是用载波的相位来携带二进制信息的调制方式。
通常用0°和180°来分别代表0和1。
其时域表达式为：
其中，2PSK的调制中an必须为双极性码。
两种方法原理图分别如图1-1和图1-2所示。
图1-1模拟调制法原理图
图1-2键控法原理图
2PSK信号的时间波形
解调原理
2PSK信号的解调方法是相干解调法，如图2-1所示。
由于PSK信号本身就是利用相位传递信息的，所以在接收端必须利用信号的相位信息来解调信号。
由于2PSK的幅度是恒定的，必须进行相干解调。
经过带通滤波的信号在相乘器中与本地载波相乘，然后用低通滤波器滤除高频分量，在进行抽样判决。
判决器是按极性来判决的。
即正抽样值判为1，负抽样值判为0。
图2-12PSK相干解调法原理图
假设相干载波的基准相位与2PSK信号的基准一致（通常默认为0相位，如图2-2所示）。
但是由于2PSK信号的载波回复过程中存在着180°的相位模糊，即恢复的本地载波与所需相干载波可能相同，也可能相反，这种相位关系的不确定性将会造成解调出的数字基带信号与发送的基带信号正好相反，即“1”变成“0”，“0”变成“1”，判决器输出数字信号全部出错。
这种现象称为2PSK方式的“倒π”现象或“反相工作”。
但在本次仿真中是直接给其同频同相的载波信号，所以不存在此问题。
图2-2相干解调中各点波形图
相关公式：
2PSK信号在一个码元的持续时间Ts内可以表示为
u1T（t）发送“1”时
ST（t）=
uoT（t）=-u1T（t）发送“0”时
其中Acosωct，0u1T（t）=
0，其他
设发送端发出的信号如上式所示，则接收端带通滤波器输出波形y（t）为
[a+nc（t）]cosωct-ns（t）sinωct，发送“1”时
y（t）=
[-a+nc（t）]cosωct-ns（t）sinωct，发送“0”时
y（t）经过想干解调（相乘—低通）后，送入抽样判决器的输入波形为
a+nc（t），发送“1”时
x（t）=
-a+nc（t），发送“0”时
由最佳判决门限分析可知，在发送“1”和“0”概率相等时，即P
（1）=P（0）时，最佳门限b*=0.此时，发“1”而错判为“0”的概率为
P（0/1）=P（x≦0）=∫0-∞f1（x）dx=1/2erfc（
）
式中：
r=a2/2σ2n
同理，发“0而错判为“1”的概率为
P（1/0）=P（x＞0）=∫0-∞f0（x）dx=1/2erfc（
）
故2PSK信号相干解调系统的总误码率为
Pe=P
（1）P（0/1）+P（0）P（0/1）=1/2erfc（
）
在大信噪比（r>>1）的条件下,上式可近似为
Pe≈e-r/2
四、源程序及相应实验结果
1、相干调解法
i=10;
j=5000;
fc=4;%载波频率
fm=i/5;%码元速率
B=2*fm;
t=linspace（0,5,j）;
a=round（rand（1,i））;%随机序列,基带信号
figure
（1）;
stem（a）;
st1=t;
forn=1:
10
ifa（n）<1;
form=j/i*（n-1）+1:
j/i*n
st1（m）=0;
end
else
form=j/i*（n-1）+1:
j/i*n
st1（m）=1;
end
end
end
figure
（2）;
subplot（411）;
plot（t,st1）;
title（'基带信号st1'）;
axis（[0,5,-1,2]）;%由于PSK中的是双极性信号，因此对上面所求单极性信号取反来与之一起构成双极性码
st2=t;
fork=1:
j;
ifst1（k）>=1;
st2（k）=0;
else
st2（k）=1;
end
end;
subplot（412）;
plot（t,st2）;
title（'基带信号反码st2'）;
axis（[0,5,-1,2]）;
st3=st1-st2;
subplot（413）;
plot（t,st3）;
title（'双极性基带信号st3'）;
axis（[0,5,-2,2]）;
s1=sin（2*pi*fc*t）;
subplot（414）;
plot（s1）;
title（'载波信号s1'）;
e_psk=st3.*s1;
figure（3）;
subplot（511）;
plot（t,e_psk）;
title（'e_2psk'）;
noise=rand（1,j）;
psk=e_psk+noise;%加入噪声
subplot（512）;
plot（t,psk）;
title（'加噪后波形'）;
psk=psk.*s1;%与载波相乘
subplot（513）;
plot（t,psk）;
title（'与载波s1相乘后波形'）;
lpf=fir1（101,[2/1000,52/1000]）;%低通滤波器设置
st=filter（lpf,1,psk）;%经过低通滤波器后的信号
subplot（514）;
plot（t,st）;
title（'低通滤波后波形'）;
form=0:
i-1;
ifst（1,m*500+250）<0;
forj=m*500+1:
（m+1）*500;
st（1,j）=0;
end
else
forj=m*500+1:
（m+1）*500;
st（1,j）=1;
end
end
end
subplot（515）;
plot（t,st）;
axis（[0,5,-1,2]）;
title（'抽样判决后波形'）
2、模拟调制法
clearall
closeall
i=10;
j=5000;
fc=4;%载波频率
fm=i/5;%码元速率
B=2*fm;
t=linspace（0,5,j）;
a=round（rand（1,i））;
figure（4）;
stem（a）;
st1=t;
forn=1:
10
ifa（n）<1;
form=j/i*（n-1）+1:
j/i*n
st1（m）=0;
end
else
form=j/i*（n-1）+1:
j/i*n
st1（m）=1;
end
end
end
figure（5）;
subplot（321）;
plot（t,st1）;
title（'绝对码'）;
axis（[0,5,-1,2]）;
b=zeros（1,i）;%全零矩阵
b
（1）=a
（1）;
forn=2:
10
ifa（n）>=1;
ifb（n-1）>=1
b（n）=0;
else
b（n）=1;
end
else
b（n）=b（n-1）;
end
end
st1=t;
forn=1:
10
ifb（n）<1;
form=j/i*（n-1）+1:
j/i*n
st1（m）=0;
end
else
form=j/i*（n-1）+1:
j/i*n
st1（m）=1;
end
end
end
subplot（323）;
plot（t,st1）;
title（'相对码st1'）;
axis（[0,5,-1,2]）;
st2=t;
fork=1:
j;
ifst1（k）>=1;
st2（k）=0;
else
st2（k）=1;
end
end;
subplot（324）;
plot（t,st2）;
title（'相对码反码st2'）;
axis（[0,5,-1,2]）;
s1=sin（2*pi*fc*t）;
subplot（325）;
plot（s1）;
title（'载波信号s1'）;
s2=sin（2*pi*fc*t+pi）;
subplot（326）;
plot（s2）;
title（'载波信号s2'）;
d1=st1.*s1;
d2=st2.*s2;
figure
（2）;
subplot（411）;
plot（t,d1）;
title（'st1*s1'）;
subplot（412）;
plot（t,d2）;
title（'st2*s2'）;
e_dpsk=d1+d2;
subplot（413）;
plot（t,e_dpsk）;
title（'调制后波形'）;
noise=rand（1,j）;
dpsk=e_dpsk+noise;%加入噪声
subplot（414）;
plot（t,dpsk）;
title（'加噪声后信号'）;
dpsk=dpsk.*s1;%与载波s1相乘
figure（3）;
subplot（411）;
plot（t,dpsk）;
title（'与载波相乘后波形'）;
lpf=fir1（101,[2/1000,52/1000]）;%低通滤波器设置
k=filter（lpf,1,dpsk）;%经过低通滤波器后的信号
subplot（412）;
plot（t,k）;
title（'低通滤波后波形'）;
st=zeros（1,i）;%全零矩阵
form=0:
i-1;
ifk（1,m*500+250）<0;
st（m+1）=0;
forj=m*500+1:
（m+1）*500;
k（1,j）=0;
end
else
forj=m*500+1:
（m+1）*500;
st（m+1）=1;
k（1,j）=1;
end
end
end
subplot（413）;
plot（t,k）;
axis（[0,5,-1,2]）;
title（'抽样判决后波形'）
dt=zeros（1,i）;%全零矩阵
dt
（1）=st
（1）;
forn=2:
10;
if（st（n）-st（n-1））<=0&&（st（n）-st（n-1））>-1;
dt（n）=0;
else
dt（n）=1;
end
end
st=t;
forn=1:
10
ifdt（n）<1;
form=j/i*（n-1）+1:
j/i*n
st（m）=0;
end
else
form=j/i*（n-1）+1:
j/i*n
st（m）=1;
end
end
end
subplot（414）;
plot（t,st）;
axis（[0,5,-1,2]）;
title（'码反变换后波形'）;
st=real（st）;
运行示例：
方法一效果图：
方法二效果图：
参考文献
[1]李晓峰.通信原理：
清华大学出版社
五、个人体会及感想
这次课程设计对我的自身能力有了进一步了解。
第一点，这进一步端正了我的学习态度，学会了实事求是、严谨的作风，提高了动手能力。
也要对自己严格要求，不能够一知半解，要力求明明白白。
浮躁的性格对于搞设计来说或者学习是致命的，一定要静下心来，踏实的做事。
第二点，我觉得动手之前，头脑里必须清楚应该怎么做，这一点是很重要的，所谓三思而后行。
第三点，在这两周的时间里，深深感觉到我对所学的知识的了解还是很局限、很肤浅的，在以后的时间里还要继续加深对这门课程的理解，加强知识储备。
在这次课程设计中我们遇到了许多的困难，由于粗心大意出了一些简单的错误，浪费了许多时间去改正。
还好在同学和老师的帮助下，给我指出了错误的原因以及改正的方法，我们组才顺利的完成了本次课程设计。
通过这次课程设计，我学到了很多书本上没有的知识。
锻炼了我们独立思考问题、分析问题、解决问题的能力。
而且本次设计有自己和本组成员共同完成。
加强了和别人沟通的能力以及团队精神，对我们走向社会是个很好的锻炼。
这个课程设计完成仓促，在编程过程中，我发现自己的程序还有很多地方可以完善，其中若有不足之前，请老师指出，我将及时改正。
通信原理课程设计评语
完
成
情
况
优秀
良好
中等
及格
不及格
出勤
方案设计
安装调试
答辩
报告
总
成
绩
升级会员 