%顯示模擬調制的波形及解調方法DSB
%信源
close all;
clear all;
dt = 0.001;          %時間采樣間隔         
fm=1;                %信源最高頻率
fc=10;               %載波中心頻率
T=5;                 %信號時長
t = 0:dt:T; 
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);
%s_dsb=s_dsb+noise;

figure(2)
subplot(211)
plot(t,mt);          %標示mt的波形
title('Input signal');
subplot(212)
plot(t,s_dsb);hold on;     %畫出DSB信號波形
plot(t,abs(mt),'r--');          %標示mt的波形
title('DSB modulated signal');

figure(1)
subplot(311)
plot(t,s_dsb);hold on;     %畫出DSB信號波形
plot(t,mt,'r--');          %標示mt的波形
%title('DSB調制信號');
title('DSB modulated signal');
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('相干解調后的信號波形與輸入信號的比較');
xlabel('t')
subplot(313)
[f,sf]=T2F(t,s_dsb);
psf = (abs(sf).^2)/T;
plot(f,psf);
axis([-2*fc 2*fc 0 max(psf)]);
title('DSB信號功率譜');
xlabel('f');