第五章参考答案(通信原理,杨建义老师,浙江大学).docVIP

运行结果 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) %滤波后波形显示 plot(t,