%% DSB(雙邊帶)調(diào)制
figure;
t0=0.15;
ts=0.001; % 抽樣間隔
fc=250;   % 載波頻率
fs=1/ts;  % 抽樣頻率
df=0.3;   % 所需的頻率分辨率
t=[0:ts:t0];
% 被調(diào)信號(hào)m(t)，載波信號(hào)c(t)，調(diào)制信號(hào)u(t)
m=[ones(1,t0/(3*ts)),-2*ones(1,t0/(3*ts)),zeros(1,t0/(3*ts)+1)];
c=cos(2*pi*fc.*t);
u=m.*c;
[M,m,df1]=fftseq(m,ts,df); % 傅立葉變換
M=M/fs;              % 縮放
[U,u,df1]=fftseq(u,ts,df);
U=U/fs;
% 繪圖
f=[0:df1:df1*(length(m)-1)]-fs/2; % 頻率向量
subplot(2,1,1);
plot(f,abs(fftshift(M)))
title('被調(diào)信號(hào)m(t)的幅度頻譜')
gtext('|M(f)|')
subplot(2,1,2);
plot(f,abs(fftshift(U)))
title('調(diào)制信號(hào)u(t)的幅度頻譜')
gtext('|U(f)|')