幅度调制(线性调制)是由调制信号去控制高频载波的幅度,使之调制信号的频谱线性变化。
载波信号:$ c(t) = A\cos\omega_ct $,基带信号为m(t),则已调信号为:(设基带信号m(t)的频谱为$M(\omega)$)
$$ s_m(t)=Am(t)\cos\omega_ct $$
$$ S_m(\omega)= \frac{A}{2}[M(\omega + \omega_c) + M(\omega - \omega_c) ] $$
可以看到,幅度调制就是把基带信号的频谱搬移到$\omega_c$处,再乘以1/2 。是线性变换。
AM调制
$$s_{AM}(t) =[A_0+m(t)]\cos\omega_ct$$
$$ S_{AM}=\pi A_0 [ \delta(\omega + \omega_c) + \delta(\omega - \omega_c) ] + \frac{1}{2}[M(\omega + \omega_c) + M(\omega - \omega_c) ] $$
为使用包络检波的方式进行解调,要求 $|m(t)|<=A_0$
clear all;
%% AM调制
fs = 800; % 采样速率,单位kHz
dt=1/fs; % 采样时间间隔,单位ms
T = 200; % 采样的总时间。频谱分辨率(df=1/T)。
t = 0 : dt : T-dt;
fm = 1; % 调制信号的频率,单位kHz
fc = 10; % 载波信号的频率,单位kHz
m = cos(2*pi*fm*t); % 调制信号
A = 3; %直流信号
s = (m+A).*cos(2*pi*fc*t); %已调信号
[f,sf] = T2F(t,s);
figure(1)
plot(t,s);
axis([0,2,-4,4]);
figure(2)
plot(f,abs(sf));
axis([-15,15,0,max(abs(sf))]);
DSB调制$s_{DSB}(t) =m(t)\cos\omega_ct$ $ S_{DSB}(\omega)=\frac{1}{2}[M(\omega + \omega_c) + M(\omega - \omega_c) ] $,只能用相干解调
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]);
其中的函数T2F是信号的傅里叶变换
%% 函数 计算信号的傅里叶变换
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
还有F2T傅里叶反变换
%% 计算信号频谱 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
来源:oschina
链接:https://my.oschina.net/u/4368300/blog/3797038