模拟幅度调制相干解调系统抗噪声性能仿真分析
一、引言
调制在通信系统中的作用至关重要。在大多数场合中,调制通常指载波调制,就是用调制信号去控制在播的参数的过程,使载波的某一个或某几个参数按照调制信号的规律而变化。本文介绍了常见的三种模拟幅度调制器及其相干解调系统,分析了它们的抗噪声性能,并对其进行了仿真分析,得出相关结论。
二、系统基础知识
1.乘法运算
运算表达式为
y
(
t
)
=
g
1
(
t
)
g
2
(
t
)
y(t)=g_1(t)g_2(t)
y(t)=g1(t)g2(t)
对于冲激函数
δ
\delta
δ(t),有
y
(
t
)
=
g
(
t
)
δ
(
t
)
=
g
(
0
)
δ
(
t
)
;
y(t)=g(t)\delta(t)=g(0)\delta(t);
y(t)=g(t)δ(t)=g(0)δ(t);
y ( t ) = g ( t ) δ ( t − t 0 ) = g ( t 0 ) δ ( t − t 0 ) ; y(t)=g(t)\delta(t-t_0)=g(t_0)\delta(t-t_0); y(t)=g(t)δ(t−t0)=g(t0)δ(t−t0);
2.卷积运算
运算表达式为:
y
(
t
)
=
g
1
(
t
)
∗
g
2
(
t
)
=
∫
∞
∞
g
1
(
t
)
g
2
(
t
−
τ
)
d
τ
y(t)=g_1(t)\ast g_2(t)=\int_ \infty^\infty g_1(t)g_2(t-\tau)d\tau
y(t)=g1(t)∗g2(t)=∫∞∞g1(t)g2(t−τ)dτ
对于冲激函数,有
y
(
t
)
=
g
(
t
)
∗
δ
(
t
)
=
g
(
t
)
;
y(t)=g(t)\ast\delta(t)=g(t);
y(t)=g(t)∗δ(t)=g(t);
y ( t ) = g ( t ) ∗ δ ( t − t 0 ) = g ( t − t 0 ) ; y(t)=g(t)\ast\delta(t-t_0)=g(t-t_0); y(t)=g(t)∗δ(t−t0)=g(t−t0);
3.确定信号通过滤波器
如图所示,确定信号通过滤波器相当于该信号x(t)与滤波器的单位冲激响应h(t)卷积,信号在时域的卷积等于信号在频域的乘积。
4.确定信号通过乘法器
如图所示,确定信号通过乘法器相当于该信号x(t)与载波信号c(t)相乘,信号在时域的乘积等于信号在频域的卷积.
三、模拟幅度调制相干解调系统
1.常规幅度调制(AM)与相干解调
调幅指数:
β
A
M
=
m
(
t
)
∣
m
a
x
\beta _{AM}=m(t)|_{max}
βAM=m(t)∣max
如图所示,图中各点信号时域表达式如下:
g
(
t
)
=
m
(
t
)
+
1
g(t)=m(t)+1
g(t)=m(t)+1
S A M ( t ) = g ( t ) c ( t ) = A c [ m ( t ) + 1 ] c o s 2 π f c t S_{AM}(t)=g(t)c(t)=A_c[m(t)+1]cos2\pi f_ct SAM(t)=g(t)c(t)=Ac[m(t)+1]cos2πfct
S ( t ) = S A M ( t ) = A c [ m ( t ) + 1 ] c o s 2 π f c t S(t)=S_{AM}(t)=A_c[m(t)+1]cos2\pi f_ct S(t)=SAM(t)=Ac[m(t)+1]cos2πfct
S d ( t ) = 1 2 A c [ m ( t ) + 1 ] [ 1 + c o s 4 π f c t ] S_d(t)=\frac{1}{2}A_c[m(t)+1][1+cos4\pi f_ct] Sd(t)=21Ac[m(t)+1][1+cos4πfct]
S
o
(
t
)
=
1
2
A
c
m
(
t
)
S_o(t)=\frac{1}{2}A_cm(t)
So(t)=21Acm(t)
对应的频域表达式分别为:
S
A
M
(
f
)
=
A
c
2
[
δ
(
f
−
f
c
)
+
δ
(
f
+
f
c
)
]
+
A
c
2
[
M
(
f
−
f
c
)
+
M
(
f
+
f
c
)
]
S_{AM}(f)=\frac{A_c}{2}[\delta(f-f_c)+\delta(f+f_c)]+\frac{A_c}{2}[M(f-f_c)+M(f+f_c)]
SAM(f)=2Ac[δ(f−fc)+δ(f+fc)]+2Ac[M(f−fc)+M(f+fc)]
S ( f ) = S A M ( f ) S(f)=S_{AM}(f) S(f)=SAM(f)
S
o
(
f
)
=
A
c
2
M
(
f
)
S_o(f)=\frac{A_c}{2}M(f)
So(f)=2AcM(f)
图中各点噪声时域表达式如下:
n
i
(
t
)
=
n
c
(
t
)
c
o
s
2
π
f
c
t
−
n
s
(
t
)
s
i
n
2
π
f
c
t
n_i(t)=n_c(t)cos2\pi f_ct-n_s(t)sin2\pi f_ct
ni(t)=nc(t)cos2πfct−ns(t)sin2πfct
n d ( t ) = 1 2 n c ( t ) + 1 2 n c ( t ) c o s 4 π f c ( t ) − 1 2 n s ( t ) s i n 4 π f c t n_d(t)=\frac{1}{2}n_c(t)+\frac{1}{2}n_c(t)cos4\pi f_c(t)-\frac{1}{2}n_s(t)sin4\pi f_ct nd(t)=21nc(t)+21nc(t)cos4πfc(t)−21ns(t)sin4πfct
n
o
(
t
)
=
1
2
n
c
(
t
)
n_o(t)=\frac{1}{2}n_c(t)
no(t)=21nc(t)
调制信号归一化平均功率:
S
A
M
2
(
t
)
‾
=
A
c
2
[
1
+
m
(
t
)
]
2
c
o
s
2
2
π
f
c
t
‾
=
1
2
A
c
2
+
1
2
A
c
2
m
2
(
t
)
‾
\overline{S_{AM}^2(t)}=\overline{A_c^2[1+m(t)]^2cos^22\pi f_ct}=\frac{1}{2}A_c^2+\frac{1}{2}A_c^2\overline{m^2(t)}
SAM2(t)=Ac2[1+m(t)]2cos22πfct=21Ac2+21Ac2m2(t)
那么就有
P
A
M
=
P
c
+
P
m
=
1
2
A
c
2
+
1
2
A
c
2
m
2
(
t
)
‾
P_{AM}=P_c+P_m=\frac{1}{2}A_c^2+\frac{1}{2}A_c^2\overline{m^2(t)}
PAM=Pc+Pm=21Ac2+21Ac2m2(t)
调制效率
η
A
M
=
P
m
P
A
M
=
m
2
(
t
)
‾
1
+
m
2
(
t
)
‾
\eta_{AM}=\frac{P_m}{P_{AM}}=\frac{\overline{m^2(t)}}{1+\overline{m^2(t)}}
ηAM=PAMPm=1+m2(t)m2(t)
AM相干解调抗噪声性能分析
假设白噪声的单边功率谱密度为
n
0
n_0
n0,基带信号
m
(
t
)
m(t)
m(t)的带宽为
B
B
B,那么有
输入信噪比
S
i
n
=
S
2
(
t
)
‾
=
S
A
M
2
(
t
)
‾
S_{in}=\overline{S^2(t)}=\overline{S_{AM}^2(t)}
Sin=S2(t)=SAM2(t)
N i n = E [ n 2 ( t ) ] = 2 n 0 B N_{in}=E[n^2(t)]=2n_0B Nin=E[n2(t)]=2n0B
(
S
N
)
i
n
=
S
i
n
N
i
n
=
A
c
2
2
[
1
+
P
m
]
2
n
0
B
(\frac{S}{N})_{in}=\frac{S_{in}}{N_{in}}=\frac{\frac{A_c^2}{2}[1+P_m]}{2n_0B}
(NS)in=NinSin=2n0B2Ac2[1+Pm]
输出信噪比
S
o
u
t
=
S
o
2
(
t
)
‾
=
A
c
2
P
m
4
S_{out}=\overline{S_o^2(t)}=\frac{A_c^2P_m}{4}
Sout=So2(t)=4Ac2Pm
N o u t = E [ n o 2 ( t ) ] = 1 2 n 0 B N_{out}=E[n_o^2(t)]=\frac{1}{2}n_0B Nout=E[no2(t)]=21n0B
( S N ) o u t = S o u t N o u t = A c 2 P m 2 n 0 B (\frac{S}{N})_{out}=\frac{S_{out}}{N_{out}}=\frac{A_c^2P_m}{2n_0B} (NS)out=NoutSout=2n0BAc2Pm
2.抑制载波双边带调制(DSB-SC)与相干解调
如图所示,图中各点信号时域表达式为
S
D
S
B
(
t
)
=
m
(
t
)
c
(
t
)
=
A
c
m
(
t
)
c
o
s
2
π
f
c
t
S_{DSB}(t)=m(t)c(t)=A_cm(t)cos2\pi f_ct
SDSB(t)=m(t)c(t)=Acm(t)cos2πfct
S ( t ) = S D S B ( t ) S(t)=S_{DSB}(t) S(t)=SDSB(t)
S d ( t ) = A c m ( t ) c o s 2 2 π f c t S_d(t)=A_cm(t)cos^22\pi f_ct Sd(t)=Acm(t)cos22πfct
S
o
(
t
)
=
A
c
2
m
(
t
)
S_o(t)=\frac{A_c}{2}m(t)
So(t)=2Acm(t)
对应的频域表达式为
S
D
S
B
(
f
)
=
A
c
2
[
M
(
f
−
f
c
)
+
M
(
f
+
f
c
)
]
S_{DSB}(f)=\frac{A_c}{2}[M(f-f_c)+M(f+f_c)]
SDSB(f)=2Ac[M(f−fc)+M(f+fc)]
S d ( f ) = A c 4 [ M ( f − 2 f c ) + M ( f + 2 f c ) ] + A c 2 M ( f ) S_d(f)=\frac{A_c}{4}[M(f-2f_c)+M(f+2f_c)]+\frac{A_c}{2}M(f) Sd(f)=4Ac[M(f−2fc)+M(f+2fc)]+2AcM(f)
S
o
(
f
)
=
A
c
2
M
(
f
)
S_o(f)=\frac{A_c}{2}M(f)
So(f)=2AcM(f)
各点噪声信号表达式为
n
(
t
)
=
n
c
(
t
)
c
o
s
2
π
f
c
t
−
n
s
(
t
)
s
i
n
2
π
f
c
t
n(t)=n_c(t)cos2\pi f_ct-n_s(t)sin2\pi f_ct
n(t)=nc(t)cos2πfct−ns(t)sin2πfct
n d ( t ) = 1 2 n c t + 1 2 n c ( t ) c o s 4 π f c t − 1 2 n s ( t ) s i n 4 π f c t n_d(t)=\frac{1}{2}n_ct+\frac{1}{2}n_c(t)cos4\pi f_ct-\frac{1}{2}n_s(t)sin4\pi f_ct nd(t)=21nct+21nc(t)cos4πfct−21ns(t)sin4πfct
n
o
(
t
)
=
1
2
n
c
t
n_o(t)=\frac{1}{2}n_ct
no(t)=21nct
调制信号归一化平均功率为
P
D
S
B
=
S
D
S
B
2
(
t
)
‾
=
A
c
2
2
m
2
(
t
)
‾
=
A
c
2
2
P
m
P_{DSB}=\overline{S_{DSB}^2(t)}=\frac{A_c^2}{2}\overline{m^2(t)}=\frac{A_c^2}{2}P_m
PDSB=SDSB2(t)=2Ac2m2(t)=2Ac2Pm
DSB-SC相干解调抗噪声性能分析
假设白噪声的单边功率谱密度为
n
0
n_0
n0,基带信号
m
(
t
)
m(t)
m(t)的带宽为
B
B
B,那么有
输入信噪比
S
i
n
=
S
2
(
t
)
‾
=
A
c
2
2
m
2
(
t
)
‾
=
A
c
2
2
P
m
S_{in}=\overline{S^2(t)}=\frac{A_c^2}{2}\overline{m^2(t)}=\frac{A_c^2}{2}P_m
Sin=S2(t)=2Ac2m2(t)=2Ac2Pm
N i n = E [ n 2 ( t ) ] = 2 n 0 B N_{in}=E[n^2(t)]=2n_0B Nin=E[n2(t)]=2n0B
(
S
N
)
i
n
=
S
i
n
N
i
n
=
A
c
2
P
m
4
n
0
B
(\frac{S}{N})_{in}=\frac{S_{in}}{N_{in}}=\frac{A_c^2P_m}{4n_0B}
(NS)in=NinSin=4n0BAc2Pm
输出信噪比
S
o
u
t
=
E
[
S
o
2
(
t
)
]
=
A
c
2
P
m
4
S_{out}=E[S_o^2(t)]=\frac{A_c^2P_m}{4}
Sout=E[So2(t)]=4Ac2Pm
N o u t = E [ n o 2 ( t ) ] = 1 2 n 0 B N_{out}=E[n_o^2(t)]=\frac{1}{2}n_0B Nout=E[no2(t)]=21n0B
( S N ) o u t = S o u t N o u t = A c 2 P m 2 n 0 B (\frac{S}{N})_{out}=\frac{S_{out}}{N_{out}}=\frac{A_c^2P_m}{2n_0B} (NS)out=NoutSout=2n0BAc2Pm
3.单边带调制(SSB)与相干解调
(1)滤波法产生SSB信号
如图所示,先产生一个双边带信号,然后让其通过一个边带滤波器,滤除不要的边带,即可得到单边带信号。
H
S
S
B
(
f
)
H_{SSB}(f)
HSSB(f)为单边带滤波器的传输函数,若它具有如下理想高通特性:
H
S
S
B
(
f
)
=
H
U
S
S
B
(
f
)
=
{
1
∣
f
∣
>
f
c
0
∣
f
∣
≤
f
c
H_{SSB}(f)=H_{USSB}(f)= \left\{ \begin{aligned} 1& \ \ \ \ |f|>f_c \\ 0& \ \ \ \ |f|\leq f_c \\ \end{aligned} \right.
HSSB(f)=HUSSB(f)={10 ∣f∣>fc ∣f∣≤fc
则可滤除下边带,保留上边带(USSB),r如图6为滤波法产生USSB信号。
若它具有如下理想低通特性:
H
S
S
B
(
f
)
=
H
L
S
S
B
(
f
)
=
{
1
∣
f
∣
<
f
c
0
∣
f
∣
≥
f
c
H_{SSB}(f)=H_{LSSB}(f)= \left\{ \begin{aligned} 1& \ \ \ \ |f|<f_c \\ 0& \ \ \ \ |f|\geq f_c \\ \end{aligned} \right.
HSSB(f)=HLSSB(f)={10 ∣f∣<fc ∣f∣≥fc
则可滤除上边带,保留下边带(LSSB),如图为滤波法产生LSSB信号。
调制信号
S
(
t
)
S(t)
S(t)的频谱表示为
S
(
f
)
=
A
c
2
[
M
(
f
−
f
c
)
+
M
(
f
+
f
c
)
]
H
S
S
B
(
f
)
S(f)=\frac{A_c}{2}[M(f-f_c)+M(f+f_c)]H_{SSB}(f)
S(f)=2Ac[M(f−fc)+M(f+fc)]HSSB(f)
(2)相移法产生SSB信号
相移法产生LSSB信号:
S
D
S
B
(
f
)
=
1
2
[
M
(
f
−
f
c
)
+
M
(
f
+
f
c
)
]
S_{DSB}(f)=\frac{1}{2}[M(f-f_c)+M(f+f_c)]
SDSB(f)=21[M(f−fc)+M(f+fc)]
H L S S B ( f ) = 1 2 [ s g n ( f + f c ) − s g n ( f − f c ) ] H_{LSSB}(f)=\frac{1}{2}[sgn(f+f_c)-sgn(f-f_c)] HLSSB(f)=21[sgn(f+fc)−sgn(f−fc)]
S L S S B ( f ) = 1 2 [ M ( f − f c ) + M ( f + f c ) ] H L S S B ( f ) S_{LSSB}(f)=\frac{1}{2}[M(f-f_c)+M(f+f_c)]H_{LSSB}(f) SLSSB(f)=21[M(f−fc)+M(f+fc)]HLSSB(f)
S
L
S
S
B
(
t
)
=
1
2
[
m
(
t
)
c
o
s
2
π
f
c
t
+
m
^
(
t
)
s
i
n
2
π
f
c
t
]
S_{LSSB}(t)=\frac{1}{2}[m(t)cos2\pi f_ct+\hat m(t)sin2\pi f_ct]
SLSSB(t)=21[m(t)cos2πfct+m^(t)sin2πfct]
相移法产生USSB信号:
H
U
S
S
B
(
f
)
=
1
−
H
L
S
S
B
(
f
)
H_{USSB}(f)=1-H_{LSSB}(f)
HUSSB(f)=1−HLSSB(f)
S U S S B ( f ) = S D S B ( f ) H U S S B ( f ) = S D S B ( f ) − S L S S B ( f ) S_{USSB}(f)=S_{DSB}(f)H_{USSB}(f)=S_{DSB}(f)-S_{LSSB}(f) SUSSB(f)=SDSB(f)HUSSB(f)=SDSB(f)−SLSSB(f)
S
U
S
S
B
(
t
)
=
S
D
S
B
(
t
)
−
S
L
S
S
B
(
t
)
=
1
2
[
m
(
t
)
c
o
s
2
π
f
c
t
−
m
^
(
t
)
s
i
n
2
π
f
c
t
]
S_{USSB}(t)=S_{DSB}(t)-S_{LSSB}(t)=\frac{1}{2}[m(t)cos2\pi f_ct-\hat m(t)sin2\pi f_ct]
SUSSB(t)=SDSB(t)−SLSSB(t)=21[m(t)cos2πfct−m^(t)sin2πfct]
如下图所示为相移法模型:
(3)SSB信号相干解调
如图所示为相干解调模型,图中各点信号时域表达式为:
S
(
t
)
=
A
c
2
m
(
t
)
c
o
s
2
π
f
c
t
±
A
c
2
m
^
(
t
)
s
i
n
2
π
f
c
t
S(t)=\frac{A_c}{2}m(t)cos2\pi f_ct\pm\frac{A_c}{2}\hat m(t)sin2\pi f_ct
S(t)=2Acm(t)cos2πfct±2Acm^(t)sin2πfct
S d ( t ) = A c 4 m ( t ) + A c 4 m ( t ) c o s 4 π f c t ± A c 4 m ^ ( t ) s i n 4 π f c t S_d(t)=\frac{A_c}{4}m(t)+\frac{A_c}{4}m(t)cos4\pi f_ct\pm\frac{A_c}{4}\hat m(t)sin4\pi f_ct Sd(t)=4Acm(t)+4Acm(t)cos4πfct±4Acm^(t)sin4πfct
S
o
(
t
)
=
A
c
4
m
(
t
)
S_o(t)=\frac{A_c}{4}m(t)
So(t)=4Acm(t)
各点噪声时域表达式为:
n
(
t
)
=
n
c
(
t
)
c
o
s
2
π
f
c
t
−
n
s
(
t
)
s
i
n
2
π
f
c
t
n(t)=n_c(t)cos2\pi f_ct-n_s(t)sin2\pi f_ct
n(t)=nc(t)cos2πfct−ns(t)sin2πfct
n d ( t ) = 1 2 n c ( t ) + 1 2 n c ( t ) c o s 4 π f c t − 1 2 n s ( t ) s i n 4 π f c t n_d(t)=\frac{1}{2}n_c(t)+\frac{1}{2}n_c(t)cos4\pi f_ct-\frac{1}{2}n_s(t)sin4\pi f_ct nd(t)=21nc(t)+21nc(t)cos4πfct−21ns(t)sin4πfct
n o ( t ) = 1 2 n c ( t ) n_o(t)=\frac{1}{2}n_c(t) no(t)=21nc(t)
SSB相干解调抗噪声性能分析
假设白噪声的单边功率谱密度为
n
0
n_0
n0,基带信号
m
(
t
)
m(t)
m(t)的带宽为
B
B
B,那么有
输入信噪比
S
i
n
=
S
2
(
t
)
‾
=
A
c
2
P
m
4
S_{in}=\overline{S^2(t)}=\frac{A_c^2P_m}{4}
Sin=S2(t)=4Ac2Pm
N i n = E [ n 2 ( t ) ] = n 0 B N_{in}=E[n^2(t)]=n_0B Nin=E[n2(t)]=n0B
(
S
N
)
i
n
=
S
i
n
N
i
n
=
A
c
2
P
m
4
n
0
B
(\frac{S}{N})_{in}=\frac{S_{in}}{N_{in}}=\frac{A^2_cP_m}{4n_0B}
(NS)in=NinSin=4n0BAc2Pm
输出信噪比
S
o
u
t
=
S
o
2
(
t
)
‾
=
A
c
2
P
m
16
S_{out}=\overline{S_o^2(t)}=\frac{A_c^2P_m}{16}
Sout=So2(t)=16Ac2Pm
N o u t = E [ n o 2 ( t ) ] = n 0 B 4 N_{out}=E[n_o^2(t)]=\frac{n_0B}{4} Nout=E[no2(t)]=4n0B
( S N ) o u t = S o u t N o u t = A c 2 P m 4 n 0 B (\frac{S}{N})_{out}=\frac{S_{out}}{N_{out}}=\frac{A^2_cP_m}{4n_0B} (NS)out=NoutSout=4n0BAc2Pm
四、抗噪声性能仿真分析
部分源代码
系统参数设置
%--------------------
%系统参数设置
%--------------------
T_start=0;%开始时间
T_stop=1;%截止时间
T=T_stop-T_start;%仿真持续时间
T_sample=1/1000;%采样间隔
f_sample=1/T_sample;%采样速率
N_sample=T/T_sample;%采样点数
基带信号、载波信号参数设置
%--------------------
%信号参数设置
%--------------------
n=0:N_sample-1;
f_m=10;%基带信号频率
f_c=100;%载波信号频率
Ac=1;%载波信号幅度
m=cos(2*pi*f_m*n*T_sample);%m(t)调制信号
c=cos(2*pi*f_c*n*T_sample);%c(t)载波信号
c_t=cos(2*pi*f_c*n*T_sample);%c_t(t)接收机端载波信号
产生白噪声信号
power_dB=-30;
noise_i=wgn(1,N_sample,power_dB);
生成低通滤波器
生成带通滤波器
h1=2*f_m*sin(2*pi*2*f_m*n*T_sample)/(2*pi*2*f_m*n*T_sample)*exp(1i*2*pi*f_c*n*T_sample)+2*f_m*sin(2*pi*2*f_m*n*T_sample)/(2*pi*2*f_m*n*T_sample)*exp(-1i*2*pi*f_c*n*T_sample);
AM相干解调
%AM信号相干解调时域波形
%-------------------
%设置带通滤波器BPF1,B_t=2B
h1=2*f_m*sin(2*pi*2*f_m*n*T_sample)/(2*pi*2*f_m*n*T_sample)*exp(1i*2*pi*f_c*n*T_sample)+2*f_m*sin(2*pi*2*f_m*n*T_sample)/(2*pi*2*f_m*n*T_sample)*exp(-1i*2*pi*f_c*n*T_sample);
%有噪声AM调幅信号信号时域波形图
figure(1)
subplot(2,2,1)
plot(n*N_sample,s_AMni);
title('有噪声AM信号'),xlabel('t'),ylabel('s_AM(t)+n_i(t)');
%通过带通滤波器BPF1
s_AMn=conv(h1,noise_i);
nn1=0:length(s_AMn)-1;
%接收端载波信号
c_t=cos(2*pi*f_c*nn1*T_sample);
subplot(2,2,2)
plot(nn1*N_sample,s_AMn);
title('有噪声AM信号通过带通滤波器的时域波形'),xlabel('t'),ylabel('s_AM(t)+n(t)');
%解调,乘以接收端载波信号
s_AMnd=s_AMn.*c_t;
subplot(2,2,3)
plot(nn1*T_sample,s_AMnd);
title('有噪声AM信号通过接收端载波信号的时域波形'),xlabel('t'),ylabel('s_AMd(t)+n_d(t)');
%设置低通滤波器
h2=Num1;%fdatool设计h(t)
m_AMno=conv(s_AMnd,h2);
mm_AMno=m_AMno(length(h2)/2:length(h2)/2+N_sample-1);%将通过滤波器之后产生的延迟去掉
nn2=0:length(m_AMno)-1;
subplot(2,2,4)
plot(nn2*T_sample,m_AMno);
title('有噪声AM信号通过相干解调器'),xlabel('t'),ylabel('m_AMo(t)+n_o(t)');
%--------------------
%AM信号相干解调频域波形
%--------------------
SS_AMni=abs(fft(s_AMni))/N_sample;%S_AM+N_i
S_AMni(1:N_sample/2)=SS_AMni(N_sample/2+1:N_sample);
S_AMni(N_sample/2+1:N_sample)=SS_AMni(1:N_sample/2);
SS_AMn=abs(fft(s_AMn))/N_sample;%S_AM+N
S_AMn(1:N_sample/2)=SS_AMn(N_sample/2+1:N_sample);
S_AMn(N_sample/2+1:N_sample)=SS_AMn(1:N_sample/2);
SS_AMnd=abs(fft(s_AMnd))/N_sample;%S_AMd+N_d(f)
S_AMnd(1:N_sample/2)=SS_AMnd(N_sample/2+1:N_sample);
S_AMnd(N_sample/2+1:N_sample)=SS_AMnd(1:N_sample/2);
MM_AMno=abs(fft(mm_AMno))/N_sample;%MM_AMo+N_o
MM_AMno(1:N_sample/2)=MM_AMno(N_sample/2+1:N_sample);
MM_AMno(N_sample/2+1:N_sample)=MM_AMno(1:N_sample/2);
figure(2)
subplot(2,2,1)
plot((n-N_sample/2)/N_sample*f_sample,S_AMni);
title('有噪声AM信号频域波形'),xlabel('f'),ylabel('S_AMni(f)');
subplot(2,2,2)
plot((n-N_sample/2)/N_sample*f_sample,S_AMn);
title('有噪声AM信号通过带通滤波器的频域波形'),xlabel('f'),ylabel('S_AMn(f)');
subplot(2,2,3)
plot((n-N_sample/2)/N_sample*f_sample,S_AMnd);
title('有噪声AM信号通过接收端载波信号频域波形'),xlabel('f'),ylabel('S_AMnd(f)');
subplot(2,2,4)
plot((n-N_sample/2)/N_sample*f_sample,MM_AMno);
title('有噪声AM信号通过相干解调器频域波形'),xlabel('f'),ylabel('MM_AMno(f)');
DSB-SC相干解调
%-------------
%DSB-SC信号信号相干解调时域波形
%-------------
figure(1)
subplot(2,2,1)
plot(n*N_sample,s_DSBni);
title('有噪声DSB-SC信号时域波形'),xlabel('t'),ylabel('s_DSB(t)+n_i(t)');
%通过带通滤波器BPF1
s_DSBn=conv(h1,noise_i);
nn1=0:length(s_DSBn)-1;
%接收端载波信号
c_t=cos(2*pi*f_c*nn1*T_sample);
subplot(2,2,2)
plot(nn1*N_sample,s_DSBn);
title('有噪声DSB信号通过带通滤波器的时域波形'),xlabel('t'),ylabel('s_DSB(t)+n(t)');
%解调,乘以接收端载波信号
s_DSBnd=s_DSBn.*c_t;
subplot(2,2,3)
plot(nn1*T_sample,s_DSBnd);
title('有噪声DSB信号通过接收端载波信号的时域波形'),xlabel('t'),ylabel('s_DSBd(t)+n_d(t)');
%设置低通滤波器
h2=Num;%fdatool设计h(t)
m_DSBno=conv(s_DSBnd,h2);
mm_DSBno=m_DSBno(length(h2)/2:length(h2)/2+N_sample-1);%将通过滤波器之后产生的延迟去掉
nn2=0:length(m_DSBno)-1;
subplot(2,2,4)
plot(nn2*T_sample,m_DSBno);
title('有噪声DSB信号通过相干解调器'),xlabel('t'),ylabel('m_DSBo(t)+n_o(t)');
%--------------------
%DSB-SC信号相干解调频域波形
%--------------------
SS_DSBni=abs(fft(s_DSBni))/N_sample;%S_DSB+N_i
S_DSBni(1:N_sample/2)=SS_DSBni(N_sample/2+1:N_sample);
S_DSBni(N_sample/2+1:N_sample)=SS_DSBni(1:N_sample/2);
SS_DSBn=abs(fft(s_DSBn))/N_sample;%S_DSB+N
S_DSBn(1:N_sample/2)=SS_DSBn(N_sample/2+1:N_sample);
S_DSBn(N_sample/2+1:N_sample)=SS_DSBn(1:N_sample/2);
SS_DSBnd=abs(fft(s_DSBnd))/N_sample;%S_DSBd+N_d(f)
S_DSBnd(1:N_sample/2)=SS_DSBnd(N_sample/2+1:N_sample);
S_DSBnd(N_sample/2+1:N_sample)=SS_DSBnd(1:N_sample/2);
MM_DSBno=abs(fft(mm_DSBno))/N_sample;%MM_DSBo+N_o
MM_DSBno(1:N_sample/2)=MM_DSBno(N_sample/2+1:N_sample);
MM_DSBno(N_sample/2+1:N_sample)=MM_DSBno(1:N_sample/2);
figure(2)
subplot(2,2,1)
plot((n-N_sample/2)/N_sample*f_sample,S_DSBni);
title('有噪声DSB信号频域波形'),xlabel('f'),ylabel('S_DSBni(f)');
subplot(2,2,2)
plot((n-N_sample/2)/N_sample*f_sample,S_DSBn);
title('有噪声DSB信号通过带通滤波器的频域波形'),xlabel('f'),ylabel('S_DSBn(f)');
subplot(2,2,3)
plot((n-N_sample/2)/N_sample*f_sample,S_DSBnd);
title('有噪声DSB信号通过接收端载波信号频域波形'),xlabel('f'),ylabel('S_DSBnd(f)');
subplot(2,2,4)
plot((n-N_sample/2)/N_sample*f_sample,MM_DSBno);
title('有噪声DSB信号通过相干解调器频域波形'),xlabel('f'),ylabel('MM_DSBno(f)');
SSB信号相干解调
%-----------------------
%有噪声SSB调幅信号信号时域波形图
%-----------------------
figure(1)
subplot(2,2,1)
plot(n*N_sample,s_SSBni);
title('有噪声SSB信号时域波形'),xlabel('t'),ylabel('s_SSB(t)+n_i(t)');
%通过带通滤波器BPF1
h1=Num1;%fdatool设置带通滤波器
s_SSBn=conv(h1,noise_i);
nn1=0:length(s_SSBn)-1;
%接收端载波信号
c_t=cos(2*pi*f_c*nn1*T_sample);
subplot(2,2,2)
plot(nn1*N_sample,s_SSBn);
title('有噪声SSB信号通过带通滤波器的时域波形'),xlabel('t'),ylabel('s_SSB(t)+n(t)');
%解调,乘以接收端载波信号
s_SSBnd=s_SSBn.*c_t;
subplot(2,2,3)
plot(nn1*T_sample,s_SSBnd);
title('有噪声SSB信号通过接收端载波信号的时域波形'),xlabel('t'),ylabel('s_SSBd(t)+n_d(t)');
%设置低通滤波器
h2=Num1;%fdatool设计h(t)
m_SSBno=conv(s_SSBnd,h2);
mm_SSBno=m_SSBno(length(h2)/2:length(h2)/2+N_sample-1);%将通过滤波器之后产生的延迟去掉
nn2=0:length(m_SSBno)-1;
subplot(2,2,4)
plot(nn2*T_sample,m_SSBno);
%--------------------
%有噪声SSB信号相干解调各点频域波形
%--------------------
SS_SSBni=abs(fft(s_SSBni))/N_sample;%S_SSB+N_i
S_SSBni(1:N_sample/2)=SS_SSBni(N_sample/2+1:N_sample);
S_SSBni(N_sample/2+1:N_sample)=SS_SSBni(1:N_sample/2);
SS_SSBn=abs(fft(s_SSBn))/N_sample;%S_SSB+N
S_SSBn(1:N_sample/2)=SS_SSBn(N_sample/2+1:N_sample);
S_SSBn(N_sample/2+1:N_sample)=SS_SSBn(1:N_sample/2);
SS_SSBnd=abs(fft(s_SSBnd))/N_sample;%S_SSBd+N_d(f)
S_SSBnd(1:N_sample/2)=SS_SSBnd(N_sample/2+1:N_sample);
S_SSBnd(N_sample/2+1:N_sample)=SS_SSBnd(1:N_sample/2);
MM_SSBno=abs(fft(mm_SSBno))/N_sample;%MM_SSBo+N_o
MM_SSBno(1:N_sample/2)=MM_SSBno(N_sample/2+1:N_sample);
MM_SSBno(N_sample/2+1:N_sample)=MM_SSBno(1:N_sample/2);
figure(2)
subplot(2,2,1)
plot((n-N_sample/2)/N_sample*f_sample,S_SSBni);
title('有噪声SSB信号频域波形'),xlabel('f'),ylabel('S_SSBni(f)');
subplot(2,2,2)
plot((n-N_sample/2)/N_sample*f_sample,S_SSBn);
title('有噪声SSB信号通过带通滤波器的频域波形'),xlabel('f'),ylabel('S_SSBn(f)');
subplot(2,2,3)
plot((n-N_sample/2)/N_sample*f_sample,S_SSBnd);
title('有噪声SSB信号通过接收端载波信号频域波形'),xlabel('f'),ylabel('S_SSBnd(f)');
subplot(2,2,4)
plot((n-N_sample/2)/N_sample*f_sample,MM_SSBno);
title('有噪声SSB信号通过相干解调器频域波形'),xlabel('f'),ylabel('MM_SSBno(f)');
仿真结果图