朝韓危機-Python模擬導彈互射

 

薩德系統是麻省理工學院的林肯實驗室牽頭開發的。林肯實驗室能夠說是美軍方手中的一張王牌。二個反導系統（末段高空區域導彈防護(THAAD)系統，和 國家導彈防護系統(NMD)）均由實驗室牽頭。它既是技術支撐方,同時又是關鍵技術的研製方,例如THAAD地基雷達(GBR),實驗室負責任務書的擬定,並負責雷達的驗收,而其中的目標識別算法等又是實驗室提供並負責試驗數據的錄取和分析。

對於在韓國部署薩德的意義和目的，著名導彈專家、美國麻省理工學院終身教授泊斯托爾曾一針見血地指出：薩德對於防範朝鮮的勞動導彈並無論用，但卻會成爲美國國家導彈防護體系（NMD）的重要組成部分。

 

發表一下政治觀點：看了戰爭之王的朋友能夠理解，和平是軍火商的噩夢。爲了賺取高額軍火利潤，美國軍火商要不停製造全球仇恨和緊張。美國在亞太不停唆使離間各個中，日，韓，朝鮮，菲律賓，製造仇恨和衝突。

 洛克希德馬丁公司是美國知名軍火商，利潤每一年上百億。薩德系統就是洛克希德馬丁的產品。

不少時候，爲了賺錢，美國政客和軍火商要不停製造全球衝突和仇恨。

（洛克希德馬丁）

切入正題，談談科學。。。。。。。

要模擬朝韓導彈互射，就要先明白彈道軌跡。咱們用python模擬最簡單的彈道軌跡。但願能幫助各位初學者學習python數學建模和matplotlib動態可視化模擬。

數學建模要用導數知識：

感謝英國大神牛頓和德國大神萊布尼茨的導數求最值方法，當導彈的瞬時速度爲0時，導彈高度達到最高值（峯值），看不懂的能夠去補補微積分知識，高中課本就能看懂。

（膜拜大神，左:牛頓，右:萊布尼茨）

 

  

Python導入math模塊，表示飛行時間t_flight：

t_flight =2*u*math.sin(theta_radians)/g

 

這是代碼運行的界面

 

運行後能夠觀察彈道數據，設置不一樣發射速度和角度能夠獲得不一樣結果。

 

生成的動態圖：

生成動態圖須要導入matplotlib模塊。

說明此語句意思animation.FuncAnimation(fig, update,generate,interval=5)

animation.FuncAnimation函數用於生成動態圖片。fig是生成的圖表對象，generate函數生成數據後傳遞給update函數更新，這樣數據不斷更新，圖形也不停變化。

interval表示時間間隔，設置的值越小，運動速度越快。

（導彈發射軌跡的動態模擬）

 

固然軍事上導彈軌跡比理論的要複雜，考慮風速等等，爲了簡化內容，文章暫時不涉及。

「薩德」系統若是進駐朝鮮半島，受影響的不只是朝鮮，中國和俄羅斯遠東地區部署的戰略戰術導彈也將面臨「威力貶值」的風險。除致使中國洲際導彈威懾力降低外，有專家分析，「薩德」部署在韓國境內，比起部署日本，能向西向北推動300千米的探測距離，從而更方便地監控中國東北、華北、東南沿海大部分地區。從理論上說，中國在東部沿海地區機動發射的「東風-21D」、「東風-26」系列反艦彈道導彈面臨「提早曝光」的危險。

薩德事件，日本人不斷煽風點火。

本做者看不下去了，決定模擬北京導彈先發制人，同時打掉薩德部署地和日本幾個大城市，東京，長崎，廣島，大阪。

咱們用python的matplotlib和basemap來完成這個模擬。matplotlib是Python經常使用的數據繪製包。它基於numpy的數組運算功能。matplotlib繪圖功能強大，能夠輕易的畫出各類統計圖形，好比散點圖，條行圖，餅圖等。matplotlib常與numpy和scipy相配合，用於許多研究領域。他們是免費工具，但其功能足能夠與科研界的大佬Matlab競爭。Basemap是Matplotlib的一個子包，負責地圖繪製。在數據可視化過程當中，咱們常須要將數據在地圖上畫出來。好比說咱們在地圖上畫出城市人口，飛機航線，軍事基地，礦藏分佈等等。這樣的地理繪圖有助於讀者理解空間相關的信息。

（basemap可用於繪製衛星地圖）

 

先爲程序收集幾個城市經緯度參數。

#東京的經緯度
tokyolat,tokyolon=35.42,139.46
#廣島座標
Hiroshima_lat, Hiroshima_lon=34.24,132.27
#大阪座標
Osaka_lat,Osaka_lon= 34.4,135.3
#長崎座標
Nagasaki_lat,Nagasaki_lon=32.46,129.52

Python程序寫好後，命名爲missle.py，而後執行腳本

 

python獲得了下面地圖

 

任務完成，最後我想說的是軍事競賽帶來的是人類種族殘殺和自我毀滅。但願數學領域高材生不要被金錢誘惑，遠離美國軍火商，珍惜和平，珍惜生命！