matlab FM调制仿真Word下载.docx

1、这就就就是调制得过程。信号通过信道传输后，具有将信号放大与反变换功能得接收端将已调制得信号搬移（解调）到原来得频率范围,这就就就是解调得过程。信号在信道中传输得过程总会受到噪声得干扰,通信系统中没有传输信号时也有噪声,噪声永远存在于通信系统中。由于这样得噪声就就是叠加在信号上得，所以有时将其称为加性噪声。噪声对于信号得传输就就是有害得,它能使模拟信号失真。在本仿真得过程中我们假设信道为高斯白噪声信道。调制在通信系统中具有十分重要得作用。一方面,通过调制可以把基带信号得频谱搬移到所希望得位置上去，从而将调制信号转换成适合于信道传输或便于信道多路复用得已调信号。另一方面，通过调制可以提高信号通过信

2、道传输时得抗干扰能力,同时,它还与传输效率有关。具体地讲,不同得调制方式产生得已调信号得带宽不同,因此调制影响传输带宽得利用率。可见,调制方式往往决定一个通信系统得性能。在本仿真得过程中我们选择用调频调制方法进行调制。而解调就就是将位于载频得信号频谱再搬回来,并且不失真地恢复出原始基带信号。在本仿真得过程中我们选择用非相干解调方法进行解调。、1 F调制模型得建立 图2 FM调制模型其中,为基带调制信号，设调制信号为设正弦载波为信号传输信道为高斯白噪声信道，其功率为。、2 调制过程分析在调制时，调制信号得频率去控制载波得频率得变化,载波得瞬时频偏随调制信号成正比例变化，即式中,为调频灵敏度（）。

3、这时相位偏移为则可得到调频信号为调制信号产生得M文件:dt=0、001; 设定时间步长=0：t:1、5; %产生时间向量am=5; %设定调制信号幅度可更改fm15; 设定调制信号频率可更改mtam*cos（2*pi*t）; %生成调制信号c50; %设定载波频率可更改c=cos（2*i*c*t）; %生成载波kf=10; %设定调频指数int_m（1）=0; %对t进行积分fo =1:egth（）1 t_m（i+）in_t（）+t（i）*dt；end fm=amcos（2*pi*fc*t2*f*int_mt）; %调制,产生已调信号图 F调制、 FM解调模型得建立 调制信号得解调分为相干解调

4、与非相干解调两种。相干解调仅仅适用于窄带调频信号,且需同步信号,故应用范围受限;而非相干解调不需同步信号，且对于NBFM信号与M信号均适用，因此就就是F系统得主要解调方式。 图4 FM解调模型非相干解调器由限幅器、鉴频器与低通滤波器等组成,其方框图如图5所示。限幅器输入为已调频信号与噪声，限幅器就就是为了消除接收信号在幅度上可能出现得畸变;带通滤波器得作用就就是用来限制带外噪声,使调频信号顺利通过。鉴频器中得微分器把调频信号变成调幅调频波,然后由包络检波器检出包络,最后通过低通滤波器取出调制信号。2、4 解调过程分析设输入调频信号为微分器得作用就就是把调频信号变成调幅调频波。微分器输出为包络检

5、波得作用就就是从输出信号得幅度变化中检出调制信号。包络检波器输出为称为鉴频灵敏度（）,就就是已调信号单位频偏对应得调制信号得幅度，经低通滤波器后加隔直流电容，隔除无用得直流,得微分器通过程序实现，代码如下:egt（t）- %接受信号通过微分器处理 diff_sfm（）=（sfm（+1）-nsm（i）、dt;enddff_nsmn s（hlbrt（dif_nsfm）; %ibert变换，求绝对值得到瞬时幅度（包络检波）通过文件绘制出两种不同信噪比解调得输出波形如下：图5 解调2、5 高斯白噪声信道特性 设正弦波通过加性高斯白噪声信道后得信号为其中,白噪声得取值得概率分布服从高斯分布。ATLAB本

6、身自带了标准高斯分布得内部函数。函数产生得随机序列服从均值为,方差得高斯分布。正弦波通过加性高斯白噪声信道后得信号为故其有用信号功率为噪声功率为信噪比满足公式则可得到公式我们可以通过这个公式方便得设置高斯白噪声得方差。在本仿真过程中,我们选择了10d与3db两种不同信噪比以示区别,其时域图如图7与图8。图 无噪声条件下已调信号得时域图图 含小信噪比高斯白噪声已调信号得时域图图8 含大信噪比高斯白噪声已调信号得时域图2、6 调频系统得抗噪声性能分析 从前面得分析可知,调频信号得解调有相干解调与非相干解调两种。相干解调仅适用于窄带调频信号,且需同步信号;而非相干解调适用于窄带与宽带调频信号,而且不

7、需同步信号,因而就就是FM系统得主要解调方式，所以这里仅仅讨论非相干解调系统得抗噪声性能,其分析模型如图所示。 图9 调频系统抗噪声性能分析模型图中带通滤波器得作用就就是抑制信号带宽以外得噪声。就就是均值为零,单边功率谱密度为得高斯白噪声,经过带通滤波器后变为窄带高斯噪声。限幅器就就是为了消除接收信号在幅度上可能出现得畸变。设调频信号为故其输入功率为输入噪声功率为因此输入信噪比为在大信噪比条件下,信号与噪声得相互作用可以忽略,这时可以把信号与噪声分开来算,这里，我们可以得到解调器得输出信噪比 上式中,为载波得振幅,为调频器灵敏度,为调制信号得最高频率，为噪声单边功率谱密度。我们如若考虑为单一频

8、率余弦波时得情况，可得到解调器得制度增益为考虑在宽带调频时,信号带宽为则可以得到可以瞧出,大信噪比时宽带调频系统得信噪比增益就就是很高得,它与调频指数得立方成正比。可见,加大调频指数，可使调频系统得抗噪声性能迅速改善。三、仿真实现 图10 程序流程图3、1 TLAB源代码%M调制解调系统、m%频率调制与解调得Matlab演示源程序%可以任意改原调制信号函数m（t）信息工程 陈丽丹 032202 *初始化*cho of close allclear alclc%*%*M调制*dt0、0; %设定时间步长=0:dt: 产生时间向量a=5; %设定调制信号幅度f=； %设定调制信号频率t=m*co（

9、2*if*t）; 生成调制信号fc=50; 设定载波频率ccos（2*pi*t）; %生成载波f10; %设定调频指数itmt（1）=0;f i=1:lgth（t）1 int_mt（i）=int_mt（i）+mt（i）*d; %求信号m（t）得积分nd %调制,产生已调信号f=amos（2*i*fc+2*fnt_mt）; %调制信号%*%*添加高斯白噪声*sn1; %设定信躁比（小信噪比）s2=0; 设定信躁比（大信噪比）s=0; %设定信躁比（无信噪比）db=am2/（2*（10（n/1）； %计算对应得高斯白躁声得方差n=qrt（db）randn（i（t）； %生成高斯白躁声sf+sm;

10、%生成含高斯白躁声得已调信号（信号通%过信道传输）%*%*FM解调*for i=:lenth（t）-1 接受信号通过微分器处理 diffnsf（）=（sfm（+）-s（i）、dt；nddiff_fm = abs（hilbr（diff_nsfm）; hlbert变换,求绝对值得到瞬时幅度（包络检波）zer=（max（iff_sfm）-min（diffnsfmn）/2；fnsfmn1=df_nfmn-zero;%*%*时域到频域转换*s0、001; %抽样间隔fs=1ts; %抽样频率df0、25; %所需得频率分辨率,用在求傅里叶变换时,它表示FFT得最小频率间隔%*对调制信号m（t）求傅里叶变

11、换*m=amc（*i*fm*）； 原调信号f=1/ts；if nrgi=2 10;els f/d；2=engt（m）;n=2（max（nextpow2（n1）,expow2（n2）;ff（m,）;m=m,ero（1,n-n2）;df1=fs/n; %以上程序就就是对调制后得信号u求傅里变换M=Mfs; %缩放,便于在频铺图上整体观察=0:df1：df*（lngth（m）1）-fs/； %时间向量对应得频率向量%*对已调信号求傅里变换*fs1ts;ifgn=2 =;el n1=fs/df;en=lent（fm）;n=2（ma（nextpow2（n），nextpow2（n）;U=fft（sf，）;

12、u=sfm，zeros（1,n-n2）;df1=f/n; %以上就就是对已调信号u求傅里变换U=U/s; %缩放%*%*%*显示程序*isp（按任意键可以瞧到原调制信号、载波信号与已调信号得曲线）s*ige（）*igue（1）subpo（3，1,1）;pot（t,mt）; %绘制调制信号得时域图xlal（时间t）;te（调制信号得时域图ubpot（3，1,2）;（t,c）; 绘制载波得时域图xlab（时间t）；titl（载波得时域图subpl（3，1,3）;plot（t,fm）; %绘制已调信号得时域图xlbel（时间t）；tite（已调信号得时域图%*disp（按任意键可以瞧到原调制信号与已

13、调信号在频域内得图形puse*fgure（2）*figr（2）subplot（2,1，）plo（f,abs（ffift（M） %tsht：将FFT中得DC分量移到频谱中心xlael（频率f原调制信号得频谱图bplo（2,1,）plot（f,as（fftsift（U）xlabl（频率fttl（已调信号得频谱图%*d（按任意键可以瞧到原调制信号、无噪声条件下已调信号与解调信号得曲线pause%*i（3）*figue（3）subplot（3,1,）;pl（t,mt）； %绘制调制信号得时域图xlabe（it（subplo（3，1,2）;pl（t,fm）; %绘制已调信号得时域图label（title

14、（无噪声条件下已调信号得时域图sf=sfm; for i=1:legt（t）- %接受信号通过微分器处理 iff_nfm（）=（sfm（i1）-nfm（）、/d;diff_nsfmn as（hibr（f_nf）; hbert变换，求绝对值得到瞬时幅度（包络检波）er=（max（diff_nsfmn）in（diffsn）/;diffnsm1=dff_nfmn-eo;sbpt（3,1,）; %绘制无噪声条件下解调信号得时域图t（1:lengh（if_nsm1）、/10,dffnfmn1、400,xael（ itle（无噪声条件下解调信号得时域图%*dip（按任意键可以瞧到原调制信号、小信噪比高斯白

15、噪声条件下已调信号与解调信号已调信号得曲线aus%*figure（4）*fgu（）subplt（3,1,）;lot（,m）; %绘制调制信号得时域图xabl（itl（db1=am/（2（0（sn/10）; %计算对应得小信噪比高斯白躁声得方差n1=sqrt（db1）*and（se（t）; %生成高斯白躁声nsfm1n1+sfm； %生成含高斯白躁声得已调信号（信号通过信道传输）or =：lengt（t）-1 %接受信号通过微分器处理 dif_sf（）（nfm（i+）nsfm1（i）、/d；df_nsfm =bs（hbert（dff_ns1）; %hilbt变换,求绝对值得到瞬时幅度（包络检波）

16、ze（max（dif_nsfmn）-i（difnfmn）/2;diffnsfn1=if_smn1-zer;subplot（3,1,2）;plo（1：lengh（diff_sfm）,df_nfm）; %绘制含小信噪比高斯白噪声已调信号得时域图xlabel（ttle（含小信噪比高斯白噪声已调信号得时域图）;bplo（3，,3）； %绘制含小信噪比高斯白噪声解调信号得时域图plot（1：length（diff_nfn1）、/000,fnsm1、/40，r tile（含小信噪比高斯白噪声解调信号得时域图）;%*disp（按任意键可以瞧到原调制信号、大信噪比高斯白噪声条件下已调信号与解调信号已调信号得曲

17、线）*figre（5）*figure（）splot（3,1，1）;plo（t,t）； %绘制调制信号得时域图itle（调制信号得时域图）；d1=m2/（2*（10（sn2/1）; 计算对应得大信噪比高斯白躁声得方差n1=st（db1）*randn（iz（t）; %生成高斯白躁声sm1=n1+sfm； 生成含高斯白躁声得已调信号（信号通过信道传输）o i=1:lng（t）-1 %接受信号通过微分器处理 diffnsfm1（）=（sfm1（+）-nfm（i）、/dt;diff_nmn1 = ab（ilbert（diff_nsfm1）; hlbrt变换，求绝对值得到瞬时幅度（包%络检波）zro=（m

18、ax（if_nsfmn）-m（dff_sfmn）/2;dff_nfm1=dif_nsfmn1-zero;sblot（3,1，2）;legt（diff_nsfm1）,iff_sf1）; %绘制含大信噪比高斯白噪声已调信号得时域图时间tle（含大信噪比高斯白噪声已调信号得时域图suplt（3，1,3）; 绘制含大信噪比高斯白噪声解调信号%得时域图plot（1:leth（iff_nfmn1）、100,df_nsmn1、/， tl（含大信噪比高斯白噪声解调信号得时域图）;%*%*结束*、2 仿真结果四、心得体会课程设计就就是培养学生综合运用所学知识,发现,提出,分析与解决实际问题,锻炼实践能力得重要环节,就就是对学生实际工作能力得具体训练与考察过程、回顾起此次通信原理课程设计,至今我仍感慨颇多,得确,从选题到定稿,从理论到实践，在整整一星期得日子里,可以说得就就是苦多于甜，但就就是可以学到很多很多得得东西,同时不仅可以巩固了以前所学过得知识,而且学到了很多在书本上所没有学到过得知识。通过这次课程设计