am.m

实现了AM的双边带调制解调功能！！！其中函数为自编

%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%AM调制%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%
%调制部分
t0=0.15;                 %信号持续时间
ts=0.0005;               %采样时间间隔
Fc=250;                  %载波频率
Fs=1/ts;                 %采用频率
df=0.3;                  %频率分辨率
a=0.8;                   %调制系数
t=[0:ts:t0];             %时间矢量
m=[ones(1,t0/(3*ts)),-2*ones(1,t0/(3*ts)),zeros(1,t0/(3*ts)+1)];
                         %定义信号序列
c=cos(2*pi*Fc.*t);       %载波信号
m_n=m/max(abs(m));       %调制信号(归一化)
[M,m,dfl]=fft_seq(m,ts,df);    %傅立叶变换
M=M/Fs;                        
u=(1+a*m_n).*c;
[U,u,dfl]=fft_seq(u,ts,df);
U=U/Fs;
f=[0:dfl:dfl*(length(m)-1)]-Fs/2;            %频率矢量
figure(1);
subplot(2,2,1);plot(t,m(1:length(t)));
title('未调制信号');
subplot(2,2,2);plot(t,u(1:length(t)));
title('已调制信号');
subplot(2,2,3);plot(f,abs(fftshift(M)));
title('未调制信号频谱');
subplot(2,2,4);plot(f,abs(fftshift(U)));
axis([0 1000 0 0.08]);
title('已调制信号频谱');


%解调部分
env=abs(Hilbert(u));                        %信号持续时间
dem1=2*(env-1)/a;                           %去掉DC分量并重新缩放
signal_power=(norm(u)^2)/length(u);         %调制信号功率
noise_std=sqrt(signal_power)/100;           %噪声标准偏差 
noise=noise_std*randn(1,length(u));         %产生噪声
r=u+noise;                                  %给调制信号加入噪声
env_r=abs(Hilbert(r));                      %找出包络
dem2=2*(env_r-1)/a;                         %去掉DC分量并重新缩放
figure(2);
subplot(3,1,1);plot(t,m(1:length(t)));
title('未调制信号');
subplot(3,1,2);plot(t,dem1(1:length(t)));
title('无噪声时的调制信号');
subplot(3,1,3);plot(t,dem2(1:length(t)));
title('有噪声时的调制信号');