第五章参考答案（通信原理，杨建义老师，浙江大学）

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) %滤波后波形显示