模拟幅度调制系统抗干扰性能仿真分析

1.引言

1.研究目的

本课程设计课题主要研究基于MATLAB的三种调幅调制和相干解调实现以及抗噪声性能的理论设计和软件仿真方法。掌握模拟系统AM、DSB-SC、SSB调制解调的原理以及叠加加性高斯白噪随后的系统的抗噪声性能分析;掌握AM、DSB-SC、SSB调制解调模拟系统的理论设计方法;掌握应用MATLAB分析系统时域、频域特性的方法,进一步锻炼应用MATLAB进行编程仿真的能力。

2.研究方法

本文通过比较解调器输出信噪比和输入信噪比的比值即信噪比增益来表示调制解调系统的抗噪声性,在三种调制解调系统对应条件相同的情况下,信噪比增益越大,说明该系统抗噪声性能越好。本文中AM、DSB-SC、SSB调制解调模拟系统的解调器均采用相干解调。

3.主要内容

本文主要内容包括AM、DSB-SC、SSB调制器,相干解调器,调制解调系统抗噪声性能分析。以及三种调制解调系统的matlab仿真实现和结果,理论结果与仿真结果的误差分析。

2.系统模型

1.AM调制器及相干解调模型

1.AM调制器

在这里插入图片描述

图2.1 AM调制模型

a图为AM调制模型。
m ( t ) \mathit m(t) m(t)为调制信号,由图可得:
s A M ( t ) = A c [ 1 + m ( t ) ] cos ⁡ 2 π f c t (2.1.1.1) \mathit s_{AM}(t) =\mathit A_c[1+\mathit m(t)]\cos2\pi\mathit f_ct \tag{2.1.1.1} sAM(t)=Ac[1+m(t)]cos2πfct(2.1.1.1)
m ( t ) \mathit m(t) m(t) s A M ( t ) \mathit s_{AM}(t) sAM(t)的时域波形图如下:
在这里插入图片描述
图2.2时域波形图

上图中: A = 1 \mathit A=1 A=1 f ( t ) = m ( t ) \mathit f(t)=\mathit m(t) f(t)=m(t) c ( t ) = A c cos ⁡ 2 π f c t \mathit c(t)=\mathit A_c\cos2\pi\mathit f_ct c(t)=Accos2πfct
为了避免包络检波发生错误,能正确解调出原信号。要求:调幅指数 β A M = m a x ∣ m ( t ) ∣ ≤ 1 \beta_{AM}=max|m(t)|\le1 βAM=maxm(t)1
m ( t ) \mathit m(t) m(t)的频谱为 M ( f ) \mathit M(f) M(f),则 s A M ( t ) \mathit s_{AM}(t) sAM(t)的频谱 s A M ( f ) \mathit s_{AM}(f) sAM(f)为:
s A M ( f ) = A c 2 [ σ ( f − f c ) + σ ( f + f c ) ] + A c 2 [ M ( f − f c ) + M ( f + f c ) ] (2.1.1.2) \mathit s_{AM}(f)={\mathit A_c\over2}[\sigma(\mathit f-\mathit f_c)+\sigma(\mathit f+\mathit f_c)]+{\mathit A_c\over2}[\mathit M(\mathit f-\mathit f_c)+\mathit M(\mathit f+\mathit f_c)]\tag{2.1.1.2} sAM(f)=2Ac[σ(ffc)+σ(f+fc)]+2Ac[M(ffc)+M(f+fc)](2.1.1.2)
m ( t ) \mathit m(t) m(t)频谱图如下,可得 s A M ( t ) \mathit s_{AM}(t) sAM(t)的频谱图:
在这里插入图片描述

图2.3 频域波形图

上图中: B \mathit B B m ( t ) \mathit m(t) m(t)的带宽

2.叠加噪声后的AM相干解调模型

在这里插入图片描述

图2.4 AM解调模型

由AM调制器得: s ( t ) = s A M ( t ) = A c [ 1 + m ( t ) ] cos ⁡ 2 π f c t (2.1.2.1) \mathit s(t)=\mathit s_{AM}(t) =\mathit A_c[1+\mathit m(t)]\cos2\pi\mathit f_ct\tag{2.1.2.1} s(t)=sAM(t)=Ac[1+m(t)]cos2πfct(2.1.2.1)

n i ( t ) \mathit n_i(t) ni(t)为加性高斯白噪声,其均值为零,双边功率谱密度为 n 0 2 {\mathit n_0\over2} 2n0。已调信号 s ( t ) \mathit s(t) s(t)与加性噪声相加后通过接收机的带通滤波器BFT,设 m ( t ) \mathit m(t) m(t)的带宽为 B \mathit B B,滤波器的等效噪声带宽为 2 B 2\mathit B 2B

宽带噪声通过BFT后成为窄带白高斯噪声 n ( t ) \mathit n(t) n(t) n ( t ) = n c ( t ) cos ⁡ 2 π f c t − n s ( t )   s i n 2 π f c t (2.1.2.2) \mathit n(t)=\mathit n_c(t) \cos2\pi\mathit f_ct- \mathit n_s(t)\ sin2\pi\mathit f_ct\tag{2.1.2.2} n(t)=nc(t)cos2πfctns(t) sin2πfct(2.1.2.2)
其中 n c ( t ) \mathit n_c(t) nc(t) n s ( t ) \mathit n_s(t) ns(t)分别是 n ( t ) \mathit n(t) n(t)的同相分量和正交分量,均值都为零,方差与 n ( t ) \mathit n(t) n(t)相同。

r ( t ) \mathit r(t) r(t)通过乘法器,与恢复载波 cos ⁡ 2 π f c t \cos2\pi\mathit f_ct cos2πfct相乘得:
输出信号: s d ( t ) = s ( t ) cos ⁡ 2 π f c t = 1 2 A c [ 1 + m ( t ) ]   [ 1 + cos ⁡ 4 π f c t ] (2.1.2.3) \mathit s_d(t)=\mathit s(t)\cos2\pi\mathit f_ct={1\over2}\mathit A_c[1+\mathit m(t)]\ [1+\cos4\pi\mathit f_ct]\tag{2.1.2.3} sd(t)=s(t)cos2πfct=21Ac[1+m(t)] [1+cos4πfct](2.1.2.3)
输出噪声: n d ( t ) = n ( t ) cos ⁡ 2 π f c t = n c ( t )   c o s 2 2 π f c t − n s ( t )   s i n 2 π f c t cos ⁡ 2 π f c t = 1 2 n c ( t ) + 1 2 n c ( t ) cos ⁡ 4 π f c t − 1 2 n s ( t ) sin ⁡ 4 π f c t (2.1.2.4) \begin{aligned}\mathit n_d(t)&=\mathit n(t)\cos2\pi\mathit f_ct\\ &=\mathit n_c(t) \ cos^22\pi\mathit f_ct- \mathit n_s(t)\ sin2\pi\mathit f_ct\cos2\pi\mathit f_ct\\ &={1\over2}\mathit n_c(t)+{1\over2}\mathit n_c(t)\cos4\pi\mathit f_ct-{1\over2}\mathit n_s(t)\sin4\pi\mathit f_ct\\ \end{aligned}\tag{2.1.2.4} nd(t)=n(t)cos2πfct=nc(t) cos22πfctns(t) sin2πfctcos2πfct=21nc(t)+21nc(t)cos4πfct21ns(t)sin4πfct(2.1.2.4)

再通过低通滤波器LPF,抑制二倍频分量,仅通过低通分量,从而实现相干解调,低通滤波器带宽为 B \mathit B B。得到相干解调后
输出信号: m o ( t ) = 1 2 A c m ( t ) (2.1.2.5) \mathit m_o(t)={1\over2}\mathit A_c\mathit m(t)\tag{2.1.2.5} mo(t)=21Acm(t)(2.1.2.5)输出噪声: n o ( t ) = 1 2 n c ( t ) (2.1.2.6) \mathit n_o(t)={1\over2}\mathit n_c(t)\tag{2.1.2.6} no(t)=21nc(t)(2.1.2.6)

2.DSB-SC调制器及相干解调模型

1.DSB-SC调制器

在这里插入图片描述

图2.5 DSB调制模型

m ( t ) \mathit m(t) m(t)为调制信号,其中令 s D S B ( t ) = s ( t ) \mathit s_{DSB}(t) =\mathit s(t) sDSB(t)=s(t),由图可得:
s D S B ( t ) = A c m ( t )   c o s 2 π f c t (2.2.1.1) \mathit s_{DSB}(t) =\mathit A_c\mathit m(t)\ cos2\pi\mathit f_ct\tag{2.2.1.1} sDSB(t)=Acm(t) cos2πfct(2.2.1.1)
m ( t ) \mathit m(t) m(t) s D S B ( t ) \mathit s_{DSB}(t) sDSB(t)的时域波形图如下:
在这里插入图片描述

图2.6 时域波形图

上图中: f ( t ) = m ( t ) \mathit f(t)=\mathit m(t) f(t)=m(t) c ( t ) = A c cos ⁡ 2 π f c t \mathit c(t)=\mathit A_c\cos2\pi\mathit f_ct c(t)=Accos2πfct
m ( t ) \mathit m(t) m(t)的频谱为 M ( f ) \mathit M(f) M(f),则 s D S B ( t ) \mathit s_{DSB}(t) sDSB(t)的频谱 s D S B ( f ) \mathit s_{DSB}(f) sDSB(f)为:
s D S B ( f ) = A c 2 [ M ( f − f c ) + M ( f + f c ) ] (2.2.1.2) \mathit s_{DSB}(f)={\mathit A_c\over2}[\mathit M(\mathit f-\mathit f_c)+\mathit M(\mathit f+\mathit f_c)]\tag{2.2.1.2} sDSB(f)=2Ac[M(ffc)+M(f+fc)](2.2.1.2)
m ( t ) \mathit m(t) m(t)频谱图如下,可得 s D S B ( t ) \mathit s_{DSB}(t) sDSB(t)的频谱图:
在这里插入图片描述

图2.7 频域波形图

2.叠加噪声后的DSB-SC相干解调模型

在这里插入图片描述

图2.8 DSB相干解调模型

由DSB-SC调制器得: s ( t ) = s D S B ( t ) = A c m ( t ) cos ⁡ 2 π f c t (2.2.2.1) \mathit s(t)=\mathit s_{DSB}(t) =\mathit A_c\mathit m(t)\cos2\pi\mathit f_ct\tag{2.2.2.1} s(t)=sDSB(t)=Acm(t)cos2πfct(2.2.2.1)

n i ( t ) \mathit n_i(t) ni(t)为宽带加性高斯白噪声,其均值为零,双边功率谱密度为 n 0 2 {\mathit n_0\over2} 2n0。已调信号 s ( t ) \mathit s(t) s(t)与加性噪声相加后通过接收机的带通滤波器BFT,设 m ( t ) \mathit m(t) m(t)的带宽为 B \mathit B B,带通滤波器的等效噪声带宽为 2 B 2\mathit B 2B

宽带噪声 n i ( t ) \mathit n_i(t) ni(t)通过BFT后成为窄带白高斯噪声 n ( t ) \mathit n(t) n(t)则, n ( t ) = n c ( t ) cos ⁡ 2 π f c t − n s ( t )   s i n 2 π f c t (2.2.2.2) \mathit n(t)=\mathit n_c(t) \cos2\pi\mathit f_ct- \mathit n_s(t)\ sin2\pi\mathit f_ct\tag{2.2.2.2} n(t)=nc(t)cos2πfctns(t) sin2πfct(2.2.2.2)
其中 n c ( t ) \mathit n_c(t) nc(t) n s ( t ) \mathit n_s(t) ns(t)分别是 n ( t ) \mathit n(t) n(t)的同相分量和正交分量,均值都为零,方差与 n ( t ) \mathit n(t) n(t)相同。

r ( t ) \mathit r(t) r(t)通过乘法器,与恢复载波 cos ⁡ 2 π f c t \cos2\pi\mathit f_ct cos2πfct相乘得:
输出信号: s d ( t ) = s ( t ) cos ⁡ 2 π f c t = 1 2 A c m ( t )   [ 1 + cos ⁡ 4 π f c t ] (2.2.2.3) \mathit s_d(t)=\mathit s(t) \cos2\pi\mathit f_ct={1\over2}\mathit A_c\mathit m(t)\ [1+\cos4\pi\mathit f_ct]\tag{2.2.2.3} sd(t)=s(t)cos2πfct=21Acm(t) [1+cos4πfct](2.2.2.3)
输出噪声:
n d ( t ) = n ( t ) cos ⁡ 2 π f c t = n c ( t )   c o s 2 2 π f c t − n s ( t )   s i n 2 π f c t cos ⁡ 2 π f c t = 1 2 n c ( t ) + 1 2 n c ( t ) cos ⁡ 4 π f c t − 1 2 n s ( t ) sin ⁡ 4 π f c t (2.2.2.4) \begin{aligned}\mathit n_d(t)&=\mathit n(t)\cos2\pi\mathit f_ct\\ =&\mathit n_c(t) \ cos^22\pi\mathit f_ct- \mathit n_s(t)\ sin2\pi\mathit f_ct\cos2\pi\mathit f_ct\\ =&{1\over2}\mathit n_c(t)+{1\over2}\mathit n_c(t)\cos4\pi\mathit f_ct-{1\over2}\mathit n_s(t)\sin4\pi\mathit f_ct\\ \end{aligned}\tag{2.2.2.4} nd(t)===n(t)cos2πfctnc(t) cos22πfctns(t) sin2πfctcos2πfct21nc(t)+21nc(t)cos4πfct21ns(t)sin4πfct(2.2.2.4)

再通过低通滤波器LPF,抑制二倍频分量,仅通过低通分量,从而实现相干解调,低通滤波器带宽为 B \mathit B B
得到相干解调后
输出信号: m o ( t ) = 1 2 A c m ( t ) (2.2.2.5) \mathit m_o(t)={1\over2}\mathit A_c\mathit m(t)\tag{2.2.2.5} mo(t)=21Acm(t)(2.2.2.5)输出噪声: n o ( t ) = 1 2 n c ( t ) (2.2.2.6) \mathit n_o(t)={1\over2}\mathit n_c(t)\tag{2.2.2.6} no(t)=21nc(t)(2.2.2.6)

3.SSB调制器及相干解调模型

1.SSB调制器

在这里插入图片描述

图2.9 SSB调制器

m ( t ) \mathit m(t) m(t) m ( t ) \mathit m(t) m(t)的希尔伯特变换,是把 m ( t ) \mathit m(t) m(t)的所有频率分量都相移 − π 2 -{\pi\over2} 2π后得到的信号,称 m ( t ) \mathit m(t) m(t)的正交信号。
由图可知: s ( t ) = s S S B ( t ) = A c 2 m ( t ) cos ⁡ 2 π f c t ∓ A c 2 m ( t ) ^ sin ⁡ 2 π f c t (2.3.1.1) \mathit s(t)=\mathit s_{SSB}(t) ={\mathit A_c\over2}\mathit m(t)\cos2\pi\mathit f_ct\mp{\mathit A_c\over2} \hat{m(t)}\sin2\pi\mathit f_ct\tag{2.3.1.1} s(t)=sSSB(t)=2Acm(t)cos2πfct2Acm(t)^sin2πfct(2.3.1.1)
式中:负对应上边带信号,正对应下边带信号

当利用上边带信号进行调幅时得到调制后信号:
s U S S B ( t ) = A c 2 m ( t ) cos ⁡ 2 π f c t − A c 2 m ( t ) ^ sin ⁡ 2 π f c t (2.3.1.2) \mathit s_{USSB}(t) ={\mathit A_c\over2}\mathit m(t)\cos2\pi\mathit f_ct-{\mathit A_c\over2}\hat {m(t)}\sin2\pi\mathit f_ct\tag{2.3.1.2} sUSSB(t)=2Acm(t)cos2πfct2Acm(t)^sin2πfct(2.3.1.2)
则调制后信号的频谱表达式为:
S U S S B ( f ) = A c 4 [ M ( f − f c ) + M ( f + f c ) ] + A c 4 [ M ( f − f c ) s g n ( f − f c ) + M ( f + f c ) s g n ( f − f c ) ] = A c 2 M ( f − f c ) u ( f − f c ) + A c 2 M ( f + f c ) u ( − f − f c ) (2.3.1.3) \begin{aligned}\mathit S_{USSB}(f) =&{\mathit A_c\over4}[\mathit M(\mathit f-\mathit f_c)+\mathit M(\mathit f+\mathit f_c)]+{\mathit A_c\over4}[\mathit M(\mathit f-\mathit f_c)\mathit sgn(\mathit f-\mathit f_c)+\mathit M(\mathit f+\mathit f_c)\mathit sgn(\mathit f-\mathit f_c)]\\ &={\mathit A_c\over2}\mathit M(\mathit f-\mathit f_c)\mathit u(\mathit f-\mathit f_c)+{\mathit A_c\over2}\mathit M(\mathit f+\mathit f_c)\mathit u(-\mathit f-\mathit f_c)\\ \end{aligned}\tag{2.3.1.3} SUSSB(f)=4Ac[M(ffc)+M(f+fc)]+4Ac[M(ffc)sgn(ffc)+M(f+fc)sgn(ffc)]=2AcM(ffc)u(ffc)+2AcM(f+fc)u(ffc)(2.3.1.3)
式中 M ( f ) \mathit M(f) M(f) m ( t ) \mathit m(t) m(t)频域表达式。
在这里插入图片描述

图2.10

当利用下边带信号进行调幅时得到调制信号:
s L S S B ( t ) = A c 2 m ( t ) cos ⁡ 2 π f c t + A c 2 m ( t ) ^ sin ⁡ 2 π f c t (2.3.1.4) \mathit s_{LSSB}(t) ={\mathit A_c\over2}\mathit m(t)\cos2\pi\mathit f_ct+{\mathit A_c\over2}\hat{m(t)}\sin2\pi\mathit f_ct\tag{2.3.1.4} sLSSB(t)=2Acm(t)cos2πfct+2Acm(t)^sin2πfct(2.3.1.4)
则调制后信号的频域表达式为:
S L S S B ( f ) = A c 4 [ M ( f − f c ) + M ( f + f c ) ] − A c 4 [ M ( f − f c ) s g n ( f − f c ) + M ( f + f c ) s g n ( f − f c ) ] = A c 2 M ( f − f c ) u ( f c − f ) + A c 2 M ( f + f c ) u ( f + f c ) (2.3.1.5) \begin{aligned}\mathit S_{LSSB}(f) =&{\mathit A_c\over4}[\mathit M(\mathit f-\mathit f_c)+\mathit M(\mathit f+\mathit f_c)]-{\mathit A_c\over4}[\mathit M(\mathit f-\mathit f_c)\mathit sgn(\mathit f-\mathit f_c)+\mathit M(\mathit f+\mathit f_c)\mathit sgn(\mathit f-\mathit f_c)]\\ &={\mathit A_c\over2}\mathit M(\mathit f-\mathit f_c)\mathit u(\mathit f_c-\mathit f)+{\mathit A_c\over2}\mathit M(\mathit f+\mathit f_c)\mathit u(\mathit f+\mathit f_c)\\ \end{aligned}\tag{2.3.1.5} SLSSB(f)=4Ac[M(ffc)+M(f+fc)]4Ac[M(ffc)sgn(ffc)+M(f+fc)sgn(ffc)]=2AcM(ffc)u(fcf)+2AcM(f+fc)u(f+fc)(2.3.1.5)
式中 M ( f ) \mathit M(f) M(f) m ( t ) \mathit m(t) m(t)频域表达式。
在这里插入图片描述

图2.11

2.叠加噪声后的SSB相干解调模型

在这里插入图片描述

图2.12 SSB相干解调模型

由SSB的相干调制器得: s ( t ) = s S S B ( t ) = A c 2 m ( t ) cos ⁡ 2 π f c t ∓ A c 2 m ( t ) ^ sin ⁡ 2 π f c t (2.3.2.1) \mathit s(t)=\mathit s_{SSB}(t) ={\mathit A_c\over2}\mathit m(t)\cos2\pi\mathit f_ct\mp{\mathit A_c\over2}\hat{ m(t)}\sin2\pi\mathit f_ct\tag{2.3.2.1} s(t)=sSSB(t)=2Acm(t)cos2πfct2Acm(t)^sin2πfct(2.3.2.1)

n i ( t ) \mathit n_i(t) ni(t)为宽带加性高斯白噪声,其均值为零,双边功率谱密度为 n 0 2 {\mathit n_0\over2} 2n0。已调信号 s ( t ) \mathit s(t) s(t)与加性噪声相加后通过接收机的带通滤波器BFT,设 m ( t ) \mathit m(t) m(t)的带宽为 B \mathit B B,则滤波器的等效噪声带宽为 B \mathit B B

宽带噪声 n i ( t ) \mathit n_i(t) ni(t)通过BFT后成为窄带白高斯噪声 n ( t ) \mathit n(t) n(t)则, n ( t ) = n c ( t ) cos ⁡ 2 π f c t − n s ( t )   s i n 2 π f c t (2.3.2.2) \mathit n(t)=\mathit n_c(t) \cos2\pi\mathit f_ct- \mathit n_s(t)\ sin2\pi\mathit f_ct\tag{2.3.2.2} n(t)=nc(t)cos2πfctns(t) sin2πfct(2.3.2.2)
其中 n c ( t ) \mathit n_c(t) nc(t) n s ( t ) \mathit n_s(t) ns(t)分别是 n ( t ) \mathit n(t) n(t)的同相分量和正交分量,均值都为零,方差与 n ( t ) \mathit n(t) n(t)相同。

r ( t ) \mathit r(t) r(t)通过乘法器,与恢复载波 cos ⁡ 2 π f c t \cos2\pi\mathit f_ct cos2πfct相乘得:
输出信号: s d ( t ) = s ( t ) cos ⁡ 2 π f c t = 1 4 A c m ( t )   + 1 4 A c m ( t ) cos ⁡ 4 π f c t ∓ 1 4 A c m ( t ) ^ sin ⁡ 4 π f c t (2.3.2.3) \mathit s_d(t)=\mathit s(t) \cos2\pi\mathit f_ct={1\over4}\mathit A_c\mathit m(t)\ +{1\over4}\mathit A_c\mathit m(t)\cos4\pi\mathit f_ct\mp{1\over4}\mathit A_c\hat{ m(t)}\sin4\pi\mathit f_ct\tag{2.3.2.3} sd(t)=s(t)cos2πfct=41Acm(t) +41Acm(t)cos4πfct41Acm(t)^sin4πfct(2.3.2.3)
输出噪声:
n d ( t ) = n ( t ) cos ⁡ 2 π f c t = n c ( t )   c o s 2 2 π f c t − n s ( t )   s i n 2 π f c t cos ⁡ 2 π f c t = 1 2 n c ( t ) + 1 2 n c ( t ) cos ⁡ 4 π f c t − 1 2 n s ( t ) sin ⁡ 4 π f c t (2.3.2.4) \begin{aligned}\mathit n_d(t)&=\mathit n(t)\cos2\pi\mathit f_ct\\ &=\mathit n_c(t) \ cos^22\pi\mathit f_ct- \mathit n_s(t)\ sin2\pi\mathit f_ct\cos2\pi\mathit f_ct\\ &={1\over2}\mathit n_c(t)+{1\over2}\mathit n_c(t)\cos4\pi\mathit f_ct-{1\over2}\mathit n_s(t)\sin4\pi\mathit f_ct\\ \end{aligned}\tag{2.3.2.4} nd(t)=n(t)cos2πfct=nc(t) cos22πfctns(t) sin2πfctcos2πfct=21nc(t)+21nc(t)cos4πfct21ns(t)sin4πfct(2.3.2.4)

再通过低通滤波器LPF,抑制二倍频和四倍频分量,仅通过低通分量,从而实现相干解调,低通滤波器带宽为 B \mathit B B
得到相干解调后
输出信号: s o ( t ) = 1 4 A c m ( t ) (2.3.2.5) \mathit s_o(t)={1\over4}\mathit A_c\mathit m(t)\tag{2.3.2.5} so(t)=41Acm(t)(2.3.2.5)输出噪声: n o ( t ) = 1 2 n c ( t ) (2.3.2.6) \mathit n_o(t)={1\over2}\mathit n_c(t)\tag{2.3.2.6} no(t)=21nc(t)(2.3.2.6)

3.抗干扰性能理论分析

因为在输入信噪比相同的条件下,输出信噪比越高,则说明解调器的抗噪声性能越好,故在此引用解调器调制制度增益 G = ( S N ) o u t ( S N ) i n \mathit G={({\mathit S\over \mathit N})_{out}\over({\mathit S\over \mathit N})_{in}} G=(NS)in(NS)out来比较不同调制解调系统的抗噪声性。 G \mathit G G越大,则表明抗噪声性越好。
式中 ( S N ) i n = S i n N i n ({\mathit S\over \mathit N})_{in}={\mathit S_{in}\over\mathit N_{in}} (NS)in=NinSin为输入信噪比, ( S N ) o u t = S o u t N o u t ({\mathit S\over \mathit N})_{out}={\mathit S_{out}\over\mathit N_{out}} (NS)out=NoutSout为输出信噪比。
---- S i n \mathit S_{in} Sin为解调器输入已调信号的平均功率。
---- N i n \mathit N_{in} Nin为解调器输入噪声的平均功率。
---- S o u t \mathit S_{out} Sout为解调器输出有用信号的平均功率。
---- N o u t \mathit N_{out} Nout为解调器输出噪声的平均功率。

在这里插入图片描述

图3.1

图中,黄色虚线框内代表解调器。
在下面三种调制解调系统的抗噪声性分析中:

---- s ( t ) \mathit s(t) s(t)为经过调制器后的信号。
----- n i ( t ) \mathit n_i(t) ni(t)为宽带高斯白噪声,均值为0,双边功率谱密度为 n 0 2 {\mathit n_0\over2} 2n0
---- P m \mathit P_m Pm表示调制信号 m ( t ) \mathit m(t) m(t)的平均功率。
---- B \mathit B B表示调制信号 m ( t ) \mathit m(t) m(t)的带宽。
---- n ( t ) \mathit n(t) n(t) n i ( t ) \mathit n_i(t) ni(t)通过带通滤波器得到的窄带高斯白噪声。
---- m o ( t ) \mathit m_o(t) mo(t)为解调器输出的有用信号。
---- n o ( t ) \mathit n_o(t) no(t)为解调器输出的噪声信号。

1.AM相干解调抗噪声性能

S i n = ⟨ s 2 ( t ) ⟩ = A c 2 2 ( 1 + P m ) (3.1.1) \mathit S_{in}=\langle \mathit s^2(t)\rangle={\mathit A_c^2\over2}(1+\mathit P_m)\tag{3.1.1} Sin=s2(t)=2Ac2(1+Pm)(3.1.1)

N i n = E [ n 2 ( t ) ] = n 0 B A M = 2 n 0 B (3.1.2) \mathit N_{in}=\mathit E[ \mathit n^2(t)]=\mathit n_0\mathit B_{AM}=2\mathit n_0\mathit B\tag{3.1.2} Nin=E[n2(t)]=n0BAM=2n0B(3.1.2)

( S N ) i n = S i n N i n = A c 2 2 ( 1 + P m ) 2 n 0 B (3.1.3) ({\mathit S\over \mathit N})_{in}={\mathit S_{in}\over\mathit N_{in}}={{\mathit A_c^2\over2}(1+\mathit P_m)\over2\mathit n_0\mathit B }\tag{3.1.3} (NS)in=NinSin=2n0B2Ac2(1+Pm)(3.1.3)

S o u t = ⟨ s o 2 ( t ) ⟩ = A c 2 P m 4 (3.1.4) \mathit S_{out}=\langle \mathit s_o^2(t)\rangle={\mathit A_c^2\mathit P_m\over4}\tag{3.1.4} Sout=so2(t)=4Ac2Pm(3.1.4)

N o u t = E [ n o 2 ( t ) ] = 1 4 E [ n c 2 ( t ) ] = 1 4 N i n = 1 2 n 0 B (3.1.5) \mathit N_{out}=\mathit E[ \mathit n_o^2(t)]={1\over4}\mathit E[ \mathit n_c^2(t)]={1\over4}\mathit N_{in}={1\over2}\mathit n_0\mathit B\tag{3.1.5} Nout=E[no2(t)]=41E[nc2(t)]=41Nin=21n0B(3.1.5)

( S N ) o u t = S o u t N o u t = A c 2 P m 2 n 0 B (3.1.6) ({\mathit S\over \mathit N})_{out}={\mathit S_{out}\over\mathit N_{out}}={\mathit A_c^2\mathit P_m\over2\mathit n_0\mathit B}\tag{3.1.6} (NS)out=NoutSout=2n0BAc2Pm(3.1.6)

G A M = ( S N ) o u t ( S N ) i n = 2 P m 1 + P m (3.1.7) \mathit G_{AM}={({\mathit S\over \mathit N})_{out}\over({\mathit S\over \mathit N})_{in}}={2\mathit P_m\over1+\mathit P_m}\tag{3.1.7} GAM=(NS)in(NS)out=1+Pm2Pm(3.1.7)

其中 B A M = 2 B \mathit B_{AM}=2\mathit B BAM=2B为窄带高斯白噪声的带宽。

2.DSB-SC相干解调抗噪声性能

S i n = ⟨ s 2 ( t ) ⟩ = A c 2 P m 2 (3.2.1) \mathit S_{in}=\langle \mathit s^2(t)\rangle={\mathit A_c^2\mathit P_m\over2}\tag{3.2.1} Sin=s2(t)=2Ac2Pm(3.2.1)

N i n = E [ n 2 ( t ) ] = n 0 B D S B − S C = 2 n 0 B (3.2.2) \mathit N_{in}=\mathit E[ \mathit n^2(t)]=\mathit n_0\mathit B_{DSB-SC}=2\mathit n_0\mathit B\tag{3.2.2} Nin=E[n2(t)]=n0BDSBSC=2n0B(3.2.2)

( S N ) i n = S i n N i n = A c 2 P m 4 n 0 B (3.2.3) ({\mathit S\over \mathit N})_{in}={\mathit S_{in}\over\mathit N_{in}}={\mathit A_c^2\mathit P_m\over4\mathit n_0\mathit B }\tag{3.2.3} (NS)in=NinSin=4n0BAc2Pm(3.2.3)

S o u t = ⟨ s o 2 ( t ) ⟩ = A c 2 P m 4 (3.2.4) \mathit S_{out}=\langle \mathit s_o^2(t)\rangle={\mathit A_c^2\mathit P_m\over4}\tag{3.2.4} Sout=so2(t)=4Ac2Pm(3.2.4)

N o u t = E [ n o 2 ( t ) ] = 1 4 E [ n c 2 ( t ) ] = 1 4 N i n = 1 2 n 0 B (3.2.5) \mathit N_{out}=\mathit E[ \mathit n_o^2(t)]={1\over4}\mathit E[ \mathit n_c^2(t)]={1\over4}\mathit N_{in}={1\over2}\mathit n_0\mathit B\tag{3.2.5} Nout=E[no2(t)]=41E[nc2(t)]=41Nin=21n0B(3.2.5)

( S N ) o u t = S o u t N o u t = A c 2 P m 2 n 0 B (3.2.6) ({\mathit S\over \mathit N})_{out}={\mathit S_{out}\over\mathit N_{out}}={\mathit A_c^2\mathit P_m\over2\mathit n_0\mathit B}\tag{3.2.6} (NS)out=NoutSout=2n0BAc2Pm(3.2.6)

G D S B − S C = ( S N ) o u t ( S N ) i n = 2 (3.2.7) \mathit G_{DSB-SC}={({\mathit S\over \mathit N})_{out}\over({\mathit S\over \mathit N})_{in}}=2\tag{3.2.7} GDSBSC=(NS)in(NS)out=2(3.2.7)

其中 B D S B = 2 B \mathit B_{DSB}=2\mathit B BDSB=2B为窄带高斯白噪声的带宽。

3.SSB相干解调抗噪声性能

S i n = ⟨ s 2 ( t ) ⟩ = A c 2 P m 4 (3.3.1) \mathit S_{in}=\langle \mathit s^2(t)\rangle={\mathit A_c^2\mathit P_m\over4}\tag{3.3.1} Sin=s2(t)=4Ac2Pm(3.3.1)

N i n = E [ n 2 ( t ) ] = n 0 B S S B = n 0 B (3.3.2) \mathit N_{in}=\mathit E[ \mathit n^2(t)]=\mathit n_0\mathit B_{SSB}=\mathit n_0\mathit B\tag{3.3.2} Nin=E[n2(t)]=n0BSSB=n0B(3.3.2)

( S N ) i n = S i n N i n = A c 2 P m 4 n 0 B (3.3.2) ({\mathit S\over \mathit N})_{in}={\mathit S_{in}\over\mathit N_{in}}={\mathit A_c^2\mathit P_m\over4\mathit n_0\mathit B }\tag{3.3.2} (NS)in=NinSin=4n0BAc2Pm(3.3.2)

S o u t = ⟨ s o 2 ( t ) ⟩ = A c 2 P m 16 (3.3.4) \mathit S_{out}=\langle \mathit s_o^2(t)\rangle={\mathit A_c^2\mathit P_m\over16}\tag{3.3.4} Sout=so2(t)=16Ac2Pm(3.3.4)

N o u t = E [ n o 2 ( t ) ] = 1 4 E [ n c 2 ( t ) ] = 1 4 N i n = 1 4 n 0 B (3.3.5) \mathit N_{out}=\mathit E[ \mathit n_o^2(t)]={1\over4}\mathit E[ \mathit n_c^2(t)]={1\over4}\mathit N_{in}={1\over4}\mathit n_0\mathit B\tag{3.3.5} Nout=E[no2(t)]=41E[nc2(t)]=41Nin=41n0B(3.3.5)

( S N ) o u t = S o u t N o u t = A c 2 P m 4 n 0 B (3.3.6) ({\mathit S\over \mathit N})_{out}={\mathit S_{out}\over\mathit N_{out}}={\mathit A_c^2\mathit P_m\over4\mathit n_0\mathit B}\tag{3.3.6} (NS)out=NoutSout=4n0BAc2Pm(3.3.6)

G S S B = ( S N ) o u t ( S N ) i n = 1 (3.3.7) \mathit G_{SSB}={({\mathit S\over \mathit N})_{out}\over({\mathit S\over \mathit N})_{in}}=1\tag{3.3.7} GSSB=(NS)in(NS)out=1(3.3.7)

其中 B S S B = B \mathit B_{SSB}=\mathit B BSSB=B为窄带高斯白噪声的带宽。

4.三种幅度调制解调系统抗噪声性能理论比较分析

DSB调制系统的信噪比增益为2,则说明DSB的解调器使信噪比改善一倍,这是因为采用相干解调把噪声中的正交分量抑制掉了。SSB调制系统的信噪比增益为1,说明SSB的解调器没有改善信噪比,因为在SSB的系统中,信号与噪声有相同的表达式,所以在相干解调时,信号与噪声的正交分量均被抑制掉了。
虽然 G S S B = 1 \mathit G_{SSB}=1 GSSB=1 G D S B − S C = 2 \mathit G_{DSB-SC}=2 GDSBSC=2但我们并不能说DSB系统的抗噪声性更好,因为在解调器输入端输入信号的带宽条件和输入噪声的带宽条件不同导致其功率不同,若能保证相应条件都相同,我们会发现两者的抗噪声性能是一样的,但是SSB所需传输带宽仅是DSB的一半。
G A M = 2 P m 1 + P m \mathit G_{AM}={2\mathit P_m\over1+\mathit P_m} GAM=1+Pm2Pm总是小于1,这说明AM解调器把输入信噪比恶化了,故其抗噪声性能差于SSB系统和DSB系统。
综上可得,抗噪声性能:SSB=DSB>AM

4.仿真实现与仿真结果

1.AM相干解调仿真及结果

echo off
close all,clear
fs=200;                 %采样频率,由采样定理决定
ts=1/fs;                %采样间隔ts由fs决定
fc=50;                  %载波频率
Ac=100;                 %载波峰值
fm=5;                   %调制信号频率
snr=20;                 %dB表示的信噪比 
snr_lin=10^(snr/10);    %线性信噪比的值
t0=0.5;
t=[0:ts:t0];            %时间向量的表达,时间范围[0,t0]
df0=0.2;                %频率分辨率的下限为0.2

%...........调制信号m(t)的时域描述...........
m=cos(2*pi*fm*t);

%........... 第一部分:调制 ...........
%正弦载波信号c(t)的时域表达
c=Ac*cos(2*pi*fc*t);

%DSB已调信号u(t)的时域表达
A0=1;
s=(A0+m).*c;

%求出消息信号和已调信号的傅立叶变换
N1=length(m);
N2=fs/df0;
N=2^(nextpow2(max(N1, N2)));    %求fft的点数N
df=fs/N;                        %最终的频率分辨率
f=linspace(-fs/2,fs/2-df,N);    %频域轴的向量,0频点在中间 

M=fftshift(fft(m,N)/fs);        %求消息信号的傅立叶变换
S=fftshift(fft(s,N)/fs);        %求已调信号的傅立叶变换  

%........... 第二部分: 加噪声 ...........
r=awgn(s,snr);           %s(t)中加高斯白噪声,噪声信噪比为20
noise=r-s;               %noise为噪声
R=fftshift(fft(r,N)/fs);        %求消息信号的傅立叶变换

在这里插入图片描述

图4.1
%........... 第三部分: 相干解调 ...........
sd=s.*cos(2*pi*fc*t)-0.5*Ac; 
noise_d=noise.*cos(2*pi*fc*t);
y=r.*cos(2*pi*fc*t)-0.5*Ac;

%求截止频率为20Hz的低通滤波器的H(f)
fL=20;                                        
H=zeros(1,N);
num1=round((fs/2-fL)/df)+1;
num2=round((fs/2+fL)/df)+1;
H(num1:num2)=ones(1,num2-num1+1);
 
%y(t)经过理想低通的输出
Y=fftshift(fft(y,N)/fs);        %y(t)的傅立叶变换
Y=Y.*H;                         %y(t)过低通滤波器后的傅立叶变换
y=ifft(ifftshift(Y)*fs);        %y(t)的时域表达
%sd(t)经过理想低通的输出
Sd=fftshift(fft(sd,N)/fs);        %sd(t)的傅立叶变换
Mo=Sd.*H;                         %sd(t)过低通滤波器后的傅立叶变换
mo=ifft(ifftshift(Mo)*fs);        %mo(t)的时域表达
 
%noise_d(t)经过理想低通的输出
Noise_d=fftshift(fft(y,N)/fs);        %noise_d(t)的傅立叶变换
Noise_o=Noise_d.*H;                         %noise_d(t)过低通滤波器后的傅立叶变换
noise_o=ifft(ifftshift(Noise_o)*fs);        %noise_o(t)的时域表达
 
 

在这里插入图片描述

图4.2
signal_o_power=(sum(mo.^2))/length(mo);       %调制器输出信号的平均功率
noise_o_power=(sum(noise_o.^2))/length(noise_o); %调制器输出噪声的平均功率
snr_in=snr_lin;
snr_out=signal_o_power/noise_o_power;
G=snr_out/snr_in;

得仿真结果:
在这里插入图片描述

图4.3

2.DSB相干解调仿真及结果

%........... 第一部分:调制 ...........
c=Ac*cos(2*pi*fc*t); 	%正弦载波信号c(t)的时域表达
s=m.*c; 		%DSB已调信号u(t)的时域表达
%........... 第二部分: 加噪声 ...........
signal_power=(sum(s.^2))/length(s);        %%%%%abs计算(已调)信号的功率或者signal_power=(norm(s)^2)/length(s);
noise_power=signal_power/snr_lin;   %计算噪声的功率
noise_std=sqrt(noise_power);        %计算噪声的标准差
noise=noise_std*randn(1,length(s));        %%%%%%得到噪声向量
r=s+noise;        %接收端的接收信号,加入了噪声
R=fftshift(fft(r,N)/fs);        %r(t)的傅立叶变换 

此次加噪声选择另一种方法。其余代码与AM相干解调代码差不多,在此不再显示。
在这里插入图片描述

图4.4

在这里插入图片描述

图4.5

在这里插入图片描述

图4.6

3.SSB相干解调仿真及结果

%........... 第一部分:调制 ...........
c=Ac*cos(2*pi*fc*t);               %正弦载波信号c(t)的时域表达,分sin,cos
si=Ac*sin(2*pi*fc*t);
m_hilbert=imag(hilbert(m));     %消息信号m(t)的希尔波特变换,产生m^(t)
s=0.5*m.*c+0.5*m_hilbert.*si;               %%%%%%USSB已调信号u(t)的时域表达,上边带,下边带

其余代码与AM相干解调代码差不多,在此不再显示。
在这里插入图片描述

图4.7

在这里插入图片描述

图4.8

在这里插入图片描述

图4.9

5.小结

本文中我们进行理论分析得出下图:
在这里插入图片描述

图5.1

但在进行仿真时,得到的DSB、SSB系统的信噪比增益都比理论值要大。因为影响信噪比增益的因素有很多,在本文仿真时,三种系统下输入的调制信号相同,即调制信号的功率谱密度相同,故误差应该是叠加的高斯白噪声产生的。由理论值和实验值以及理论分析可得,DSB与SSB的抗噪声性相同,均优于AM的抗噪声性,但是SSB的传输带宽只有DSB的一半,故SSB得到普遍应用。

6.参考文献

[1] 第三章模拟调制系统-抗噪声性能
[2] 模拟幅度调制系统仿真
[3] 模拟幅度调制系统抗干扰性能仿真分析[模板]
[4]幅度调制系统的抗噪声性能

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值