基于matlab的数字调制.doc
2ASK、2FSK、2PSK数字调制系统的MATLAB实现及性能分析与比较引言数字带通传输系统为了进行长距离传输,克服传输失真,传输损耗,同时保证带内特性。必须对数字信号进行载波调制,将信号频谱搬移到高频段才能在信道中传输,因而现代通信系统采取数字调制技术。通过数字基带信号对载波某些参量进行控制,使之随机带信号的变化而变化。根据控制载波参量大的不同,数字调制有调幅(ASK),调频FSK,调相PSK三种基本形式。MATLAB用于仿真,分析和修改,还可以应用图形界面功能GUI能为仿真系统生成一个人机交互界面,便于仿真系统的操作,因此采用MATLAB对数字系统进行仿真。通过对系统的仿真,我们可以更加直观的了解数字调制系统的性能()及影响性能的因素,从而便于改进系统,获得更佳的传输性能。关键词数字系统性能ASKFSKPSKMATLAB仿真一数字调制与解调原理112ASK(1)2ASK2ASK就是把频率、相位作为常量,而把振幅作为变量,信息比特是通过载波的幅度来传递的。由于调制信号只有0或1两个电平,相乘的结果相当于将载频或者关断,或者接通,它的实际意义是当调制的数字信号“1时,传输载波;当调制的数字信号为“0“时,不传输载波。公式为00,COS2KASKATTS当,当122FSK2FSK可以看做是2个不同频率的2ASK的叠加,其调制与解调方法与2ASK差不多,主要频率F1和F2,不同的组合产生所要求的2FSK调制信号。公式如下132PSK2PSK以载波的相位变化为基准,载波的相位随数字基带序列信号的1或者0而改变,通常用已经调制完的载波的0或者Π表示数据1或者0,每种相位与之0COS1,212KFSKATATS当,当一一对应。二.数字调制技术的仿真实现本课程设计需要借助MATLAB的M文件编程功能,对2ASK2PSK2FSK进行调制与解调的设计,并绘制出调制与解调后的波形,误码率的情况分析,软件仿真可在已有平台上实现。12ASK代码主函数CLOSEALLCLEARALLN16FC1000000BITRATE1000000N50NOISETINOISE10SIGNALSOURCEN,N生成二进制代码TRANSMITTEDSIGNALASKMODUSIGNAL,BITRATE,FC,N调制后信号SIGNAL1GUSSIANTRANSMITTEDSIGNAL,NOISE加噪声CONFIGUESIGNALDEMOASKSIGNAL1,BITRATE,FC,N,NSOURCE代码FUNCTIONSENDSIGNALSOURCEN,NSENDSIGNALRANDINT1,NBITFORI1LENGTHSENDSIGNALIFSENDSIGNALI0BIT1ZEROS1,NELSEBIT1ONES1,NENDBITBIT,BIT1ENDFIGURE1PLOT1LENGTHBIT,BIT,TITLE TRANSMITTINGOFBINARY ,GRIDONAXIS0,NLENGTHSENDSIGNAL,2,2ENDASKMODU代码FUNCTIONTRANSMITTEDSIGNALASKMODUSIGNAL,BITRATE,FC,NSIGNAL为输入信号,BITRATE为BIT速率,FC调制信号频率,NSIGNAL00101101BITRATE1000000FC1000000N32TLINSPACE0,1/BITRATE,NCSIN2PITFCTRANSMITTEDSIGNALFORI1LENGTHSIGNALTRANSMITTEDSIGNALTRANSMITTEDSIGNAL,SIGNALICENDFIGURE2画调制图PLOT1LENGTHTRANSMITTEDSIGNAL,TRANSMITTEDSIGNALTITLE MODULATIONOFASK GRIDONFIGURE3画频谱实部M0LENGTHTRANSMITTEDSIGNAL1FFFTTRANSMITTEDSIGNALPLOTM,ABSREALF,TITLE ASK_FREQUENCYDOMAINANALYSISREAL GRIDONFIGURE4画频谱虚部PLOTM,IMAGFTITLE ASK_FREQUENCYDOMAINANALYSISIMAG GRIDONENDCHECKRATEPE代码FUNCTIONPEWRONGCHECKRATEPESIGNAL1,SIGNAL2,SRIGHTS0WRONGS0FORKI1S2IFSIGNAL1KISIGNAL2KIRIGHTSRIGHTS1ELSEWRONGSWRONGS1ENDENDPEWRONGWRONGS/WRONGSRIGHTSENDDEMOASK代码FUNCTIONBITSTREAMDEMOASKRECEIVEDSIGNAL,BITRATE,FC,N,NLOADNUMSIGNAL1RECEIVEDSIGNALSIGNAL2ABSSIGNAL1ÕÛÁ÷SIGNAL3FILTERNUM1,1,SIGNAL2LPF,°ÜÂǼ̲¨INFIXLENGTHNUM1/2ÑÓ³Ùʱ¼ÄBITSTREAMLLFC/BITRATENIINLL/2WHILEI05IILLENDFIGURE6SUBPLOT3,1,1接收波形PLOT1LENGTHSIGNAL1,SIGNAL1TITLE WAVEOFRECEIVINGTERMINALINCLUDINGNOISE GRIDONSUBPLOT3,1,2接收整流后波形PLOT1LENGTHSIGNAL2,SIGNAL2TITLE WAVEOFCOMMUTATE GRIDONSUBPLOT3,1,3包络检波波形PLOT1LENGTHSIGNAL3,SIGNAL3TITLE WAVEOFLPF GRIDONBITFORI1LENGTHBITSTREAMIFBITSTREAMI0BIT1ZEROS1,NELSEBIT1ONES1,NENDBITBIT,BIT1ENDFIGURE7解调后的二进制波形PLOTBIT,TITLE BINARYOFRECEIVINGTERMINAL ,GRIDONAXIS0,NLENGTHBITSTREAM,25,25ENDGUSSIAN代码加高斯白噪声FUNCTIONSIGNALGUSSIANTRANSMITTEDSIGNAL,NOISESIGNALSQRT2TRANSMITTEDSIGNALSIGNALAWGNSIGNAL,NOISEFIGURE5PLOT1LENGTHSIGNAL,SIGNALTITLE WAVEINCLUDINGNOISE ,GRIDONEND/FSK主函数代码CLOSEALLCLEARALLN16二进制代码长度F118000000频率1F26000000频率2BITRATE1000000BIT速率N50码元宽度NOISETINOISE10家性噪声大小SIGNALSO