主頁 > 後端開發 > 超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

2022-03-04 06:12:23 後端開發

超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

 

大家好,我是J哥,

這個沒有點數學基礎是很難算出來的,但是我們有了計算機就不一樣了,依靠計算機極快速的運算速度,我們利用微分的思想,加上一點簡單的三角學知識,就可以實作它,

好,話不多說,我們來看看它的演算法原理,看圖:

超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

 

由于待會要用pygame演示,它的坐標系是y軸向下,所以這里我們也用y向下的坐標系,

干貨主要有:

① 200 多本 Python 電子書(和經典的書籍)應該有

② Python標準庫資料(最全中文版)

③ 專案原始碼(四五十個有趣且可靠的練手專案及原始碼)

④ Python基礎入門、爬蟲、網路開發、大資料分析方面的視頻(適合小白學習)

⑤ Python學習路線圖(告別不入流的學習)

Python學習交流Q群101677771

 

演算法總的思想就是根據上圖,把時間t分割成足夠小的片段(比如1/1000,這個時間片越小越精確),每一個片段分別構造如上三角形,計算出導彈下一個時間片走的方向(即∠a)和走的路程(即vt=|AC|),這時候目標再在第二個時間片移動了位置,這時剛才計算的C點又變成了第二個時間片的初始點,這時再在第二個時間片上在C點和新的目標點構造三角形計算新的vt,然后進入第三個時間片,如此反復即可,

假定導彈和目標的初始狀態下坐標分別是(x1,y1),(x,y),構造出直角三角形ABE,這個三角形用來求∠a的正弦和余弦值,因為vt是自己設定的,我們需要計算A到C點x和y坐標分別移動了多少,移動的值就是AD和CD的長度,這兩個分別用vt乘cos(a)和sin(a)即可,

計算sin(a)和cos(a),正弦對比斜,余弦鄰比斜,斜邊可以利用兩點距離公式計算出,即:

超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

 

于是

超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

 

AC的長度就是導彈的速度乘以時間即 |AC|=vt,然后即可計算出AD和CD的長度,于是這一個時間片過去后,導彈應該出現在新的位置C點,他的坐標就是老的點A的x增加AD和y減去CD,

于是,新的C點坐標就是:

超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

 

只要一直反復回圈執行這個操作即可,好吧,為了更形象,把第一個時間片和第二個時間片放在一起看看:

超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

 

第一個是時間片構造出的三角形是ABE,經過一個時間片后,目標從B點走到了D點,導彈此時在C點,于是構造新的三角形CDF,重復剛才的計算程序即可,圖中的角∠b就是導彈需要旋轉的角度,現實中只需要每個時間片修正導彈的方向就可以了,具體怎么讓導彈改變方向,這就不是我們需要研究的問題了

好,由于最近在用Python的pygame庫制作小游戲玩,接下來我們就用pygame來演示一下這個效果,效果如下圖:

超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

 

很簡單的代碼如下:

import pygame,sys
from math import *
pygame.init()
screen=pygame.display.set_mode((800,700),0,32)
missile=pygame.image.load('element/red_pointer.png').convert_alpha()
x1,y1=100,600           #導彈的初始發射位置
velocity=800            #導彈速度
time=1/1000             #每個時間片的長度
clock=pygame.time.Clock()
old_angle=0
while True:
    for event in pygame.event.get():
        if event.type==pygame.QUIT:
            sys.exit()
    clock.tick(300)
    x,y=pygame.mouse.get_pos()          #獲取滑鼠位置,滑鼠就是需要打擊的目標
    distance=sqrt(pow(x1-x,2)+pow(y1-y,2))      #兩點距離公式
    section=velocity*time               #每個時間片需要移動的距離
    sina=(y1-y)/distance
    cosa=(x-x1)/distance
    angle=atan2(y-y1,x-x1)              #兩點線段的弧度值
    x1,y1=(x1+section*cosa,y1-section*sina)
    d_angle = degrees(angle)        #弧度轉角度
    screen.blit(missile, (x1-missile.get_width(), y1-missile.get_height()/2))
    dis_angle=d_angle-old_angle          #dis_angle就是到下一個位置需要改變的角度
    old_angle=d_angle                    #更新初始角度
    pygame.display.update()

如果僅把導彈考慮為一個質點的話,那么以上演算法就已經足矣,我沒有做導彈的旋轉,因為一個質點也不分頭尾不需要旋轉,當然這前提得是你加載的導彈圖片很小的時候不旋轉看起來也沒什么問題,但是在pygame里面做旋轉并不是一件容易的事情,我們先把圖片替換成一張矩形的,再加入旋轉函式看看效果如何

超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

 

missiled = pygame.transform.rotate(missile, -(d_angle))
screen.blit(missiled, (x1-missile.get_width(), y1-
missile.get_height()/2))

因為圖片的坐標點是它的左上角的點,所以如果我們想讓圖片的坐標固定在箭頭尖點,那么把圖片實際列印位置x減少圖片長度,y減少一半寬度就行,但是實際運行效果并不好:

超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

 

大致方向相同,但是圖片箭頭的尖點并沒有一直跟隨滑鼠,這是為什么呢,經過一番研究,我發現原來是這個圖旋轉的機制問題,我們看看旋轉后的圖片變成什么樣了:

超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

 

旋轉后的圖片變成了藍色的那個范圍,根據旋轉角度的不同,所變成的圖片大小也不一樣,我們看旋轉90的情況

超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

 

超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

 

我們發現,旋轉后的圖片不僅面積變大了,導彈頭的位置也變了,那應該怎么解決這個問題呢?思路是,每一次旋轉圖片以后,求出旋轉圖的頭位置(圖中的綠色箭頭點),然后把綠圖的列印位置移動一下,下,x,y分別移動兩個頭的距離,就可以讓旋轉后的導彈頭對準實際我們參與運算的那個導彈頭的位置,移動后應該是這樣的:

超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

 

這樣,兩個導彈頭的點就一致了,接下來我們分析求旋轉后的導彈頭的演算法,根據旋轉角度的不同,旋轉角在不同象限引數不一樣,所以我們分為這四種情況

1,2象限

超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

 

3,4象限,它的旋轉只有正負0—180,所以3,4象限就是負角

超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

 

顯示圖片的時候我們將他移動

screen.blit(missiled, (x1-width+(x1-C[0]),y1-height/2+(y1-C[1])))

這里的 (x1-width, y1-height/2) 其實才是上圖中的 (x1, y1)

所以最后我們加入相關演算法代碼,效果就比較完美了

超燃!利用 Python 實作導彈自動追蹤,這個演算法幾人能看懂

 

大功告成,最后附上全部的演算法代碼

import pygame,sys
from math import *
pygame.init()
font1=pygame.font.SysFont('microsoftyaheimicrosoftyaheiui',23)
textc=font1.render('*',True,(250,0,0))
screen=pygame.display.set_mode((800,700),0,32)
missile=pygame.image.load('element/rect1.png').convert_alpha()
height=missile.get_height()
width=missile.get_width()
pygame.mouse.set_visible(0)
x1,y1=100,600           #導彈的初始發射位置
velocity=800            #導彈速度
time=1/1000             #每個時間片的長度
clock=pygame.time.Clock()
A=()
B=()
C=()
while True:
    for event in pygame.event.get():
        if event.type==pygame.QUIT:
            sys.exit()
    clock.tick(300)
    x,y=pygame.mouse.get_pos()          #獲取滑鼠位置,滑鼠就是需要打擊的目標
    distance=sqrt(pow(x1-x,2)+pow(y1-y,2))      #兩點距離公式
    section=velocity*time               #每個時間片需要移動的距離
    sina=(y1-y)/distance
    cosa=(x-x1)/distance
    angle=atan2(y-y1,x-x1)              #兩點間線段的弧度值
    fangle=degrees(angle)               #弧度轉角度
    x1,y1=(x1+section*cosa,y1-section*sina)
    missiled=pygame.transform.rotate(missile,-(fangle))
    if 0<=-fangle<=90:
        A=(width*cosa+x1-width,y1-height/2)
        B=(A[0]+height*sina,A[1]+height*cosa)
    if 90<-fangle<=180:
        A = (x1 - width, y1 - height/2+height*(-cosa))
        B = (x1 - width+height*sina, y1 - height/2)
    if -90<=-fangle<0:
        A = (x1 - width+missiled.get_width(), y1 - height/2+missiled.get_height()-height*cosa)
        B = (A[0]+height*sina, y1 - height/2+missiled.get_height())
    if -180<-fangle<-90:
        A = (x1-width-height*sina, y1 - height/2+missiled.get_height())
        B = (x1 - width,A[1]+height*cosa )
    C = ((A[0] + B[0]) / 2, (A[1] + B[1]) / 2)
    screen.fill((0,0,0))
    screen.blit(missiled, (x1-width+(x1-C[0]),y1-height/2+(y1-C[1])))
    screen.blit(textc, (x,y)) #滑鼠用一個紅色*代替
    pygame.display.update()

以上便是用Python模擬導彈自動追蹤的代碼實體,

 

 

轉載請註明出處,本文鏈接:https://www.uj5u.com/houduan/436950.html

標籤:Python

上一篇:利用Python搶回在螞蟻森林把逝去的能量

下一篇:Python基礎入門:正則re.sub使用自定義替換方法

標籤雲
其他(157675) Python(38076) JavaScript(25376) Java(17977) C(15215) 區塊鏈(8255) C#(7972) AI(7469) 爪哇(7425) MySQL(7132) html(6777) 基礎類(6313) sql(6102) 熊猫(6058) PHP(5869) 数组(5741) R(5409) Linux(5327) 反应(5209) 腳本語言(PerlPython)(5129) 非技術區(4971) Android(4554) 数据框(4311) css(4259) 节点.js(4032) C語言(3288) json(3245) 列表(3129) 扑(3119) C++語言(3117) 安卓(2998) 打字稿(2995) VBA(2789) Java相關(2746) 疑難問題(2699) 细绳(2522) 單片機工控(2479) iOS(2429) ASP.NET(2402) MongoDB(2323) 麻木的(2285) 正则表达式(2254) 字典(2211) 循环(2198) 迅速(2185) 擅长(2169) 镖(2155) 功能(1967) .NET技术(1958) Web開發(1951) python-3.x(1918) HtmlCss(1915) 弹簧靴(1913) C++(1909) xml(1889) PostgreSQL(1872) .NETCore(1853) 谷歌表格(1846) Unity3D(1843) for循环(1842)

熱門瀏覽
  • 【C++】Microsoft C++、C 和匯編程式檔案

    ......

    uj5u.com 2020-09-10 00:57:23 more
  • 例外宣告

    相比于斷言適用于排除邏輯上不可能存在的狀態,例外通常是用于邏輯上可能發生的錯誤。 例外宣告 Item 1:當函式不可能拋出例外或不能接受拋出例外時,使用noexcept 理由 如果不打算拋出例外的話,程式就會認為無法處理這種錯誤,并且應當盡早終止,如此可以有效地阻止例外的傳播與擴散。 示例 //不可 ......

    uj5u.com 2020-09-10 00:57:27 more
  • Codeforces 1400E Clear the Multiset(貪心 + 分治)

    鏈接:https://codeforces.com/problemset/problem/1400/E 來源:Codeforces 思路:給你一個陣列,現在你可以進行兩種操作,操作1:將一段沒有 0 的區間進行減一的操作,操作2:將 i 位置上的元素歸零。最終問:將這個陣列的全部元素歸零后操作的最少 ......

    uj5u.com 2020-09-10 00:57:30 more
  • UVA11610 【Reverse Prime】

    本人看到此題沒有翻譯,就附帶了一個自己的翻譯版本 思考 這一題,它的第一個要求是找出所有 $7$ 位反向質數及其質因數的個數。 我們應該需要質數篩篩選1~$10^{7}$的所有數,這里就不慢慢介紹了。但是,重讀題,我們突然發現反向質數都是 $7$ 位,而將它反過來后的數字卻是 $6$ 位數,這就說明 ......

    uj5u.com 2020-09-10 00:57:36 more
  • 統計區間素數數量

    1 #pragma GCC optimize(2) 2 #include <bits/stdc++.h> 3 using namespace std; 4 bool isprime[1000000010]; 5 vector<int> prime; 6 inline int getlist(int ......

    uj5u.com 2020-09-10 00:57:47 more
  • C/C++編程筆記:C++中的 const 變數詳解,教你正確認識const用法

    1、C中的const 1、區域const變數存放在堆疊區中,會分配記憶體(也就是說可以通過地址間接修改變數的值)。測驗代碼如下: 運行結果: 2、全域const變數存放在只讀資料段(不能通過地址修改,會發生寫入錯誤), 默認為外部聯編,可以給其他源檔案使用(需要用extern關鍵字修飾) 運行結果: ......

    uj5u.com 2020-09-10 00:58:04 more
  • 【C++犯錯記錄】VS2019 MFC添加資源不懂如何修改資源宏ID

    1. 首先在資源視圖中,添加資源 2. 點擊新添加的資源,復制自動生成的ID 3. 在解決方案資源管理器中找到Resource.h檔案,編輯,使用整個專案搜索和替換的方式快速替換 宏宣告 4. Ctrl+Shift+F 全域搜索,點擊查找全部,然后逐個替換 5. 為什么使用搜索替換而不使用屬性視窗直 ......

    uj5u.com 2020-09-10 00:59:11 more
  • 【C++犯錯記錄】VS2019 MFC不懂的批量添加資源

    1. 打開資源頭檔案Resource.h,在其中預先定義好宏 ID(不清楚其實ID值應該設定多少,可以先新建一個相同的資源項,再在這個資源的ID值的基礎上遞增即可) 2. 在資源視圖中選中專案資源,按F7編輯資源檔案,按 ID 型別 相對路徑的形式添加 資源。(別忘了先把檔案拷貝到專案中的res檔案 ......

    uj5u.com 2020-09-10 01:00:19 more
  • C/C++編程筆記:關于C++的參考型別,專供新手入門使用

    今天要講的是C++中我最喜歡的一個用法——參考,也叫別名。 參考就是給一個變數名取一個變數名,方便我們間接地使用這個變數。我們可以給一個變數創建N個參考,這N + 1個變數共享了同一塊記憶體區域。(參考型別的變數會占用記憶體空間,占用的記憶體空間的大小和指標型別的大小是相同的。雖然參考是一個物件的別名,但 ......

    uj5u.com 2020-09-10 01:00:22 more
  • 【C/C++編程筆記】從頭開始學習C ++:初學者完整指南

    眾所周知,C ++的學習曲線陡峭,但是花時間學習這種語言將為您的職業帶來奇跡,并使您與其他開發人員區分開。您會更輕松地學習新語言,形成真正的解決問題的技能,并在編程的基礎上打下堅實的基礎。 C ++將幫助您養成良好的編程習慣(即清晰一致的編碼風格,在撰寫代碼時注釋代碼,并限制類內部的可見性),并且由 ......

    uj5u.com 2020-09-10 01:00:41 more
最新发布
  • Rust中的智能指標:Box<T> Rc<T> Arc<T> Cell<T> RefCell<T> Weak

    Rust中的智能指標是什么 智能指標(smart pointers)是一類資料結構,是擁有資料所有權和額外功能的指標。是指標的進一步發展 指標(pointer)是一個包含記憶體地址的變數的通用概念。這個地址參考,或 ” 指向”(points at)一些其 他資料 。參考以 & 符號為標志并借用了他們所 ......

    uj5u.com 2023-04-20 07:24:10 more
  • Java的值傳遞和參考傳遞

    值傳遞不會改變本身,參考傳遞(如果傳遞的值需要實體化到堆里)如果發生修改了會改變本身。 1.基本資料型別都是值傳遞 package com.example.basic; public class Test { public static void main(String[] args) { int ......

    uj5u.com 2023-04-20 07:24:04 more
  • [2]SpinalHDL教程——Scala簡單入門

    第一個 Scala 程式 shell里面輸入 $ scala scala> 1 + 1 res0: Int = 2 scala> println("Hello World!") Hello World! 檔案形式 object HelloWorld { /* 這是我的第一個 Scala 程式 * 以 ......

    uj5u.com 2023-04-20 07:23:58 more
  • 理解函式指標和回呼函式

    理解 函式指標 指向函式的指標。比如: 理解函式指標的偽代碼 void (*p)(int type, char *data); // 定義一個函式指標p void func(int type, char *data); // 宣告一個函式func p = func; // 將指標p指向函式func ......

    uj5u.com 2023-04-20 07:23:52 more
  • Django筆記二十五之資料庫函式之日期函式

    本文首發于公眾號:Hunter后端 原文鏈接:Django筆記二十五之資料庫函式之日期函式 日期函式主要介紹兩個大類,Extract() 和 Trunc() Extract() 函式作用是提取日期,比如我們可以提取一個日期欄位的年份,月份,日等資料 Trunc() 的作用則是截取,比如 2022-0 ......

    uj5u.com 2023-04-20 07:23:45 more
  • 一天吃透JVM面試八股文

    什么是JVM? JVM,全稱Java Virtual Machine(Java虛擬機),是通過在實際的計算機上仿真模擬各種計算機功能來實作的。由一套位元組碼指令集、一組暫存器、一個堆疊、一個垃圾回收堆和一個存盤方法域等組成。JVM屏蔽了與作業系統平臺相關的資訊,使得Java程式只需要生成在Java虛擬機 ......

    uj5u.com 2023-04-20 07:23:31 more
  • 使用Java接入小程式訂閱訊息!

    更新完微信服務號的模板訊息之后,我又趕緊把微信小程式的訂閱訊息給實作了!之前我一直以為微信小程式也是要企業才能申請,沒想到小程式個人就能申請。 訊息推送平臺🔥推送下發【郵件】【短信】【微信服務號】【微信小程式】【企業微信】【釘釘】等訊息型別。 https://gitee.com/zhongfuch ......

    uj5u.com 2023-04-20 07:22:59 more
  • java -- 緩沖流、轉換流、序列化流

    緩沖流 緩沖流, 也叫高效流, 按照資料型別分類: 位元組緩沖流:BufferedInputStream,BufferedOutputStream 字符緩沖流:BufferedReader,BufferedWriter 緩沖流的基本原理,是在創建流物件時,會創建一個內置的默認大小的緩沖區陣列,通過緩沖 ......

    uj5u.com 2023-04-20 07:22:49 more
  • Java-SpringBoot-Range請求頭設定實作視頻分段傳輸

    老實說,人太懶了,現在基本都不喜歡寫筆記了,但是網上有關Range請求頭的文章都太水了 下面是抄的一段StackOverflow的代碼...自己大修改過的,寫的注釋挺全的,應該直接看得懂,就不解釋了 寫的不好...只是希望能給視頻網站開發的新手一點點幫助吧. 業務場景:視頻分段傳輸、視頻多段傳輸(理 ......

    uj5u.com 2023-04-20 07:22:42 more
  • Windows 10開發教程_編程入門自學教程_菜鳥教程-免費教程分享

    教程簡介 Windows 10開發入門教程 - 從簡單的步驟了解Windows 10開發,從基本到高級概念,包括簡介,UWP,第一個應用程式,商店,XAML控制元件,資料系結,XAML性能,自適應設計,自適應UI,自適應代碼,檔案管理,SQLite資料庫,應用程式到應用程式通信,應用程式本地化,應用程式 ......

    uj5u.com 2023-04-20 07:22:35 more