am.m

这是关于模拟调制的MATLAB源代码

%% 准带限信号的AM(双边带)调制
t0=0.15;
ts=0.001; 
fc=250;   
fs=1/ts; 
a=0.85;
df=0.2;   
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));  % 归一化后的被调信号
u=(1+a*m_n).*c;
[M,m,df1]=fftseq(m,ts,df); 
M=M/fs;             
[U,u,df1]=fftseq(u,ts,df);
U=U/fs;
f=[0:df1:df1*(length(m)-1)]-fs/2; 
subplot(4,2,1);
plot(f,abs(fftshift(M)))
title('被调信号m(t)的幅度频谱')
gtext('|M(f)|')
subplot(4,2,2);
plot(f,abs(fftshift(U)))
title('调制信号u(t)的幅度频谱')
gtext('|U(f)|')
subplot(4,2,3);
plot(t,m(1:length(t)));
title('被调信号m(t)')
gtext('m(t)')
subplot(4,2,4);
plot(t,u(1:length(t)));
title('调制信号u(t)')
gtext('u(t)')