(完整word版)基于MATLAB的FM调制实现.docVIP

FM MATLAB建模与仿真 一、 FM调制原理： FM 属于角度调制，角度调制与线性调制不同，已调信号频谱不再是原调制信号频谱的线性搬移， 而是频谱的非线性变换， 会产生与频谱搬移不同的新的频率成分，故又称为非线性调制。 在本实验中使用正弦信号作为基带信号进行调制的分析 . [1] 频率调制的一般表达式 为： FM 调制是相位偏移随 ｍ（ｔ）的积分呈线性变化。 FM 调制模型的建立 图 1 FM 调制模型 其中， m(t ) 为基带调制信号，设调制信号为 mt()  Acos(2 fmt) 设正弦载波为 c(t )  cos(2  fct ) 信号传输信道为高斯白噪声信道，其功率为  2 。 二  发送 m(t) SFM (t ) SFM (t ) 信息源 FM 调制器 信道 FM 解调器 m0 (t) 接收 图 2 总体模型 调制过程的分析 ： 在调制时，调制信号的频率去控制载波的频率的变化， 载波的瞬时频偏随调 制信号 m(t) 成正比例变化，即sFM FM MATLAB建模与仿真 d ( t ) K f m ( t ) dt 式中， Kf 为调频灵敏度（ rad ）。 ( s V ) 这时相位偏移为 (t ) K f m( )d 则可得到调频信号为 sFM(t) Acos ct Kf m( )d FM 调制 对 FM 调制信号的频谱分析 clear all ts=0.00125;  %信号抽样时间间隔 t=0:ts:10-ts;  %时间向量 am=10; fs=1/ts; df=fs/length(t);  % 抽样频率 %fft 的频率分辨率 msg=am*cos(2*pi*10*[0:0.01:0.99]); msg1=msg*ones(1,fs/10);  %扩展成取样信号形式 msg2=reshape(msg1.,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);  %对 msg进行傅利叶变换 N=(0:length(m)-1)*fs/length(m)-fs/2; subplot(3,1,2) plot(N,abs(m));  % 调制信号频谱图 title( 调制信号频谱 ) int_msg(1)=0; for ii=1:length(t)-1 int_msg(ii+1)=int_msg(ii)+msg2(ii)*ts;  %消息信号积分 end kf=50; fc=250;  %载波频率 Sfm=am*cos(2*pi*fc*t+2*pi*kf*int_msg); % 调频信号 FM MATLAB建模与仿真 Pfm=fft(Sfm)/fs;  % FM 信号频谱 subplot(3,1,3); plot(f,fftshift(abs(Pfm)))  % 画出已调信号频谱 title(FM  信号频谱  ) Pc=sum(abs(Sfm).^2)/length(Sfm)  %已调信号功率 Ps=sum(abs(msg2).^2)/length(msg2)  %消息信号功率 fm=50; betaf=kf*max(msg)/fm W=2*(betaf+1)*fm  % 调制指数 % 调制信号带宽 用 FFT函数进行傅利叶变换， 进行傅立叶变化便分别得到调制信号与调制之后的 FM信号的频谱图如下： x 10 5 消息信号频谱 2 0 -2 0 1 2 3 4 5 6 7 8 9 10 调制信号频谱 600 400 200 0 -500 -400 -300 -200 -100 0 100 200 300 400 500 FM信号频谱 10 5 0 -400 -300 -200 -100 0 100 200 300 400 500 -500 图2-5 通过频谱图的对照比较我们可以看出 FM调制并不是使原正弦信号的频谱在原来位置上通过移动得到调制波形， 调制后的波形与调制前的完全不同， 这证明 FM调制并不是线性的，而是非线性的。由图像可知 fft 变换出来的图形是左右 FM MATLAB建模与仿真 对称的。 改变采样频率对调制信号进行分析 采样频率分别为 100HZ，200HZ， 400HZ，800HZ时的时域信号波形 采样频率为100HZ时，时域信号的波形 1 度 0.5 幅 号0 信 域 时 -0.5 -1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0 t(s) 采样频率为200HZ时，时域信号