FM调制解调系统设计与仿真汇总

1.1FM调制模型的建立1.2调制过程分析制信号m(t)成正比例变化，即=Kfm(t)则可得到调频信号为%设定时间步长mt=am*cos(2*pi*fm*t);ct=cos(2*pi*fc*t);int_mt(i+1)=int_mt(i)+mt(i)*dt;sfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt);%产生时间向量%设定调制信号幅度←可更改%设定调制信号频率←可更改%生成调制信号%设定载波频率←可更改%设定调频指数%调制，产生已调信号01000.511.5010-100.5001000.501.3FM解调模型的建立1.4解调过程分析dHzK称为鉴频灵敏度(VdHzt(包络检波)5000.505000.505000.501.5高斯白噪声信道特性设正弦波通过加性高斯白噪声信道后的信号为A22N=2N信噪比S满足公式N则可得到公式2=B5000.505000.505000.50图6无噪声条件下已调信号的时域图5000.50声已调信号的时域图0050010001500声解调信号的时域图000.50图7含小信噪比高斯白噪声已调信号的时域图5000.511.5声已调信号的时域图0050010001500声解调信号的时域图5000.50图8含大信噪比高斯白噪声已调信号的时域图初始化程序初始化程序FM调制程序添加高斯白噪声程序FM解调程序绘制频谱图程序3.1MATLAB源代码FM统.m%*****************初始化******************offlearall%*****************************************%****************FM调制*******************%设定时间步长生时间向量%设定调制信号幅度mt=am*cos(2*pi*fm*t);ct=cos(2*pi*fc*t);intmti+1)=int_mt(i)+mt(i)*dt;sfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt);%设定调制信号频率%生成调制信号%设定载波频率载波定调频指数%调制，产生已调信号调制信号%*****************************************%*************添加高斯白噪声**************%设定信躁比(小信噪比)%设定信躁比(大信噪比)%设定信躁比(无信噪比)%计算对应的高斯白躁声的方差%生成高斯白躁声%生成含高斯白躁声的已调信号(信号通%*****************************************%****************FM解调*******************forilengtht%接受信号通过微分器处理diff_nsfmn=abs(hilbert(diff_nsfm));%hilbert变换，求绝对值得到瞬时幅度(包络检波)%*****************************************%**************时域到频域转换**************%*****对调制信号m(t)求傅里叶变换*****m=am*cos(2*pi*fm*t);抽样间隔信号M=fft(m,n);换M=M/fs;%以上程序是对调制后的信号u求傅里变%缩放，便于在频铺图上整体观察%时间向量对应的频率向量%************对已调信号u求傅里变换**********%******************************************%*****************************************%***************显示程序******************disp('按任意键可以看到原调制信号、载波信号和已调信号的曲线')%**************figure(1)******************title('调制信号的时域图');title('载波的时域图');title('已调信号的时域图');%绘制调制信号的时域图%绘制载波的时域图%绘制已调信号的时域图%******************************************disp('按任意键可以看到原调制信号和已调信号在频域内的图形')%************figure(2)*********************title('原调制信号的频谱图')title('已调信号的频谱图')%fftshift:将FFT中的DC分量移到频谱中心%******************************************disp('按任意键可以看到原调制信号、无噪声条件下已调信号和解调信号的曲线')%**************figure(3)******************title('调制信号的时域图');title('无噪声条件下已调信号的时域图');%绘制调制信号的时域图%绘制已调信号的时域图%接受信号通过微分器处理diff_nsfmn=abs(hilbert(diff_nsfm));%hilbert变换，求绝对值得到瞬时幅度(包络检波)subplot噪声条件下解调信号的时域图title('无噪声条件下解调信号的时域图');%*****************************************disp('按任意键可以看到原调制信号、小信噪比高斯白噪声条件下已调信号和解调信号已调信号的曲线')%**************figure(4)******************title('调制信号的时域图');%绘制调制信号的时域图%计算对应的小信噪比高斯白躁声的方差%生成高斯白躁声%生成含高斯白躁声的已调信号(信号通forilengtht)-1%接受信号通过微分器处理diff_nsfmn1=abs(hilbert(diff_nsfm1));%hilbert变换，求绝对值得到瞬时幅度(包络检波)plotlengthdiffnsfmdiffnsfm噪比高斯白噪声已调信号的时域图title('含小信噪比高斯白噪声已调信号的时域图');subplot绘制含小信噪比高斯白噪声解调信号的时域图title('含小信噪比高斯白噪声解调信号的时域图');%*****************************************disp('按任意键可以看到原调制信号、大信噪比高斯白噪声条件下已调信号和解调信号已调信号的曲线')%**************figure(5)******************title('调制信号的时域图');输)%绘制调制信号的时域图%计算对应的大信噪比高斯白躁声的方差%生成高斯白躁声%生成含高斯白躁声的已调信号(信号通%接受信号通过微分器处理%绘制含大信噪比高斯白噪声已调信号title('含大信噪比高斯白噪声已调信号的时域