基于matlab的fm调制实现。.doc

MATLAB建模与仿真题  目基于MATLAB的FM调制实现学  院   信息工程学院   专  业    通信工程    班  级  11通信工程 学  号  11031    姓 名         指导教师             一、FM调制原理:FM属于角度调制,角度调制与线性调制不同,已调信号频谱不再是原调制信号频谱的线性搬移,而是频谱的非线性变换,会产生与频谱搬移不同的新的频率成分,故又称为非线性调制。[1]为:                                                                     FM调制是相位偏移随m(t)的积分呈线性变化。FM调制模型的建立  图1 FM调制模型其中,为基带调制信号,设调制信号为设正弦载波为信号传输信道为高斯白噪声信道,其功率为。图2总体模型二 调制过程的分析:在调制时,调制信号的频率去控制载波的频率的变化,载波的瞬时频偏随调制信号成正比例变化,即式中,为调频灵敏度()。  clearallts=;             %信号抽样时间间隔t=0:ts:10-ts;           %时间向量am=10;fs=1/ts;              %抽样频率df=fs/length(t);          %fft的频率分辨率msg=am*cos(2*pi*10*[0::]);   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;           %消息信号积分forii=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);%调频信号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*(bet