FM的MATLAB调制解调.doc

通信原理仿真实验报告 实验名称： FM调制及解调 姓 名： 专 业： 年 级： 学 号： 201X年 X 月X日  FM调制及解调 设输入信号为 ，载波中心频率为 ，VCO的压控振荡系数为 ，载波平均功率为1W。试画出： 已调信号的时域波形； 已调信号的振幅谱； 用鉴频器解调该信号，并与输入信号比较。 程序代码 clear all t0=2; tz=0.0001; %时间向量精度 fs=1/tz; %设定抽样频率 t=[-t0:tz:t0]; %产生时间向量 kf=5; %设定压控振荡器系数 fc=10; %设定载波频率 kd=0.8; %设定鉴频增益/鉴频器灵敏度 m_fun=cos(2*pi*t); int_m(1)=0; %对m_fun积分 for i=1:length(t)-1 int_m(i+1)=int_m(i)+m_fun(i)*tz; end x=sqrt(2)*cos(2*pi*fc*t+kf*int_m); %调制信号 y=m_fun.*kd*kf; %解调信号 z=-sqrt(2)*(2*pi*fc+kf*m_fun).*sin(2*pi*fc*t+kf*int_m); Nf=4096*32; M=fft(m_fun,Nf); %对原始信号快速傅里叶变换 f=[0:1:Nf-1]./Nf.*fs; X=fft(x,Nf); %对已调信号快速傅里叶变换 Y=fft(y,Nf); %对解调信号快速傅里叶变换 figure(1); %生成原始信号的时域图形 plot(t,m_fun(1:length(t)),linewidth,2); title(原始信号的时域图形); xlabel(时间/s); legend(m(t)) figure(2); %生成原始信号的频域图形 h1=plot(f,abs(fftshift(M))/max(abs(M)),linewidth,1); title(原始信号的频域图形); xlabel(频率/Hz); legend(M(f)); figure(3); %生成已调信号的时域图形 plot(t,x(1:length(t)),linewidth,2); title(已调信号的时域图形); xlabel(时间/s); legend(x(t)); figure(4); %生成已调信号的频域图形 plot(f,abs(fftshift(X))/max(abs(X)),linewidth,1); title(已调信号的频域图形); xlabel(频率/Hz); legend(X(f)); figure(5); %鉴频微分电路输出 plot(t,z(1:length(t)),linewidth,2); title(鉴频微分电路输出的时域图形); xlabel(时间/s); legend(z(t)); figure(6); %生成解调信号的时域图形 plot(t,y(1:length(t)),linewidth,2); title(解调信号的时域图形); xlabel(时间/s); legend(y(t)); figure(7); %生成解调信号的频域图形 plot(f,abs(fftshift(Y))/max(abs(Y)),linewidth,1); title(解调信号的频域图形); xlabel(频率/Hz); legend(Y(f)); 二、实验结果与分析 信号的时域波形 调频信号的一般表达式为：，已知，，,,由于载波平均为故由于的频率变化不明显，对比第、三图可知，时张图，现象更加明显 （2）已调信号的振幅谱 原信号频域表达式由程序的加窗截断了频谱泄露，