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

plot(t,c2(1:length(t))); axis([-0.05 0.15 -3 3]); title('滤波后的波形');

subplot(2,2,4) %滤波波形频域显示 plot(f,abs(fftshift(C2))); axis([-600 600 0 1500]); title('滤波后的频谱'); figure(3);

subplot(2,2,1); plot(t,m(1:length(t)));; title('原始信号')

axis([-0.05 0.15 -3 3]);

subplot(2,2,2) plot(f,abs(fftshift(M))); title('原始信号的频谱') axis([-600 600 0 1500]);

subplot(2,2,3) plot(t,c2(1:length(t))); axis([-0.05 0.15 -3 3]); title('滤波后的波形')

subplot(2,2,4) plot(f,abs(fftshift(C2))); axis([-600 600 0 1500]); title('滤波后的频谱')

figure(4); subplot(4,2,1) plot(t,m(1:length(t))); title('原始信号')

axis([-0.05 0.15 -3 3]);

subplot(4,2,2) plot(f,abs(fftshift(M))); title('原始信号的频谱') axis([-600 600 0 1500]);

subplot(4,2,3) plot(t,c(1:length(t))); title('载波信号')

axis([-0.05 0.15 -3 3]);

%消息信号时域显示%消息信号频域显示 %滤波后波形显示 %滤波波形频域显示

%消息信号时域显示 %消息信号频域显示 %调制载波时域显示

subplot(4,2,4) %调制载波频域显示 plot(f,abs(fftshift(C))); title('载波的频谱')

axis([-600 600 0 1500]);

subplot(4,2,5) %调制信号时域显示 plot(t,u(1:length(t))); title('双边带调制信号的波形') axis([-0.05 0.15 -3 3]);

subplot(4,2,6) %调制信号频域显示 plot(f,abs(fftshift(U))); title('双边带调制信号的频谱') axis([-600 600 0 1500]);

Hl=heaviside(f+fc)-heaviside(f-fc); %低通滤波器（截止频率为fc） C1= U.*abs(fftshift(Hl)); %通过该低通滤波器

c1=ifft(C1) ; %傅里叶反变换得SSB下边带调制的最终波形

subplot(4,2,7) ; %滤波后波形显示 plot(t,c1(1:length(t))); axis([-0.05 0.15 -3 3]); title('下边带调制信号')

subplot(4,2,8) %滤波波形频域显示 plot(f,abs(fftshift(C1))); axis([-600 600 0 1500]); title('下边带调制的频谱')

figure(5);

subplot(3,2,1) %消息信号时域显示 plot(t,m(1:length(t))); title('原始信号')

axis([-0.05 0.15 -3 3]);

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]);

v = 0.85*m.*c+c; %调制

[U,u,df1]=fftseq(u,ts,df); %调制信号频谱

subplot(3,2,5) %调制信号时域显示 plot(t,v(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]);

附上自定义函数fftseq()

function [M,m,df]=fftseq(m,ts,df) fs=1/ts; if nargin == 2 n1=0; else

n1=fs/df; end

n2=length(m);

n=2^(max(nextpow2(n1),nextpow2(n2))); M=fft(m,n);

m=[m,zeros(1,n-n2)]; df=fs/n; end

2. 维弗法：滤波法和相移法相结合，仅对载波移相，不需要宽带移相。 2至8点处的信号光谱变化如下

3. 发送信号是s?t??Acm?t?cos?2?fct????Ac?t?m?t?sin?2?fct???

^??j?2?fct???? ?Re?Ac?m?t??jm?t??e? ?????^???j?2?fct? ?Re?Ac?m?t??jm?t??ee? ?????^???j?其复包络是sL?t??Acm?t??jm?t?e

????^