【调制解调】ISB 独立边带调幅

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

说明

学习数字信号处理算法时整理的学习笔记。同系列文章目录可见 《DSP 学习之路》目录,代码已上传到 Github - ModulationAndDemodulation。本篇介绍 ISB 独立边带调幅信号的调制与解调,内附全套 MATLAB 代码。

目录
  • 说明
  • 1. ISB 调制算法
    • 1.1 算法描述
    • 1.2 ISB 信号调制示例
  • 2. ISB 解调算法
    • 2.1 数字正交解调
  • 参考资料
  • 附录代码
    • 附.1 文件 lpf_filter.m
    • 附.2 文件 mod_isb.m
    • 附.3 文件 demod_isb.m
    • 附.4 文件 main_modISB_example.m
    • 附.5 文件 main_demodISB_example.m

1. ISB 调制算法

1.1 算法描述

如果发射机仍然发射两个边带,但是和双边带 DSB 信号不同,两个边带中含有两种不同的信息,这种调制方式叫独立边带调制(ISB, Independent Side Band)。ISB 信号的带宽等于两个基带信号(调制信号)带宽之和,即 \(B_{ISB}={f_{HU}}+{f_{HL}}\)。ISB 信号的时域表达式为:

\[s_{ISB}(t)=\left[m_U(t)+m_L(t)\right]cos(\omega_ct) - \left[\hat{m}_U(t)-\hat{m}_L(t)\right]sin(\omega_ct) \tag{1} \]

式中,\(m_U(t)\)\(m_L(t)\) 分别为上、下边带信号,\(\hat{m}_U(t)\)\(\hat{m}_L(t)\) 分别是上、下边带信号的 Hilbert 变换,推导过程可参考本人同系列文章 【调制解调】SSB 单边带调幅。可仿照 SSB 移相法的思路获得 ISB 信号。

1.2 ISB 信号调制示例

调制信号 \(m(t)\) 可以是确知信号,也可以是随机信号。当 \(m(t)\) 是确知信号时,不妨假设上边带信号 \(m_U(t)\) 的时域表达式如下:

\[m_U(t) = sin(2{\pi}{f_m}t)+cos({\pi}{f_m}t) \tag{2} \]

下边带信号 \(m_L(t)\) 的时域表达式如下:

\[m_L(t) = sin(3{\pi}{f_m}t)+cos(4{\pi}{f_m}t) \tag{3} \]

各调制参数取值:\(f_m=2500Hz\)\(f_c=20000Hz\)。信号采样率 \(f_s=8{f_c}\),仿真总时长为 \(2s\)。ISB 独立边带调制效果如下图所示(为了美观,时域只显示前 500 个点),上边带调制信号 \(m_U(t)\) 双边幅度谱有四根离散谱线(\({\pm}2500Hz\)\({\pm}1250Hz\)),下边带调制信号 \(m_L(t)\) 双边幅度谱有四根离散谱线(\({\pm}3750Hz\)\({\pm}5000Hz\)),ISB 独立边带信号有八根离散谱线(\(\pm15000Hz\)\(\pm16250Hz\)\(\pm21250Hz\)\(\pm22500Hz\))。

【调制解调】ISB 独立边带调幅

代码详见 mod_isb.mmain_modISB_example.m


2. ISB 解调算法

ISB 信号解调的相关资料较少,这里仿照 SSB 信号数字正交解调的方法做了一个 ISB 信号数字正交解调。

2.1 数字正交解调

ISB 数字正交解调一般有以下三个步骤:

  1. 第一步:乘以正交相干载波得到 \({s_I}(t)\)\({s_Q}(t)\),即 \({s_I}(t)=s(t)cos({\omega_ct}+{\phi_0})\)\({s_Q}(t)=-s(t)sin({\omega_ct}+{\phi_0})\)
  2. 第二步:低通滤波器滤除 \({s_I}(t)\)\({s_Q}(t)\) 中的高频分量。
  3. 第三步:计算 \({s_Q}(t)\) 的希尔伯特变换 \(\hat{s}_Q(t)\),进一步可得上边带信号 \(m_U(t)\) 的解调结果 \(m_{Uo}=s_{I}(t)-\hat{s}_{Q}(t)\),下边带信号 \(m_L(t)\) 的解调结果 \(m_{Lo}=s_{I}(t)+\hat{s}_{Q}(t)\)

对 1.2 节中的 ISB 信号,设定信噪比 \(SNR=50dB\),上边带解调效果如下,计算误差,有:\(\sqrt{\sum{{\lvert}m(t_i)-\hat{m}(t_i){\rvert}^2}}/\sqrt{\sum{{\lvert}m(t_i){\rvert}^2}}\approx0.0022\)

【调制解调】ISB 独立边带调幅

下边带解调效果如下,计算误差,有:\(\sqrt{\sum{{\lvert}m(t_i)-\hat{m}(t_i){\rvert}^2}}/\sqrt{\sum{{\lvert}m(t_i){\rvert}^2}}\approx0.0022\)

【调制解调】ISB 独立边带调幅

代码详见 demod_isb.mmain_demodISB_example.m。更改相干载波的初始相位为 \({\phi_0}=\pi/4,\pi/2\),或者更改相干载波的中心频率为 \(0.8f_c,1.2f_c\) 后,解调效果变差,说明这种方法对相干载波同频同相的要求较高。


参考资料

[1] 楼才义,徐建良,杨小牛.软件无线电原理与应用[M].电子工业出版社,2014.文章来源地址https://www.toymoban.com/news/detail-558050.html


附录代码

附.1 文件 lpf_filter.m

function sig_lpf = lpf_filter(sig_data, cutfre)
% LPF_FILTER    自定义理想低通滤波器
% 输入参数:
%       sig_data        待滤波数据
%       cutfre          截止频率,范围 (0,1)
% 输出参数:
%       sig_lpf         低通滤波结果
% @author 木三百川

nfft = length(sig_data);
lidx = round(nfft/2-cutfre*nfft/2);
ridx = nfft - lidx;
sig_fft_lpf = fftshift(fft(sig_data));
sig_fft_lpf([1:lidx,ridx:nfft]) = 0;
sig_lpf = real(ifft(fftshift(sig_fft_lpf)));

end

附.2 文件 mod_isb.m

function [ sig_isb ] = mod_isb(fc, fs, mut, mlt, t)
% MOD_ISB        ISB 独立边带调制
% 输入参数:
%       fc      载波中心频率
%       fs      信号采样率
%       mut     上边带调制信号
%       mlt     下边带调制信号
%       t       采样时间
% 输出参数:
%       sig_isb ISB 独立边带调幅实信号
% @author 木三百川

% 计算 mu(t) 与 ml(t) 的希尔伯特变换(相移)
hmut = imag(hilbert(mut));
hmlt = imag(hilbert(mlt));

% 与正交载波相合成
sig_isb = (mut+mlt).*cos(2*pi*fc*t)-(hmut-hmlt).*sin(2*pi*fc*t);

% 绘图
nfft = length(sig_isb);
freq = (-nfft/2:nfft/2-1).'*(fs/nfft);
figure;set(gcf,'color','w');
plot_length = min(500, length(sig_isb));
subplot(3,2,1);
plot(t(1:plot_length), mut(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('上边带调制信号mu(t)');
subplot(3,2,2);
plot(freq, 10*log10(fftshift(abs(fft(mut,nfft)/nfft))+eps));xlim([freq(1),freq(end)]);
xlabel('频率/hz');ylabel('幅度/dB');title('上边带调制信号mu(t)双边幅度谱');

subplot(3,2,3);
plot(t(1:plot_length), mlt(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('下边带调制信号ml(t)');
subplot(3,2,4);
plot(freq, 10*log10(fftshift(abs(fft(mlt,nfft)/nfft))+eps));xlim([freq(1),freq(end)]);
xlabel('频率/hz');ylabel('幅度/dB');title('下边带调制信号ml(t)双边幅度谱');

subplot(3,2,5);
plot(t(1:plot_length), sig_isb(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('ISB独立边带调幅信号s(t)');
subplot(3,2,6);
plot(freq, 10*log10(fftshift(abs(fft(sig_isb,nfft)/nfft))+eps));xlim([freq(1),freq(end)]);
xlabel('频率/hz');ylabel('幅度/dB');title('ISB独立边带调幅信号s(t)双边幅度谱');

end

附.3 文件 demod_isb.m

function [ sig_isbu_demod,sig_isbl_demod ] = demod_isb(sig_isb_receive, fc, fs, t, phi0)
% DEMOD_ISB                 ISB 数字正交解调
% 输入参数:
%       sig_isb_receive     SSB 接收信号,行向量
%       fc                  载波中心频率
%       fs                  信号采样率
%       t                   采样时间
%       phi0                载波初始相位
% 输出参数:
%       sig_isbu_demod      上边带解调结果,与 sig_isb_receive 等长
%       sig_isbl_demod      下边带解调结果,与 sig_isb_receive 等长
% @author 木三百川

% 第一步:乘以正交相干载波
sig_isb_i = sig_isb_receive.*cos(2*pi*fc*t+phi0);
sig_isb_q = -sig_isb_receive.*sin(2*pi*fc*t+phi0);

% 第二步:低通滤波
sig_isb_i_lpf = lpf_filter(sig_isb_i, fc/(fs/2));
sig_isb_q_lpf = lpf_filter(sig_isb_q, fc/(fs/2));

% 第三步:计算希尔伯特变换
sig_isb_q_lpf = imag(hilbert(sig_isb_q_lpf));
sig_isbu_demod = sig_isb_i_lpf-sig_isb_q_lpf;
sig_isbl_demod = sig_isb_i_lpf+sig_isb_q_lpf;

end

附.4 文件 main_modISB_example.m

clc;
clear;
close all;
% ISB 调制仿真(调制信号为确知信号,相移法)
% @author 木三百川

% 调制参数
fm = 2500;              % 调制信号参数
fc = 20000;             % 载波频率
fs = 8*fc;              % 采样率
total_time = 2;         % 仿真时长,单位:秒

% 采样时间
t = 0:1/fs:total_time-1/fs;

% 调制信号为确知信号
mut = sin(2*pi*fm*t)+cos(pi*fm*t);
mlt = sin(3*pi*fm*t)+cos(4*pi*fm*t);

% ISB 调制
[ sig_isb ] = mod_isb(fc, fs, mut, mlt, t);

附.5 文件 main_demodISB_example.m

clc;
clear;
close all;
% ISB 解调仿真(调制信号为确知信号,数字正交解调)
% @author 木三百川

% 调制参数
fm = 2500;              % 调制信号参数
fc = 20000;             % 载波频率
fs = 8*fc;              % 采样率
total_time = 2;         % 仿真时长,单位:秒

% 采样时间
t = 0:1/fs:total_time-1/fs;

% 调制信号为确知信号
mut = sin(2*pi*fm*t)+cos(pi*fm*t);
mlt = sin(3*pi*fm*t)+cos(4*pi*fm*t);

% ISB 调制
[ sig_isb_send ] = mod_isb(fc, fs, mut, mlt, t);

% 加噪声
snr = 50;               % 信噪比
sig_isb_receive = awgn(sig_isb_send, snr, 'measured');

% 数字正交解调
phi0 = 0;
[ sig_isbu_demod,sig_isbl_demod ] = demod_isb(sig_isb_receive, fc, fs, t, phi0);

% 绘图
nfft = length(sig_isb_receive);
freq = (-nfft/2:nfft/2-1).'*(fs/nfft);
figure;set(gcf,'color','w');
plot_length = min(500, length(sig_isb_receive));
subplot(1,2,1);
plot(t(1:plot_length), sig_isb_receive(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('ISB接收信号');
subplot(1,2,2);
plot(freq, 10*log10(fftshift(abs(fft(sig_isb_receive,nfft)/nfft))+eps));xlim([freq(1),freq(end)]);
xlabel('频率/hz');ylabel('幅度/dB');title('ISB接收信号双边幅度谱');

figure;set(gcf,'color','w');
plot(t(1:plot_length), mut(1:plot_length));xlim([t(1),t(plot_length)]);
hold on;
plot(t(1:plot_length), sig_isbu_demod(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('上边带解调效果');
legend('上边带调制信号','上边带解调信号');

figure;set(gcf,'color','w');
plot(t(1:plot_length), mlt(1:plot_length));xlim([t(1),t(plot_length)]);
hold on;
plot(t(1:plot_length), sig_isbl_demod(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('下边带解调效果');
legend('下边带调制信号','下边带解调信号');

coefu = mean(abs(mut))/mean(abs(sig_isbu_demod));
fprintf('norm(上边带调制信号 - %.2f * 上边带解调信号)/norm(上边带调制信号) = %.4f.\n', coefu, norm(mut-coefu*sig_isbu_demod)/norm(mut));

coefl = mean(abs(mlt))/mean(abs(sig_isbl_demod));
fprintf('norm(下边带调制信号 - %.2f * 下边带解调信号)/norm(下边带调制信号) = %.4f.\n', coefl, norm(mlt-coefl*sig_isbl_demod)/norm(mlt));

到了这里,关于【调制解调】ISB 独立边带调幅的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

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

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

相关文章

  • 数字带通调制系统实验-BPSK的调制与解调

    1、实验原理 BPSK的调制原理 在二进制数字调制中,当正弦载波的相位随二进制数字基带信号离散变化 时,则产生二进制移相键控(2PSK)信号。通常用已调信号载波的0度和180 度分别表示二进制数字基带信号的1和0.二进制移相键控信号的时域表达式为 这种以载波的不同相位直接

    2024年02月02日
    浏览(65)
  • GMSK调制与解调学习记录

    GMSK(高斯最小频移键控)信号是在MSK(最小频移键控)信号的基础上得到的。而MSK信号本质上为连续相位调制(CPM)信号,是一种特殊的连续相位的频移键控 (CPFSK)。其最大频移为比特速率的1/4,即MSK是调制系数为0.5的连续相位的FSK。在FSK调制方式中,根据原始的信息序列,相

    2024年02月03日
    浏览(36)
  • 2ASK调制解调实验

    一,数字信号的键控调制 在数字调制技术中,由于数字信号具有离散值的特点,因此数字调制技术通常有两种 方法:一是利用模拟调制的方法去实现数字调制,即把数字调制看成是模拟调制的一个特 例,把数字基带信号当做模拟信号的特殊情况来处理;二是利用数字信号的

    2023年04月08日
    浏览(43)
  • 【调制解调】PM 调相

    学习数字信号处理算法时整理的学习笔记。同系列文章目录可见 《DSP 学习之路》目录,代码已上传到 Github - ModulationAndDemodulation。本篇介绍 PM 调相信号的调制与解调,内附全套 MATLAB 代码。 目录 说明 1. PM 调制算法 1.1 PM 信号描述 1.2 PM 信号调制示例 2. PM 解调算法 2.1 FM 解调

    2024年02月16日
    浏览(41)
  • Verilog正交调制解调

    了解正交调制解调的原理和实现方法 学会 I P IP I P 核的使用 学会利用 m o d e s i m modesim m o d es im 进行仿真 相关参数: (1)直线阵通道数: 96 96 96 (2)信号频率: 10 k H z 10kHz 10 k Hz (3)采样率: 400 k H z 400kHz 400 k Hz (4)低通滤波器阶数: 64 64 64 (5)低通滤波器截止频率:

    2023年04月15日
    浏览(68)
  • FM调制解调---FPGA

            实验通过编写一个DMA读模块获取FM调制的数据源,DMA模块的实现是基于AXI协议。因为在数据的传输中,Xilinx提供的官方DMA IP核在传输完一次突发数据后需要在PS 端重新启动一次都或者写操作,如此的话,在进行大量数据的传输工作时,尤其是对DDR 不同地址区域同时进

    2024年02月11日
    浏览(47)
  • 【调制解调】FM 调频

    学习数字信号处理算法时整理的学习笔记。同系列文章目录可见 《DSP 学习之路》目录,代码已上传到 Github - ModulationAndDemodulation。本篇介绍 FM 调频信号的调制与解调,内附全套 MATLAB 代码。 目录 说明 1. FM 调制算法 1.1 FM 信号描述 1.2 FM 信号的带宽与功率分配 1.3 FM 信号的调制

    2024年02月16日
    浏览(41)
  • 利用matlab实现AM调制解调

    基本原理 幅度调制(AM)是用调制信号去控制高频载波的振幅,使之随调制信号作线性变化的过程。     上图给出了AM调制解调的原理模型,从图中可知发送信号m(t)和直流分量A 0 叠 加 后乘以高频载波cos(w c t)后即可形成AM调制信号。 具体时域波形为   对应的频谱波形为   在解

    2023年04月08日
    浏览(49)
  • BLE调制与解调的一些东西

    BLE是GFSK的IQ调制 所谓IQ调制,就是利用IQ两个分量序列去控制两路正交信号,I和Q两个序列可以是任意数字,也可以是符合某些规律的序列。 总的原理公式就是: 看这篇博客的图片,很清楚。 https://blog.csdn.net/qq_41019681/article/details/111305603 想要更深入的理解,就看这篇,别人写

    2024年02月02日
    浏览(39)
  • 关于64QAM调制,软硬解调详解

    本文主要梳理记录一下64QAM的调制原理,以及软解调和硬解调的区别。 64QAM调制是M-QAM调制的一种,属于正交振幅调制,即采用幅度相位相结合的调制方式使得一个码片包含更多bit的信息。 64QAM调制可得到64个不同的波形,分别代表000000,000001…这也意味着一共有64种符号,一

    2023年04月22日
    浏览(48)

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

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

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

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

二维码1

领取红包

二维码2

领红包