通信干扰信号识别
抗干扰通信编程作业——通信干扰信号的识别
一、 题目要求
基本部分:1)生成单音干扰、多音干扰、宽带噪声干扰、部分频带噪声干扰、宽带梳
状谱干扰、线性调频干扰等 6 种通信干扰信号;2)选择合适的特征参数,采用决策树法
实现对上述干扰信号的识别,高斯白噪声信道,干噪比(JNR)为 0~15dB,识别正确率大于
95%。扩展部分:选择合适的特征参数,采用 NN 或者 SVM 机器学习实现对上述干扰信
号的识别,高斯白噪声信道,干噪比(JNR)为 0~15dB,识别正确率大于 95%。
二、 通信干扰信号的生成
(一) 单音干扰
单音干扰信号在某个频点上发射,是一个单频连续正弦波形。时域表达式如下
J
( )
t
JP
cos(2
J
f
t
J
)
J
~ (0
,2 )
(1)
其中, JP 为干扰信号功率, Jf 为干扰信号频率, J 为干扰信号相位且服从均匀分布。
设定 1w
J ,得到如 XX 所示时域波形和单边频谱。
JP ,
50Hz
Jf
,
0
single-tone jam
1.5
1
0.5
0
-0.5
-1
-1.5
0
10
20
30
50
60
40
time (milliseconds)
70
80
90
100
Figure 1 单音干扰时域波形
第 1页
通信干扰信号识别
Single-Sided Amplitude Spectrum of single-tone jam
100
10-1
10-2
|
)
f
(
Y
|
10-3
0
20
40
60
100
80
Frequency (Hz)
120
140
160
180
200
Figure 2 单音干扰单边频谱
(二) 多音干扰
多音干扰信号在多个离散的频点上发射,时域表达式如下
( )
J t
Q
l
1
P
l
cos(2
l
f t
l
) ~ (0,2 )
U
l
(2)
其中, lP 为干扰信号功率, lf 为干扰信号频率, l为干扰信号相位且服从均匀分布。
,设定
,
100
+
H
z
l
10Hz,
l
1,2,
..,
.
Q
Q 表示多音干扰的频点个数。设定 ~ (0.9,1.1)
10Q 。得到如 XX 所示时域波形和单边频谱。
lP U
lf
6
4
2
0
-2
-4
-6
0
multiple-tone jam
20
40
60
100
120
80
time (milliseconds)
140
160
180
200
Figure 3 多音干扰时域波形
第 2页
|
)
f
(
Y
|
100
10-1
10-2
10-3
0
通信干扰信号识别
Single-Sided Amplitude Spectrum of multiple-tone jam
50
100
150
200
Frequency (Hz)
250
300
350
Figure 4 多音干扰单边频谱
由图 X 可以看出,多音干扰信号从 100Hz 的频点开始,每间隔 10Hz 有一个冲激,代
表一个干扰频点,一个有 10 个冲激。
(三) 宽带噪声干扰
宽带噪声干扰是将高斯白噪声通过一个宽带滤波器,得到在每一个宽带范围内的阻塞
/ 2J
噪声。设白噪声 ~ (0,
则宽带噪声干扰可以表示为 ( )
n t
N P
(
t
)
W
I
)
,其中 JP 是干扰功率, IW 是低通滤波器的截止频率。
( )
t
,滤波器的频域表达式为
( )
h t
( 2
H j
f
)
1, |
0,
|
f W
I
otherwise
(3)
设定 1w
JP ,
IW
600Hz
,得到时域宽带噪声干扰信号和其单边频谱如图所示。
3
2
1
0
-1
-2
-3
0
Broadband Noise Jam
100
200
300
500
600
400
time (milliseconds)
700
800
900
1000
Figure 5 宽带噪声干扰时域波形
第 3页
)
B
d
(
|
)
f
(
Y
|
-20
-40
-60
-80
-100
-120
0
通信干扰信号识别
Single-Sided Amplitude Spectrum of Broadband Noise Jam
100
200
300
400
500
600
700
800
900
1000
Frequency (Hz)
Figure 6 多音干扰单边频谱
(四) 窄带噪声干扰
与宽带噪声干扰产生方法类似,将白噪声通过一个窄带滤波器,得到瞄准某一个窄带
范围的阻塞信号。窄带滤波器的频域表达式为
( 2
H j
f
)
1, |
0,
f
f
J
|
otherwise
W
I
2
(4)
Jf 表示瞄准频率。设定
Jf
600Hz
,
IW
50Hz
,得到时域窄带噪声干扰信号和其单
边频谱如图所示。
0.8
0.6
0.4
0.2
0
-0.2
-0.4
-0.6
-0.8
Narrowband Noise jam
0
100
200
300
500
400
600
time (milliseconds)
700
800
900
1000
Figure 7 窄带噪声干扰信号时域波形
第 4页
)
B
d
(
|
)
f
(
Y
|
-20
-40
-60
-80
-100
-120
-140
0
通信干扰信号识别
Single-Sided Amplitude Spectrum of Narrowband Noise jam
100
200
300
600
700
800
900
1000
400
500
Frequency (Hz)
Figure 8 窄带噪声干扰信号单边频谱
(五) 宽带梳状谱干扰
将白噪声通过多个窄带滤波器,得到瞄准多个窄带范围的阻塞信号。宽带梳状谱干扰
的滤波器频域表达式为
( 2
H j
f
)
1, |
0,
|
f
f
,
J i
otherwise
W
,
I i
2
(5)
i
其中 1,2,...,
100Hz+(
f
频率 ,
J i
i
Q
,Q 表示窄带滤波器的个数。设定 4Q ,第 i 个窄带滤波器的瞄准
1)
。得到时域宽带梳状谱干扰信号
1,2,...,
50Hz
Hz,
i
100
, ,
I iW
Q
Partialband Noise Jam
和其单边频谱如图所示。
1.5
1
0.5
0
-0.5
-1
-1.5
0
100
200
300
500
400
600
time (milliseconds)
700
800
900
1000
Figure 9 宽带梳状谱干扰信号时域波形
第 5页
通信干扰信号识别
Single-Sided Amplitude Spectrum of Partialband Noise Jam
-20
-30
-40
-50
-60
-70
-80
-90
)
B
d
(
|
)
f
(
Y
|
-100
0
100
200
300
400
500
600
700
800
900
1000
Frequency (Hz)
Figure 10 宽带梳状谱干扰信号时域波形
(六) 扫频干扰
扫频干扰在部分频谱内呈现出随时间线性扫描的特性。其时域表达式为
( )
J t
( )exp(
A t
2
j
t
i
/
2
)
j
j
t
i
( )
A t
Arect t T
( /
)
, |
A
0,
|
t T
else
/ 2
(6)
(7)
其中 i为扫频速率, i为初始角频率,为初始相位,T 为扫频持续时间。由时域表
达式可以看出,扫频干扰的瞬时频率随时间呈不断的线性变化。设定 100 rad/s
20 r
i
T ,得到扫频干扰时域信号和其单边频谱如图所示。
, 0 , 2s
ad
i
,
Sweeping Jam
1
0.5
0
-0.5
-1
0
100
200
300
500
400
600
time (milliseconds)
700
800
900
1000
Figure 11 扫频干扰信号时域波形
第 6页
通信干扰信号识别
Single-Sided Amplitude Spectrum of Sweeping Jam
50
100
150
300
350
400
450
500
200
250
Frequency (Hz)
Figure 12 扫频干扰信号时域波形
)
B
d
(
|
)
f
(
Y
|
-20
-40
-60
-80
-100
-120
0
三、 特征参数提取
(一) 信号预处理
在对信号进行特征参数提取前,需要首先对信号进行预处理,包括归一化、中心化。
1. 归一化
设样本数据为 ( ),
x n n
1,2,...,
N
,采用中间值法将样本数据归一化到[-1,1]内,公式表
示为
( )
y k
( )
x k
1 (
x
max
2
x
mid
x
min
)
,
k
1,2,...,
N
x
mid
x
max
x
min
2
其中 maxx 表示最大值, minx 表示最小值。
2. 中心化
采用 Z 分法
( )
y k
( )
x k
,
k
1,2,...,
N
其中为随机变量 X 的均值,为其标准差。
(二) 时域参数
1. R 参数
第 7页
(8)
(9)
(10)
通信干扰信号识别
R
2
2
(11)
其中为随机变量 X 的均值,为其标准差,反映了变量 X 包络的变化程度。R 越大,
说明信号包络的变化程度越大;当 R 较小时,说明信号包络的变化程度较小。对 6 种类型
的通信干扰信号计算其 R 参数,得到如 Figure 13 所示结果。
Figure 13 不同干扰信号 R 参数随 JNR 变化趋势
可以看出,单音干扰和扫频干扰由于其信号包络恒定,因此其 R 参数相比其他 4 种干
扰信号要小,并且二者的 R 参数值几乎相同,随着 JNR 增大,R 参数减小。宽带噪声干扰
和窄带噪声干扰其信号包络与 JNR 大小无关,因此其 R 参数几乎不随 JNR 变化。梳状谱
干扰随着 JNR 增大,R 参数增大,说明其 JNR 越大时梳状谱信号的包络变化程度越大。多
音干扰的 R 参数没有明显的特征。
2. 时域矩偏度系数 a3
a
3
3
)
(
E x
3
(11)
为时域信号包络 x 的均值,为其标准差。 3b 描述了频谱偏离正态分布的程度。 3b 越
大,说明频谱偏离正态分布的程度越大。
(三) 频域参数
第 8页