%顯示模擬調(diào)制的波形及解調(diào)方法DSB，文件mdsb.m
%信源
close all;
clear all;
dt=0.001;                      %采樣時間間隔
fm=1;                          %信源的最高頻率
fc=10;                         %載波中心頻率
T=5;                           %信號時長
t=0:dt:T-dt;                      
mt=sqrt(2)*cos(2*pi*fm*t);     %信源
%N0=0.01;                      %白噪單邊功率譜密度
%DSB modulation
s_dsb=mt.*cos(2*pi*fc*t);
B=2*fm;
%noise=noise_nb(fc,B,N0,t);
%zs_dsb=s_dsb+noise;
figure(1);
subplot(311)
plot(t,s_dsb);hold on;         %畫出DSB信號的波形
plot(t,mt,'r--');              %標(biāo)出mt的波形
title('DSB調(diào)制信號');
xlabel('t');
%DSB demodulation
rt=s_dsb.*cos(2*pi*fc*t);
rt=rt-mean(rt);                                                            
[f,rf]=T2F(t,rt); 
[t,rt]=lpf(f,rf,2*fm);                                                                                                                                 
subplot(312)
plot(t,rt);hold on;
plot(t,mt/2,'r--');
title('相干解調(diào)后的信號波形與輸入信號的比較');
xlabel('t')
subplot(313)
[f,sf]=T2F(t,s_dsb);          %求調(diào)制信號的頻譜
psf=(abs(sf).^2)/T;           %求調(diào)制信號的功率譜密度
plot(f,psf);
axis([-2*fc 2*fc 0 max(psf)]);
title('DSB的信號功率譜');
xlabel('f');