基于MATLAB的FM调制实现

1、FMMATLAB1模与仿真一、FM调制原理：FM属于角度调制，角度调制与线性调制不同，已调信号频谱不再是原调制信号频谱的线性搬移，而是频谱的非线性变换，会产生与频谱搬移不同的新的频率成分，故又称为非线性调制。在本实验中使用正弦信号作为基带信号进行调制的分析.频率调制的一般表达式1为：FM调制是相位偏移随m(t)的积分呈线性变化FM调制模型的建立切FMSfm")图1FM调制模型其中，m(t)为基带调制信号，设调制信号为m(t)=Acos(2fmt)设正弦载波为c(t)=cos(27fct)信号传输信道为高斯白噪声信道，其功率为二-2信息源发送m(t)Sfm(t),I信道调制过程的分析：

2、在调制时，调制信制信号m(t)成正比例变化，即Sfm(t)FM解调器(0图2总体JFMm0(t).接收FMMATLAB1模与仿真d(t)dt-Kfm(t)式中，Kf为调频灵敏度（ra%s.V）这时相位偏移为(t)=Kfm()d则可得到调频信号为SFM(t)=AcosctKfm()dFM调制1.对FM调制信号的频谱分析clearallts=0.00125;t=0:ts:10-ts;am=10;fs=1/ts;df=fs/length(t);msg=am*cos(2*pi*10*0:0.01:0.99);msg1=msg'*ones(1,fs/10);msg2=reshape(msg1.&

3、#39;,1,length(t);Pm=fft(msg2);f=-fs/2:df:fs/2-df;subplot(3,1,1)plot(t,fft(abs(Pm)title('消息信号频谱')m=fft(msg,1024);N=(0:length(m)-1)*fs/length(m)-fs/2;subplot(3,1,2)plot(N,abs(m);title('调制信号频谱')%信号抽样时间间隔%时间向量%抽样频率%fft的频率分辨率%扩展成取样信号形式%求消息信号的频谱%又msg进行傅利叶变换%调制信号频谱图int_msg(1)=0;%消息信号积分forii

4、=1:length(t)-1int_msg(ii+1)=int_msg(ii)+msg2(ii)*ts;endkf=50;fc=250;%载波频率Sfm=am*cos(2*pi*fc*t+2*pi*kf*int_msg);%调频信号FMMATLAB1模与仿真Pfm=fft(Sfm)/fs;subplot(3,1,3);plot(f,fftshift(abs(Pfm)title('FM信号频谱')%FM信号频谱%画出已调信号频谱Pc=sum(abs(Sfm).A2)/length(Sfm)Ps=sum(abs(msg2).A2)/length(msg2)%已调信号功率%消息信号功

5、率fm=50;betaf=kf*max(msg)/fmW=2*(betaf+1)*fm%调制指数%调制信号带宽用FFT®数进行傅利叶变换,的FMJ号的频谱图如下：进行傅立叶变化便分别得到调制信号与调制之后5消息信号频谱0¥r-26106004002000x10-500-400-300-200-1000100200300400500调制信号频谱FM信号频谱图2-5通过频谱图的对照比较我们可以看出FMS制并不是使原正弦信号的频谱在原来位置上通过移动得到调制波形，调制后的波形与调制前的完全不同，这证明FMS制并不是线性的，而是非线性的。由图像可知fft变换出来的图形是左右FMMA

6、TLAB1模与仿真00.10.20.30.40.50.60.70.80.91t(s)度幅号信域时度幅号信域时改变采样频率对调制信号进行分析采样频率分别为100HZ；200HZ400HZ800Hz时的时域信号波形采样频率为100Hz寸，时域信号的波形10.50-0.5-1采样频率为200Hz时，时域信号的波形10.50-0.5-100.10.20.30.40.50.60.70.80.91t(s)采样频率为400HZM,时域信号的波形15051-oO-度幅号信域时,4O.6OFMMATLAB1模与仿真采样频率分别为100Hz200Hz400HZ；800Hz时的频域信号波形100Hz采样频率FM信号

7、的频域频谱50r(TI(T-30-20-1001020304050频率分量1O。O432200Hz采样频率FM信号的频域频谱-60-40-20020406080100频率分量1。O。O5432-8谱率功号信制调P过进谱率功号信制调MF过进FMMATLAB1模与仿真50400Hz采样频率FM信号的频域频谱-200-150-100-50050100150200频率分量oooO4321谱率功号信制调笆过进谱率功号信制调MF过进800Hz采样频率FM信号的频域频谱50403020100-400-300-200-1000100200300400频率分量载波频率不变，改变采样频率后，如果将采样频率取的太小

8、，到小于两倍的载波频率，MATLA颤件将不能输出波形，说明采样频率只有在大于两倍载波在频域内，信号的能量主要集中在采样频率为|100|HZ的范围内，可以看出采用等位置采样方法进行周期采样，频率随机波动的大小不会产生谱分析误差，当频率随机波动水平较大时，不仅会导致功率谱峰值的显著降低，。3 .加入噪声，绘制出加入噪声后的FM信号的时频图，并对其进行分析噪声对信号的影响。设其载波频率fc=150;采样频率fs=800;调制信号x=sin(2*pi*30*t);FM信号y=modulate(x,fc,fs,'FM');利用y1=y+awgn(y,1,0);加入噪声。得到其时域波形图如

9、下：FMMATLAB1模与仿真FM信号加入噪声的时域波形图通过b=fft(y1,1024)；f=(0:length(b)-1)*fs/length(b)傅利叶变换，得到加入噪声后的FM信号频谱为：-fs对y进行-300-N口口-1OO口11ooIFreciugmu白11nJ00n005&31三fids警dFM信号加入噪声的频谱图通过对函数y1=y+awgn(y,10,0)改变信噪比后，信噪比由1变为10,得到信噪比为10的时域图为:信噪比改为10的时域图频谱图为:FMMATLAB1模与仿真00000e-s032三ls-s-c-Ann-s>nri信噪比改为10的频域图通过对函数y1

10、=y+awgn(y,20,0)改变信噪比后，信噪比由10变为20,得到信噪比为20的时域图为:信噪比改为20的时域图频域图为:-kin-rm-innninn?nnirnWLJ0suaOIn)：-inn-annqno信噪比改为20的频域图加入噪声后时域波形与原来的时域波形相比，波形明显失真，波形不仅不如FMMATLAB1模与仿真原本波形般规则，而且曲线之间还出现了为链接在一起的断裂，但随着信噪比的增大，与原有的波形的相似度也增大了，说信噪比越大，噪声对信号的影响也变小了。从加入噪声的图形与未加入噪声的对比中我们还可以看出噪声对时域图的变化明显比频域图的变化更为突出，白噪声在整个频谱内每个频点的能

11、量为常数，且基本恒定，所以他对于时域的影响更大。通过对噪声的理解我们可以知道对于调频系统来说，增加传输带宽就可以改善抗噪声性能。调频方式的这种以带宽取信噪比的特性是十分有益的4 .分析信噪声比大于多少的情况下，噪声会完全淹没在信号中，在接收端,解调出来的信号和没有加入噪声时传输到接收端的信号几乎完全一致设置信噪比在10db,15db,25db,30db,40db,50db,80db原始消息信号200-20-1000-800-600-400-20002004006008001000无噪声时的解调信号200-20111*hiLctrrr-1000-800-600-400-2000200400600

12、8001000信噪比为10dB时的解调信号原始消息信号-1000-800-600-400-20002004006008001000无噪声时的解调信号200-20卜NWfF1I-A,liJrrT2L.rr卜200-20-1000-800-600-400-20002004006008001000-1000-800-600-400-20002004006008001000信噪比为15dB时的解调信号FMMATLAB1模与仿真原始消息信号012345678910无噪声时的解调信号012345678910信噪比为30dB时的解调信号-20012345678910无噪声时的解调信号20-201J-1000

13、-800-600-400-20002004006008001000信噪比为40dB时的解调信号200-20-1000-800-600-400-20002004006008001000FMMATLAB1模与仿真20原始消息信号0无噪声时的解调信号2004006008001000-20-1000-800-600-400-20020(1(1E"IIIIII1J1IbI0- 20J11c1-1000-800-600-400-20002004006008001000信噪比为50dB时的解调信号201111匚(；'|''- 2011L1-1000-800-600-400-2000200