通信原理与MATLAB(一):AM的调制解调

这篇具有很好参考价值的文章主要介绍了通信原理与MATLAB(一):AM的调制解调。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

1.AM的调制原理

AM的调制原理如下图所示,基带信号m(t)和直流量A0相加,然后和高频载波相乘实现AM信号的调制。
通信原理与MATLAB(一):AM的调制解调

2.AM的解调原理

AM的解调原理如下图所示,AM信号经过信道传输之后,再和载波相乘,然后经过低通滤波,隔直流之后恢复出原始基带信号。
通信原理与MATLAB(一):AM的调制解调

3. AM调制解调的代码

AM.m文件,主文件

% AM调制解调过程
%% 基本参数
clear all;                  % 清除变量
close all;                  % 关闭所有窗口图像
fm = 100;                   % 基带信号频率
T = 2;                      % 信号时长
fs = 20000;                 % 采样频率 奈奎斯特采样定理为最大频率的两倍,这里取20倍为了绘制更多的细节,让时域信号更平滑
dt=1/fs;                    % 时间采样间隔,采样频率的倒数
N=T/dt;                     % 采样点个数,总时长除以采样间隔
t=[0:N-1]*dt;               % 采样点的时间序列,作为横坐标

%% ******************调制信号时域波形******************
Am=1;                       % 基带信号幅度 
mt=Am*cos(2*pi*fm*t);       % 基带信号
figure(1);                  % 绘制第一幅图
subplot(221);               % 窗口分割,将一幅图分割成2*2plot(t,mt,'Linewidth',2);   % 时间t为横坐标,基带信号mt为纵坐标绘图,线宽为2
xlabel('t/时间');           % 横坐标标注
ylabel('幅度');             % 纵坐标标注
title('基带信号');          % 图标题标注
axis([0,0.1,-1.1,1.1]);     % 横纵坐标范围设置
line([0,0.1],[0,0],'color','b','Linewidth',2);% 绘制一条从(0,0)(0.1,0)的蓝色实线,线宽为2

%% ******************调制信号频域波形******************
subplot(222);
[mf,msf]=T2F(t,mt);         % 傅里叶变换,得到纵坐标频谱和横坐标频率
plot(mf,abs(msf),'Linewidth',2);    % 画出基带信号频谱,线宽为2
title('基带信号的频谱');            % 图标题标注
xlabel('f/Hz');                     % 横坐标标注
ylabel('幅度/H(f)');                % 纵坐标标注
axis([-150 150 -inf inf]);          % 横纵坐标范围设置

%% ******************载波信号时域波形******************
subplot(223);
fc=1000;                            % 载波频率
zaibo=cos(2*pi*fc*t);               % 载波时域信号
plot(t,zaibo,'r','Linewidth',2);    % 时间t为横坐标,载波信号zaibo为纵坐标绘图,线宽为2,红色
xlabel('t/时间');                   % 横坐标标注
ylabel('幅度');                     % 纵坐标标注
title('载波信号');                  % 图标题标注
axis([0,0.01,-1.1,1.1]);            % 横纵坐标范围设置
line([0,0.01],[0,0],'color','b','Linewidth',2);% 绘制一条从(0,0)(0.01,0)的蓝色实线,线宽为2

%% ******************载波信号频域波形******************
subplot(224);
[mf1,msf1]=T2F(t,zaibo);            % 傅里叶变换,得到纵坐标频谱和横坐标频率
plot(mf1,abs(msf1),'r','Linewidth',2);  % 载波信号频谱,线宽为2,红色
title('载波信号的频谱');            % 图标题标注
xlabel('f/Hz');                     % 横坐标标注
ylabel('幅度/H(f)');                % 纵坐标标注
axis([-1200 1200 -inf inf]);        % 横纵坐标范围设置

%% ******************AM波信号时域波形******************
Ma=0.5;                             % 调制指数
A0=Am/Ma;                           % 直流信号幅度
sam=(A0+mt).*zaibo;                 % 画出AM信号波形
%% 加噪声
SNR=20;                              %信噪比为  dB
sam=awgn(sam,SNR,'measured');
figure(2);                          % 绘制第二幅图
subplot(211);                       % 窗口分割,将一幅图分割成2*1plot(t,sam,'Linewidth',2);          % 画出AM信号波形,线宽为2
hold on;                            % 当前轴及图像保持而不被刷新,准备接受此后将绘制的图形,多图共存
plot(t,A0+mt,'r-');                 % 画出AM信号包络,红色实线
title('调制指数为0.5的AM调制信号及其包络');% 图标题标注
xlabel('t/时间');                   % 横坐标标注
ylabel('幅度');                     % 纵坐标标注
axis([0,0.02,-3.1,3.1]);            % 横纵坐标范围设置
line([0,0.02],[0,0],'color','b','Linewidth',2);% 绘制一条从(0,0)(0.02,0)的蓝色实线,线宽为2

%% ******************AM波信号频域波形******************
[mf2,msf2]=T2F(t,sam);              % 傅里叶变换,得到纵坐标频谱和横坐标频率
subplot(212);
plot(mf2,abs(msf2),'Linewidth',2);  % 画出AM信号频谱
title('AM波信号的频谱');            % 图标题标注
xlabel('f/Hz');                     % 横坐标标注
ylabel('幅度/H(f)');                % 纵坐标标注
axis([-1500 1500 -inf inf]);        % 横纵坐标范围设置

%% ******************相干解调******************
%% ******************已调信号与载波信号相乘******************
st=sam.*zaibo;                      % 已调信号与载波信号相乘
figure(3);                          % 绘制第三幅图
subplot(211);                       % 窗口分割,将一幅图分割成2*1plot(t,st,'Linewidth',2);           % 时间t为横坐标,相乘信号st为纵坐标绘图,线宽为2
title('已调信号与载波信号相乘');    % 图标题标注
xlabel('t/时间');                   % 横坐标标注
ylabel('幅度');                     % 纵坐标标注
axis([0 0.02 -4 4]);                % 横纵坐标范围设置
line([0,0.02],[0,0],'color','b','Linewidth',2);% 绘制一条从(0,0)(0.02,0)的蓝色实线,线宽为2

%% ******************已调信号与载波信号相乘的频谱******************
[f,sf]=T2F(t,st);                   % 傅里叶变换,得到纵坐标频谱和横坐标频率
subplot(212);                       % 窗口分割,将一幅图分割成2*1plot(f,sf,'Linewidth',2);           % 绘制相乘信号st的频谱,线宽为2
title('已调信号与载波信号相乘的频谱');% 图标题标注
xlabel('f/Hz');                     % 横坐标标注
ylabel('幅度/H(f)');                % 纵坐标标注
axis([-2200 2200 -inf inf]);        % 横纵坐标范围设置

%% ******************解调出来的信号******************
[t,st]=lpf(f,sf,2*fm);              % 频域低通滤波
figure(4);                          % 绘制第三幅图
subplot(211);                       % 窗口分割,将一幅图分割成2*1plot(t,2*st-A0,'Linewidth',2);      % 绘制解调波形
title('经过低通滤波的相干解调信号波形');% 图标题标注
xlabel('t/时间');                   % 横坐标标注
ylabel('幅度');                     % 纵坐标标注
axis([0 0.1 -1.1 1.1]);                 % 横纵坐标范围设置 
line([0,0.1],[0,0],'color','b','Linewidth',2);% 绘制一条从(0,0)(0.1,0)的蓝色实线,线宽为2

%% ******************原始信号******************
subplot(212);
plot(t,mt,'r-','Linewidth',2);      % 绘制原始基带信号
title('原始基带信号');                    % 图标题标注
xlabel('t/时间');                   % 横坐标标注
ylabel('幅度');                     % 纵坐标标注
axis([0 0.1 -1.1 1.1]);                 % 横纵坐标范围设置 
line([0,0.1],[0,0],'color','b','Linewidth',2);% 绘制一条从(0,0)(0.1,0)的蓝色实线,线宽为2

T2F.m文件,FFT功能

function [f,sf]= T2F(t,st)      % FFT
% dt = t(2)-t(1);
T=t(end);                       % 输入信号的时间最大值为T
df = 1/T;                       % dt=1/fs; 时间采样间隔,采样频率的倒数;
                                % N=T/dt;  采样点个数,总时长除以采样间隔
                                % 两式联合推导 df = 1/T 
N = length(st);                 % 输入信号时间的长度为采样点数
f=-N/2*df : df : N/2 * df-df;   % 频率分布
sf = fft(st);                   % 做FFT
sf = T/N * fftshift(sf);        % 最后输出,将0-fs频谱搬移到-fs/2-fs/2频谱

F2T.m文件,IFFT功能

function[t,st]=F2T(f,Sf)        % IFFT
df=f(2)-f(1);                   % 频率间隔
fmax=(f(end)-f(1)+df);          % 最大频率减最低频率加上频率间隔为带宽
dt=1/fmax;                      % 采样间隔
N=length(f);                    % 采样点数
t=[0:N-1] * dt;                 % 时间分布
Sf=fftshift(Sf);                %0-fs频谱搬移到-fs/2-fs/2频谱
st=fmax * ifft(Sf);             % 做IFFT
st=real(st);                    % 取实部

lpf.m文件,低通滤波功能

function[t,st]=lpf(f,sf,B)          % 频率LPF
df=f(2)-f(1);                       % 频率间隔
fN=length(f);                       % 采样点数
ym=zeros(1,fN);                     % 生成1行fN列的0向量
xm=floor(B/df);                     % 低频带宽频率除以间隔后的点数向下取整
xm_shift=[-xm:xm-1]+floor(fN/2);    % 因为前面做FFT将0频率搬移到中心处,
                                    % 因此,低通低频频率相应地搬移fN/2,才是对应的频率点
ym(xm_shift)=1;                     % 低通通过频率处幅度为1,其余为0,相当于理想低通
yf=ym.* sf;                         % FFT信号的频谱和对应低频带宽处频率值为1的行向量相乘
[t,st]=F2T(f,yf);                   % IFFT

注意:这四个文件要放在同一文件夹中,第一个是主文件,另外三个是函数文件

4.AM调制解调结果图

结果图中AM信号是经过信道,加了高斯白噪声的。
如果不想加噪声,把下面这行代码注释即可。

sam=awgn(sam,SNR,'measured');

通信原理与MATLAB(一):AM的调制解调
通信原理与MATLAB(一):AM的调制解调
通信原理与MATLAB(一):AM的调制解调
通信原理与MATLAB(一):AM的调制解调

5.AM的优缺点

优点:实现简单,具有直流量便于载波同步;
缺点:直流量占用较大能量,但不传递信息,功率利用率低。文章来源地址https://www.toymoban.com/news/detail-403989.html

到了这里,关于通信原理与MATLAB(一):AM的调制解调的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • 通信原理与MATLAB(十二):MSK的调制解调

    MSK调制原理如下图所示,基带码元先差分编码,然后经过串并转换分成I、Q两路,再与对应的载波相乘,然后再相加完成MSK的调制。 其中注意:I、Q两路码元分别是差分编码后的相对码的奇数和偶数位置上的码元,I路对应于奇数,Q路对应于偶数。而pk是I路码元,但是其码元宽

    2024年02月12日
    浏览(42)
  • 通信原理与MATLAB(八):2PSK的调制解调

    2PSK调制原理如下图所示,和2ASK调制原理相似,只不过基带码元是双极性不归零码,基带码元d(t)和高频载波相乘实现2PSK信号的调制。 波形图如下图所示 2PSK的解调原理如下图所示,2PSK信号经过信道传输之后,再和载波相乘,然后经过低通滤波后抽样判决恢复出原始基带码元

    2024年02月07日
    浏览(39)
  • 通信原理与MATLAB(六):2ASK的调制解调

    2ASK调制原理如下图所示,基带码元d(t)和高频载波相乘实现2ASK信号的调制。 波形图如下图所示 2ASK的解调原理如下图所示,2ASK信号经过信道传输之后,再和载波相乘,然后经过低通滤波后抽样判决恢复出原始基带码元信号。 结果图中2ASK信号是经过信道,加了高斯白噪声的。

    2024年02月10日
    浏览(40)
  • 通信原理与MATLAB(七):2FSK的调制解调

    2FSK调制原理如下图所示,基带码元d(t)中码元为1时,波形为频率为f1的高频载波;基带码元d(t)中码元为0时,波形为频率为f2的高频载波实现2FSK信号的调制,即基带码元和f1的高频正弦波相乘生成2ASK,基带码元的反码和f2的高频正弦波相乘生成第二个2ASK,两个2ASK相加得到2FSK。

    2024年02月09日
    浏览(35)
  • 【通信原理(含matlab程序)】实验二:FM的调制和解调

    💥💥💞💞欢迎来到本博客❤️❤️💥💥 本人持续分享更多关于电子通信专业内容以及嵌入式和单片机的知识,如果大家喜欢,别忘点个赞加个关注哦,让我们一起共同进步~ 理解FM基本原理、FM信号时域和频域主要特点; 理解FM信号调制和解调的实现方法; 掌握matlab程序

    2024年02月08日
    浏览(40)
  • 基于matlab实现双路音频信号的AM调制与解调

    使用audioread读取音频信号,将采样率设置为100kHz,并利用awgn函数为信号添加高斯白噪声。 分别绘制两路音频信号的时频域图。 结果如下: 设置两路不同频率的载波,这里设置频率分别为35kHz和25kHz,载波频率可以根据需要进行更改,但是注意需要满足奈奎斯特采样定理,即

    2023年04月08日
    浏览(51)
  • 【调制解调】AM 调幅

    学习数字信号处理算法时整理的学习笔记。同系列文章目录可见 《DSP 学习之路》目录,代码已上传到 Github - ModulationAndDemodulation。本篇介绍 AM 调幅信号的调制与解调,内附全套 MATLAB 代码。 目录 说明 1. AM 调制算法 1.1 算法描述 1.2 调制信号 m(t) 为确知信号时 1.3 调制信号 m

    2024年02月13日
    浏览(57)
  • 通信原理板块——幅度调制(线性调制)原理(AM、DSB、SSB、VSB)

    微信公众号上线,搜索公众号 小灰灰的FPGA ,关注可获取相关源码,定期更新有关FPGA的项目以及开源项目源码,包括但不限于各类检测芯片驱动、低速接口驱动、高速接口驱动、数据信号处理、图像处理以及AXI总线等 1、调制的定义、目的及分类 (1)调制——将信号形式转换成

    2024年02月04日
    浏览(41)
  • 《移动通信原理与应用》——QPSK调制解调仿真

    目录 一、QPSK调制与解调流程图: 二、仿真运行结果:  三、MATLAB仿真代码:  QPSK调制流程图: QPSK解调流程图:    1、Figure1:为发送端比特流情况图:             从Figure1看出发送端发送的比特流信息…[ak,bk]…情况:奇数进入I路,偶数进入Q路。比特进入I路与Q路情况如

    2024年01月23日
    浏览(39)
  • QPSK调制解调原理及MATLAB仿真

    ⭐️作者简介: 小瑞同学 ,主要学习 FPGA、信号处理、通信 等。 🍎个人主页:小瑞同学的博客主页 🌻个人信条:越努力,越幸运! ⏰日期:2023.11.16 📖文章内容概述:简单介绍了 QPSK 的基本原理和具体流程,并使用 MATLAB 进行了仿真分析。 👻 QPSK(正交相移键控) ,是

    2024年04月15日
    浏览(40)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包