双边带调幅建模与仿真（matlab仿真）

.求

已于 2023-07-16 17:01:04 修改

阅读量1k

点赞数 2

分类专栏： matla通信建模仿真文章标签： matlab 开发语言信息与通信信号处理

于 2023-07-15 22:52:09 首次发布

本文链接：https://blog.csdn.net/m0_56324525/article/details/131744943

版权

matla通信建模仿真专栏收录该内容

6 篇文章 7 订阅

订阅专栏

一、实验目的

了解双边带调幅调制解调工作原理，掌握AM调制解调系统的构成，基带信号、载波和已调信号之间的关系；掌握利用matlab对双边带调幅进行建模与仿真的分析方法。

二、实验任务

某消息信号m(t) ，其中m(t) 为单音频信号，频率为6420Hz；用AM方式调制载波c(t)=cos⁡(2πfct) ，用包络检波器或相干解调法解调该信号，分别给出：

（1）消息信号、已调信号和解调信号的波形，并绘制它们的频谱；

（2）已调信号的功率和调制效率；

（3）假设已调信号通过AWGN信道，信噪比为20dB，给出解调后的信号与原始信号。

简述AM调制解调系统的工作原理和程序设计思路，分别给出（1）无噪声时，消息信号、已调信号和解调信号的波形及其频谱，已调信号的功率和调制效率；（2）加入AWGN时，消息信号和解调信号的波形，并对仿真结果进行分析获得有效结论。

三、实验代码

Fs=1000000;              %采样频率
fc=100000;               %载波中心频率
A0=1.5;                  %直流分量
T=2;                    %信号时长
N=T*Fs;                 %采样点个数
t=[0:N-1]/Fs;            	%采样点的时间序列
w=2*pi*642*10;           %信号频率
wc=2*pi*fc;

mt=cos(w*t);             %信源
figure(1);
subplot(321);plot(t,mt);
title('基带调制信号');
axis([0 0.002 -1.5 1.5]);
line([0,0.004],[0,0],'color','k');
Mk=fft(mt);              %信号频谱：1024点fft
fmk=(0:length(Mk)-1)*Fs/length(Mk);
subplot(322);plot(fmk,abs(Mk));
axis([0 10000 0 1.1*max(abs(Mk))]);
title('基带调制信号频谱');
%mt的最大值是2

Fc=cos(wc * t);          %载波信号
subplot(323);plot(t,Fc);
title('载波信号');axis([0 0.001 -1.5 1.5]);
Fck=fft(Fc);            %载波频谱
fck=(0:length(Fck)-1)*Fs/length(Fck);
subplot(324);plot(fck,abs(Fck));
axis([0 1.1*fc 0 1.1*max(abs(Fck))]);
title('载波频谱');
am=(mt+A0).* Fc;      %AM信号
subplot(325);plot(t,am);
hold on;                %画出AM信号波形
plot(t,mt+A0,'r--');title('调制信号及其包络');
axis([0 0.001 -4 4]);
Amk=fft(am);            %AM频谱：
fam=(0:length(Amk)-1)*Fs/length(Amk);
subplot(326);plot(fam,abs(Amk));
axis([0 120000 0 1.1*max(abs(Amk))]);
title('AM信号频谱');


%无噪声时相干解调
figure(2);
st=am.* cos(wc*t);
subplot(411);
plot(t,am);title('无噪声时AM信号与载波信号相乘');
axis([0 0.002 -4 4]);

st_fft=fft(st);         %傅里叶变换
%巴特沃斯低通滤波器
wp=0.2*pi;wst=0.3*pi;ap=1;as=15;
T0=1;fs=1/T0;
[N,Wc]=buttord(wp,wst,ap,as,'s');
[B,A]=butter(N,Wc,'s');
[Bz,Az]=bilinear(B,A,fs);
mt0=filter(Bz,Az,st);

subplot(412);plot(t,mt0);
title('无噪声时经过低通滤波的相干解调信号波形');
axis([0 0.002 0 1.5]);
Mk0=fft(mt0);              
fmk0=(0:length(Mk0)-1)*Fs/length(Mk0);
subplot(413);plot(fmk0,abs(Mk0));
axis([0 10000 0 1.1*max(abs(Mk0))]);
title('无噪声时解调出调制信号频谱');
subplot(414);plot(t,mt);
title('基带调制信号');
axis([0 0.002 -1.5 1.5]);

% 通过AWGN信道
figure(3);
snr=20;
am_awgn=awgn(am,snr);
st=am_awgn.* cos(wc*t);
subplot(411);
plot(t,am_awgn);
title('过AWGN后AM信号与载波信号相乘');
axis([0 0.002 -4 4]);
mt1=filter(Bz,Az,st);
subplot(412);plot(t,mt1);
title('过AWGN后经过低通滤波的相干解调信号波形');
axis([0 0.002 0 1.5]);
Mk1=fft(mt1);           
fmk1=(0:length(Mk1)-1)*Fs/length(Mk1);
subplot(413);plot(fmk1,abs(Mk1));
axis([0 10000 0 1.1*max(abs(Mk1))]);
title('有噪声时解调出调制信号频谱');

subplot(414);plot(t,mt);
title('基带调制信号');
axis([0 0.002 -1.5 1.5]);