DSB（双边带）调制

clear all;
%% DSB调制
% DSB（双边带）只需将调制信号m(t)与载波信号cos(wt)直接相乘即可
dt=1/800;
T = 200;  % 采样的总时间。频谱分辨率（df=1/T）。
t = 0 : dt : T-dt;
fm = 2; % 调制信号的频率，单位kHz
fc = 20; % 载波信号的频率，单位kHz
m = cos(2*pi*fm*t); % 调制信号
s = m.*cos(2*pi*fc*t); %DSB已调信号
[f,sf] = T2F(t,s);
figure(1)
plot(t,s);
axis([0,1,-1,1]);
figure(2)
plot(f,abs(sf));
axis([-30,30,0,55]);

%% 函数  计算信号的傅里叶变换
function[f, sf] = T2F(t,st) % t为时域采样点；st为采样的时域信号
dt = t(2) -t(1);
% T = t(end);
T = t(end)-t(1)+dt;
df = 1/T;
N = length(st);
f = -N/2*df : df: N/2*df - df;
sf = fft(st);
sf = T/N * fftshift(sf);
end

%% 计算信号频谱 sf 的傅里叶反变换
function [ t, st ] = F2T( f, sf )
df = f(2)-f(1);
Fmx = f(end)-f(1)+df;
dt = 1/Fmx;
N = length(sf);
T = dt * N;
t = 0 : dt : T-dt; % 或 t = -T/2 : dt : T/2-dt;
sff = fftshift(sf);
st = Fmx * ifft(sff);
end