FM的MATLAB调制解调

通信原理仿真实验报告实验名称：FM调制及解调姓 名：专 业：年 级：学 号：201X年 X 月X日FM调制及解调设输入信号为 ，载波中心频率为 ，VCO的压控振荡系数为 ，载波平均功率为1W。试画出：l 已调信号的时域波形；l 已调信号的振幅谱；l 用鉴频器解调该信号，并与输入信号比较。一、 程序代码clear allt0=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;endx=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);二、实验结果与分析（1）已调信号的时域波形调频信号的一般表达式为：，已知，,由于载波平均功率为1W，故。由于，调制波形的频率变化不明显，对比第二、三图可知，越大时（第三张图），调制现象更加明显。（2）已调信号的振幅谱原信号频域表达式，由于Matlab程序的加窗截断效应产生了频谱泄露现象，导致冲激函数变为脉冲，可以观察出在处有两个冲激。已调信号载频为10Hz，从图像可以看出在处有两个冲激。（3）解调输出的时域、频域图像1、解调时域波形第一幅图是鉴频电路微分电路的输出波形，该电路使得载波信号的振幅和频率都随原信号m(t)发生变化。第二幅图是鉴频电路的包络检波器输出，此时输出信号