FM的MATLAB调制解调(word文档良心出品).docxVIP

PAGE PAGE # 通信原理仿真实验报告 实验名称： FM调制及解调 姓 名： 专 业： 年 级： 学 号： 201X年X月X日 PAGE PAGE # FM调制及解调 设输入信号为m(t)二 设输入信号为m(t)二cos2二t，载波中心频率为 fc OHz ，VCO的 压控振荡系数为5Hz 压控振荡系数为5Hz/V ,载波平均功率为1W试画出: 已调信号的时域波形； 已调信号的振幅谱； 用鉴频器解调该信号，并与输入信号比较 、程序代码 clear allt0=2; tz=O.OOO1; fs=1/tz; t=[-tO:tz:tO];kf=5;fc=10; kd=0.8;%时间向量精度%设定抽样频率%产生时间向量%设定压控振荡器系数%设定载波频率 clear all t0=2; tz=O.OOO1; fs=1/tz; t=[-tO:tz:tO]; kf=5; fc=10; kd=0.8; %时间向量精度 %设定抽样频率 %产生时间向量 %设定压控振荡器系数 %设定载波频率 %设定鉴频增益/鉴频器灵敏度 m_fun=cos(2*pi*t); in t_m(1)=0; for i=1:le ngth(t)-1 in t_m(i+1)=i nt_m(i)+m_fu n( i)*tz; %对m_fun积分 end x=sqrt(2)*cos(2*pi*fc*t+kf*i nt_m); % 调制信号 y=m_fu n. *kd*kf; %解调信号 z=-sqrt(2)*(2*pi*fc+kf*m_fu n).*si n(2*pi*fc*t+kf*i nt_m); Nf=4096*32; M=fft(m_fu n,Nf); f=[0:1:Nf-1]./Nf.*fs; X=fft(x,Nf); Y=fft(y,Nf); %对原始信号快速傅里叶变换 %对已调信号快速傅里叶变换 %对解调信号快速傅里叶变换 figure(1); plot(t,m_fu n(1:le ngth(t)),li newidth,2); title(原始信号的时域图形’)； xlabel(时间 /s); lege nd(m(t)) %生成原始信号的时域图形 figure(2); %生成原始信号的频域图形 h1=plot(f,abs(fftshift(M))/max(abs(M)),l in ewidth,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)); 二、实验结果与分析 已调信号的时域波形 已调信号的时域图形(Kf=5 )原始信号的时域图形时间/S1.50.5-0.5-1x(t)hHllLblmlli-1.5 已调信号的时域图形(Kf=5 ) 原始信号的时域图形 时间/S 1.5 0.5 -0.5 -1 x(t) hHllLblmlli -1.5 I 1 1 1 ■ 1 1 1 1 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 时间/s -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 时间/s 调频信号的一般表达式为：x(t) = Acos[ ?‘ct Kf nm )d ]，已知 m(t)二cos2二t，Wc=2二仁，仁=10Hz, Kf =5Hz/V,由于载波平