https://ieeexplore.ieee.org/document/7812508/references#references
弹道模型
一个超音速子弹在沿着轨迹飞行时,它受一个重力和一个拖曳力。对于水平发射而言,重力远小于拉力,可被忽略。根据这一假定,超音速子弹的路径是近似线性的,让
v
(
x
)
v(x)
v(x)表示子弹在与开火点距离为x的地方的超音速度,它所受的拉力可以表示为:
F
D
=
−
m
d
v
d
t
=
1
2
ρ
C
D
A
v
2
(
1
)
F_D=-m \frac{dv}{dt}=\frac 12\rho C_DAv^2 (1)
FD=−mdtdv=21ρCDAv2(1)
其中
m
m
m和
A
A
A分别是子弹的质量和横截面积,
ρ
\rho
ρ是空气密度,
C
D
C_D
CD是拉力系数。在很多实际应用中,如果开火的仰角
<
1
0
。
<10^。
<10。,(1)可以被用来得到一个准确的子弹速度表达式。
拉力系数被给出为:
C
D
=
κ
(
c
v
)
η
(
2
)
C_D=\kappa(\frac cv)^\eta (2)
CD=κ(vc)η(2)
其中
c
c
c是空气中的声速,
κ
\kappa
κ是一个由子弹的真实形状决定的因子,并且因子
η
≥
0
\eta\ge0
η≥0,将式(2)代入(1)并且化简为
d
v
d
t
=
−
2
C
b
−
1
v
2
−
η
(
3
)
\frac {dv}{dt}=-2C_b^{-1}v^{2-\eta} (3)
dtdv=−2Cb−1v2−η(3)
其中
C
b
C_b
Cb作为弹道常量被给出为:
C
b
=
4
m
ρ
κ
A
c
η
(
4
)
C_b=\frac{4m}{\rho\kappa Ac^\eta}(4)
Cb=ρκAcη4m(4)
将式子(3)等号左边写作
v
d
v
/
d
x
vdv/dx
vdv/dx,然后式(3)等号两端以
x
x
x进行积分,设定初始状态为:当
x
=
0
x=0
x=0时
v
=
V
0
v=V_0
v=V0,可以得到一个子弹速度的表达式,也是
x
x
x的函数(即速度曲线):
v
(
x
)
=
(
V
0
η
−
2
η
C
b
−
1
x
)
1
/
η
,
0
≤
x
≤
1
2
η
C
b
(
V
0
η
−
c
η
)
(
5
)
v(x)=(V_0^\eta-2\eta C_b^{-1}x)^{1/\eta},0\le x\le \frac 1{2\eta}C_b(V_0^\eta-c^\eta) (5)
v(x)=(V0η−2ηCb−1x)1/η,0≤x≤2η1Cb(V0η−cη)(5)
其中
V
0
V_0
V0是子弹的枪口速度(或初始速度),在式(5)中,
x
x
x的上限决定了超过子弹速度进入亚声速的点的最大距离。在此界限中,当
C
b
→
∞
C_b\to\infin
Cb→∞时,对所有的
x
x
x都有
v
(
x
)
→
V
0
v(x)\to V_0
v(x)→V0,这等价于子弹以定长速度
V
0
V_0
V0飞行。
本文中所考量的5种不同的弹道模型对应着5个不同的 η \eta η值:0,0.5,0.514,1和2。这5种模型分别被称作指数速度,二次速度,Duckworth,线性速度,和常拉力模型。
二次速度模型
当
η
=
0.5
\eta=0.5
η=0.5时,式(5)变成
v
(
x
)
=
(
V
0
1
/
2
−
C
b
−
1
x
)
2
,
0
≤
x
≤
C
b
(
V
0
(
1
/
2
)
−
c
1
/
2
)
.
(
6
)
v(x)=(V_0^{1/2}-C_b^{-1}x)^2,0\le x\le C_b(V_0^(1/2)-c^{1/2}).(6)
v(x)=(V01/2−Cb−1x)2,0≤x≤Cb(V0(1/2)−c1/2).(6)
子弹速度的衰减呈飞行距离的二次函数。Kozick参照此模型作为拉力模型并用它为他们的研究生成仿真数据,该模型的结果与开火点的准确定位并不匹配。
Duckworth 模型
通过假定声音辐射(冲击波SW)决定了超音速子弹的拉力,Duckworth发现拉力系数 C D C_D CD 表达式中的指数 η \eta η的数值为0.514。将该值代入式(5)即可得到Duckworth模型。逐渐降低的子弹速度是飞行距离的近似二次函数。
指数速度模型
当取边界值
η
→
0
\eta\to0
η→0时,此模型可从式(5)中推导出来
v
(
x
)
=
V
0
e
−
2
C
b
−
1
x
,
0
≤
x
≤
1
2
C
b
(
ln
V
0
−
ln
c
)
.
(
7
)
v(x)=V_0e^{-2C_b^{-1}x},0\le x\le\frac 12C_b(\ln V_0-\ln c).(7)
v(x)=V0e−2Cb−1x,0≤x≤21Cb(lnV0−lnc).(7)
线性速度模型
当
η
=
1
\eta =1
η=1,(5)变为
v
(
x
)
=
V
0
−
2
C
b
−
1
x
,
0
≤
x
≤
1
2
C
b
(
V
0
−
c
)
.
(
8
)
v(x)=V_0-2C_b^{-1}x,0\le x\le\frac12C_b(V_0-c).(8)
v(x)=V0−2Cb−1x,0≤x≤21Cb(V0−c).(8)
子弹的速度随着飞行距离的增加线性降低(即
d
v
/
d
x
=
c
o
n
s
t
a
n
t
dv/dx=constant
dv/dx=constant)。
常拉力模型
当
η
=
2
\eta=2
η=2时,可以发现拉力
F
D
F_D
FD是独立于子弹速度的(即
F
D
=
c
o
n
s
t
a
n
t
F_D=constant
FD=constant)。将
η
=
2
\eta=2
η=2代入式(5)可得到常拉力模型:
v
(
x
)
=
(
V
0
2
−
4
C
b
−
1
x
)
1
/
2
,
0
≤
x
≤
1
4
C
b
(
V
0
2
−
c
2
)
.
(
9
)
v(x)=(V_0^2-4C_b^{-1}x)^{1/2},0\le x\le \frac14C_b(V_0^2-c^2).(9)
v(x)=(V02−4Cb−1x)1/2,0≤x≤41Cb(V02−c2).(9)
模型评估
使用36种用于不同活动的不同种类的子弹的真实数据来对五种不同的模型的预测子弹速度曲线的准确性进行评估。这36种子弹的速度曲线可以从生产商提供的数据表中获得。子弹从水平方向上射出并且它的速度曲线通过使用多普勒雷达测得,多普勒雷达提供了最准确的子弹飞行特征测量方式。通过多普勒雷达技术,每一次射击的枪口速度通常可被确定在0.1m/s之内。多次射击被用来测量不同类型子弹的速度曲线。总共有36个数据集,每个代表了对某一特定的子弹雷达所测量到的速度曲线。速度数据的有效距离是0,100,200,300,600,800,且对某些子弹类型来说,可达到900或1000m(从开火点)。速度数据降落到330m/s以下即被丢弃。此36种子弹的口径范围从5.56mm到12.95mm,且步枪枪管的长度从510mm到960mm。
一个弹道模型(5)(针对
η
\eta
η的某一特定值)下的针对一个给定子弹类型的雷达测量速度数据的最小二乘拟合被给出如下:
v
^
L
S
(
x
k
)
=
(
V
^
0
η
−
2
η
C
^
b
−
1
x
k
)
1
/
η
,
1
≤
k
≤
K
.
(
10
)
\hat v_{LS}(x_k)=(\hat V_0^{\eta}-2\eta\hat C_b^{-1}x_k)^{1/\eta},1\le k\le K.(10)
v^LS(xk)=(V^0η−2ηC^b−1xk)1/η,1≤k≤K.(10)
其中,
K
K
K是数据点的个数,且
(
V
^
0
,
C
^
b
)
(\hat V_0,\hat C_b)
(V^0,C^b)是
(
V
0
,
C
b
)
(V_0,C_b)
(V0,Cb)经过LS拟合后的值,或者更精确地说,最小化观测到的子弹速度与它门的预测值之间的方差和
(
V
^
0
,
C
^
b
)
=
a
r
g
m
i
n
(
V
0
,
C
b
)
∑
k
=
1
K
[
v
^
(
x
k
)
−
(
V
0
η
−
2
η
C
b
−
1
x
k
)
1
/
η
]
2
.
(
11
)
(\hat V_0,\hat C_b)=\mathop{argmin}_{(V_0,C_b)}\sum_{k=1}^K[\hat v(x_k)-(V_0^\eta-2\eta C_b^{-1}x_k)^{1/\eta}]^2.(11)
(V^0,C^b)=argmin(V0,Cb)k=1∑K[v^(xk)−(V0η−2ηCb−1xk)1/η]2.(11)
其中
v
^
(
x
k
)
\hat v(x_k)
v^(xk)是在距离开火点
x
k
x_k
xk处观测到的子弹速度。从定义的速度数据得到的弹道模型的LS拟合的均方根偏差
ϵ
r
m
s
\epsilon_{rms}
ϵrms和最大绝对误差
ϵ
m
a
x
\epsilon_{max}
ϵmax分别如下:
ϵ
r
m
s
=
1
K
∑
k
=
1
K
[
v
^
(
x
k
)
−
v
^
L
S
(
x
k
)
]
2
.
(
12
)
\epsilon_{rms}=\sqrt {\frac1K\sum_{k=1}^K[\hat v(x_k)-\hat v_{LS}(x_k)]^2}.(12)
ϵrms=K1k=1∑K[v^(xk)−v^LS(xk)]2.(12)
ϵ m a x = max k ∣ v ^ ( x k ) − v ^ L S ( x k ) ∣ . ( 13 ) \epsilon_{max}=\max_k|\hat v(x_k)-\hat v_{LS}(x_k)|.(13) ϵmax=kmax∣v^(xk)−v^LS(xk)∣.(13)
依次用式(10)-(13)对这36个数据集进行处理,并分别使用本文中5个特定的弹道模型的 η \eta η值。非线性LS最小化(11)在MATLAB中实现使用最优化工具箱函数 l s q n o n l i n lsqnonlin lsqnonlin。
对所有的36中子弹类型,模型2和模型1间在计算RMS误差和MA误差的数值时都只存在一个较小的差值,这表明二阶速度模型和达克沃斯模型有相近的准确性。对于33和36类型子弹来说,二阶速度模型叫指数速度模型更为精准,且对所有的子弹来说,都比线性速度模型和常拉力模型更加精准。对于大约75%的子弹类型来说,线性速度模型和指数速度模型有相似的准确率。常拉力模型的准确率最低,其次是线性速度模型和指数速度模型。
每个弹道模型的准确率可以用全部36种子弹的 ϵ r m s \epsilon_{rms} ϵrms和 ϵ m a x \epsilon_{max} ϵmax的平均值和标准差来进行数量上的描述。如下图所示:[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-g45x7JZB-1594794664658)(C:\Users\weijiahao123\AppData\Roaming\Typora\typora-user-images\image-20200711172848282.png)]
推荐使用二阶速度模型或着达克沃斯模型来描述超音速子弹沿轨迹随着飞行距离变化,速度所产生的变化。然而,与达克沃斯模型相比,二阶速度模型具有的优势是SW信号的分离点可以通过解6阶多项式的相关根来有效地计算出来。
使用二阶速度模型和单传感器节点进行定位
已知的弹道模型参数
二阶速度模型被用来定位开火点时通常使用MB和SW信号,或只需要SW信号,假定这两个模型参数都是已知的。然而MB-SW方法只需要一个单独的传感器节点(其中包含一个麦克风阵列),SW方法要求一个传感器节点网或者分布式节点。使用MB-SW方法,开火点的距离通过测量MB和SW的到传感器节点到达时间差(TDOA)和到达角度(AOA);开火点的方向估计通过在节点上的每个传感器对测量MB的TDOA(因此方向的估计是不受应用的弹道模型影响的)。下图展示了超音速子弹的弹道、开火点、传感器节点。在图中,P是分离点(到达传感器节点的SW的产生点);
x
P
x_P
xP是从开火点到P点的距离;
θ
\theta
θ是P点的马赫锥角;
γ
\gamma
γ是MB和SW到节点处的AOA差;
α
\alpha
α和
d
d
d分别是节点的错过角度和错过距离。运用子弹的二阶速度模型,从节点到开火点的距离给出为:
r
=
x
p
sin
γ
(
V
0
1
/
2
−
C
b
−
1
x
p
)
4
−
c
2
(
V
0
1
/
2
−
C
b
−
1
x
p
)
2
.
(
14
)
r=\frac {x_p}{\sin \gamma}\frac{\sqrt{(V_0^{1/2}-C_b^{-1}x_p)^4-c^2}}{(V_0^{1/2}-C_b^{-1}x_p)^2}.(14)
r=sinγxp(V01/2−Cb−1xp)2(V01/2−Cb−1xp)4−c2.(14)
其中
x
p
x_p
xp可以通过解六阶多项式的相关根获得,此多项式的系数是
Δ
τ
\Delta \tau
Δτ,
γ
\gamma
γ,
V
0
V_0
V0和
C
b
C_b
Cb。如果常速度模型被应用到子弹上(即假设$C_b\to \infin $),那么上述该距离为:
r
=
c
Δ
τ
1
−
cos
γ
.
(
15
)
r=\frac {c\Delta\tau}{1-\cos\gamma}.(15)
r=1−cosγcΔτ.(15)
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-kFJQzCIq-1594794664664)(C:\Users\weijiahao123\AppData\Roaming\Typora\typora-user-images\image-20200711185320091.png)]
通过计算从一次战场实验记录的真实数据,表明当使用二阶速度模型而不是常速度模型时,MB-SW定位法的距离准确性可以明显提高。然而,二阶速度模型需要枪口速度和子弹的弹道常数的先验知识;这两个参数的任何不确定性都将导致定距准确性降低。
未知的弹道模型常数
一个子弹类型的弹道常数与另一个类型相比来说变化很大,即使子弹的口径相同。因此,子弹的弹道常数和它的口径之间没有一个一对一的关系。另一方面,子弹的枪口速度并不仅仅由子弹类型决定,它还与枪管长度有关。枪口速度随枪管长度增加而增加。因此,枪口速度和弹道常数的先验知识需要关于武器类型的先验知识作为基础。
在实际情况中,子弹的武器类型或二次速度模型参数未知时,可以使用传感器节点网络同时估计子弹的炮口速度和弹道常数以及开火点的距离。
长度增加而增加。因此,枪口速度和弹道常数的先验知识需要关于武器类型的先验知识作为基础。
在实际情况中,子弹的武器类型或二次速度模型参数未知时,可以使用传感器节点网络同时估计子弹的炮口速度和弹道常数以及开火点的距离。