m06_amod.m

MATLAB仿真技术与实例应用教程

t0=0.15;    %信号保持时间
ts=0.001;   %采样时间间隔
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;        %频率矢量

clf
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))); %已调制信号频谱
title('已调制信号频谱');