基于matlab的AM—DSB调制系统仿真
title('滤波前AM解调信号频谱'); axis([187960,188040,0,200000]);
滤波前AM解调信号波形6420-1x 10-0.85-0.6-0.4-0.200.20.40.60.81滤波前AM解调信号频谱21.510.501.87961.87971.87981.87991.881.88011.88021.88031.8804x 105
如图6.7 AM解调信号
%=======================AM解调信号FIR滤波=================
t=-1:0.00001:1;
A0=10; %载波信号振幅 A1=5; %调制信号振幅
A2=3; %已调信号振幅 ************************************* f=6000; %载波信号频率 w0=f*pi;
m=0.15; %调制度 **************************************
Uc=A0.*cos(w0*t); %载波信号 subplot(5,2,1); plot(t,Uc);
title('载波信号'); axis([0,0.01,-15,15]);
T1=fft(Uc); %傅里叶变换 subplot(5,2,2); plot(abs(T1));
title('载波信号频谱');
axis([5800,6200,0,1000000]);
10
基于matlab的AM—DSB调制系统仿真
mes=A1*cos(0.001*w0*t); %调制信号 subplot(5,2,3); plot(t,mes);
title('调制信号'); T2=fft(mes); subplot(5,2,4);
plot(abs(T2)); title('调制信号频谱');
axis([198000,202000,0,2000000]);
Uam=A2*(1+m*mes).*cos((w0).*t); %AM 已调信号 ***************** subplot(5,2,5); plot(t,Uam); title('已调信号'); T3=fft(Uam); subplot(5,2,6); plot(abs(T3));
title('已调信号频谱');
axis([5950,6050,0,500000]);
Dam=Uam.*cos(w0*t); %对AM已调信号进行解调 subplot(5,2,7);
plot(t,Dam);
title('滤波前的AM解调信号波形');
T4=fft(Dam); %求AM信号的频谱 subplot(5,2,8); plot(abs(T4));
title('滤波前的AM解调信号频谱'); axis([187960,188040,0,200000]);
z21=fftfilt(b21,Dam); %FIR低通滤波 subplot(5,2,9);
plot(t,z21,'r');
title('滤波后的AM解调信号波形');
T5=fft(z21); %求AM信号的频谱 subplot(5,2,10); plot(abs(T5),'r');
title('滤波后的AM解调信号频谱'); axis([198000,202000,0,200000]);
11
基于matlab的AM—DSB调制系统仿真
载波信号10100-10050-5-11040-10-11050-12-0.8-0.6-0.4-0.200.20.40.60.810.0010.0020.0030.0040.0050.0060.0070.0080.0090.01调制信号215x 105载波信号频谱05800x 10658505900595060006050610061506200调制信号频谱01.98-0.8-0.6-0.4-0.20已调信号0.20.40.60.81x 1051.9851.991.99522.0052.012.0152.02x 105已调信号频谱0595059605970598059906000601060206030604060502101.87961.87971.87981.8799x 105滤波前的AM解调信号波形x 105滤波前的AM解调信号频谱-0.8-0.6-0.4-0.200.20.40.60.811.881.88011.88021.88031.8804x 105滤波后的AM解调信号波形50-5-121滤波后的AM解调信号频谱-0.8-0.6-0.4-0.200.20.40.60.8101.981.9851.991.99522.0052.012.0152.02x 105
如图6.9 AM解调信号FIR低通滤波
%=======================AM解调信号FIR滤波======加噪==========
t=-1:0.00001:1;
A0=10; %载波信号振幅 A1=5; %调制信号振幅
A2=3; %已调信号振幅 ************************************* f=6000; %载波信号频率 w0=f*pi;
m=0.15; %调制度 **************************************
12
基于matlab的AM—DSB调制系统仿真
Uc=A0.*cos(w0*t); %载波信号 subplot(5,2,1); plot(t,Uc);
title('载波信号'); axis([0,0.01,-15,15]);
T1=fft(Uc); %傅里叶变换 subplot(5,2,2); plot(abs(T1));
title('载波信号频谱');
axis([5800,6200,0,1000000]);
mes=A1*cos(0.001*w0*t); %调制信号 subplot(5,2,3); plot(t,mes);
title('调制信号'); T2=fft(mes); subplot(5,2,4);
plot(abs(T2)); title('调制信号频谱');
axis([198000,202000,0,2000000]);
Uam1=A2*(1+m*mes).*cos((w0).*t); %AM 已调信号 ***************** subplot(5,2,5); plot(t,Uam1); title('已调信号'); T3=fft(Uam1); subplot(5,2,6); plot(abs(T3));
title('已调信号频谱');
axis([5950,6050,0,500000]);
sn1=20; %信噪比
db1=A1^2/(2*(10^(sn1/10))); %计算对应噪声方差 n1=sqrt(db1)*randn(size(t)); %生成高斯白噪声 Uam=n1+Uam1;
Dam=Uam.*cos(w0*t); %对AM已调信号进行解调 subplot(5,2,7);
plot(t,Dam);
title('滤波前的AM解调信号波形');
T4=fft(Dam); %求AM信号的频谱 subplot(5,2,8);
13
基于matlab的AM—DSB调制系统仿真
plot(abs(T4));
title('滤波前的AM解调信号频谱'); axis([187960,188040,0,200000]);
z21=fftfilt(b21,Dam); %FIR低通滤波 subplot(5,2,9);
plot(t,z21,'r');
title('滤波后的AM解调信号波形');
T5=fft(z21); %求AM信号的频谱 subplot(5,2,10); plot(abs(T5),'r');
title('滤波后的AM解调信号频谱'); axis([198000,202000,0,200000]);
载波信号00.005调制信号50-5-1-0.50已调信号100-10-1-0.500.510.510.0110505800x 105载波信号频谱59006000610062002101.98x 106调制信号频谱1.9922.012.02x 10已调信号频谱600060505505950x 105滤波前的AM解调信号波形100-10-1-0.500.512101.87961.87981.881.88021.8804x 10滤波后的AM解调信号频谱x 10555滤波前的AM解调信号频谱x 10滤波后的AM解调信号波形50-5-1-0.500.512101.981.9922.012.02x 105
14