jmam.m

通信中一些基本实验的程序AM、FM、ASK、FSK、PSK等等

%显示模拟调制的波形及解调AM，文件mam.m
%信源
close all;
clear all;
dt=0.001;                    %采样时间间隔
fm=1;                        %信源最高频率
fc=10;                       %载波中心频率
T=5;                         %信号时长
t=0:dt:T-dt;
mt=sqrt(2)*cos(2*pi*fm*t);   %信源
N0=5;                    %白噪声单边功率谱密度

%AM modulation
A=2;
s_am=(A+mt).*cos(2*pi*fc*t);
figure(1)
subplot(411)
plot(t,s_am);hold on;         %画出AM信号波形
plot(t,A+mt,'r--');           %标示AM的包络
title('AM调制信号及其包络');
xlabel('t');

B=2*fm;                      %带通滤波器带宽
noise=noise_nb(fc,B,N0,t);  % 窄带高斯白噪声
J=(sqrt(2)+noise).*cos(2*pi*t+pi);
s_am=s_am+J;


%AM demodulation
rt=s_am.*cos(2*pi*fc*t);      %相干解调
rt=rt-mean(rt);
[f,rf]=T2F(t,rt);
[t,rt]=lpf(f,rf,2*fm);        %低通滤波
subplot(412)
plot(t,rt);hold on;
plot(t,mt/2,'r--');
title('相干解调后的信号波形与输入信号的比较');
xlabel('t')

[f,mf]=T2F(t,mt);
[f,rf]=T2F(t,rt);            
prf=(abs(rf).^2)/T;           
pmf=(abs(mf).^2)/T;
subplot(413)
plot(f,pmf);                  %画出原信号的频谱
axis([-2*fc 2*fc 0 max(pmf)]);
title('原信号的功率谱');
subplot(414)
plot(f,prf,'r--');           %画出干扰后的信号频谱
axis([-2*fc 2*fc 0 max(prf)]);
title('噪声干扰后的功率谱');
xlabel('f');