调制与解调（一）——模拟调制系统调制解调方法

---

前言

根据天线理论，天线所能够发送的波长与天线长度存在着一定的联系，天线长度一般要大于波长的四分之一。那么在此规律下，若想传播低频信号则需要架设高度为数千米甚至数十前面的天线，这种方式显然是无法实现的。那么如何才能收发低频信号呢，调制的方法应运而生，所谓调制即将低频信号（调制信号或基带信号）加到高频信号（载波）上，得到已调信号。而解调是调制的逆过程，其目的是将调制信号从已调信号中恢复出来。
调制的方式有很多，本篇文章主要讨论模拟调制和解调的实现方法，以AM,DSB,SSB,NBFM,WBFM为例，进行matlab调制仿真，并使用相干解调和非相干解调对上述信号进行解调仿真。

一、幅度调制

幅度调制是用调制信号去控制高频载波的幅度，使之随调制信号做线性变化的过程。

1、调幅（AM）

常规的双边带幅度调制，简称调幅（AM）,其时域表达式为：$$\begin{gathered} s_{AM}=\left[A_0+m\left(t\right)\right]\cos {\omega }_{c}t \\ =A_0\cos {\omega }_{c}t+m\left(t\right)\cos {\omega }_{c}t \end{gathered}$$
其第一项为载波项，第二项为边带项，在调制过程中，携带调制信号信息的只有边带项，因此载波项为无用项，由此就造成了功率的浪费。因此要寻求对这一问题的改进。

2、抑制载波的双边带调制（DSB）

由上一部分对AM的分析，为了避免功率的浪费，考虑到如果在调制过程中直流部分去掉，即去掉$A_0$，既可以抑制载波，其时域表达式如下：
$$s_{DSB}=m\left(t\right)\cos {\omega }_{c}t$$
此时比较AM和DSB调制信号的频谱，AM信号的频偏由上边带，下边带，载波分量组成，而DSB由上边带和下边带组成，抑制了载波分量，避免了功率的浪费，但是上下边带都各自携带了全部的调制信号信息，就相当于同样的信息出现了两次，这就造成了带宽的浪费（结论为DSB和AM信号的带宽都为2倍的调制信号带宽），进一步寻求改进，即想办法消除一个边带以减小带宽，增大频带利用率。

3、抑制载波的单边带调制（SSB）

产生SSB信号的方法有两种，滤波法和相移法。

3.1滤波法

首先滤波法比较容易理解：首先产生一个双边带信号，再让其通过一个边带滤波器滤除不要的边带即可得到单边带信号。器边带滤波器选取理想高通或低通滤波器，如采用高通滤波器：$$H(\omega )=\{\begin{array}{ll}1,& |\omega |>{\omega }_c\\ 0,& |\omega |<{\omega }_c\end{array}$$
则可滤除下边带，保留上边带，若采用低通滤波器：
$$H(\omega )=\{\begin{array}{ll}1,& |\omega |<{\omega }_c\\ 0,& |\omega |>{\omega }_c\end{array}$$
则滤除上边带保留下边带，此种方法的难点是边带滤波器的制作，实际上的滤波器都不具备如此理想的边缘特性。

3.2相移法

相移法的原理为：利用相依网络，对载波信号和调制信号适当相移以便在何处过程中将其中一个边带抵消掉而获得SSB，由此就不需要实现滤波法中边带滤波器陡峭的边带条件。将相移法的过程简单推导如下：
以单频信号为例，调制信号为:$m\left(t\right)=A_m\cos {\omega }_{m}t$，载波信号为：$c\left(t\right)=\cos {\omega }_{c}t$,
则其DSB的时域表达式为：$$\begin{array}{rl}{s}_{DSB}& =m(t)\cos {\omega }_{c}t\\ & =A_m\cos {\omega }_{m}t\cos {\omega }_{c}t\\ & =\frac{1}{2}{A}_m\cos \left({\omega }_c+{\omega }_m\right)t+\frac{1}{2}{A}_m\cos \left({\omega }_c-{\omega }_m\right)t\end{array}$$
保留上边带为：$$\begin{array}{rl}& \frac{1}{2}{A}_m\cos \left({\omega }_c+{\omega }_m\right)t\\ =& \frac{1}{2}{A}_m\cos {\omega }_{m}t\cos {\omega }_{c}t-\frac{1}{2}{A}_m\sin {\omega }_{m}t\sin {\omega }_{c}t\end{array}$$
保留下边带为：
$$\begin{array}{rl}& \frac{1}{2}{A}_m\cos \left({\omega }_c-{\omega }_m\right)t\\ =& \frac{1}{2}{A}_m\cos {\omega }_{m}t\cos {\omega }_{c}t+\frac{1}{2}{A}_m\sin {\omega }_{m}t\sin {\omega }_{c}t\end{array}$$

因此可以得出:$$\begin{array}{rl}{s}_{SSB}\left(t\right)& \\ =& \frac{1}{2}{A}_m\cos {\omega }_{m}t\cos {\omega }_{c}t\mp \frac{1}{2}{A}_m\sin {\omega }_{m}t\sin {\omega }_{c}t\end{array}$$
或者写成：$$s_{SSB}\left(t\right)=\frac{1}{2}m\left(t\right)\cos {\omega }_{c}t\mp \frac{1}{2}\hat{m}\left(t\right)\sin {\omega }_{c}t$$
其中，$\hat{m}\left(t\right)$为$m\left(t\right)$的希尔伯特变换。
而相移法的技术难点在于相移网络的制作，很难将$m\left(t\right)$的所有频率分量精准相移$\frac{\pi }{2}$。
综上所述，要实现SSB无论是滤波法和相移法都有其技术难点，因此需要找出一种折衷的方法，能滤除部分边带且较容易实现。

4、残留边带调制

残留边带调制（VSB）就是这样一种折衷方法，既克服了DSB占用频带宽的特点，又解决了SSB实现困难的特点。VSB不用陡峭的滤波器，而是用比较平缓的滤波器切除无用边带。

二、角度调制

1、FM和PM的联系

对于一个信号而言，由于其频率和相位存在微积分的关系，因此调频和调相是可以相互转换的，由此可引出：直接调频，间接调频，直接调相，间接调相，其系统框图如图所示。由于调频与调相有这样的关系，在下面讨论时就以FM为例。

调频信号的时域表达式为：
$$s_{FM}\left(t\right)=A\cos \left[{\omega }_{c}t+K_f\int m\left(\tau \right)d\tau \right]$$
其中$K_f$称为调频灵敏度。
调相信号的时域表达式为：$$s_{PM}\left(t\right)=A\cos \left[{\omega }_{c}t+K_{P}m\left(\tau \right)\right]$$
其中$K_p$称为调相灵敏度。
两种调频方法相比较，直接调频和间接调频各有利弊，直接调频法在实现线性调频的前提下可获得较大频偏，但频率稳定度不高；间接调频法稳定性好但电路复杂，且频偏相对于直接跳频要小。

2、窄带调频NBFM与窄带调频WBFM

若FM信号满足：$$\left|K_f\int m\left(\tau \right)d\tau \right|\ll 0.5$$
则称其为窄带调频，否则则称为宽带调频。
窄带调频信号的时域表达式为：$$S_{NBFM}\left(t\right)=A\cos \left({\omega }_{c}t\right)-\left[A{K}_f\int m\left(\tau \right)d\tau \right]\sin {\omega }_{c}t$$
其频谱为：$$\begin{gathered} S_{NBFM}\left(\omega \right)=\pi A\left[\delta \left(\omega +{\omega }_c\right)+\delta \left(\omega -{\omega }_c\right)\right]+ \\ \frac{A{K}_f}{2}\left[\frac{M\left(\omega -{\omega }_c\right)}{\omega -{\omega }_c}-\frac{M\left(\omega +{\omega }_c\right)}{\omega +{\omega }_c}\right] \end{gathered}$$

宽带调频时，频谱不再是调制信号的线性搬移，而是一个非线性过程。

三、相干解调

相干解调的原理与调制相似，都是频谱搬移，需要在接收端提供与已调信号完全同步的本地载波，与已调信号相乘后，经低通滤波器滤出低频分量。适用于对所有线性调制系统的解调。

四、非相干解调（包络检波）

适用于当${\left|m\left(t\right)\right|}_{\max }⩽A_0$的情况，包络与调制信号$m(t)$的形状一样，可以直接从包络中提取载波信息，在两种方法都适用的情况下，非相干解调较为简单，但存在门限效应。
门限效应：当输入信噪比低于某一门限后，输出信噪比不再按比例下降，而是急剧恶化，这是由包络检波器的非线性解调作用引起的。

五、调制解调的matlab仿真

1、线性调制：

首先定义载波信号和调制信号：

mt=cos(2*pi*fmax*t); %基带信号
s_am=A*cos(2*pi*fc*t);%载波信号

再把调制信号加到载波信号上实现调幅，这里用相移法实现单边调制。

s_am=(A+mt).*cos(2*pi*fc*t);%AM调制信号
s_dsb=mt.*cos(2*pi*fc*t);%DSB调制
s_ssb=real(hilbert(mt).*exp(j*2*pi*fc*t));%SSB调制
[f,xf_am]=FFT_SHIFT(t,s_am); %AM调制信号频谱
[f,xf_dsb]=FFT_SHIFT(t,s_dsb); %DSB调制信号频谱
[f,xf_ssb]=FFT_SHIFT(t,s_ssb); %SSB调制信号频谱

通过以上步骤得出线性调制后的时域结果及其频谱：



其频谱如图所示：

因此所得结果符合预计，DSB信号较AM信号而言滤除了载波分量，而SSB较DSB信号而言只保留了一个边带。

2、非线性调制

若要仿真宽带调频和窄带调频，需要注意其判别条件即：$$\left|K_f\int m\left(\tau \right)d\tau \right|\ll 0.5$$

Kf_1=150; %宽带调频系数
Kf_2=90;  %窄带调频系数
mt=cos(2*pi*fm*t) ; %载波信号
mti=1/2/pi/fm*sin(2*pi*fm*t);%mt的积分函数
FMt=A0*cos(2*pi*fc*t+2*pi*Kf_1*mti);
figure(1)
subplot(211);
plot(t,FMt);hold on;
xlabel('t');ylabel('WBFM调制信号')
subplot(212)
[f,sf]=FFT_SHIFT(t,FMt);
plot(f,abs(sf));
axis([-1000 1000 0 1])
xlabel('f');xlabel('f /hz');ylabel('WBFM调频信号的幅度谱');
%NBFM调制
FMt_2=A0*cos(2*pi*fc*t+2*pi*Kf_2*mti);
figure(2)
subplot(211);
plot(t,FMt_2);hold on;
xlabel('t');ylabel('NBFM调制信号')
subplot(212)
[f,sf]=FFT_SHIFT(t,FMt_2);
plot(f,abs(sf));
axis([-1000 1000 0 1])
xlabel('f');xlabel('f /hz');ylabel('NBFM调频信号的幅度谱');

3、相干解调

[t,am_t]=RECT_LPF(f,AMf,B);
figure()
subplot(311);
plot(t,am_t);
xlabel('AM加噪后解调信号');
n_dsb=n_dsb-mean(n_dsb);
[f,DSBF]=FFT_SHIFT(t,n_dsb);
[t,dsb_t]=RECT_LPF(f,DSBF,B);
subplot(312)
plot(t,dsb_t);
xlabel('DSB加噪后解调信号');

其中低通滤波器函数RECT_LPF定义为：

function [t,st]=RECT_LPF(f,Sf,B)
df=f(2)-f(1);
fN=length(f);
RectH=zeros (1,fN);
BN=f1oor(B/df);
BN_SHIFT=[-BN :BN-1]+floor(fN/2);
RectH(BN_SHIFT)=1;
Yf=RectH.*Sf;
[t,st]=IFFT_SHIFT(f,Yf);