Page 138 - 《应用声学》2023年第4期
P. 138
800 2023 年 7 月
计结果、经验概率分布统计结果和分布检验结果 表 1 参数估计结果统计
放在一起显示,如图 10、图 11 所示。从图中可以看 Table 1 Statistics of parameter estimation results
出,统计分析结果,包括均值、方差、偏度、峰度、
近似背景噪声 冲击噪声占优 船舶航行噪声占优
α 稳定分布的 4 个参数估计,可以用来刻画环境噪
α = 1.52 α = 1.49 α = 1.87
声模型,表征不同时间或地点的变异性和比较声学 2 ∼ 30 Hz 2.48/11.89 2.09/9.02 1.34/6.59
环境。从表 1 看出,这段数据中的冲击噪声主要来 上 0.46/3.07 0.59/2.99 −0.37/2.43
表 α = 1.63 α = 1.60 α = 1.81
源是在 2∼30 Hz 频段内很明显的冰源噪声,上表
面 30 ∼ 100 Hz 1.92/10.41 1.87/27.29 1.09/13.35
面声道中 α、β、γ、δ 参数估计值分别为 1.60、0.36、 声 0.32/3.12 0.40/3.06 −0.08/2.94
0.08、0.05,下表面声道中 α、β、γ、δ 参数估计值分 道 α = 1.93 α = 1.83 α = 1.86
100 ∼ 500 Hz 0.77/27.83 0.85/27.29 4.02/101.16
别为 1.98、0.82、0.08、0.17。图中给出的峰度结果
0.01/3.47 0.01/3.58 0.01/3.33
对于分析噪声是否存在强冲击噪声及线谱信号具
α = 2.0 α = 2.0 α = 2.0
有重要意义,各处理结果的峰度图中黄色标记线 2 ∼ 30 Hz 1.38/4.63 1.15/5.63 0.58/3.42
为 2,红色标记线为 3,峰度大于 3,说明噪声的分 下 −0.04/2.36 −0.18/1.95 −0.02/1.86
表 α = 1.98 α = 1.98 α = 1.56
布与高斯噪声相比较为陡峭,具有尖顶峰,而峰度
面 30 ∼ 100 Hz 0.91/14.35 1.04/14.74 0.84/27.15
小于 3,则说明噪声的分布与高斯噪声相比较为平 声 −0.01/2.98 −0.03/2.92 −0.02/3.12
坦,具有平峰值。偏度与峰度类似,它也是描述噪 道 α = 1.95 α = 1.85 α = 1.92
100 ∼ 500 Hz 2.48/11.89 1.29/19.45 2.91/71.51
声分布形态的统计量,其描述的是某总体取值分布
0.01/3.28 0.01/3.50 0.01/3.20
的对称性,与 α 稳定分布中的 β 值存在一定的对应
注:表格中,上为 α 值,中为偏度最大值/峰度最大值,下为偏度均
关系。 值/峰度均值。
4000 ᮠဋᄰவڏ
1.0
ࣨए/(10 -3 V) -0.5 0 3500 ࠛएͥᝠڏ
0.5
গѬ࣋ࠛएڏ
3000
-1.0
50 100 150 200 250 300 350 400 450 500 550 2500
f↼x↽ 2000
100 -40 1500
ᮠဋ/Hz 80 -60 1000
60
-80
40 500
-100
50 100 150 200 250 300 350 400 450 500 550 600 0
T/s -1.5 -1.0 -0.5 0 0.5 1.0 1.5
x⊳ ֓
(a) 15:28, 30~100 HzηՂᮠڏ (b) 15:28, 30~100 HzηՂѬ࣋
1.0 6
4 ᰴளѬ࣋
0.5 2 ࠄ
کϙ 0 Ϡए 0 0.999
0.997
-2 0.990
0.980
-0.5 0.950
-4 0.900
-1.0 -6 0.750
0 200 400 600 0 200 400 600 ഐဋ 0.500
0.2 8 0.250
0.100
0.050
0.020
வࣀ 0.1 ए 6 0.010
0.003
0.001
4 α=1.60, β=0.36,
γ=0.08, δ=-0.05
0 2
0 200 400 600 0 200 400 600 -1.0 -0.5 0 0.5 1.0
T/s T/s ࣨϙ⊳ ֓
(c) 15:28, 30~100 HzηՂѬ࣋ፒᝠ (d) 15:28, 30~100 HzѬ࣋ೝᰎፇ౧
图 10 冲击噪声占优时上表面声道噪声分布特性 (时间:0150)
Fig. 10 Noise distribution characteristics of upper duct when impact noise dominates(Time: 0150)