利用matlab实现AM调制解调-Toy模板网

这篇具有很好参考价值的文章主要介绍了利用matlab实现AM调制解调。希望对大家有所帮助。如果存在错误或未考虑完全的地方，请大家不吝赐教，您也可以点击"举报违法"按钮提交疑问。

基本原理

幅度调制(AM)是用调制信号去控制高频载波的振幅，使之随调制信号作线性变化的过程。

利用matlab实现AM调制解调

文章来源地址https://www.toymoban.com/news/detail-400189.html

上图给出了AM调制解调的原理模型，从图中可知发送信号m(t)和直流分量A0叠加后乘以高频载波cos(wct)后即可形成AM调制信号。

具体时域波形为

利用matlab实现AM调制解调

对应的频谱波形为

利用matlab实现AM调制解调

在解调端信道输出信号SAM(t)乘以跟发送端同频同相的高频载波cos(wct)后，经过一个低通滤波器提取低频分量之后，得到的信号再扩大为原来的两倍在减去一个直流分量A0(因为经过低通滤波器之后，幅度会变为原来的一半)，最后即可得到原始的基带信号。理论推导如下：

送入解调器的AM表达式如下：

利用matlab实现AM调制解调

与同频同相的相干载波相乘后得到的信号：

利用matlab实现AM调制解调

经历低通滤波器滤除高频信号后得

利用matlab实现AM调制解调

将上述经过低通滤波器之后信号乘以2再减去一个直流分量A0即可得到解调信号Sd(t)

利用matlab实现AM调制解调

利用代码实现AM调制解调

clc
close
A0=1;%直流分量
fc=1000;%载波频率
fs=5000;
t0=1;
t=0:1/fs:t0-1/fs;
x=2*cos(2*pi*100*t)+3*cos(2*pi*500*t);%原始信号
y1=cos(2*pi*fc*t);%载波
N=length(t);
X=fft(x);
Y1=fft(y1);
w=(-N/2:1:N/2-1)./t0;
figure(1)
subplot(211);plot(t,x);title('原始信号的时域波形');
subplot(212);plot(w,(2/N)*abs(fftshift(fft(x))));title('原始信号的频域波形');
figure(2)
subplot(211);plot(t,y1);title('载波信号的时域波形');
subplot(212);plot(w,(2/N)*abs(fftshift(fft(y1))));title('载波信号的频域波形');
y=(A0+x).*y1;%调制信号
Y=fft(y);
figure(3)
subplot(211);plot(t,y);title('调制信号的时域波形');
subplot(212);plot(w,(2/N)*abs(fftshift(Y)));title('调制信号的频域波形');

%解调部分
y2=y.*y1;%相干解调
Y2=fft(y2);
figure(4)
subplot(211);plot(t,y2);title('相干解调之后的时域波形');
subplot(212);plot(w,(2/N)*abs(fftshift(Y2)));title('相干解调之后的频域波形');
y2_filter=filter(hm,y2);%经过低通滤波器
y3=2*(y2_filter)-A0;
y3_fft=fft(y3);
figure(5)
subplot(211);plot(t,y3);title('解调信号的时域波形');
subplot(212);plot(w,(2/N)*abs(fftshift(y3_fft)));title('解调信号的频域波形');

hm.m

function Hd = hm
%HM Returns a discrete-time filter object.

% MATLAB Code
% Generated by MATLAB(R) 9.6 and Signal Processing Toolbox 8.2.
% Generated on: 12-May-2022 10:30:18

% FIR Window Lowpass filter designed using the FIR1 function.

% All frequency values are in Hz.
Fs = 2000; % Sampling Frequency

N = 200; % Order
Fc = 550; % Cutoff Frequency
flag = 'scale'; % Sampling Flag

% Create the window vector for the design algorithm.
win = hamming(N+1);

% Calculate the coefficients using the FIR1 function.
b = fir1(N, Fc/(Fs/2), 'low', win, flag);
Hd = dfilt.dffir(b);

% [EOF]
结果如下：

利用matlab实现AM调制解调