2PSK调制解调系统的设计与仿真

《2PSK调制解调系统的设计与仿真》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、郑州航空工业管理学院《电子信息系统仿真》课程设计13级电子信息工程专业81班级题目2PSK调制解调系统设计与仿真姓名韩啟典学号131308109指导教师王赴二O—五年12月10日一,MAT3B软件简介MATLAB（矩阵实验室）是一种专业的计算机程序，它是美国MathWorks公司出品的商业数学软件，用于算法开发，数据可视化，数据分析以及工程科学的矩阵数学运算。在以后几年里，逐渐发展为一种极其灵活的计算体系，用于解决各种重要的技术问题。它Mathematica以及Maple并称为三大数学软件。Matlab程序

2、执行MATLAB语言，并提供了一个极其广泛的预定义函数库，这样就使得技术工作变的简单高效。MATLAB是一个庞大的程序，拥有难以置信的各种丰富的函数,基本的MATLAB语言已经拥有了超过1000多个函数，而它的工具包带有更多的函数，由此扩展了它在许多专业领域的能力。二，理论分析2.1,2PSK调制解调系统设计与仿真的原理调制原理：2PSK调制器可以采用相乘器，也可以采用相位选择器。开关电路2PSK信号的调制原理框图解调原理：2PSK信号的解调方法是相干解调法。由于PSK信号本身就是利用相位传递信息的，所以在

3、接收端必须利用信号的相位信息来解调信号。经过带通滤波的信号在相乘器中与本地载波相乘，然后用低通滤波器滤除高频分量，再进行抽样判决。2psk信号的解调原理框图2.2,程序清单clearall;closeall;fs=7e4;%抽样频率fm=14e3;%基带频率n二3*(7*fs/fm);final=(l/fs)*(n-l);fc=3e4;%载波频率t=O:l/fs:(final);Fn=fs/2;%奈奎斯特频率%用正弦波产纶方波%=================================twopi_f

4、c_t=2*pi*fm*t;A=2;phi=O;x=A*cos(twopi_fc_t+phi);%方波am=3;x(x>O)=am;x(x<0)=-3;figure(l)subplot(321);plot(t,x);axis([O2e-4-55]);title。基带信号J;gridoncar二cos(2*pi*fc*t);%载波ask=x.*car;%载波调制subplot(322);plot(t,ask);axis([O20e-5-33]);title(PSK信号');gridon;%==========

5、=======vn=0.1;noise=vn*(randn(size(t)));%产生噪音subplot(323);plot(t,noise);gridon;title(噪音信号J;axis([00.1e-2-0.30.3]);askn=(ask+noise);%调制后加噪subplot(324);plot(t,askn);axis([020e-5-33]);titled加噪后信号)；gridon;%带通滤波%========================================ffiW=40e3

6、;f=[0:3e3:4e5];w=2*pi*f/fs;z=exp(w*j);BW=2*pi*fBW/fs;a=.8547;%BW=2(l-a)/sqrt(a)p=(jA2*aA2);gain=135;Hz=gain*(z+1).*(z-1)./(z.A2-(pi));subplot(325);plot(f,abs(Hz));titled带通滤波器J;axis([O25e4070]);gridon;Hz(Hz==0)=10A(8);%avoidlog(0)subplot(326);plot(f,20*log1

7、0(abs(Hz)));gridon;titleC'Receiver-3dBFilterResponse*);axis([025e41038]);%滤波器系数a=[l00.7305];%[1Op]b=[0.1350-0.135];%gain*[l0-1]faskn=filter(b,a,askn);figure(2)subplot(321);plot(t,faskn);axis([0200e-5-1.51.5]);titled通过带通滤波后输出*);gridon;cm=faskn.*car;%解调subpl

8、ot(322);plot(t,cm);axis([0200e-5-1.51.5]);gridon;title(通过相乘器后输出*);%低通滤波器%=====二=========二========p=0.72;gain1=0.14;%gain=(1-p)/2Hz1=gain1*(z+1)./(z-(p));subplot(323);Hzl(Hz1==0)=10A(-8);%avoidlog(O)plot(f,20*