1.运行结果
DSB调制
DSB解调
原始信号与解调后的波形对比
LSSB调制
AM调制
程序源代码
t0 = 0.15; %信号持续时间 ts = 0.0001; %信号采样率 fc = 250; %载波信号 fs = 1/ts; %采样频率 df = 0.3; %频谱分辩率 t = [-t0/3:ts:t0]; %时间矢量
m=heaviside(t)-3*heaviside(t-t0/3)+2*heaviside(t-2*t0/3); %消息信号 [M,m,df1]=fftseq(m,ts,df); %消息信号的频谱 c = cos(2*pi*fc.*t); %载波信号 [C,c,df]=fftseq(c,ts,df); %载波信号频谱 u = m.*c; %调制信号
[U,u,df1]=fftseq(u,ts,df); %调制信号频谱
figure(1);
subplot(3,2,1) %消息信号时域显示 plot(t,m(1:length(t))); title('原始信号')
axis([-0.05 0.15 -3 3]);
f = [0:df:df*(length(m)-1)]-fs/2; %频谱矢量
subplot(3,2,2) %消息信号频域显示 plot(f,abs(fftshift(M))); title('原始信号的频谱') axis([-600 600 0 1500]);
subplot(3,2,3) %调制载波时域显示 plot(t,c(1:length(t))); title('载波信号')
axis([-0.05 0.15 -3 3]);
subplot(3,2,4) %调制载波频域显示 plot(f,abs(fftshift(C))); title('载波的频谱')
axis([-600 600 0 1500]);
subplot(3,2,5) %调制信号时域显示 plot(t,u(1:length(t))); title('调制信号')
axis([-0.05 0.15 -3 3]);
subplot(3,2,6) %调制信号频域显示 plot(f,abs(fftshift(U))); title('调制信号的频谱') axis([-600 600 0 1500]);
c1=u.*c; %与同频同相的调制载波相乘 [C1,c1,df1]=fftseq(c1,ts,df); %解调信号信号频谱
figure(2);
subplot(2,2,1) %解调信号时域显示 plot(t,c1(1:length(t))); title('解调信号')
axis([-0.05 0.15 -3 3]);
subplot(2,2,2) %解调信号频域显示 plot(f,abs(fftshift(C1))); title('解调信号的频谱') axis([-600 600 0 1500]);
H=heaviside(f+fc)-heaviside(f-fc); %低通滤波器（截止频率为fc） C2= C1.*abs(fftshift(H)); %通过该低通滤波器
c2=ifft(2*C2); %傅里叶反变换得出相干解调最终的波形
subplot(2,2,3) %滤波后波形显示
