【調製解調】ISB 獨立邊帶調幅

来源:https://www.cnblogs.com/young520/archive/2023/07/14/17552681.html
-Advertisement-
Play Games

學習數字信號處理演算法時整理的學習筆記。本篇介紹 ISB 獨立邊帶調幅信號的調製與解調，內附全套 MATLAB 代碼。 ...


說明

學習數字信號處理演算法時整理的學習筆記。同系列文章目錄可見 《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\)）。

Oh Shit!-圖片走丟了-打個廣告-歡迎來博客園關註“木三百川”

代碼詳見 mod_isb.m、main_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\)。

Oh Shit!-圖片走丟了-打個廣告-歡迎來博客園關註“木三百川”

下邊帶解調效果如下，計算誤差，有：\(\sqrt{\sum{{\lvert}m(t_i)-\hat{m}(t_i){\rvert}^2}}/\sqrt{\sum{{\lvert}m(t_i){\rvert}^2}}\approx0.0022\)。

Oh Shit!-圖片走丟了-打個廣告-歡迎來博客園關註“木三百川”

代碼詳見 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

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));

本文作者：木三百川

本文鏈接：https://www.cnblogs.com/young520/p/17552681.html

版權聲明：本文系博主原創文章，著作權歸作者所有。商業轉載請聯繫作者獲得授權，非商業轉載請附上出處鏈接。遵循 署名-非商業性使用-相同方式共用 4.0 國際版 (CC BY-NC-SA 4.0) 版權協議。


您的分享是我們最大的動力!

-Advertisement-
Play Games
更多相關文章
  • Java入門13（socket）
    ## Socket編程（網路通信） #### 伺服器端Demo（ServreSocket） ​ 創建服務端時，如果不提供IP地址，則預設為本地連接（127.0.0.1），但是一定需要手動配置監聽埠！ ```java public static void main(String[] args) { ...
  • C++將WSAGetLastError轉換成字元串信息
    #include <iostream> #include <Windows.h> #include <WinSock2.h> std::string GetLastErrorMessage() { DWORD errorCode = WSAGetLastError(); LPSTR errorMes ...
  • 爬蟲突破驗證碼技術 - 2Captcha
    在互聯網世界中，驗證碼作為一種防止機器人訪問的工具，是爬蟲最常遇到的阻礙。驗證碼的類型眾多，從簡單的數字、字母驗證碼，到複雜的圖像識別驗證碼，再到更為高級的互動式驗證碼，每一種都有其獨特的識別方法和應對策略。在這篇文章中，我們將一一介紹各種驗證碼的工作原理和使用[2Captcha](https://... ...
  • Excel巨集教程_編程入門自學教程_菜鳥教程-免費教程分享
    ## 教程簡介 巨集語言Visual Basic for Application(VBA).Visual Basic是windows環境下開發應用軟體的一種通用程式設計語言，功能強大，簡便易用。 Excel巨集是Excel中的一種編程功能，它可以讓用戶錄製一系列的操作，以便在需要時自動執行這些操作。用戶 ...
  • jmeter教程_編程入門自學教程_菜鳥教程-免費教程分享
    ## 教程簡介 Apache JMeter 是 Apache 組織基於 Java 開發的壓力測試工具，用於對軟體做壓力測試。JMeter 最初被設計用於 Web 應用測試，但後來擴展到了其他測試領域，可用於測試靜態和動態資源，如靜態文件、Java 小服務程式、CGI 腳本、Java 對象、資料庫和 ...
  • Miller_rabin 素數測試 學習筆記
    ## Miller_rabin 素數測試 一種用來判斷素數的演算法。 ### 前置芝士 #### 威爾遜定理 若 $p$ 為素數，$(p-1)! \equiv -1 (\mod p)$。 證明： 充分性證明： 如果 $p$ 不是素數，那麼他的因數必定存在於$ 1,2,3,\dots,p−1$ 之中，所 ...
  • 網路流學習筆記
    網路流 何為網路流 想要弄清楚網路流，首先要知道網路的概念，通常在運籌學中，網路是指一個有向圖$G\ =\ (V,E)$ 。其每條邊$(u,v)\in E$都有一個權值$c(u,v)$,稱為這條邊的流量(Capacity)，還有兩個特殊的點，一個是源點(Source),一個是匯點(Sink)在圖論中 ...
  • 【pandas小技巧】--創建測試數據
    學習`pandas`的過程中，為了嘗試`pandas`提供的各類功能強大的函數，常常需要花費很多時間去創造測試數據。 在`pandas`中，快速創建測試數據可以更快的評估 `pandas` 函數。通過生成一組測試數據，可以評估例如 `read_csv`、`read_excel`、`groupby`等 ...
一周排行
    -Advertisement-
    Play Games
  • 移動開發（一）：使用.NET MAUI開發第一個安卓APP
    移動開發（一）：使用.NET MAUI開發第一個安卓APP 對於工作多年的C#程式員來說，近來想嘗試開發一款安卓APP，考慮了很久最終選擇使用.NET MAUI這個微軟官方的框架來嘗試體驗開發安卓APP，畢竟是使用Visual Studio開發工具，使用起來也比較的順手，結合微軟官方的教程進行了安卓 ...
  • wpf ToggleButton選中效果和一個登錄界面
    前言 QuestPDF 是一個開源 .NET 庫，用於生成 PDF 文檔。使用了C# Fluent API方式可簡化開發、減少錯誤並提高工作效率。利用它可以輕鬆生成 PDF 報告、發票、導出文件等。 項目介紹 QuestPDF 是一個革命性的開源 .NET 庫，它徹底改變了我們生成 PDF 文檔的方 ...
  • 彈幕樹洞項目功能新增篇
    項目地址 項目後端地址： https://github.com/ZyPLJ/ZYTteeHole 項目前端頁面地址： ZyPLJ/TreeHoleVue (github.com) https://github.com/ZyPLJ/TreeHoleVue 目前項目測試訪問地址： http://tree ...
  • 第27篇 sqlserver2022詳細安裝步驟
    話不多說，直接開乾 一.下載 1.官方鏈接下載： https://www.microsoft.com/zh-cn/sql-server/sql-server-downloads 2.在下載目錄中找到下麵這個小的安裝包 SQL2022-SSEI-Dev.exe，運行開始下載SQL server； 二. ...
  • .NET 開源高性能 MQTT 類庫
    前言 隨著物聯網（IoT）技術的迅猛發展，MQTT（消息隊列遙測傳輸）協議憑藉其輕量級和高效性，已成為眾多物聯網應用的首選通信標準。 MQTTnet 作為一個高性能的 .NET 開源庫，為 .NET 平臺上的 MQTT 客戶端與伺服器開發提供了強大的支持。 本文將全面介紹 MQTTnet 的核心功能 ...
  • Serilog文檔翻譯系列（六） - 可用的接收器、增強器、格式化輸出
    Serilog支持多種接收器用於日誌存儲，增強器用於添加屬性，LogContext管理動態屬性，支持多種輸出格式包括純文本、JSON及ExpressionTemplate。還提供了自定義格式化選項，適用於不同需求。 ...
  • 警惕 Visual Studio 屬性求值副作用導致邏輯不符合預期
    目錄簡介獲取 HTML 文檔解析 HTML 文檔測試參考文章 簡介 動態內容網站使用 JavaScript 腳本動態檢索和渲染數據，爬取信息時需要模擬瀏覽器行為，否則獲取到的源碼基本是空的。 本文使用的爬取步驟如下： 使用 Selenium 獲取渲染後的 HTML 文檔 使用 HtmlAgility ...
  • [使用目前最新版]HybridCLR6.9.0+YooAsset2.2.4實現純C# Unity熱更新方案 (一)
    1.前言 什麼是熱更新 游戲或者軟體更新時，無需重新下載客戶端進行安裝，而是在應用程式啟動的情況下，在內部進行資源或者代碼更新 Unity目前常用熱更新解決方案 HybridCLR，Xlua，ILRuntime等 Unity目前常用資源管理解決方案 AssetBundles，Addressable， ...
  • 在 ASP.NET Core Web API 中使用操作篩選器統一處理通用操作
    本文章主要是在C# ASP.NET Core Web API框架實現向手機發送驗證碼簡訊功能。這裡我選擇是一個互億無線簡訊驗證碼平臺，其實像阿裡雲，騰訊雲上面也可以。 首先我們先去 互億無線 https://www.ihuyi.com/api/sms.html 去註冊一個賬號 註冊完成賬號後，它會送 ...
  • 第28篇 如何.net中實現高效可靠數據同步api
    通過以下方式可以高效，並保證數據同步的可靠性 1.API設計 使用RESTful設計，確保API端點明確，並使用適當的HTTP方法（如POST用於創建，PUT用於更新）。 設計清晰的請求和響應模型，以確保客戶端能夠理解預期格式。 2.數據驗證 在伺服器端進行嚴格的數據驗證，確保接收到的數據符合預期格 ...
所有分類