模拟幅度调制系统抗干扰性能仿真分析
文章目录
1.引言
使载波的振幅按照所需传送信号的变化规律而变化,但频率保持不变的调制方法称之为幅度调制。调幅在有线电或无线电通信和广播中应用甚广。模拟基带信号 m ( t ) m(t) m(t)对载波 c ( t ) c(t) c(t)进行调制是按 m ( t ) m(t) m(t)的变化规律去控制载波的幅度。振幅调制主要包括模拟常规调幅(AM)、抑制载波双边带调幅(DSB-SC)、单边带调幅(SSB)等。本文利用MATLAB仿真模拟传输系统,分析各种系统对信号时域以及频域特性的变化关系,得到其输入输出信噪比与相关的增益,进行对比分析其抗噪声性能。
2.系统模型
2.1常规调幅(AM)
图
1
图1
图1
AM调制器模型如图1所示。其时域表达式为:
s
A
M
(
t
)
=
A
c
[
1
+
m
(
t
)
]
c
o
s
2
π
f
c
t
(1)
s_{AM}(t)=A_c[1+m(t)]cos2πf_ct\tag1
sAM(t)=Ac[1+m(t)]cos2πfct(1)
进一步,我们定义调幅指数:
β
A
M
=
m
a
x
∣
m
(
t
)
∣
(2)
β_{AM} =max∣m(t)∣\tag2
βAM=max∣m(t)∣(2)
根据(1),进行傅里叶变换,可以得到其频谱表达式:
M
(
f
)
=
A
c
2
[
δ
(
f
−
f
c
)
+
δ
(
f
+
f
c
)
]
+
A
c
2
[
M
(
f
−
f
c
)
+
M
(
f
+
f
c
)
]
(3)
M(f)=\frac{A_c}{2}[δ(f-f_c)+δ(f+f_c)]+\frac{A_c}{2}[M(f-f_c)+M(f+f_c)]\tag3
M(f)=2Ac[δ(f−fc)+δ(f+fc)]+2Ac[M(f−fc)+M(f+fc)](3)
图
2
图2
图2
AM信号的解调除了可以采用想干解调外,还可以采用非相干的包络检波。AM非相干解调如图2所示,理想包络检波器的数输出是输入带通信号
r
(
t
)
r(t)
r(t)的包络
∣
r
L
(
t
)
∣
|r_L(t)|
∣rL(t)∣。隔直流电路去除直流分量
A
c
A_c
Ac,输出为
y
0
(
t
)
=
A
m
(
t
)
y_0(t)=Am(t)
y0(t)=Am(t)。
2.2抑制载波双边带调幅(DSB-SC)
图
3
图3
图3
双边带抑制载波调幅信号
s
(
t
)
s(t)
s(t)是利用均值为零的模拟基带信号
m
(
t
)
m(t)
m(t)与正弦载波相乘得到的,如图3所示。
此调幅信号的时域表达式为:
s
D
S
B
(
t
)
=
m
(
t
)
A
c
c
o
s
2
π
f
c
t
(4)
s_{DSB}(t)=m(t)A_ccos2πf_ct\tag4
sDSB(t)=m(t)Accos2πfct(4)
进行傅里叶变换,得到其频域表达式为:
s
D
S
B
(
f
)
=
A
c
2
[
M
(
f
−
f
c
)
+
M
(
f
+
f
c
)
]
(5)
s_{DSB}(f)=\frac{A_c}{2}[M(f-f_c)+M(f+f_c)]\tag5
sDSB(f)=2Ac[M(f−fc)+M(f+fc)](5)
图
4
图4
图4
由图4可以看出,若基带信号的带宽是B,则已调信号的带宽为2B,其中
∣
f
∣
>
f
c
|f|>f_c
∣f∣>fc的部分叫上边带,
∣
f
∣
<
f
c
|f|<f_c
∣f∣<fc的部分叫下边带。
图
5
图5
图5
DSB调制可以采用图5的相干解调器来解调,该解调器适用于多种模拟及数字调制。
s
d
(
t
)
=
m
(
t
)
c
o
s
2
π
f
c
t
c
o
s
[
2
π
(
f
c
+
Δ
f
)
t
+
θ
0
]
=
1
2
m
(
t
)
c
o
s
(
2
π
Δ
f
t
+
Δ
θ
)
+
c
o
s
[
2
π
(
2
f
c
+
Δ
f
)
t
+
Δ
θ
]
(6)
s_d(t)=m(t)cos2πf_ctcos[2π(f_c+Δf)t+θ_0]=\frac{1}{2}m(t){cos(2πΔft+Δθ)+cos[2π(2f_c+Δf)t+Δθ]}\tag6
sd(t)=m(t)cos2πfctcos[2π(fc+Δf)t+θ0]=21m(t)cos(2πΔft+Δθ)+cos[2π(2fc+Δf)t+Δθ](6)
s
0
(
t
)
=
1
2
m
(
t
)
c
o
s
(
2
π
Δ
f
t
+
Δ
θ
)
(7)
s_0(t)=\frac{1}{2}m(t)cos(2πΔft+Δθ)\tag7
s0(t)=21m(t)cos(2πΔft+Δθ)(7)
2. 3单边带调幅(SSB)
双边带调制(包括DSB-SC以及AM)的两个频谱都有两个对称的边带,使其戴宽师基带信号带宽的2倍。为了节约信号传输中对信道带宽的占用,可以去除其中一个边带来传输,这种方式就是单边带(SSB)调制,如图6所示:
图
6
图6
图6
我们可以采用滤波的方法,从双边带信号中将下边带以及上边带信号分离出来。图6中(a)为SSB调制器,其中
h
S
S
B
(
t
)
h_{SSB}(t)
hSSB(t)可以为下边带滤波器(b)或者上边带滤波器(c),因此单边带信号的波形表达式为:
s
S
S
B
(
t
)
=
A
c
m
(
t
)
c
o
s
2
π
f
c
t
∗
h
S
S
B
(
t
)
(8)
s_{SSB}(t)=A_ cm(t)cos2πf_ct∗h_{SSB}(t)\tag8
sSSB(t)=Acm(t)cos2πfct∗hSSB(t)(8)
傅立叶变换得到其频域表达式为:
S
S
S
B
(
f
)
=
A
c
2
[
M
(
f
−
f
c
)
+
M
(
f
+
f
c
)
]
H
(
S
S
B
)
(
f
)
(9)
S_{SSB}(f)=\frac{A_c}{2}[M(f-f_c)+M(f+f_c)]H_(SSB)(f)\tag9
SSSB(f)=2Ac[M(f−fc)+M(f+fc)]H(SSB)(f)(9)
图
7
图7
图7
图7中,(a)为下边带
S
L
S
S
B
(
f
)
S_{LSSB}(f)
SLSSB(f),(b)为上边带
S
H
S
S
B
(
f
)
S_{HSSB}(f)
SHSSB(f)。
单边带信号的解调可以采用相干解调器,恢复载波理想同步时,解调器输出正比于输入信号复包络的实部:
y
0
(
t
)
=
A
c
2
m
(
t
)
(10)
y_0(t)=\frac{A_c}{2}m(t)\tag{10}
y0(t)=2Acm(t)(10)。
3.抗干扰性能理论分析
3.1AM相干解调抗噪声性能
图
8
图8
图8
Input SNR:
S
i
n
=
<
s
2
(
t
)
>
=
A
c
2
2
[
1
+
P
m
]
S_{in}=<s^2(t)>=\frac{A^2_c}{2}[1+P_m]
Sin=<s2(t)>=2Ac2[1+Pm]
N
i
n
=
E
[
n
2
(
t
)
]
=
n
0
B
A
M
=
2
n
0
B
N_{in}=E[n^2(t)]=n_0B_{AM}=2n_0B
Nin=E[n2(t)]=n0BAM=2n0B
(
S
/
N
)
i
n
=
S
i
n
N
i
n
=
A
c
2
2
[
1
+
P
m
]
2
n
0
B
(11)
(S/N)_{in}=\frac{S_{in}}{N_{in}}=\frac{\frac{A^2_c}{2}[1+P_m]}{2n_0B}\tag{11}
(S/N)in=NinSin=2n0B2Ac2[1+Pm](11)
Output SNR:
S
o
u
t
=
<
s
0
2
(
t
)
>
=
A
c
2
P
m
4
S_{out}=<s^2_0(t)>=\frac{A^2_cP_m}{4}
Sout=<s02(t)>=4Ac2Pm
N
o
u
t
=
E
[
n
0
2
(
t
)
]
=
1
4
E
[
n
c
2
(
t
)
]
=
1
4
N
i
n
N_{out}=E[n^2_0(t)]=\frac{1}{4}E[n^2_c(t)]=\frac{1}{4}N_{in}
Nout=E[n02(t)]=41E[nc2(t)]=41Nin
(
S
/
N
)
o
u
t
=
S
o
u
t
N
o
u
t
=
A
c
2
P
m
2
n
0
B
(12)
(S/N)_{out}=\frac{S_{out}}{N_{out}}=\frac{A^2_cP_m}{2n_0B}\tag{12}
(S/N)out=NoutSout=2n0BAc2Pm(12)
G
A
M
=
2
P
M
1
+
P
M
(13)
G_{AM}=\frac{2P_M}{1+P_M}\tag{13}
GAM=1+PM2PM(13)
3.2DSB-SC相干解调抗噪声性能
图
9
图9
图9
Input SNR:
S
i
n
=
s
2
(
t
)
‾
=
A
c
2
2
m
(
t
)
2
‾
=
A
c
2
P
m
2
S_{in}=\overline{s^2(t)}=\frac{A^2_c}{2}\overline{m(t)^2}=\frac{A^2_cP_m}{2}
Sin=s2(t)=2Ac2m(t)2=2Ac2Pm
N
i
n
=
E
[
n
2
(
t
)
]
=
n
0
B
D
S
B
−
S
C
=
2
n
0
B
N_{in}=E[n^2(t)]=n_0B_{DSB-SC}=2n_0B
Nin=E[n2(t)]=n0BDSB−SC=2n0B
(
S
/
N
)
i
n
=
S
i
n
N
i
n
=
1
2
A
c
2
P
m
2
n
0
B
(14)
(S/N)_{in}=\frac{S_{in}}{N_{in}}=\frac{1}{2}\frac{A^2_cP_m}{2n_0B}\tag{14}
(S/N)in=NinSin=212n0BAc2Pm(14)
Output SNR:
S
o
u
t
=
E
[
s
0
2
(
t
)
]
=
A
c
2
P
M
4
=
1
2
S
i
n
S_{out}=E[s^2_0(t)]=\frac{A^2_cP_M}{4}=\frac{1}{2}S_{in}
Sout=E[s02(t)]=4Ac2PM=21Sin
N
o
u
t
=
E
[
n
0
2
(
t
)
]
=
1
4
E
[
n
c
2
(
t
)
]
=
1
4
N
i
n
N_{out}=E[n^2_0(t)]=\frac{1}{4}E[n^2_c(t)]=\frac{1}{4}N_{in}
Nout=E[n02(t)]=41E[nc2(t)]=41Nin
(
S
/
N
)
o
u
t
=
S
o
u
t
N
o
u
t
=
2
S
i
n
N
i
n
=
A
c
2
P
m
2
n
0
B
(15)
(S/N)_{out}=\frac{S_{out}}{N_{out}}=2\frac{S_{in}}{N_{in}}=\frac{A^2_cP_m}{2n_0B}\tag{15}
(S/N)out=NoutSout=2NinSin=2n0BAc2Pm(15)
G
D
S
B
−
S
C
=
2
(16)
G_{DSB-SC}=2\tag{16}
GDSB−SC=2(16)
3.3SSB相干解调抗噪声性能
图
10
图10
图10
Input SNR:
S
i
n
=
<
s
2
(
t
)
>
=
A
c
2
P
m
4
S_{in}=<s^2(t)>=\frac{A^2_cP_m}{4}
Sin=<s2(t)>=4Ac2Pm
N
i
n
=
E
[
n
2
(
t
)
]
=
n
0
B
S
S
B
=
n
0
B
N_{in}=E[n^2(t)]=n_0B_{SSB}=n_0B
Nin=E[n2(t)]=n0BSSB=n0B
(
S
/
N
)
i
n
=
S
i
n
N
i
n
=
1
4
A
c
2
P
m
n
0
B
(17)
(S/N)_{in}=\frac{S_{in}}{N_{in}}=\frac{1}{4}\frac{A^2_cP_m}{n_0B}\tag{17}
(S/N)in=NinSin=41n0BAc2Pm(17)
Output SNR:
S
o
u
t
=
<
s
0
2
(
t
)
>
=
<
m
2
(
t
)
>
16
=
1
4
S
i
n
S_{out}=<s^2_0(t)>=\frac{<m^2(t)>}{16}=\frac{1}{4}S_{in}
Sout=<s02(t)>=16<m2(t)>=41Sin
N
o
u
t
=
E
[
n
0
2
(
t
)
]
=
1
4
E
[
n
c
2
(
t
)
]
=
1
4
N
i
n
N_{out}=E[n^2_0(t)]=\frac{1}{4}E[n^2_c(t)]=\frac{1}{4}N_{in}
Nout=E[n02(t)]=41E[nc2(t)]=41Nin
(
S
/
N
)
o
u
t
=
S
o
u
t
N
o
u
t
=
A
c
2
P
m
4
n
0
B
(18)
(S/N)_{out}=\frac{S_{out}}{N_{out}}=\frac{A^2_cP_m}{4n_0B}\tag{18}
(S/N)out=NoutSout=4n0BAc2Pm(18)
G
S
S
B
=
1
(19)
G_{SSB}=1\tag{19}
GSSB=1(19)
4.仿真实现与仿真结果
4.1AM常规调幅仿真
系统参数设置
%------------------
%系统参数设置
%-----------------
T_start=0;%开始时间
T_stop=1;%截止时间
T=T_stop-T_start;%仿真持续时间
T_sample=0.001;%采样间隔
f_sample=1/T_sample; % 采样速率
N_sample=T/T_sample;% 采样点数
%-----------------
%参数设置
%-----------------
fm=5;
fc=200;
Ac=1;
N0=30
f_res=f_sample/N_sample;%频率分辨率
f_max=f_res*N_sample/2;%最大频率
m=cos(2*pi*fm*n*T_sample);%基带信号
c=Ac*cos(2*pi*fc*n*T_sample);%载波信号
AM调幅
n=0:N_sample;
nt = wgn(1,length(n),N0/2);
m=1*cos(2*pi*fm*n*T_sample);
am=(1+m).*cos(2*pi*fc*n*T_sample);
subplot(3,1,1);
plot(n*T_sample,m);ylabel('m(t)');xlabel('t');
subplot(3,1,2);
plot(n*T_sample,am);ylabel('Sam(t)');xlabel('t');
subplot(3,1,3);
plot(n*T_sample,nt+am);ylabel('S(t)');xlabel('t');title('加入噪声后信号');
figure;
F=abs(fft(m));
F_rearrange=[F(N_sample/2+1:N_sample-1),F(1:N_sample/2)];
subplot(3,1,1);
plot((-N_sample/2+1:N_sample/2-1)*f_res,F_rearrange(1:N_sample-1));
xlabel('f');
ylabel('M(f)');
F=abs(fft(am));
F_rearrange=[F(N_sample/2+1:N_sample-1),F(1:N_sample/2)];
subplot(3,1,2);
plot((-N_sample/2+1:N_sample/2-1)*f_res,F_rearrange(1:N_sample-1));
xlabel('f');
ylabel('Sam(f)');
F=abs(fft(nt+am));
F_rearrange=[F(N_sample/2+1:N_sample-1),F(1:N_sample/2)];
subplot(3,1,3);
plot((-N_sample/2+1:N_sample/2-1)*f_res,F_rearrange(1:N_sample-1));
xlabel('f');ylabel('S(f)');
title('加入噪声后信号的傅里叶变换');
fct = cos(2*pi*fc*n*T_sample);
amt = (am).*fct;
figure
subplot(2,1,1);
plot(n*T_sample,m);ylabel('m(t)');xlabel('t');title('消息信号');
subplot(2,1,2);
n = (0:length(amt)-1)*T_sample;
[amt,~]=bw(n,amt);
plot(n,amt);
ylabel('AM(t)');xlabel('t');title('解调信号');
仿真截图:
4.2抑制载波双边带调幅(DSB-SC)仿真
DSB-SC
n=0:N_sample;
nt = wgn(1,length(n),N0/2);
m=1*cos(2*pi*fm*n*T_sample);
dsb=m.*cos(2*pi*fc*n*T_sample);
subplot(3,1,1);
plot(n*T_sample,m);ylabel('m(t)');xlabel('t');
subplot(3,1,2);
plot(n*T_sample,dsb);ylabel('Sdsb(t)');xlabel('t');
subplot(3,1,3);
plot(n*T_sample,nt+dsb);ylabel('S(t)');xlabel('t');title('加入噪声后信号');
figure;
F=abs(fft(m));
F_rearrange=[F(N_sample/2+1:N_sample-1),F(1:N_sample/2)];
subplot(3,1,1);
plot((-N_sample/2+1:N_sample/2-1)*f_res,F_rearrange(1:N_sample-1));
xlabel('f');
ylabel('M(f)');
F=abs(fft(dsb));
F_rearrange=[F(N_sample/2+1:N_sample-1),F(1:N_sample/2)];
subplot(3,1,2);
plot((-N_sample/2+1:N_sample/2-1)*f_res,F_rearrange(1:N_sample-1));
xlabel('f');
ylabel('Sdsb(f)');
F=abs(fft(nt+dsb));
F_rearrange=[F(N_sample/2+1:N_sample-1),F(1:N_sample/2)];
subplot(3,1,3);
plot((-N_sample/2+1:N_sample/2-1)*f_res,F_rearrange(1:N_sample-1));
xlabel('f');ylabel('S(f)');
title('加入噪声后信号的傅里叶变换');
figure
subplot(3,1,1);
plot(n*T_sample,m);ylabel('m(t)');xlabel('t');title('消息信号');
subplot(3,1,2);
mm = 2*dsb.*cos(2*pi*fc*n*T_sample);
plot(n,mm);
ylabel('DSB(t)');xlabel('t');title('未通过滤波器');
subplot(3,1,3);
plot(n*T_sample,m);
ylabel('DSB(t)');xlabel('t');title('解调信号');
仿真截图:
4.3 单边带调幅(SSB)仿真
SSB调幅
n=0:N_sample;
m=0.5*cos(2*pi*fm*n*T_sample);
ssc=m.*cos(2*pi*fc*n*T_sample);
ssc = conv(ssc,num);
subplot(3,1,1);
plot((0:length(m)-1)*T_sample,m);ylabel('m(t)');xlabel('t');
subplot(3,1,2);
plot((0:length(ssc)-1)*T_sample,ssc);ylabel('Sssc(t)');xlabel('t');
subplot(3,1,3);
nt = wgn(1,length(ssc),N0/2);
plot((0:length(ssc)-1)*T_sample,nt+ssc);ylabel('S(t)');xlabel('t');title('加入噪声后信号');
figure;
f_res=f_sample/N_sample;%频率分辨率
f_max=f_res*N_sample/2;%最大频率
F=abs(fft(m));
F_rearrange=[F(N_sample/2+1:N_sample-1),F(1:N_sample/2)];
subplot(3,1,1);
plot((-N_sample/2+1:N_sample/2-1)*f_res,F_rearrange(1:N_sample-1));
xlabel('f');
ylabel('M(f)');
F=abs(fft(ssc));
F_rearrange=[F(N_sample/2+1:N_sample-1),F(1:N_sample/2)];
subplot(3,1,2);
plot((-N_sample/2+1:N_sample/2-1)*f_res,F_rearrange(1:N_sample-1));
xlabel('f');
ylabel('Sssc(f)');
F=abs(fft(nt+ssc));
F_rearrange=[F(N_sample/2+1:N_sample-1),F(1:N_sample/2)];
subplot(3,1,3);
plot((-N_sample/2+1:N_sample/2-1)*f_res,F_rearrange(1:N_sample-1));
xlabel('f');ylabel('S(f)');
title('加入噪声后信号的傅里叶变换');
fct = cos(2*pi*fc*(0:length(ssc)-1)*T_sample);
amt = 4*(ssc).*fct;
figure;
subplot(3,1,1);
n = (0:length(m)-1)*T_sample;
plot(n,m);ylabel('m(t)');xlabel('t');title('消息信号');
subplot(3,1,2);
n = (0:length(amt)-1)*T_sample;
plot(n,amt);
ylabel('ssc(t)');xlabel('t');title('未滤波');
subplot(3,1,3);
amt = conv(amt,num);
n = (0:length(amt)-1)*T_sample;
plot(n(100:1100)-0.1,amt(100:1100));
ylabel('ssc(t)');xlabel('t');title('解调信号');
amt = 4*(ssc).*fct;
figure;
仿真截图:
5.小结
通过本次的仿真实验可以得出在AM、DSB-SC、SSB这三个系统中,DSB-SC的解调增益为2,SSB系统的解调增益为1,而AM相干解调系统的解调增益不为定值,且其大小随着基带信号功率的变化而变化但小于2。从表格可看出SSB信号中有两个DSB信号,但只有前一项对解调输出有贡献,这一部分的功率占SSB信号总功率
P
R
P_R
PR的一半,后一项虽对解调输出无贡献,但他能消掉DSB的一个边带,使信号带宽变窄,进而使噪声功率减半。即DSB-SC系统与SSB系统的抗干扰能力相同,AM系统的抗干扰能力与基带信号的功率相关,且功率越大抗干扰能力越强。
6.参考文献
[1]通信原理(第四版).周炯槃.北京邮电大学出版社``
[2]现代通信原理6.1 常规调幅调制(AM)与抑制载波双边带(DSB-SC)调制
[3]现代通信原理6.2:单边带(SSB)调制