通信原理实验（一）模拟线性调制

最新推荐文章于 2024-08-22 20:41:15 发布
最新推荐文章于 2024-08-22 20:41:15 发布
阅读量327 收藏 6
点赞数 5
文章标签： matlab 开发语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Awanyyy/article/details/141305834
版权

余弦波产生

function signal = cosx(A, fc, t1, t2, fs)
    t = t1 : 1/fs : t2;
    signal = A * cos(2 * pi * fc * t);
end

A为振幅，fc为频率，t1、t2分别为截至时间，fs为采样率。

AM信号

AM信号产生

AM调制程序

function signal = AM(mt, ct)
    A = abs(min(mt));
    signal = (mt + 1.5 * A) .* ct;
end

mt为调制信号，ct为载波信号。

实验主程序

clc;clear;
%%参数设定
A1 = sqrt(2); A2 = 2;
fc1 = 10; fc2 = 1000;
fs = 10000;
t1 = -0.5; t2 = 0.5;
t = t1 : 1/fs : t2;
%% 调制信号与载波设定
mt = cosx(A1, fc1, t1, t2, fs);
ct = cosx(A2, fc2, t1, t2, fs);
%% 调制
yam = AM(mt, ct);
%% 解调
xam = coherent(yam, ct, 5 * fc1, fc2 / 3, fs);
%% 频谱分析
HM = fft(mt); %调制信号FFT
HS = fft(ct); %载波信号FFT
HAM = fft(yam); %已调信号FFT
% Hxam = fft(xam); %解调信号FFT
NF = length(HAM);
nf = 0 : 1/NF : 1 - 1/NF; f = nf * fs;

%% 画图
figure(1);
subplot(2, 3, 1);plot(t, mt,'LineWidth', 1.2);title('调制信号时域波形')
subplot(2, 3, 2);plot(t, ct,'LineWidth', 1.2);title('载波信号时域波形')
subplot(2, 3, 3);plot(t, yam,'LineWidth', 1.2);title('AM信时域波形')
subplot(2, 3, 4);plot(f, abs(HM),'LineWidth', 1.2);title('调制信号频谱')
subplot(2, 3, 5);plot(f, abs(HS),'LineWidth', 1.2);title('载波信号频谱')
subplot(2, 3, 6);plot(f, abs(HAM),'LineWidth', 1.2);title('AM信号频谱')
% figure(2);
% subplot(2, 1, 1);plot(t, xam);title('AM解调信号时域波形')
% subplot(2, 1, 2);plot(f, abs(Hxam));title('AM解调信号频谱')

实验结果

AM信号频谱含有纯载波分量，而只有两个边带分量携带信息，所以AM调制方式的功率利用率较低。

AM信号解调

滤波器设计

function Hd = LPF(fpass, fstop, fs)
%LPF 返回离散时间滤波器对象。

% MATLAB Code
% Generated by MATLAB(R) 23.2 and Signal Processing Toolbox 23.2.
% Generated on: 14-Aug-2024 14:31:40

% Butterworth Lowpass filter designed using FDESIGN.LOWPASS.

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

Fpass = fpass;         % Passband Frequency
Fstop = fstop;         % Stopband Frequency
Apass = 1;           % Passband Ripple (dB)
Astop = 80;          % Stopband Attenuation (dB)
match = 'stopband';  % Band to match exactly

% Construct an FDESIGN object and call its BUTTER method.
h  = fdesign.lowpass(Fpass, Fstop, Apass, Astop, Fs);
Hd = design(h, 'butter', 'MatchExactly', match);

% [EOF]

相干解调程序

function signal = coherent(y, s, fpass, fstop, fs)
    lowpass = LPF(fpass, fstop, fs);
    signal = y .* s;
    signal = filter(lowpass, signal);
end

实验主程序

clc;clear;
%%参数设定
A1 = sqrt(2); A2 = 2;
fc1 = 10; fc2 = 1000;
fs = 10000;
t1 = -0.5; t2 = 0.5;
t = t1 : 1/fs : t2;
%% 调制信号与载波设定
mt = cosx(A1, fc1, t1, t2, fs);
ct = cosx(A2, fc2, t1, t2, fs);
%% 调制
yam = AM(mt, ct);
%% 解调
xam = coherent(yam, ct, 5 * fc1, fc2 / 3, fs);
%% 频谱分析
HM = fft(mt); %调制信号FFT
HS = fft(ct); %载波信号FFT
HAM = fft(yam); %已调信号FFT
Hxam = fft(xam); %解调信号FFT
NF = length(HAM);
nf = 0 : 1/NF : 1 - 1/NF; f = nf * fs;

%% 画图
figure(1);
subplot(2, 3, 1);plot(t, mt,'LineWidth', 1.2);title('调制信号时域波形')
subplot(2, 3, 2);plot(t, ct,'LineWidth', 1.2);title('载波信号时域波形')
subplot(2, 3, 3);plot(t, yam,'LineWidth', 1.2);title('AM信时域波形')
subplot(2, 3, 4);plot(f, abs(HM),'LineWidth', 1.2);title('调制信号频谱')
subplot(2, 3, 5);plot(f, abs(HS),'LineWidth', 1.2);title('载波信号频谱')
subplot(2, 3, 6);plot(f, abs(HAM),'LineWidth', 1.2);title('AM信号频谱')
figure(2);
subplot(2, 1, 1);plot(t, xam);title('AM解调信号时域波形')
subplot(2, 1, 2);plot(f, abs(Hxam));title('AM解调信号频谱')

实验结果

可看出解调信号中含有直流分量，可以理解为AM信号的调制信号必须恒为正值，所以AM的调制信号一定含有直流分量，所以解调后的信号也含有直流分量。

DSB信号

DSB信号产生

DSB调制程序

function signal = DSB(mt, ct)
    signal = mt .* ct;
end

实验主程序


clc;clear;
%%参数设定
A1 = sqrt(2); A2 = 2;
fc1 = 10; fc2 = 1000;
fs = 10000;
t1 = -0.5; t2 = 0.5;
t = t1 : 1/fs : t2;
%% 调制信号与载波设定
mt = cosx(A1, fc1, t1, t2, fs);
st = cosx(A2, fc2, t1, t2, fs);
%% 调制
y = DSB(mt, st);
%% 解调
x = coherent(y, st, 5 * fc1, fc2 / 3, fs);
%% 频谱分析
HM = fft(mt); %调制信号FFT
HS = fft(st); %载波信号FFT
Hy = fft(y); %已调信号FFT
% Hx = fft(x); %解调信号FFT
NF = length(Hy);
nf = 0 : 1/NF : 1 - 1/NF; f = nf * fs;

%% 画图
figure(1);
subplot(2, 3, 1);plot(t, mt);title('调制信号时域波形')
subplot(2, 3, 2);plot(t, st);title('载波信号时域波形')
subplot(2, 3, 3);plot(t, y);title('DSB信时域波形')
subplot(2, 3, 4);plot(f, abs(HM));title('调制信号频谱')
subplot(2, 3, 5);plot(f, abs(HS));title('载波信号频谱')
subplot(2, 3, 6);plot(f, abs(Hy));title('DSB信号频谱')
% figure(2);
% subplot(2, 1, 1);plot(t, x);title('DSB解调信号时域波形')
% subplot(2, 1, 2);plot(f, abs(Hx));title('DSB解调信号频谱')

实验结果

DSB信号的频谱不含载波成分，所以功率利用率为1，所以DSB也称作抑制载波的调制方式，频谱成分由上下边带组成。

DSB信号解调

相干解调程序相同。

实验主程序

clc;clear;
%%参数设定
A1 = sqrt(2); A2 = 2;
fc1 = 10; fc2 = 1000;
fs = 10000;
t1 = -0.5; t2 = 0.5;
t = t1 : 1/fs : t2;
%% 调制信号与载波设定
mt = cosx(A1, fc1, t1, t2, fs);
st = cosx(A2, fc2, t1, t2, fs);
%% 调制
y = DSB(mt, st);
%% 解调
x = coherent(y, st, 5 * fc1, fc2 / 3, fs);
%% 频谱分析
HM = fft(mt); %调制信号FFT
HS = fft(st); %载波信号FFT
Hy = fft(y); %已调信号FFT
Hx = fft(x); %解调信号FFT
NF = length(Hy);
nf = 0 : 1/NF : 1 - 1/NF; f = nf * fs;

%% 画图
figure(1);
subplot(2, 3, 1);plot(t, mt);title('调制信号时域波形')
subplot(2, 3, 2);plot(t, st);title('载波信号时域波形')
subplot(2, 3, 3);plot(t, y);title('DSB信时域波形')
subplot(2, 3, 4);plot(f, abs(HM));title('调制信号频谱')
subplot(2, 3, 5);plot(f, abs(HS));title('载波信号频谱')
subplot(2, 3, 6);plot(f, abs(Hy));title('DSB信号频谱')
figure(2);
subplot(2, 1, 1);plot(t, x);title('DSB解调信号时域波形')
subplot(2, 1, 2);plot(f, abs(Hx));title('DSB解调信号频谱')

实验结果

DSB解调后的信号不含直流成分，因为实验中DSB的调制信号为纯交流信号，不含直流成分。

SSB信号

SSB信号产生

产生原理

SSB信号产生可以使用滤波法，也可以使用相移法，本实验使用相移法。

希尔伯特变换器：就是一个宽带相移器，将每个频率成分都相移pi/2。

具体操作，将调制信号和载波通过希尔伯特变换器后相乘得到y1，将调制信号和载波相乘得到y2，

y = y1 + y2得到保留下边带的SSB信号，y = y1 - y2得到保留上边带的SSB信号。

调制程序

function signal = SSB(mt, ct, mod)
     sth = imag(hilbert(st));
     mth = imag(hilbert(mt));
     if mod == 0
         signal = (1/2) * mt .* ct - (1/2) * mth .* sth;
     end

     if mod == 1
         signal = (1/2) * mt .* ct + (1/2) * mth .* sth;
     end
         
end

注意：matlab中的hilbert(x) 是得到信号x的解析信号，解析信号的虚部才是数学意义上信号x的希尔伯特变换。

实验主程序

%截图保存路径 H:\OneDrive\learn\experimental\principle of communication
clc;clear;
%%参数设定
A1 = sqrt(2); A2 = 2;
fc1 = 10; fc2 = 1000;
fs = 10000;
t1 = -0.5; t2 = 0.5;
t = t1 : 1/fs : t2;
%% 调制信号与载波设定
mt = cosx(A1, fc1, t1, t2, fs);
st = cosx(A2, fc2, t1, t2, fs);
%% 调制
mod = 0; %保留上边带
%mod = 1; %保留下边带
y = SSB(mt, st, mod);
%% 解调
x = coherent(y, st, 5 * fc1, fc2 / 3, fs);
%% 频谱分析
HM = fft(mt); %调制信号FFT
HS = fft(st); %载波信号FFT
Hy = fft(y); %已调信号FFT
% Hx = fft(x); %解调信号FFT
NF = length(Hy);
nf = 0 : 1/NF : 1 - 1/NF; f = nf * fs;

%% 画图
figure(1);
subplot(2, 3, 1);plot(t, mt);title('调制信号时域波形')
subplot(2, 3, 2);plot(t, st);title('载波信号时域波形')
subplot(2, 3, 3);plot(t, y);title('SSB信时域波形')
subplot(2, 3, 4);plot(f, abs(HM));title('调制信号频谱')
subplot(2, 3, 5);plot(f, abs(HS));title('载波信号频谱')
subplot(2, 3, 6);plot(f, abs(Hy));title('SSB信号频谱')
% figure(2);
% subplot(2, 1, 1);plot(t, x);title('SSB解调信号时域波形')
% subplot(2, 1, 2);plot(f, abs(Hx));title('SSB解调信号频谱')

实验结果

 

该实验保留上边带，所以已调信号只有 10 + 1000 Hz这一个中心频率而没有 1000 - 10 Hz这个中心频率。

SSB解调

采用相干解调

实验主程序

%截图保存路径 H:\OneDrive\learn\experimental\principle of communication
clc;clear;
%%参数设定
A1 = sqrt(2); A2 = 2;
fc1 = 10; fc2 = 1000;
fs = 10000;
t1 = -0.5; t2 = 0.5;
t = t1 : 1/fs : t2;
%% 调制信号与载波设定
mt = cosx(A1, fc1, t1, t2, fs);
st = cosx(A2, fc2, t1, t2, fs);
%% 调制
mod = 0; %保留上边带
%mod = 1; %保留下边带
y = SSB(mt, st, mod);
%% 解调
x = coherent(y, st, 5 * fc1, fc2 / 3, fs);
%% 频谱分析
HM = fft(mt); %调制信号FFT
HS = fft(st); %载波信号FFT
Hy = fft(y); %已调信号FFT
Hx = fft(x); %解调信号FFT
NF = length(Hy);
nf = 0 : 1/NF : 1 - 1/NF; f = nf * fs;

%% 画图
figure(1);
subplot(2, 3, 1);plot(t, mt);title('调制信号时域波形')
subplot(2, 3, 2);plot(t, st);title('载波信号时域波形')
subplot(2, 3, 3);plot(t, y);title('SSB信时域波形')
subplot(2, 3, 4);plot(f, abs(HM));title('调制信号频谱')
subplot(2, 3, 5);plot(f, abs(HS));title('载波信号频谱')
subplot(2, 3, 6);plot(f, abs(Hy));title('SSB信号频谱')
figure(2);
subplot(2, 1, 1);plot(t, x);title('SSB解调信号时域波形')
subplot(2, 1, 2);plot(f, abs(Hx));title('SSB解调信号频谱')

实验结果

解调信号同样没有直流成分，原因和DSB相同，但SSB信号占用的带宽为AM和DSB的0.5倍。

Awanyyy
关注
  • 5
    点赞
  • 踩
  • 6
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
基于Matlab实现常见线性调制方法（通信原理）【100011810】
04-11
在实际应用中，为了满足辐射效率，实现对信道的多路复用，提高系统抗干扰能力调制，要对信号进行调制，把原始信号转换成适合在信道中传输的信号。 模拟调制是指用来自信源的基带模拟信号去调制某个载波，常见的模拟调制方式有幅度（线性）调制和角度（非线性）调制，其中调幅（AM）、双边带（DSB）、单边带（SSB）和残留边带（VSB）等调制方式属于幅度调制。
模拟调制仿真_通信原理实验_通信原理MATLAB_
10-01
基于通信原理的基本原理对调制过程进行仿真实验。
通信原理-模拟线性调制
CCC_B的博客
02-27 1428
一些概念 连续波调制：以正弦波为载波的调制方式 线性调制：调制后信号的频谱为调制信号的平移和线性变换  1. 常规双边带调幅 AM  2. 抑制载波双边带调幅DSB  3. 单边带调幅SSB  4. 残留边带调幅VSB 非线性调制：调制后信号与基带信号并非平移和线性变换的关系。  1. 频率调幅FM  2. 相位调幅PM 双边带调幅AM 常规双边带调幅 SAM(t)=[A0+f(t)]cos(ωct+θc)S_{AM}(t)=[A_0+f(t)]cos(\omega_ct+\theta_c)SAM​(t
通信原理之模拟角度调制(非线性调制)详解
starmier的博客
03-12 7848
角度调制（非线性调制）详解1、FM和PM2、FM参数与带宽2.1、调频指数和最大频偏2.2、FM的频谱2.3、FM带宽3、FM的产生与解调3.1、FM的产生3.2、FM的解调（鉴频）4、FM的特点与应用    正弦载波有三个参量: 幅度、频率、相位。       c(t)=Acos[wct+ψ]c(t)=Acos[w_ct+ψ]c(t)=Acos[wc​t+ψ] 这三个参量都可以用来携载消息信息 m(t) 。若 m(t) 被载在幅度上，则为幅度调制,AM(Amplitude Modulation)；若在频率
通信原理之模拟幅度调制（线性调制）详解
热门推荐
starmier的博客
03-12 1万+
幅度调制（线性调制）详解1、基本概念1.1、调制与解调1.2、调制的目的1.3、调制的分类2、幅度调制2.1、基础知识2.2、幅度调制的一般模型2.3、常规调幅 AM(Amplitude Modulation)2.3.1、AM-调幅2.3.2、AM-包络检波器2.3.3、调幅系数（调幅度）m2.3.4、AM调制效率（功率利用率）2.4、双边带调制（DSB-SC）2.4.1、双边带调制（DSB-SC）2.4.2、DSB相干解调2.5、单边带调制 SSB2.6、残留边带调制 VSB 1、基本概念 1.1、调制与
通信原理02：模拟调制
m0_46521579的博客
09-12 3585
模拟调制：线性调制（AM，DSB，SSB，VSB）；非线性调制（FM，PM）；抗噪声性能
matlab通信原理实验低频高频,通信原理实验一 模拟调制已调信号的波形和频谱.doc...
weixin_39890517的博客
03-20 1446
通信原理实验一 模拟调制已调信号的波形和频谱.doc (17页) 本资源提供全文预览，点击全文预览即可全文预览,如果喜欢文档就下载吧，查找使用更方便哦！14.90 积分《通信原理》课程实验报告 HUNAN UNIVERSITY通信原理课程实验报告 题 目： 模拟调制已调信号的波形和频谱 学生姓名： 学生学号： 专业班级： 　 指导老师： 杜青松 目 录1.实验目的 12.实验原理 12.1AM...
通信原理模拟线性调制系统(共55张PPT)精选.pptx
11-15
通信原理模拟线性调制系统(共55张PPT)精选.pptx
通信原理课件：第三章 模拟线性调制.ppt
07-07
通信原理课件：第三章 模拟线性调制.ppt
模拟线性调制系统的仿真
01-31
对通信中模拟调制系统理解，DSB，SSB，AM调制的理解，仿真，滤波器，星座图等 1. 编写AM、DSB、SSB的调制程序，并画出时域波形和频谱图。 2. 完成DSB调制的相干解调，并绘制解调波形和频谱图。
通信原理实验一.txt
11-06
1. 加深对模拟调制系统基本原理和基本组成的理解。 2. 熟悉调制和解调过程中信号时域和频域特性的变化过程。 3. 编写AM、DSB、SSB的调制程序，并画出时域波形和频谱图。 4. 完成DSB调制的相干解调，并绘制解调波形和...
Matlab simulink建模与仿真 第一章（simulink入门）
Zevalin的博客
08-18 1079
（1）MATLAB是一种以矩阵为来处理数据的的计算软件，其应用范围十分广泛，该产品由若干模块组成，simulink是其中一个，它可利用图形化的工具来进行建模与仿真。
MATLAB 手动实现投影密度法分割建筑物立面 （73）
weixin_44329757的博客
08-18 142
从原始点云中，自动分割提取建筑物立面点云用于立面绘图，可以减少人为操作流程。这里从0开始，手动实现一种基于投影密度法的建筑物立面点云提取操作，其原理简单，效果优秀，对于城市内的规则建筑物具有极好的提取立面效果，下面是提取的效果预览，和具体的实现代码，直接复制粘贴即可使用。（红色点位提取的立面点，顶面和树木，地面都被有效过滤）这是一种稳定，常用，易于实现的立面点云提取算法，效果较好。
3.5、matlab打开显示保存点云文件(.ply/.pcd)以及经典点云模型数据
XU157303764的博客
08-22 1433
matlab打开显示保存点云文件(.ply/.pcd)以及经典点云模型数据 (点云模型连接附下)
路径规划 | 灰狼算法+B样条曲线优化无人机三维路径规划（Matlab）
关注并私信文章链接，获取对应文章源码和数据，机器学习之心的博客。
08-18 442
路径规划 | 灰狼算法+B样条曲线优化无人机三维路径规划（Matlab）
MATLAB CSF布料模拟滤波分类地面点和地物点（71）
weixin_44329757的博客
08-18 195
CSF（Cloth Simulation Filter）地面滤波算法是一种用于点云数据处理的技术，特别是在自动驾驶、机器人视觉和地理信息系统等领域中，用于从测量数据中提取出地面点和非地面点。CSF地面滤波方法，始终是开源，可快速上手使用，滤波效果较好，能适应大多数场景的地面与非地面分离方法，作者提供了不同语言的版本使用，这里的MATLAB版之前写过C++版的配置和使用，也可作为参考PCL点云处理之CSF布料模拟滤波。
MATLAB根据数值画直方图
u013600306的博客
08-18 185
一直很纠结MATLAB为什么不提供根据数值（或统计值）画直方图的函数，只给一个不专业的bar，原来histogram支持。
Ceres&PCL 最小二乘插值（曲线拟合）
dayuhaitang1的博客
08-18 56
Ceres&PCL 最小二乘插值（曲线拟合）
3.3、matlab彩色图和灰度图的二值化算法汇总
最新发布
XU157303764的博客
08-22 1477
matlab彩色图和灰度图的二值化算法汇总
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值