matlab中fftseq怎么用,[转载]AM调制解调的matlab仿真

用matlab产生一个频率为1Hz，振幅为1的余弦信源，设载波频率为10Hz,A=2.想干解调器输入信噪比为15dB

主程序：

%AM调制

fc=10;

f0=1;

A=2;

t0=5;%如果频率为1，那至少要信号时长5个周期

snr=15;

dt=0.003;

fs=1/dt;

B=2;

df=0.0003;

t=0:dt:t0;

Lt=length(t);

snr_lin=10^(snr/10);

figure(1);

subplot(321)%调制信号波形

m1=cos(2*pi*f0*t);

plot(t,m1(1:length(t)))

R=2*max(m1);

axis([0 t0 -R/2

R/2])

ylabel('调制信号');

subplot(322)%调制信号频谱

[M1,m1,df1,f]=T2F(m1,dt,df,fs);

plot(f,fftshift(abs(M1)));

axis([-2 2 0

3]);

xlabel('f');

ylabel('调制信号频谱');

subplot(323)%载波波形

m2=cos(2*pi*fc*t);

plot(t,m2(1:length(t)));

ylabel('载波信号波形');

axis([0 0.4 -R/2

R/2])

subplot(324);%载波频谱

[M2,m2,df1,f]=T2F(m2,dt,df,fs);

plot(f,fftshift(abs(M2)));

xlabel('f');

ylabel('载波信号频谱');

axis([-20 20 0

3])

subplot(325)%已调波波形

m=(m1(1:length(t))+A).*m2(1:length(t));

plot(t,m(1:length(t)));

R=2*max(m)+A;

axis([0 3 -R/2

R/2])

ylabel('已调信号波形');

subplot(326)

[M,m,df1,f]=T2F(m,dt,df,fs);