【调制解调】PM调相-6ren

【调制解调】PM调相

转载作者：我是一只小鸟更新时间：2023-07-18 14:34:16

说明

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

说明
1. PM 调制算法
- 1.1 PM 信号描述
- 1.2 PM 信号调制示例
2. PM 解调算法
- 2.1 FM 解调积分法
- 2.2 希尔伯特求瞬时相位法
- 2.3 数字正交解调
3. PM 仿真（MATLAB Communications Toolbox）
参考资料
附录代码
- 附.1 文件 mod_pm.m
- 附.2 文件 main_modPM_example.m
- 附.3 文件 lpf_filter.m
- 附.4 文件 demod_pm_method1.m
- 附.5 文件 main_demodPM_example1.m
- 附.6 文件 demod_pm_method2.m
- 附.7 文件 main_demodPM_example2.m
- 附.8 文件 demod_pm_method3.m
- 附.9 文件 main_demodPM_example3.m
- 附.10 文件 main_CommPM_example.m

1. PM 调制算法

1.1 PM 信号描述

用调制信号去控制载波的瞬时相位，使其按照调制信号的规律变化，当调制信号是模拟信号时，这个过程就被称为调相（PM）。在 FM 调频与 PM 调相的过程中，载波的幅度都保持恒定不变，而频率和相位的变化都表现为载波瞬时相位的变化，故把调频和调相统称为角度调制，与幅度调制技术相比，角度调制最突出的优势是其较高的抗噪声性能，然而获得这种优势的代价是角度调制占用比幅度调制信号更宽的带宽。PM 信号的时域表达式为:

\[s_{PM}(t)=Acos{\left[\omega_ct+{K_p}m(t)\right]} \tag{1} \]

式中： $A$ 为载波恒定振幅， $K_p$ 为调相灵敏度（单位 $rad/V$ ）， $m(t)$ 是调制信号（携带要发出去的信息）， $cos{\omega_ct}$ 是载波， $\omega_c$ 是载波角频率，与载波频率 $f_c$ 之间的关系为 $\omega_c=2{\pi}f_c$ 。由式 $(1)$ 可得 PM 信号相对于载波相位 ${\omega}_ct$ 的瞬时相位偏移为:

\[\varphi(t)={K_p}m(t) \tag{2} \]

由 $(1)$ 式可知 PM 信号相对于载波相位 ${\omega}_ct$ 的瞬时相位偏移随 $m(t)$ 呈线性变化，比例系数为 $K_p$ 。PM 信号与 FM 信号有很多相同之处，如果预先不知道调制信号 $m(t)$ 的具体形式，则无法判断已调信号是 PM 信号还是 FM 信号。PM 信号的调相指数（调制指数） $\beta$ 被定义为最大的相位偏移:

\[\beta={K_p}{\lvert}{m(t)}{\rvert}_{max} \tag{3} \]

若 $m(t)$ 为单一频率的正弦波（即 $m(t)={A_m}cos({2{\pi}{f_m}t})$ ），则调制指数的表达式如下.

\[{\beta}={K_p}{\lvert}{m(t)}{\rvert}_{max}={K_p}{A_m} \tag{4} \]

1.2 PM 信号调制示例

与 FM 信号一样，PM 信号的调制方法也可分为 3 种，一种是直接调相法，一种是间接调相法，第三种是正交调制法.

（1）直接调相法。

直接调相法是用调制信号 $m(t)$ 直接改变谐振回路的参数，使载波信号通过谐振回路时产生相移而形成调相波，最常用的直接调相电路是变容二极管调相器.

（2）间接调相法。

间接调相电路比直接调相电路复杂，由于频率与相位之间存在微分与积分的关系，所以 FM 与 PM 之间是可以相互转换的。先对调制信号进行微分，然后对微分结果进行调频，就得到了调相信号，这种方式叫间接调相。先对调制信号进行积分，然后对积分结果进行调相，就得到了调频信号，这种方式叫间接调频.

（3）正交调制法。

将 $(1)$ 式进行三角展开，可以得到:

\[\begin{aligned} s_{PM}(t)&=Acos{\left[\omega_ct+{K_p}m(t)\right]}\\[1em] &=Acos{(\omega_ct)}cos{\left[{K_p}m(t)\right]}-Asin{(\omega_ct)}sin{\left[{K_p}m(t)\right]}\\[1em] \end{aligned} \tag{5} \]

正交调制流程如下，其中第 2 步也可以先将 $I(t)$ 和 $Q(t)$ 组成一个复信号 $Z(t)=I(t)+iQ(t)$ ，然后乘以复载波 $exp(i\omega_ct)$ ，最后取实部，即 $s_{PM}(t)=Real\left[Z(t)exp(i\omega_ct)\right]$ ，两种方法是等价的。:

对调制信号 $m(t)$ 分别取余弦和正弦，得到 $I$ 路数据 $I(t)=cos{\left[{K_p}m(t)\right]}$ 与 $Q$ 路数据 $Q(t)=sin{\left[{K_p}m(t)\right]}$ 。
分别乘以载波 $Acos(\omega_ct)$ 与 $-Asin(\omega_ct)$ 后相加，得到 PM 信号 $s_{PM}(t)=I(t)Acos(\omega_ct)-Q(t)Asin(\omega_ct)$ 。

调制信号 $m(t)$ 可以是确知信号，也可以是随机信号。当 $m(t)$ 是确知信号时，不妨假设 $m(t)$ 的时域表达式如下:

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

各调制参数取值： $A=1$ ， $f_m=2500Hz$ ， ${\beta}=4$ ， $f_c=20000Hz$ 。信号采样率 $f_s=8{f_c}$ ，仿真总时长为 $2s$ 。PM 调制效果如下图所示（为了美观，时域只显示前 500 个点），调制信号 $m(t)$ 双边幅度谱有四根离散谱线（ ${\pm}2500Hz$ 、 ${\pm}1250Hz$ ），载波 $c(t)$ 的双边幅度谱有两根离散谱线（ ${\pm}20000Hz$ ），代码详见附录 main_modPM_example.m 与 mod_pm.m ，设置 beta = 4 即可.

Oh Shit!-图片走丢了-打个广告-欢迎来博客园关注“木三百川”

2. PM 解调算法

解调是调制的逆过程，其作用是从接收的已调信号中恢复原基带信号（即调制信号）。PM 解调时，可以先使用 FM 解调器，然后对解调结果进行积分。另一种思路是直接提取 PM 信号的相位信息，有希尔伯特变换法以及数字正交解调法。下面分别用几种不同方法对 1.2 节中确知信号的 PM 调制结果进行解调.

2.1 FM 解调积分法

FM 信号的解调详见本人同系列博客【调制解调】FM 调频，内有 FM 解调器相关代码。这种方法的流程如下:

第一步 ：进行 FM 解调。
第二步 ：对 FM 解调结果进行积分。
第三步 ：去除直流分量（减去自身均值）。

对 1.2 节中的 PM 信号，设定信噪比 $SNR=50dB$ ，解调效果如下，解调后幅度放大系数 $k=\overline{{\lvert}m(t){\rvert}}/\overline{{\lvert}\hat{m}(t){\rvert}}\approx70439.28$ ，使用这个系数放大解调信号幅值，然后计算误差，有： $\sqrt{\sum{{\lvert}m(t_i)-k\hat{m}(t_i){\rvert}^2}}/\sqrt{\sum{{\lvert}m(t_i){\rvert}^2}}\approx0.0388$ .

Oh Shit!-图片走丢了-打个广告-欢迎来博客园关注“木三百川”

代码详见 demod_pm_method1.m 和 main_demodPM_example1.m ，调频解调器的代码详见本人同系列博客【调制解调】FM 调频 .

2.2 希尔伯特求瞬时相位法

可以通过希尔伯特变换的方法求得 PM 信号的瞬时相位，然后减去由载波导致的线性分量 ${\omega_ct}$ ，获得非线性瞬时相位，也就是相对于载波相位 ${\omega}_ct$ 的瞬时相位偏移。这种方法的流程如下:

第一步 ：计算 $s(t)$ 的希尔伯特变换，得到一个复信号（实部为 $s(t)$ ，虚部为其希尔伯特变换结果），对所得复信号求相位角并去卷叠，得到 PM 信号的瞬时相位 ${\phi}(t)$ 。
第二步 ：减去由载波导致的线性分量 $\omega_ct$ 以及载波初相 $\varphi_0$ ，得到解调结果 $m_o(t)={\phi}(t)-\omega_ct-\varphi_0$ 。

对 1.2 节中的 PM 信号，设定信噪比 $SNR=50dB$ ，解调效果如下，解调后幅度放大系数 $k=\overline{{\lvert}m(t){\rvert}}/\overline{{\lvert}\hat{m}(t){\rvert}}\approx0.44$ ，使用这个系数放大解调信号幅值，然后计算误差，有： $\sqrt{\sum{{\lvert}m(t_i)-k\hat{m}(t_i){\rvert}^2}}/\sqrt{\sum{{\lvert}m(t_i){\rvert}^2}}\approx0.0010$ .

Oh Shit!-图片走丢了-打个广告-欢迎来博客园关注“木三百川”

代码详见 demod_pm_method2.m 和 main_demodPM_example2.m .

2.3 数字正交解调

数字正交解调属于相干解调的一种，但这种方法具有较强的抗载频失配能力，不要求相干载波严格的同相。PM 数字正交解调一般有以下四个步骤:

第一步 ：乘以正交相干载波得到 ${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})$ 。
第二步 ：低通滤波器滤除 ${s_I}(t)$ 与 ${s_Q}(t)$ 中的高频分量。
第三步 ：通过反正切函数计算相位，得到解调结果 $m_o(t)=atan\left[\frac{s_Q(t)}{s_I(t)}\right]=km(t)+\Delta\phi$ 。
第四步 ：去除直流分量（减去自身均值）。

Oh Shit!-图片走丢了-打个广告-欢迎来博客园关注“木三百川”

代码详见 demod_pm_method3.m 和 main_demodPM_example3.m 。更改相干载波的初始相位为 ${\phi_0}=\pi/4,\pi/2$ 后，解调效果依然很好；但更改相干载波的中心频率为 $0.8f_c,1.2f_c$ 后，解调效果变得很差，说明这种方法具有较好的抗载频相位失配能力.

3. PM 仿真（MATLAB Communications Toolbox）

MATLAB 的 Communications Toolbox 中提供了 PM 调制函数 pmmod ，高斯白噪声函数 awgn ，以及 PM 解调函数 pmdemod ，可以很方便地完成 PM 信号仿真。使用这三个函数实现上面 1.2 节中确知信号 $m(t)$ 的 PM 调制解调（将调相指数 $\beta$ 改为 $1$ ），调制后加噪声的效果如下:

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.0022$ 。代码详见附录 main_CommPM_example.m .

参考资料

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

[2] 樊昌信,曹丽娜.通信原理.第7版[M].国防工业出版社,2012. 。

[3] 刘学勇.详解MATLAB/Simulink通信系统建模与仿真[M].电子工业出版社,2011. 。

[4] 王丽娜,王兵.卫星通信系统.第2版[M].国防工业出版社,2014. 。

附录代码

附.1 文件 mod_pm.m

                        
                          function [ sig_pm ] = mod_pm(fc, beta, fs, mt, t, A)
% MOD_PM        PM 调相
% 输入参数：
%       fc      载波中心频率
%       beta    调频指数/调制指数
%       fs      信号采样率
%       mt      调制信号
%       t       采样时间
%       A       载波恒定振幅
% 输出参数：
%       sig_pm  调相(PM)实信号
% @author 木三百川

% 计算调相灵敏度
Kp = beta/max(abs(mt));

% 生成信号
ct = A*cos(2*pi*fc*t);           % 载波 
sig_pm = A*cos(2*pi*fc*t+Kp*mt); % PM调相信号

% 绘图
nfft = length(sig_pm);
freq = (-nfft/2:nfft/2-1).'*(fs/nfft);
figure;set(gcf,'color','w');
plot_length = min(500, length(sig_pm));
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_pm(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('PM调相信号s(t)');
subplot(3,2,6);
plot(freq, 10*log10(fftshift(abs(fft(sig_pm,nfft)/nfft))+eps));xlim([freq(1),freq(end)]);
xlabel('频率/hz');ylabel('幅度/dB');title('PM调相信号s(t)双边幅度谱');

end

附.2 文件 main_modPM_example.m

                        
                          clc;
clear;
close all;
% PM 调制仿真(调制信号为确知信号)
% @author 木三百川

% 调制参数
A = 1;                  % 载波恒定振幅
fm = 2500;              % 调制信号参数
beta = 4;               % 调相指数/调制指数
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);

% PM 调制
[ sig_pm ] = mod_pm(fc, beta, fs, mt, t, A);

附.3 文件 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

附.4 文件 demod_pm_method1.m

                        
                          function [ sig_pm_demod ] = demod_pm_method1(sig_pm_receive, fc, fs, t, phi0)
% DEMOD_PM_METHOD1        PM 解调(FM 解调积分法)
% 输入参数：
%       sig_pm_receive      PM 接收信号，行向量
%       fc                  载波中心频率
%       fs                  信号采样率
%       t                   采样时间
%       phi0                相干载波初始相位
% 输出参数：
%       sig_pm_demod        解调结果，与 sig_pm_receive 等长
% @author 木三百川

% 第一步：进行 FM 解调
[ sig_pm_demod ] = demod_fm_method4(sig_pm_receive, fc, fs, t, phi0);

% 第二步：对 FM 解调结果进行积分
sig_pm_demod = cumtrapz(t, sig_pm_demod);

% 第三步：去直流
sig_pm_demod = sig_pm_demod - mean(sig_pm_demod);

end

附.5 文件 main_demodPM_example1.m

                        
                          clc;
clear;
close all;
% PM 解调仿真(调制信号为确知信号，FM 解调积分法)
% @author 木三百川

% 调制参数
A = 1;                  % 载波恒定振幅
fm = 2500;              % 调制信号参数
beta = 4;               % 调相指数/调制指数
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);

% PM 调制
[ sig_pm_send ] = mod_pm(fc, beta, fs, mt, t, A);

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

% FM 解调积分法
ini_phase = 0;
[ sig_pm_demod ] = demod_pm_method1(sig_pm_receive, fc, fs, t, ini_phase);

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

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

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), coef*sig_pm_demod(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('解调效果');
legend('调制信号','解调信号(放大后)');

附.6 文件 demod_pm_method2.m

                        
                          function [ sig_pm_demod ] = demod_pm_method2(sig_pm_receive, fc, fs, t, phi0)
% DEMOD_PM_METHOD2        PM 解调(希尔伯特求瞬时相位法)
% 输入参数：
%       sig_pm_receive      PM 接收信号，行向量
%       fc                  载波中心频率
%       fs                  信号采样率
%       t                   采样时间
%       phi0                相干载波初始相位
% 输出参数：
%       sig_pm_demod        解调结果，与 sig_pm_receive 等长
% @author 木三百川

% 第一步：使用希尔伯特变换计算瞬时相位
inst_phase = unwrap(angle(hilbert(sig_pm_receive)));

% 第二步：去除载波分量
sig_pm_demod = inst_phase - 2*pi*fc*t - phi0;

end

附.7 文件 main_demodPM_example2.m

                        
                          clc;
clear;
close all;
% PM 解调仿真(调制信号为确知信号，希尔伯特求瞬时相位法)
% @author 木三百川

% 调制参数
A = 1;                  % 载波恒定振幅
fm = 2500;              % 调制信号参数
beta = 4;               % 调相指数/调制指数
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);

% PM 调制
[ sig_pm_send ] = mod_pm(fc, beta, fs, mt, t, A);

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

% 希尔伯特求瞬时相位法
ini_phase = 0;
[ sig_pm_demod ] = demod_pm_method2(sig_pm_receive, fc, fs, t, ini_phase);

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

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

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), coef*sig_pm_demod(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('解调效果');
legend('调制信号','解调信号(放大后)');

附.8 文件 demod_pm_method3.m

                        
                          function [ sig_pm_demod ] = demod_pm_method3(sig_pm_receive, fc, fs, t, phi0)
% DEMOD_PM_METHOD4        PM 数字正交解调/相干解调
% 输入参数：
%       sig_pm_receive      PM 接收信号，行向量
%       fc                  载波中心频率
%       fs                  信号采样率
%       t                   采样时间
%       phi0                相干载波初始相位
% 输出参数：
%       sig_pm_demod        解调结果，与 sig_pm_receive 等长
% @author 木三百川

% 第一步：乘以正交相干载波
sig_pm_i = sig_pm_receive.*cos(2*pi*fc*t+phi0);
sig_pm_q = -sig_pm_receive.*sin(2*pi*fc*t+phi0);

% 第二步：低通滤波
sig_pm_i_lpf = lpf_filter(sig_pm_i, fc/(fs/2));
sig_pm_q_lpf = lpf_filter(sig_pm_q, fc/(fs/2));

% 第三步：计算相位
sig_pm_demod = unwrap(atan2(sig_pm_q_lpf, sig_pm_i_lpf));

% 第四步：去直流
sig_pm_demod = sig_pm_demod - mean(sig_pm_demod);

end

附.9 文件 main_demodPM_example3.m

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

% 调制参数
A = 1;                  % 载波恒定振幅
fm = 2500;              % 调制信号参数
beta = 4;               % 调相指数/调制指数
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);

% PM 调制
[ sig_pm_send ] = mod_pm(fc, beta, fs, mt, t, A);

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

% 数字正交解调
ini_phase = 0;
[ sig_pm_demod ] = demod_pm_method3(sig_pm_receive, fc, fs, t, ini_phase);

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

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

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), coef*sig_pm_demod(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('解调效果');
legend('调制信号','解调信号(放大后)');

附.10 文件 main_CommPM_example.m

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

% 调制参数
A = 1;                  % 载波恒定振幅
fm = 2500;              % 调制信号参数
beta = 1;               % 调相指数/调制指数
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);

% FM 调制
ini_phase = 0;
sig_pm_send = A*pmmod(mt,fc,fs,beta,ini_phase);

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

% PM 解调
[ sig_pm_demod ] = pmdemod(sig_pm_receive,fc,fs,beta,ini_phase);

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

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_pm_demod(1:plot_length));xlim([t(1),t(plot_length)]);
xlabel('t/s');ylabel('幅度');title('解调效果');
legend('调制信号','解调信号');

fprintf('norm(调制信号 - 解调信号)/norm(调制信号) = %.4f.\n', norm(mt-sig_pm_demod)/norm(mt));

最后此篇关于【调制解调】PM调相的文章就讲到这里了,如果你想了解更多关于【调制解调】PM调相的内容请搜索CFSDN的文章或继续浏览相关文章，希望大家以后支持我的博客！。

文章推荐：我是如何组织Go代码的（目录结构依赖注入wire）

文章推荐： SAS-PROCFCMP

文章推荐：极速安装kubernetes-1.22.0（三台CentOS7服务器）

文章推荐：并发编程---信号量线程同步

PHP 设置 [pm.max_children、pm.start_servers、pm.min_spare_servers 和 pm.max_spare_servers]
在将网站从 GoDaddy 共享服务器传输到 EC2 实例的过程中。处理流量，在典型的一天的高峰时段，大约有 300 名活跃访客，至少可以说是有问题的。我的 CPU 使用率缓慢上升，最终达到 100%
php - PM 系统 - 预览以前的 PM
我在一个 PM 系统上工作，我希望在最后收到的 PM 上方列出一次对话的先前发送的 PM。但我的问题是:如何在数据库中设置这样的表？我玩了一会儿关于为每个特定对话使用一个 ID，但是该 ID 的来源是
excel - 如何使用 excel 公式从文本字符串 "7.00 pm to 8.00 pm"或 "11.00 am to 12.00 pm"计算小时差
我想从给定的文本字符串中计算小时数和分钟数，目前我使用一些 excel 内置函数创建了以下公式，示例文本字符串: 7.00pm to 8.00pm or 9.00am to 11.45pm or 1
c# - 日期时间自定义格式 - 将 AM/PM 显示为大写而不是 am/pm
我有一个字符串，其中包含:14 Dec 2011 9:45 am(注意“AM”没有大写) 现在我想用它创建一个 datetime 变量。我试着查看 this ，但这与我想要的相反。我也试过this
java - 如何在 Java 中更改时间格式 - "am/pm"与 "AM/PM"
之前创建的测试使用 DateTimeFormatter.ofPattern("ha"); 并返回 "10AM"(对于 '2017-04-09T10:00-06 :00[美国/山区]'). 在我的 Ma
python - 时间格式正则表达式 HH :MM AM/am/PM/pm in python
我编写了正则表达式来捕获 HH:MM AM/PM/am/pm 但它无法提取精确的模式正则表达式代码: import re def replace_entities(example): res
Mysql - 智能查找 "-"并附加 "AM"和 "PM"- 如果不附加 AM、PM
我的 Sql 查询如下所示: SELECT monday FROM `restaurantslive` 结果是我想做的是在分钟后附加“AM”和“PM”，例如其中一条记录是 6:30–10:30A
javascript - jquery datetimepicker 插件 am/pm - 数据库中的所有 pm 字符串
我正在使用 codeigniter 和 jquery datetimepicker 插件。在我的 View 页面上，我有类似的日期时间选择器格式 $(".traveler").datetimepic
java - 安卓时间0 :00 pm should be displaying 12:00 pm instead
这个问题在这里已经有了答案: 关闭 10 年前。 Possible Duplicate: Android date/time displaying 0 instead of 12 下面是我使用的代码
c# - 如何检查 6 :00 PM to 9:00 PM 等字符串是否存在时间冲突
我正在构建类似考试日期表的内容。我目前在查找时间之间的冲突时遇到问题.. 我有一个存储时间间隔的字符串列表，例如- List times = new List(); times.Add("6:00 P
java - 如何将 am-pm-of-day 格式化为 "AM/PM"而不是 "a.m./p.m."？
DateTimeFormatter's API reference似乎错过了正确格式化 am-pm-of-day 的细节，我总是可以使用 String.replace() 但我觉得将选项更改为 AM
android - 将 12 (PM/AM)Hour AM PM 时间中的字符串转换为 24 小时时间 android
我在转换来自服务器的时间时遇到问题，我想将其转换为 24 小时制。我正在使用以下代码: String timeComeFromServer = "3:30 PM"; SimpleDateFormat
android - Android 上 pm clear 和 pm uninstall -k 的区别
从 pm 帮助，我得到这个: pm uninstall: removes a package from the system. Options: -k: keep the data and c
php - 将 AM PM 时间字符串存储到 MySQL 中的 TIME 数据类型中并在显示时使用 AM PM 检索？
我在前端输入日期为 10:00 AM 、12:00 PM 等...(表示 12 小时格式)。现在我想将该值保存在数据库中的 time 数据类型列中。我如何将该 AM PM 值保存到 MySQL 中的
perl - 执行使用 Elasticsearch.pm 的 perl 脚本时，Module/Runtime.pm 处的版本格式无效(需要版本)
我在执行使用 Elasticsearch.pm(新版本，小写)模块的 perl 脚本时遇到问题。该脚本是正确的(我还使用 perl -c 选项检查了语法)，但是当我尝试执行它时，我收到了这个错误:
java - 为什么在 JDO 中是 "pm.newQuery(Employee.class)"而不是 "pm.newQuery(Employee)"？
在这个JDO中，为什么这里需要.class？ Query averageSalaryQuery = pm.newQuery(Employee.class); 如果可能，我更愿意编写这种更简洁的语法？
SQL - 将 24 小时 ("military") 时间 (2145) 转换为 "AM/PM time"(9 :45 pm)
我正在使用 2 个字段，这些字段存储为 smallint 军事结构化时间。编辑我在 IBM Informix Dynamic Server 版本 10.00.FC9 上运行 beg_tm 和 end_
java - Cron 表达式 : Run on specific time (1:30 pm , 15:24 pm) 仅适用于周一和周三
我刚刚查看了 http://www.cronmaker.com/并尝试为以下场景创建 cron 表达式。在周一和周二的两个特定时间(即 1:24 和 3:34)运行作业。我是为此生成以下表达式的。
c# - 转换 dd/MM/yyyy hh :mm am/pm to MM/dd/yyyy hh:mm am/pm
我想将字符串:24/11/2016 04:30 pm 转换为日期时间值:11/24/2016 04:30 pm。我的代码为: DateTime date = DateTime.ParseExact(
ruby - 在 ruby 的 strftime 中，为什么 %P 是 "am/pm"而 %p 是 "AM/PM"
这不是更困惑吗？难道不应该反过来这样更容易记住吗？最佳答案查看此备忘单:http://cheat.errtheblog.com/s/strftime/ Ruby 1.8 之前好像没有%P 选项，所

我是一只小鸟

个人简介

我是一名优秀的程序员,十分优秀！

作者热门文章

滴滴打车优惠券免费领取

全站热门文章

首页

博学

6Ren·AI

商城