火星19

弾道ミサイルの軌道（５）からの続き

ロイターが伝えた朝鮮中央通信のニュースによれば，10月31日に発射したのは大陸間弾道ミサイル「火星19」だとのこと。ICBMの最終完結形態で世界最強のミサイルをうたっている。知らんけど。飛行距離は1001.2km，飛行時間は5156秒，最高高度は7687.5km だ。

無意味なアラートを出すよりも，こういった基本的な情報がほしいところだ。しかし，防衛省が韓国軍や米軍からこれらの情報を入手してもことさらにマスコミに流す必要はないと考えているのだろう。日本政府にしてみれば，アラートを出すだけで危機意識の醸成という目的は達成される。

2年前の「火星17」については，同じロフテッド軌道で，飛行距離 999.2 km ，飛行時間 4135 秒，最高高度 6040.9 km だった。これに比べると，今回は飛行時間が1000秒伸びて，最高高度も1600kmほど高くなっている。

そこで，忘れかけていたMathematicaのコードを引っ張り出して再計算してみた。目の子でパラメータサーチをした。単純なモデルなので，そんなにピッタリ再現できるわけではないが，ほぼほぼ雰囲気は得られる。飛行時間を5156秒=86分あたりに固定すると，最高高度は170kmほど高くなって7850km（+2.1%）になってしまう。最高高度を7600km台にするためには，飛行時間を2分ほど減らして84分にする（-2.3%）必要がある。燃焼時加速度a，燃焼時間τ，燃料重量比pなどを変えても，飛行時間で制約がかかる範囲では，飛行高度を再現するには至らなかった。

これらのパラメタのままで，投射角度だけを45度にすれば（飛行時間Tは，地表に到達する点で調整する），大円到達距離は16500kmを越えてくる。

g = 0.0098; R = 6350; τ = 90.8; p = 0.75; a = 0.045; s = 86.7 Degree; T = 5160;

g = 0.0098; R = 6350; τ = 90.8; p = 0.75; a = 0.045; s = 45 Degree; T = 7050; 

(* g=0.0098; R=6350; τ=90.4; p=0.75; a=0.045;s=86.7Degree;T=5160-120; *)

(* g=0.0098; R=6350; τ=90.4; p=0.75; a=0.045;s=45 Degree;T=6850; *)

fr[t_, τ_] := a*Sin[s]*HeavisideTheta[τ - t]

ft[t_, τ_] := a*Cos[s]*r[t]*HeavisideTheta[τ - t]

fm[t_, τ_] := -p/(τ - p*t)*HeavisideTheta[τ - t]

sol = NDSolve[{r''[t] == -fm[t, τ]*r'[t] + h[t]^2/r[t]^3 - g R^2/r[t]^2 + fr[t, τ], 

r[0] == R, r'[0] == 0,  h'[t] == -fm[t, τ]*h[t] + ft[t, τ], h[0] == 0}, {r,  h}, {t, 0, T}]

f[t_] := r[t] /. sol[[1, 1]]

d[t_] := h[t] /. sol[[1, 2]]

Plot[{6350, f[t]}, {t, 0, T}]

Plot[{f[t + 1] - f[t], d[t]*R/f[t]^2, d[t]/f[t]}, {t, 0, T}, 

 PlotRange -> {-8, 8}]

{f[T/2] - 6350, tyx[T], v[T]/0.340}

→ {7851.42, {5160, -46.3658, 1007.58}, 24.5585}

g0 = ParametricPlot[{NIntegrate[R d[t]/f[t]^2 , {t, 0, tt}], f[tt] - R}, {tt, 0, T}]

→　{8231.43, {7050, 121.626, 16657.9}, 25.6423}

g9 = ParametricPlot[{NIntegrate[R d[t]/f[t]^2 , {t, 0, tt}], f[tt] - R}, {tt, 0, T}]

図：北朝鮮のミサイルの軌道（ロイター記事から引用）