FM的MATLAB调制解调

1 通信原理仿真实验报告 实验名称： FM 调制及解调 姓 名： 专 业： 年 级： 学 号： 201X 年 X 月 X 日 2 FM 调制及解调 设输入信号为( )cos 2m tt ，载波中心频率为 10czfH，VCO 的压控振荡系数为 5/zHV ，载波平均功率为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); 3 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'); ...