振幅调制和解调的原理及MATLAB编程实现(6)

jjj

n=0:N-1; t=n/Fs

ft=Am*cos(2*pi*fm*t+Qm);%调制信号 subplot(421) plot(t,ft) title('调制信号') y1=fft(ft,N)%傅利叶变换 mag1=abs(y1)%取复变函数幅值

f1=(0:length(y1)-1)*Fs/length(y1);%确定频率范围 subplot(423)%子窗口位置

plot(f1(1:N/2/10),mag1(1:N/2/10))%图形输出 title('调制信号频谱')%标题

SAM=(Ac+ft).*cos(2*pi*fc*t+Qc);%已调波信号AM subplot(422) plot(t,SAM)

title('已调波信号AM') y2=fft(SAM,N); mag2=abs(y2);

f2=(0:length(y2)-1)*Fs/length(y2); subplot(424)

plot(f2(1:N/2/5),mag2(1:N/2/5))

title('已调波信号AM频谱') %以上为调制部分 spt=SAM.*cos(2*pi*fcP*t+QcP)%相干解调 subplot(425)

plot(t,spt) title('相干解调波形') y3=fft(spt,N); mag3=abs(y3);

f3=(0:length(y3)-1)*Fs/length(y3); subplot(426)

plot(f3(1:N/2/5),mag3(1:N/2/5))

搜索“diyifanwen.net”或“第一范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，第一范文网，提供最新资格考试认证振幅调制和解调的原理及MATLAB编程实现(6)全文阅读和word下载服务。