主頁 > 後端開發 > 【調制解調】DSB 雙邊帶調幅

【調制解調】DSB 雙邊帶調幅

2023-07-11 08:21:46 後端開發

說明

學習數字信號處理演算法時整理的學習筆記，同系列文章目錄可見 《DSP 學習之路》目錄，代碼已上傳到 Github - ModulationAndDemodulation，本篇介紹 DSB 雙邊帶調幅信號的調制與解調，內附全套 MATLAB 代碼，

目錄
  • 說明
  • 1. DSB 調制演算法
    • 1.1 演算法描述
    • 1.2 DSB 信號調制示例
  • 2. DSB 解調演算法
    • 2.1 插入載波包絡檢波法
    • 2.2 相干解調（同步檢測）
    • 2.3 數字正交解調
  • 3. DSB 仿真（MATLAB Communications Toolbox）
  • 參考資料
  • 附錄代碼
    • 附.1 檔案 mod_dsb.m
    • 附.2 檔案 main_modDSB_example.m
    • 附.3 檔案 demod_dsb_method1.m
    • 附.4 檔案 main_demodDSB_example1.m
    • 附.5 檔案 lpf_filter.m
    • 附.6 檔案 demod_dsb_method2.m
    • 附.7 檔案 main_demodDSB_example2.m
    • 附.8 檔案 demod_dsb_method3.m
    • 附.9 檔案 main_demodDSB_example3.m
    • 附.10 檔案 main_CommDSB_example.m


1. DSB 調制演算法

1.1 演算法描述

在 AM 調幅信號中，載波分量并不攜帶資訊，資訊完全由邊帶傳送，如果在 AM 調制模型中將直流 \(A_0\) 去掉，即可得到一種高調制效率的調制方式——抑制載波雙邊帶信號（DSB - SC, Double Side Band with Suppressed Carrier），簡稱雙邊帶信號（DSB），其時域運算式為：

\[s_{DSB}(t)=m(t)cos{\omega_ct} \tag{1} \]

式中：\(m(t)\) 是調制信號（攜帶要發出去的資訊），它可以是確知信號，也可以是隨機信號，其均值通常為 0；\(cos{\omega_ct}\) 是載波，\(\omega_c\) 是載波角頻率，與載波頻率 \(f_c\) 之間的關系為 \(\omega_c=2{\pi}f_c\)，DSB 的頻譜與 AM 頻譜相近，只是沒有了在 \(\pm\omega_c\) 處的 \(\delta\) 函式，對式 \((1)\) 進行傅里葉變換，得到 DSB 信號的頻譜（幅度譜）運算式：

\[S_{DSB}(\omega)=\frac{1}{2}\left[M(\omega+\omega_c)+M(\omega-\omega_c)\right] \tag{2} \]

式中，\(M(\omega)\) 是調制信號 \(m(t)\) 的頻譜，DSB 信號的特性如下：

  • DSB 信號的頻譜由上邊帶與下邊帶兩部分組成，不存在載波分量，它的帶寬仍是基帶信號（調制信號）帶寬 \(f_H\) 的 2 倍，即 \(B_{DSB}=2f_{H}\)，與 AM 信號帶寬相同，

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

  • 由于不存在載波分量，有用功率 \(P_s\) 就是信號總功率 \(P_{DSB}\)，即 \(P_s=P_{DSB}\)，全部功率都用于資訊傳輸，調制效率 \({\eta_{DSB}}=100\%\)，

1.2 DSB 信號調制示例

調制信號 \(m(t)\) 可以是確知信號，也可以是隨機信號，當 \(m(t)\) 是確知信號時，不妨假設 \(m(t)\) 的時域運算式如下：

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

各調制引數取值：\(f_m=2500Hz\)，\(f_c=20000Hz\)，信號采樣率 \(f_s=8{f_c}\)，仿真總時長為 \(2s\)，DSB 調制效果如下圖所示（為了美觀，時域只顯示前 500 個點），調制信號 \(m(t)\) 雙邊幅度譜有四根離散譜線（\({\pm}2500Hz\)、\({\pm}1250Hz\)），載波 \(c(t)\) 的雙邊幅度譜有兩根離散譜線（\({\pm}20000Hz\)），DSB 信號有八根離散譜線（\(\pm17500Hz\)、\(\pm18750Hz\)、\(\pm21250Hz\)、\(\pm22500Hz\)），代碼詳見附錄 main_modDSB_example.m 與 mod_dsb.m，

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


2. DSB 解調演算法

解調是調制的逆程序，其作用是從接收的已調信號中恢復原基帶信號（即調制信號），DSB 信號的包絡不再與調制信號 \(m(t)\) 的變化規律一致，因而不能采用簡單的包絡檢波來恢復調制信號，通常采用相干解調的方法來進行解調，另一種方法是，插入很強的載波，使其成為或近似為 AM 信號，則可利用包絡檢波器恢復調制信號，這種方法被稱為插入載波包絡檢波法，為了保證檢波質量，插入的載波振幅應遠大于信號的振幅，同時也要求插入的載波與調制載波同頻同相，下面介紹三種解調方法并對 1.2 節中的 DSB 信號進行解調，

2.1 插入載波包絡檢波法

插入幅值為 \(A_0\) 的載波，得到：

\[s_{DSB}(t)+{A_0}cos{\omega_ct}=\left[A_0+m(t)\right]cos{\omega_ct} \tag{4} \]

其中 \(A_0 \geq {\lvert}{m(t)}{\rvert}_{max}\)，這樣就得到了一個 AM 信號，使用 AM 解調器進行解調即可，步驟如下：

  1. 第一步：加上載波 \({A_0}cos{\omega_ct}\)，其中 \(A_0 \geq {\lvert}{m(t)}{\rvert}_{max}\)，獲得 AM 信號，
  2. 第二步：使用 AM 解調器進行解調，

對 1.2 節中的 DSB 信號，設定信噪比 \(SNR=50dB\)，解調效果如下，計算誤差，有：\(\sqrt{\sum{{\lvert}m(t_i)-\hat{m}(t_i){\rvert}^2}}/\sqrt{\sum{{\lvert}m(t_i){\rvert}^2}}\approx0.0022\)，更改插入載波的初始相位為 \({\phi_0}=\pi/4,\pi/2\)，或者更改插入載波的中心頻率為 \(0.8f_c,1.2f_c\) 后，解調效果變差，說明這種方法對插入載波同頻同相的要求較高，

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

代碼詳見 demod_dsb_method1.m 和 main_demodDSB_example1.m，AM 解調器詳見本人同系列博客 【調制解調】AM 調幅，

2.2 相干解調（同步檢測）

將 DSB 信號與同頻同相的相干載波相乘，得到：

\[\begin{aligned} s_{DSB}(t){\cdot}cos{(\omega_ct)}&=m(t)cos{(\omega_ct)}{\cdot}cos{(\omega_ct)}\\[1em] &=\frac{1}{2}m(t)+\frac{1}{2}m(t)cos(2\omega_ct) \end{aligned} \tag{5} \]

然后通過一個低通濾波器即可獲得解調結果，步驟如下：

  1. 第一步：乘以相干載波（即乘以 \(2cos({\omega_ct}+{\phi_0})\)，前面的 2 被用來做幅度補償，
  2. 第二步：低通濾波器濾除高頻載波，濾除 \(2{\omega}_c\)，

對 1.2 節中的 DSB 信號，設定信噪比 \(SNR=50dB\)，解調效果如下，計算誤差，有：\(\sqrt{\sum{{\lvert}m(t_i)-\hat{m}(t_i){\rvert}^2}}/\sqrt{\sum{{\lvert}m(t_i){\rvert}^2}}\approx0.0016\)，更改相干載波的初始相位為 \({\phi_0}=\pi/4,\pi/2\) 后，解調幅值發生失真，當與真實相位相差 \(\pi/2\) 時幅值失真最大；但更改相干載波的中心頻率為 \(0.8f_c,1.2f_c\) 后，解調效果變得很差，波形完全失真，說明這種方法對相干載波同頻同相的要求也較高，

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

代碼詳見 lpf_filter.m、demod_dsb_method2.m 和 main_demodDSB_example2.m，

2.3 數字正交解調

DSB 數字正交解調一般有以下兩個步驟，它與相干解調（同步檢測）法是等效的：

  1. 第一步：乘以正交相干載波得到 \({s_I}(t)\) 與 \({s_Q}(t)\)，即 \({s_I}(t)=2s(t)cos({\omega_ct}+{\phi_0})\)，\({s_Q}(t)=-2s(t)sin({\omega_ct}+{\phi_0})\)，前面的 2 被用來做幅度補償，
  2. 第二步：低通濾波器濾除 \({s_I}(t)\) 與 \({s_Q}(t)\) 中的高頻分量，所得的 \(s_I(t)\) 即為解調結果，

對 1.2 節中的 DSB 信號，設定信噪比 \(SNR=50dB\)，解調效果如下，計算誤差，有：\(\sqrt{\sum{{\lvert}m(t_i)-\hat{m}(t_i){\rvert}^2}}/\sqrt{\sum{{\lvert}m(t_i){\rvert}^2}}\approx0.0016\)，與相干解調（同步檢測）一樣，這種方法對相干載波同頻同相的要求較高，

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

代碼詳見 lpf_filter.m、demod_dsb_method3.m 和 main_demodDSB_example3.m，


3. DSB 仿真（MATLAB Communications Toolbox）

MATLAB 的 Communications Toolbox 中提供了 AM 調制函式 ammod，高斯白噪聲函式 awgn，以及 AM 解調函式 amdemod，可以很方便地完成 DSB 信號仿真，設定 ammod 與 amdemod 的輸入引數 carramp = 0 即為 DSB 的調制與解調（carramp 引數的默認值就是 0，不顯式設定這個引數也可以），使用這三個函式實作上面 1.2 節中確知信號 \(m(t)\) 的 DSB 調制解調，調制后加噪聲的效果如下：

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

解調效果如下：

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

解調信號與調制信號波形基本重回，計算誤差，有：\(\sqrt{\sum{{\lvert}m(t_i)-\hat{m}(t_i){\rvert}^2}}/\sqrt{\sum{{\lvert}m(t_i){\rvert}^2}}\approx0.0025\)，代碼詳見附錄 main_CommDSB_example.m，


參考資料

[1] 樓才義,徐建良,楊小牛.軟體無線電原理與應用[M].電子工業出版社,2014.

[2] 樊昌信,曹麗娜.通信原理.第7版[M].國防工業出版社,2012.

[3] CSDN - 通信原理之模擬幅度調制（線性調制）詳解，


附錄代碼

附.1 檔案 mod_dsb.m

function [ sig_dsb ] = mod_dsb(fc, fs, mt, t)
% MOD_DSB        DSB 雙邊帶調幅
% 輸入引數：
%       fc      載波中心頻率
%       fs      信號采樣率
%       mt      調制信號
%       t       采樣時間
% 輸出引數：
%       sig_dsb DSB 雙邊帶調幅實信號
% @author 木三百川

% 生成信號
ct = cos(2*pi*fc*t);  
sig_dsb = mt.*ct;   % DSB 雙邊帶調幅信號

% 繪圖
nfft = length(sig_dsb);
freq = (-nfft/2:nfft/2-1).'*(fs/nfft);
figure;set(gcf,'color','w');
plot_length = min(500, length(sig_dsb));
subplot(3,2,1);
plot(t(1:plot_length), mt(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('調制信號m(t)');
subplot(3,2,2);
plot(freq, 10*log10(fftshift(abs(fft(mt,nfft)/nfft))+eps));xlim([freq(1),freq(end)]);
xlabel('頻率/hz');ylabel('幅度/dB');title('調制信號m(t)雙邊幅度譜');

subplot(3,2,3);
plot(t(1:plot_length), ct(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('載波c(t)');
subplot(3,2,4);
plot(freq, 10*log10(fftshift(abs(fft(ct,nfft)/nfft))+eps));xlim([freq(1),freq(end)]);
xlabel('頻率/hz');ylabel('幅度/dB');title('載波c(t)雙邊幅度譜');

subplot(3,2,5);
plot(t(1:plot_length), sig_dsb(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('DSB雙邊帶調幅信號s(t)');
subplot(3,2,6);
plot(freq, 10*log10(fftshift(abs(fft(sig_dsb,nfft)/nfft))+eps));xlim([freq(1),freq(end)]);
xlabel('頻率/hz');ylabel('幅度/dB');title('DSB雙邊帶調幅信號s(t)雙邊幅度譜');

end

附.2 檔案 main_modDSB_example.m

clc;
clear;
close all;
% DSB 調制仿真(調制信號為確知信號)
% @author 木三百川

% 調制引數
fm = 2500;              % 調制信號引數
fc = 20000;             % 載波頻率
fs = 8*fc;              % 采樣率
total_time = 2;         % 仿真時長，單位：秒

% 采樣時間
t = 0:1/fs:total_time-1/fs;

% 調制信號為確知信號
mt = sin(2*pi*fm*t)+cos(pi*fm*t);

% DSB 調制
[ sig_dsb ] = mod_dsb(fc, fs, mt, t);

附.3 檔案 demod_dsb_method1.m

function [ sig_dsb_demod ] = demod_dsb_method1(sig_dsb_receive, fc, fs, t, phi0)
% DEMOD_DSB_METHOD1        DSB 插入載波包絡檢波法
% 輸入引數：
%       sig_dsb_receive     DSB 接收信號，行向量
%       fc                  載波中心頻率
%       fs                  信號采樣率
%       t                   采樣時間
%       phi0                載波初始相位
% 輸出引數：
%       sig_dsb_demod       解調結果，與 sig_dsb_receive 等長
% @author 木三百川

% 第一步：插入載波
A0 = max(abs(sig_dsb_receive))/0.8;
sig_dsb2am = sig_dsb_receive + A0*cos(2*pi*fc*t+phi0);

% 第二步：使用 AM 解調器進行解調
[ sig_dsb_demod ] = demod_am_method4(sig_dsb2am, fs, t);

end

附.4 檔案 main_demodDSB_example1.m

clc;
clear;
close all;
% DSB 解調仿真(調制信號為確知信號，插入載波包絡檢波法)
% @author 木三百川

% 調制引數
fm = 2500;              % 調制信號引數
fc = 20000;             % 載波頻率
fs = 8*fc;              % 采樣率
total_time = 2;         % 仿真時長，單位：秒

% 采樣時間
t = 0:1/fs:total_time-1/fs;

% 調制信號為確知信號
mt = sin(2*pi*fm*t)+cos(pi*fm*t);

% DSB 調制
[ sig_dsb_send ] = mod_dsb(fc, fs, mt, t);

% 加噪聲
snr = 50;               % 信噪比
sig_dsb_receive = awgn(sig_dsb_send, snr, 'measured');

% 插入載波包絡檢波法
phi0 = 0;
[ sig_dsb_demod ] = demod_dsb_method1(sig_dsb_receive, fc, fs, t, phi0);

% 繪圖
nfft = length(sig_dsb_receive);
freq = (-nfft/2:nfft/2-1).'*(fs/nfft);
figure;set(gcf,'color','w');
plot_length = min(500, length(sig_dsb_receive));
subplot(1,2,1);
plot(t(1:plot_length), sig_dsb_receive(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('DSB接收信號');
subplot(1,2,2);
plot(freq, 10*log10(fftshift(abs(fft(sig_dsb_receive,nfft)/nfft))+eps));xlim([freq(1),freq(end)]);
xlabel('頻率/hz');ylabel('幅度/dB');title('DSB接收信號雙邊幅度譜');

figure;set(gcf,'color','w');
plot(t(1:plot_length), mt(1:plot_length));xlim([t(1),t(plot_length)]);
hold on;
plot(t(1:plot_length), sig_dsb_demod(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('解調效果');
legend('調制信號','解調信號');

coef = mean(abs(mt))/mean(abs(sig_dsb_demod));
fprintf('norm(調制信號 - %.2f * 解調信號)/norm(調制信號) = %.4f.\n', coef, norm(mt-coef*sig_dsb_demod)/norm(mt));

附.5 檔案 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

附.6 檔案 demod_dsb_method2.m

function [ sig_dsb_demod ] = demod_dsb_method2(sig_dsb_receive, fc, fs, t, phi0)
% DEMOD_DSB_METHOD2        DSB 相干解調（同步檢測）
% 輸入引數：
%       sig_dsb_receive     DSB 接收信號，行向量
%       fc                  載波中心頻率
%       fs                  信號采樣率
%       t                   采樣時間
%       phi0                載波初始相位
% 輸出引數：
%       sig_dsb_demod       解調結果，與 sig_dsb_receive 等長
% @author 木三百川

% 第一步：乘以相干載波
sig_dsbct = 2*sig_dsb_receive.*cos(2*pi*fc*t+phi0);

% 第二步：低通濾波
sig_dsb_demod = lpf_filter(sig_dsbct, fc/(fs/2));

end

附.7 檔案 main_demodDSB_example2.m

clc;
clear;
close all;
% DSB 解調仿真(調制信號為確知信號，相干解調（同步檢測）)
% @author 木三百川

% 調制引數
fm = 2500;              % 調制信號引數
fc = 20000;             % 載波頻率
fs = 8*fc;              % 采樣率
total_time = 2;         % 仿真時長，單位：秒

% 采樣時間
t = 0:1/fs:total_time-1/fs;

% 調制信號為確知信號
mt = sin(2*pi*fm*t)+cos(pi*fm*t);

% DSB 調制
[ sig_dsb_send ] = mod_dsb(fc, fs, mt, t);

% 加噪聲
snr = 50;               % 信噪比
sig_dsb_receive = awgn(sig_dsb_send, snr, 'measured');

% 相干解調（同步檢測）
phi0 = 0;
[ sig_dsb_demod ] = demod_dsb_method2(sig_dsb_receive, fc, fs, t, phi0);

% 繪圖
nfft = length(sig_dsb_receive);
freq = (-nfft/2:nfft/2-1).'*(fs/nfft);
figure;set(gcf,'color','w');
plot_length = min(500, length(sig_dsb_receive));
subplot(1,2,1);
plot(t(1:plot_length), sig_dsb_receive(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('DSB接收信號');
subplot(1,2,2);
plot(freq, 10*log10(fftshift(abs(fft(sig_dsb_receive,nfft)/nfft))+eps));xlim([freq(1),freq(end)]);
xlabel('頻率/hz');ylabel('幅度/dB');title('DSB接收信號雙邊幅度譜');

figure;set(gcf,'color','w');
plot(t(1:plot_length), mt(1:plot_length));xlim([t(1),t(plot_length)]);
hold on;
plot(t(1:plot_length), sig_dsb_demod(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('解調效果');
legend('調制信號','解調信號');

coef = mean(abs(mt))/mean(abs(sig_dsb_demod));
fprintf('norm(調制信號 - %.2f * 解調信號)/norm(調制信號) = %.4f.\n', coef, norm(mt-coef*sig_dsb_demod)/norm(mt));

附.8 檔案 demod_dsb_method3.m

function [ sig_dsb_demod ] = demod_dsb_method3(sig_dsb_receive, fc, fs, t, phi0)
% DEMOD_DSB_METHOD3        DSB 數字正交解調，與相干解調（同步檢測）是等效的
% 輸入引數：
%       sig_dsb_receive     DSB 接收信號，行向量
%       fc                  載波中心頻率
%       fs                  信號采樣率
%       t                   采樣時間
%       phi0                載波初始相位
% 輸出引數：
%       sig_dsb_demod       解調結果，與 sig_dsb_receive 等長
% @author 木三百川

% 第一步：乘以正交相干載波
sig_dsb_i = 2*sig_dsb_receive.*cos(2*pi*fc*t+phi0);
sig_dsb_q = -2*sig_dsb_receive.*sin(2*pi*fc*t+phi0);

% 第二步：低通濾波
sig_dsb_i_lpf = lpf_filter(sig_dsb_i, fc/(fs/2));
sig_dsb_q_lpf = lpf_filter(sig_dsb_q, fc/(fs/2));
sig_dsb_demod = sig_dsb_i_lpf;

end

附.9 檔案 main_demodDSB_example3.m

clc;
clear;
close all;
% DSB 解調仿真(調制信號為確知信號，數字正交解調)
% @author 木三百川

% 調制引數
fm = 2500;              % 調制信號引數
fc = 20000;             % 載波頻率
fs = 8*fc;              % 采樣率
total_time = 2;         % 仿真時長，單位：秒

% 采樣時間
t = 0:1/fs:total_time-1/fs;

% 調制信號為確知信號
mt = sin(2*pi*fm*t)+cos(pi*fm*t);

% DSB 調制
[ sig_dsb_send ] = mod_dsb(fc, fs, mt, t);

% 加噪聲
snr = 50;               % 信噪比
sig_dsb_receive = awgn(sig_dsb_send, snr, 'measured');

% 數字正交解調
phi0 = 0;
[ sig_dsb_demod ] = demod_dsb_method3(sig_dsb_receive, fc, fs, t, phi0);

% 繪圖
nfft = length(sig_dsb_receive);
freq = (-nfft/2:nfft/2-1).'*(fs/nfft);
figure;set(gcf,'color','w');
plot_length = min(500, length(sig_dsb_receive));
subplot(1,2,1);
plot(t(1:plot_length), sig_dsb_receive(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('DSB接收信號');
subplot(1,2,2);
plot(freq, 10*log10(fftshift(abs(fft(sig_dsb_receive,nfft)/nfft))+eps));xlim([freq(1),freq(end)]);
xlabel('頻率/hz');ylabel('幅度/dB');title('DSB接收信號雙邊幅度譜');

figure;set(gcf,'color','w');
plot(t(1:plot_length), mt(1:plot_length));xlim([t(1),t(plot_length)]);
hold on;
plot(t(1:plot_length), sig_dsb_demod(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('解調效果');
legend('調制信號','解調信號');

coef = mean(abs(mt))/mean(abs(sig_dsb_demod));
fprintf('norm(調制信號 - %.2f * 解調信號)/norm(調制信號) = %.4f.\n', coef, norm(mt-coef*sig_dsb_demod)/norm(mt));

附.10 檔案 main_CommDSB_example.m

clc;
clear;
close all;
% DSB 調制解調仿真(使用Communications Toolbox工具箱)
% @author 木三百川

% 調制引數
fm = 2500;              % 調制信號引數
fc = 20000;             % 載波頻率
fs = 8*fc;              % 采樣率
total_time = 2;         % 仿真時長，單位：秒

% 采樣時間
t = 0:1/fs:total_time-1/fs;

% 調制信號為確知信號
mt = sin(2*pi*fm*t)+cos(pi*fm*t);

% DSB 調制
ini_phase = 0;
sig_dsb_send = ammod(mt, fc, fs, ini_phase);

% 加噪聲
snr = 50;               % 信噪比
sig_dsb_receive = awgn(sig_dsb_send, snr, 'measured');

% DSB 解調
[ sig_dsb_demod ] = amdemod(sig_dsb_receive, fc, fs, ini_phase);

% 繪圖
nfft = length(sig_dsb_receive);
freq = (-nfft/2:nfft/2-1).'*(fs/nfft);
figure;set(gcf,'color','w');
plot_length = min(500, length(sig_dsb_receive));
subplot(1,2,1);
plot(t(1:plot_length), sig_dsb_receive(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('DSB接收信號');
subplot(1,2,2);
plot(freq, 10*log10(fftshift(abs(fft(sig_dsb_receive,nfft)/nfft))+eps));xlim([freq(1),freq(end)]);
xlabel('頻率/hz');ylabel('幅度/dB');title('DSB接收信號雙邊幅度譜');

figure;set(gcf,'color','w');
plot(t(1:plot_length), mt(1:plot_length));xlim([t(1),t(plot_length)]);
hold on;
plot(t(1:plot_length), sig_dsb_demod(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('解調效果');
legend('調制信號','解調信號');

coef = mean(abs(mt))/mean(abs(sig_dsb_demod));
fprintf('norm(調制信號 - %.2f * 解調信號)/norm(調制信號) = %.4f.\n', coef, norm(mt-coef*sig_dsb_demod)/norm(mt));

本文作者：木三百川

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

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

轉載請註明出處，本文鏈接：https://www.uj5u.com/houduan/556998.html

標籤：其他

上一篇：為什么使用ioutil.ReadAll 函式需要注意

下一篇：返回列表

標籤雲
其他(162370) Python(38273) JavaScript(25530) Java(18294) C(15239) 區塊鏈(8275) C＃(7972) AI(7469) 爪哇(7425) MySQL(7294) html(6777) 基礎類(6313) sql(6102) 熊猫(6058) PHP(5876) 数组(5741) R(5409) Linux(5347) 反应(5209) 腳本語言(PerlPython)(5129) 非技術區(4971) Android(4615) 数据框(4311) css(4259) 节点.js(4032) C語言(3288) json(3245) 列表(3129) 扑(3119) C++語言(3117) 安卓(2998) 打字稿(2995) VBA(2789) Java相關(2746) 疑難問題(2699) 细绳(2522) 單片機工控(2479) iOS(2438) ASP.NET(2404) MongoDB(2323) 麻木的(2285) 正则表达式(2254) 字典(2211) 循环(2198) 迅速(2185) 擅长(2169) 镖(2155) HtmlCss(1995) .NET技术(1986) 功能(1967) Web開發(1951) C++(1942) python-3.x(1918) 弹簧靴(1913) xml(1889) PostgreSQL(1882) .NETCore(1863) 谷歌表格(1846) Unity3D(1843) for循环(1842)

熱門瀏覽
  • 【C++】Microsoft C++、C 和匯編程式檔案

    ......

    uj5u.com 2020-09-10 00:57:23 more
  • 例外宣告

    相比于斷言適用于排除邏輯上不可能存在的狀態，例外通常是用于邏輯上可能發生的錯誤。 例外宣告 Item 1：當函式不可能拋出例外或不能接受拋出例外時，使用noexcept 理由 如果不打算拋出例外的話，程式就會認為無法處理這種錯誤，并且應當盡早終止，如此可以有效地阻止例外的傳播與擴散。 示例 //不可 ......

    uj5u.com 2020-09-10 00:57:27 more
  • Codeforces 1400E Clear the Multiset(貪心 + 分治)

    鏈接：https://codeforces.com/problemset/problem/1400/E 來源：Codeforces 思路：給你一個陣列，現在你可以進行兩種操作，操作1：將一段沒有 0 的區間進行減一的操作，操作2：將 i 位置上的元素歸零。最終問：將這個陣列的全部元素歸零后操作的最少 ......

    uj5u.com 2020-09-10 00:57:30 more
  • UVA11610 【Reverse Prime】

    本人看到此題沒有翻譯，就附帶了一個自己的翻譯版本 思考 這一題，它的第一個要求是找出所有 $7$ 位反向質數及其質因數的個數。 我們應該需要質數篩篩選1~$10^{7}$的所有數，這里就不慢慢介紹了。但是，重讀題，我們突然發現反向質數都是 $7$ 位，而將它反過來后的數字卻是 $6$ 位數，這就說明 ......

    uj5u.com 2020-09-10 00:57:36 more
  • 統計區間素數數量

    1 #pragma GCC optimize(2) 2 #include <bits/stdc++.h> 3 using namespace std; 4 bool isprime[1000000010]; 5 vector<int> prime; 6 inline int getlist(int ......

    uj5u.com 2020-09-10 00:57:47 more
  • C/C++編程筆記：C++中的 const 變數詳解，教你正確認識const用法

    1、C中的const 1、區域const變數存放在堆疊區中，會分配記憶體（也就是說可以通過地址間接修改變數的值）。測驗代碼如下： 運行結果： 2、全域const變數存放在只讀資料段（不能通過地址修改，會發生寫入錯誤）， 默認為外部聯編，可以給其他源檔案使用（需要用extern關鍵字修飾） 運行結果： ......

    uj5u.com 2020-09-10 00:58:04 more
  • 【C++犯錯記錄】VS2019 MFC添加資源不懂如何修改資源宏ID

    1. 首先在資源視圖中，添加資源 2. 點擊新添加的資源，復制自動生成的ID 3. 在解決方案資源管理器中找到Resource.h檔案，編輯，使用整個專案搜索和替換的方式快速替換 宏宣告 4. Ctrl+Shift+F 全域搜索，點擊查找全部，然后逐個替換 5. 為什么使用搜索替換而不使用屬性視窗直 ......

    uj5u.com 2020-09-10 00:59:11 more
  • 【C++犯錯記錄】VS2019 MFC不懂的批量添加資源

    1. 打開資源頭檔案Resource.h，在其中預先定義好宏 ID（不清楚其實ID值應該設定多少，可以先新建一個相同的資源項，再在這個資源的ID值的基礎上遞增即可） 2. 在資源視圖中選中專案資源，按F7編輯資源檔案，按 ID 型別 相對路徑的形式添加 資源。（別忘了先把檔案拷貝到專案中的res檔案 ......

    uj5u.com 2020-09-10 01:00:19 more
  • C/C++編程筆記：關于C++的參考型別，專供新手入門使用

    今天要講的是C++中我最喜歡的一個用法——參考，也叫別名。 參考就是給一個變數名取一個變數名，方便我們間接地使用這個變數。我們可以給一個變數創建N個參考，這N + 1個變數共享了同一塊記憶體區域。(參考型別的變數會占用記憶體空間，占用的記憶體空間的大小和指標型別的大小是相同的。雖然參考是一個物件的別名，但 ......

    uj5u.com 2020-09-10 01:00:22 more
  • 【C/C++編程筆記】從頭開始學習C ++：初學者完整指南

    眾所周知，C ++的學習曲線陡峭，但是花時間學習這種語言將為您的職業帶來奇跡，并使您與其他開發人員區分開。您會更輕松地學習新語言，形成真正的解決問題的技能，并在編程的基礎上打下堅實的基礎。 C ++將幫助您養成良好的編程習慣（即清晰一致的編碼風格，在撰寫代碼時注釋代碼，并限制類內部的可見性），并且由 ......

    uj5u.com 2020-09-10 01:00:41 more
最新发布
  • 【調制解調】DSB 雙邊帶調幅

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

    uj5u.com 2023-07-11 08:21:46 more
  • 為什么使用ioutil.ReadAll 函式需要注意

    # 1. 引言 當我們需要將資料一次性加載到記憶體中，`ioutil.ReadAll` 函式是一個方便的選擇，但是`ioutil.ReadAll` 的使用是需要注意的。 在這篇文章中，我們將首先對`ioutil.ReadAll`函式進行基本介紹，之后會介紹其存在的問題，以及引起該問題的原因，最后給出了 ......

    uj5u.com 2023-07-11 08:21:41 more
  • 為什么使用ioutil.ReadAll 函式需要注意

    # 1. 引言 當我們需要將資料一次性加載到記憶體中，`ioutil.ReadAll` 函式是一個方便的選擇，但是`ioutil.ReadAll` 的使用是需要注意的。 在這篇文章中，我們將首先對`ioutil.ReadAll`函式進行基本介紹，之后會介紹其存在的問題，以及引起該問題的原因，最后給出了 ......

    uj5u.com 2023-07-11 08:20:28 more
  • Rust 使用egui創建一個簡單的下載器demo

    倉庫連接: https://github.com/GaN601/egui-demo-download-util 這是我第一個rust gui demo, 學習rust有挺長時間了, 但是一直沒有落實到實踐中, 本著對桌面應用的興趣, 考察了slint、egui兩種框架, 最后還是選擇了egui. 這 ......

    uj5u.com 2023-07-11 07:48:40 more
  • python筆記：第六章函式&方法

    # 1.系統函式 由系統提供，直接拿來用或是匯入模塊后使用 ``` a = 1.12386 result = round(a,2) print(result) > 1.12 ``` # 2.自定義函式 * 函式是結構化編程的核心 * 使用關鍵詞`def`來定義函式 ``` #函式定義 def fun ......

    uj5u.com 2023-07-11 07:48:37 more
  • == 與 equals 的區別？

    一. 介紹： Java中的 "==" 是一個運算子，是用于比較兩個物件地址值或基本資料型別之間的值是否相等。它的來源可以追溯到C語言，以及受C語言影響的許多其他編程語言。 Java中的equals() 是一個方法，可重寫該方法用于比較兩個物件屬性內容是否相等的方法。該方法繼承自Object類，在Ja ......

    uj5u.com 2023-07-11 07:48:32 more
  • Java 包、訪問修飾符

    # Java 包、訪問修飾符 # 1. 包 ## 包可以理解為創建不同的目錄來分別存放類，類似計算機當中檔案夾 > ## 通過包可以讓相同的類在不同的目錄下使用，防止重名的問題 > > ## 通過包可以很好的管理我們撰寫的類 > > ## 通過包可以控制訪問范圍 ## 使用 idea 工具創建包通過 ......

    uj5u.com 2023-07-11 07:48:27 more
  • Java入門11（JDBC）

    ## JDBC 驅動加載 => 連接創建 => 創建編譯 / 預編譯陳述句 => 獲取結果集 => 遍歷結果集 => 回傳結果集 | 介面 | | | | | | Driver | 驅動 | | Connection | 連接 | | Statement | 操作 | | ResultSet | 結果 ......

    uj5u.com 2023-07-11 07:46:29 more
  • C++類模板實作工廠模式（優化if else/switch case）

    引自：https://blog.csdn.net/weixin_43795921/article/details/127224633 template <typename IdentifierType, class AbstractProduct, class ProductCreator = Ab ......

    uj5u.com 2023-07-11 07:46:24 more
  • Java 中怎樣將 bytes 轉換為 long 型別?

    將bytes 轉換為long型別： 第一種方式： String 接收 bytes 的構造器轉成 String，再 Long.parseLong； 但此種情況需要注意：位元組陣列中的每個位元組都必須是有效的數字字符。如果位元組陣列包含非數字字符，則會引發NumberFormatException例外。確保在 ......

    uj5u.com 2023-07-11 07:46:20 more

Copyright © 2010-2020 有解無憂