ToB企服应用市场:ToB评测及商务社交产业平台

标题: 【调制解调】ISB 独立边带调幅 [打印本页]
作者: 欢乐狗    时间: 2023-7-14 08:55
标题: 【调制解调】ISB 独立边带调幅
说明

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

目录


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\)）。

代码详见 mod_isb.m、main_modISB_example.m。

2. ISB 解调算法

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

ISB 数字正交解调一般有以下三个步骤：
对 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\)。

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

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

参考资料

[1] 楼才义,徐建良,杨小牛.软件无线电原理与应用[M].电子工业出版社,2014.

附录代码

附.1 文件 lpf_filter.m
  1. function sig_lpf = lpf_filter(sig_data, cutfre)
  2. % LPF_FILTER    自定义理想低通滤波器
  3. % 输入参数：
  4. %       sig_data        待滤波数据
  5. %       cutfre          截止频率，范围 (0,1)
  6. % 输出参数：
  7. %       sig_lpf         低通滤波结果
  8. % @author 木三百川
  10. nfft = length(sig_data);
  11. lidx = round(nfft/2-cutfre*nfft/2);
  12. ridx = nfft - lidx;
  13. sig_fft_lpf = fftshift(fft(sig_data));
  14. sig_fft_lpf([1:lidx,ridx:nfft]) = 0;
  15. sig_lpf = real(ifft(fftshift(sig_fft_lpf)));
  17. end
复制代码
附.2 文件 mod_isb.m
  1. function [ sig_isb ] = mod_isb(fc, fs, mut, mlt, t)
  2. % MOD_ISB        ISB 独立边带调制
  3. % 输入参数：
  4. %       fc      载波中心频率
  5. %       fs      信号采样率
  6. %       mut     上边带调制信号
  7. %       mlt     下边带调制信号
  8. %       t       采样时间
  9. % 输出参数：
  10. %       sig_isb ISB 独立边带调幅实信号
  11. % @author 木三百川
  13. % 计算 mu(t) 与 ml(t) 的希尔伯特变换（相移）
  14. hmut = imag(hilbert(mut));
  15. hmlt = imag(hilbert(mlt));
  17. % 与正交载波相合成
  18. sig_isb = (mut+mlt).*cos(2*pi*fc*t)-(hmut-hmlt).*sin(2*pi*fc*t);
  20. % 绘图
  21. nfft = length(sig_isb);
  22. freq = (-nfft/2:nfft/2-1).'*(fs/nfft);
  23. figure;set(gcf,'color','w');
  24. plot_length = min(500, length(sig_isb));
  25. subplot(3,2,1);
  26. plot(t(1:plot_length), mut(1:plot_length));xlim([t(1),t(plot_length)]);
  27. xlabel('t/s');ylabel('幅度');title('上边带调制信号mu(t)');
  28. subplot(3,2,2);
  29. plot(freq, 10*log10(fftshift(abs(fft(mut,nfft)/nfft))+eps));xlim([freq(1),freq(end)]);
  30. xlabel('频率/hz');ylabel('幅度/dB');title('上边带调制信号mu(t)双边幅度谱');
  32. subplot(3,2,3);
  33. plot(t(1:plot_length), mlt(1:plot_length));xlim([t(1),t(plot_length)]);
  34. xlabel('t/s');ylabel('幅度');title('下边带调制信号ml(t)');
  35. subplot(3,2,4);
  36. plot(freq, 10*log10(fftshift(abs(fft(mlt,nfft)/nfft))+eps));xlim([freq(1),freq(end)]);
  37. xlabel('频率/hz');ylabel('幅度/dB');title('下边带调制信号ml(t)双边幅度谱');
  39. subplot(3,2,5);
  40. plot(t(1:plot_length), sig_isb(1:plot_length));xlim([t(1),t(plot_length)]);
  41. xlabel('t/s');ylabel('幅度');title('ISB独立边带调幅信号s(t)');
  42. subplot(3,2,6);
  43. plot(freq, 10*log10(fftshift(abs(fft(sig_isb,nfft)/nfft))+eps));xlim([freq(1),freq(end)]);
  44. xlabel('频率/hz');ylabel('幅度/dB');title('ISB独立边带调幅信号s(t)双边幅度谱');
  46. end
复制代码
附.3 文件 demod_isb.m
  1. function [ sig_isbu_demod,sig_isbl_demod ] = demod_isb(sig_isb_receive, fc, fs, t, phi0)
  2. % DEMOD_ISB                 ISB 数字正交解调
  3. % 输入参数：
  4. %       sig_isb_receive     SSB 接收信号，行向量
  5. %       fc                  载波中心频率
  6. %       fs                  信号采样率
  7. %       t                   采样时间
  8. %       phi0                载波初始相位
  9. % 输出参数：
  10. %       sig_isbu_demod      上边带解调结果，与 sig_isb_receive 等长
  11. %       sig_isbl_demod      下边带解调结果，与 sig_isb_receive 等长
  12. % @author 木三百川
  14. % 第一步：乘以正交相干载波
  15. sig_isb_i = sig_isb_receive.*cos(2*pi*fc*t+phi0);
  16. sig_isb_q = -sig_isb_receive.*sin(2*pi*fc*t+phi0);
  18. % 第二步：低通滤波
  19. sig_isb_i_lpf = lpf_filter(sig_isb_i, fc/(fs/2));
  20. sig_isb_q_lpf = lpf_filter(sig_isb_q, fc/(fs/2));
  22. % 第三步：计算希尔伯特变换
  23. sig_isb_q_lpf = imag(hilbert(sig_isb_q_lpf));
  24. sig_isbu_demod = sig_isb_i_lpf-sig_isb_q_lpf;
  25. sig_isbl_demod = sig_isb_i_lpf+sig_isb_q_lpf;
  27. end
复制代码
附.4 文件 main_modISB_example.m
  1. clc;
  2. clear;
  3. close all;
  4. % ISB 调制仿真(调制信号为确知信号，相移法)
  5. % @author 木三百川
  7. % 调制参数
  8. fm = 2500;              % 调制信号参数
  9. fc = 20000;             % 载波频率
  10. fs = 8*fc;              % 采样率
  11. total_time = 2;         % 仿真时长，单位：秒
  13. % 采样时间
  14. t = 0:1/fs:total_time-1/fs;
  16. % 调制信号为确知信号
  17. mut = sin(2*pi*fm*t)+cos(pi*fm*t);
  18. mlt = sin(3*pi*fm*t)+cos(4*pi*fm*t);
  20. % ISB 调制
  21. [ sig_isb ] = mod_isb(fc, fs, mut, mlt, t);
复制代码
附.5 文件 main_demodISB_example.m
  1. clc;
  2. clear;
  3. close all;
  4. % ISB 解调仿真(调制信号为确知信号，数字正交解调)
  5. % @author 木三百川
  7. % 调制参数
  8. fm = 2500;              % 调制信号参数
  9. fc = 20000;             % 载波频率
  10. fs = 8*fc;              % 采样率
  11. total_time = 2;         % 仿真时长，单位：秒
  13. % 采样时间
  14. t = 0:1/fs:total_time-1/fs;
  16. % 调制信号为确知信号
  17. mut = sin(2*pi*fm*t)+cos(pi*fm*t);
  18. mlt = sin(3*pi*fm*t)+cos(4*pi*fm*t);
  20. % ISB 调制
  21. [ sig_isb_send ] = mod_isb(fc, fs, mut, mlt, t);
  23. % 加噪声
  24. snr = 50;               % 信噪比
  25. sig_isb_receive = awgn(sig_isb_send, snr, 'measured');
  27. % 数字正交解调
  28. phi0 = 0;
  29. [ sig_isbu_demod,sig_isbl_demod ] = demod_isb(sig_isb_receive, fc, fs, t, phi0);
  31. % 绘图
  32. nfft = length(sig_isb_receive);
  33. freq = (-nfft/2:nfft/2-1).'*(fs/nfft);
  34. figure;set(gcf,'color','w');
  35. plot_length = min(500, length(sig_isb_receive));
  36. subplot(1,2,1);
  37. plot(t(1:plot_length), sig_isb_receive(1:plot_length));xlim([t(1),t(plot_length)]);
  38. xlabel('t/s');ylabel('幅度');title('ISB接收信号');
  39. subplot(1,2,2);
  40. plot(freq, 10*log10(fftshift(abs(fft(sig_isb_receive,nfft)/nfft))+eps));xlim([freq(1),freq(end)]);
  41. xlabel('频率/hz');ylabel('幅度/dB');title('ISB接收信号双边幅度谱');
  43. figure;set(gcf,'color','w');
  44. plot(t(1:plot_length), mut(1:plot_length));xlim([t(1),t(plot_length)]);
  45. hold on;
  46. plot(t(1:plot_length), sig_isbu_demod(1:plot_length));xlim([t(1),t(plot_length)]);
  47. xlabel('t/s');ylabel('幅度');title('上边带解调效果');
  48. legend('上边带调制信号','上边带解调信号');
  50. figure;set(gcf,'color','w');
  51. plot(t(1:plot_length), mlt(1:plot_length));xlim([t(1),t(plot_length)]);
  52. hold on;
  53. plot(t(1:plot_length), sig_isbl_demod(1:plot_length));xlim([t(1),t(plot_length)]);
  54. xlabel('t/s');ylabel('幅度');title('下边带解调效果');
  55. legend('下边带调制信号','下边带解调信号');
  57. coefu = mean(abs(mut))/mean(abs(sig_isbu_demod));
  58. fprintf('norm(上边带调制信号 - %.2f * 上边带解调信号)/norm(上边带调制信号) = %.4f.\n', coefu, norm(mut-coefu*sig_isbu_demod)/norm(mut));
  60. coefl = mean(abs(mlt))/mean(abs(sig_isbl_demod));
  61. fprintf('norm(下边带调制信号 - %.2f * 下边带解调信号)/norm(下边带调制信号) = %.4f.\n', coefl, norm(mlt-coefl*sig_isbl_demod)/norm(mlt));
复制代码
免责声明：如果侵犯了您的权益，请联系站长，我们会及时删除侵权内容，谢谢合作！



欢迎光临 ToB企服应用市场:ToB评测及商务社交产业平台 (https://dis.qidao123.com/) Powered by Discuz! X3.4