用 Python 实现导弹自动追踪,超燃!

69 阅读9分钟

由于待会要用pygame演示,它的坐标系是y轴向下,所以这里我们也用y向下的坐标系。

算法总的思想就是根据上图,把时间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),正弦对比斜,余弦邻比斜,斜边可以利用两点距离公式计算出,即:

b763c9c34e48c59c94b3157efa202f6b.png

于是

4eee42f78ab0f57e6c4b5f04b58e3633.png

AC的长度就是导弹的速度乘以时间即 |AC|=vt,然后即可计算出AD和CD的长度,于是这一个时间片过去后,导弹应该出现在新的位置C点,他的坐标就是老的点A的x增加AD和y减去CD。

于是,新的C点坐标就是:

775f83d16c71cc0a766544ba575db343.png

只要一直反复循环执行这个操作即可,好吧,为了更形象,把第一个时间片和第二个时间片放在一起看看:

faaa334dfee54dd1db75e7d294837570.png

第一个是时间片构造出的三角形是ABE,经过一个时间片后,目标从B点走到了D点,导弹此时在C点,于是构造新的三角形CDF,重复刚才的计算过程即可,图中的角∠b就是导弹需要旋转的角度,现实中只需要每个时间片修正导弹的方向就可以了,具体怎么让导弹改变方向,这就不是我们需要研究的问题了

好,由于最近在用Python的pygame库制作小游戏玩,接下来我们就用pygame来演示一下这个效果,效果如下图:

6d4412ceb817ee85008427d74b84ddb5.gif

很简单的代码如下:

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+sectioncosa,y1-sectionsina)

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里面做旋转并不是一件容易的事情,我们先把图片替换成一张矩形的,再加入旋转函数看看效果如何

0951e33c2274b185de24efb98f8657c2.png

missiled = pygame.transform.rotate(missile, -(d_angle))

screen.blit(missiled, (x1-missile.get_width(), y1-

missile.get_height()/2))

因为图片的坐标点是它的左上角的点,所以如果我们想让图片的坐标固定在箭头尖点,那么把图片实际打印位置x减少图片长度,y减少一半宽度就行。但是实际运行效果并不好:

4ee63eb678c3eb07a47e7355e3b2830a.gif

大致方向相同,但是图片箭头的尖点并没有一直跟随鼠标,这是为什么呢。经过一番研究,我发现原来是这个图旋转的机制问题,我们看看旋转后的图片变成什么样了:

049213e6ecfe7a70065c0be40b165144.png

旋转后的图片变成了蓝色的那个范围,根据旋转角度的不同,所变成的图片大小也不一样,我们看旋转90的情况

5989a58ffa9eff1ee49d41bd27905fd3.png

8b5d3f0faee8447a812477021dd0e611.png

我们发现,旋转后的图片不仅面积变大了,导弹头的位置也变了。那应该怎么解决这个问题呢?思路是,每一次旋转图片以后,求出旋转图的头位置(图中的绿色箭头点),然后把绿图的打印位置移动一下,下,x,y分别移动两个头的距离,就可以让旋转后的导弹头对准实际我们参与运算的那个导弹头的位置,移动后应该是这样的:

da5c7625a19fc4542b70b3488be4e505.png

这样,两个导弹头的点就一致了。接下来我们分析求旋转后的导弹头的算法。根据旋转角度的不同,旋转角在不同象限参数不一样,所以我们分为这四种情况

1,2象限

4f3c4eced8d74979e413de8ea47b4386.png

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

d858510eebda86d46d6ac603fc872f3b.png

显示图片的时候我们将他移动

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

这里的 (x1-width, y1-height/2) 其实才是上图中的 (x1, y1) 

所以最后我们加入相关算法代码,效果就比较完美了

1c303e0c36fe3b60dc3a277ca457f7be.gif

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

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+sectioncosa,y1-sectionsina)

missiled=pygame.transform.rotate(missile,-(fangle))

if 0<=-fangle<=90:

A=(width*cosa+x1-width,y1-height/2)

B=(A[0]+heightsina,A[1]+heightcosa)

if 90<-fangle<=180:

A = (x1 - width, y1 - height/2+height*(-cosa))

B = (x1 - width+height*sina, y1 - height/2)

最后

Python崛起并且风靡,因为优点多、应用领域广、被大牛们认可。学习 Python 门槛很低,但它的晋级路线很多,通过它你能进入机器学习、数据挖掘、大数据,CS等更加高级的领域。Python可以做网络应用,可以做科学计算,数据分析,可以做网络爬虫,可以做机器学习、自然语言处理、可以写游戏、可以做桌面应用…Python可以做的很多,你需要学好基础,再选择明确的方向。这里给大家分享一份全套的 Python 学习资料,给那些想学习 Python 的小伙伴们一点帮助!

👉Python所有方向的学习路线👈

Python所有方向的技术点做的整理,形成各个领域的知识点汇总,它的用处就在于,你可以按照上面的知识点去找对应的学习资源,保证自己学得较为全面。

👉Python必备开发工具👈

工欲善其事必先利其器。学习Python常用的开发软件都在这里了,给大家节省了很多时间。

👉Python全套学习视频👈

我们在看视频学习的时候,不能光动眼动脑不动手,比较科学的学习方法是在理解之后运用它们,这时候练手项目就很适合了。

👉实战案例👈

学python就与学数学一样,是不能只看书不做题的,直接看步骤和答案会让人误以为自己全都掌握了,但是碰到生题的时候还是会一筹莫展。

因此在学习python的过程中一定要记得多动手写代码,教程只需要看一两遍即可。

👉大厂面试真题👈

我们学习Python必然是为了找到高薪的工作,下面这些面试题是来自阿里、腾讯、字节等一线互联网大厂最新的面试资料,并且有阿里大佬给出了权威的解答,刷完这一套面试资料相信大家都能找到满意的工作。

了解详情:docs.qq.com/doc/DSnl3ZG…