基于势能的平面运动模拟

运动模拟

接前面

递归调用------单向汉诺塔

里面说要讲Python解方程。Python还可以做数值计算进行物理现象的模拟。

七边形最小分割弦

等多篇都用到了公式。

F = m a F=ma F=ma

物体加速度的大小跟作用力成正比,跟物体的质量成反比,且与物体质量的倒数成正比;加速度的方向跟作用力的方向相同。 F 和 a F和a F和a的方向相同。

F = m a = m d 2 r d t 2 F = ma = m \frac{{{d^2}\ r}}{{d{t^2}}} F=ma=mdt2d2 r

a = Δ v / Δ t a = Δv/Δt a=Δv/Δt得出上面的公式。模拟时用0.01秒距离的变化表示加速度。

F = k Q q r 2 r ^ F= k \frac{{Qq}}{{{r^2}}} \hat{r} F=kr2Qqr^

库仑定律(Coulomb's law)是静止点电荷相互作用力的规律。1785年由法国科学家C,-A.de库仑根据实验得出,真空中两个静止的点电荷之间的相互作用力与它们的电荷量的乘积成正比,与它们的距离的二次方成反比,作用力的方向在它们的连线上,同名电荷相斥,异名电荷相吸。

代码

python 复制代码
import numpy as np
import matplotlib.pyplot as plt
import math
#常数
Q = (((0.0,0.0),10.0),((5.0,-5.0),5.0))#电荷的位置和值
TIMELIMIT = 20.0 #模拟终止的时刻
RLIMIT = 0.1     #距离r的最小值
H = 0.01         #时刻的跨度
t = 0.0
#画出电荷的位置
for qi in Q:
    plt.plot(qi[0][0],qi[0][1],".")

vx,vy,x,y = -2,-1,2,2
#当前时刻与位置
#print("{:.7f}{:.7f}{:.7f}{:.7f}{:.7f}".format(t,x,y,vx,vy))
#追加现在位置
xlist = [x]
ylist = [y]
#平面运动的计算
while t < TIMELIMIT:
    t = t + H
    rmin = float("inf") #将距离最小值初始化
    for qi in Q:
        rx = qi[0][0] - x #计算距离rx
        ry = qi[0][1] - y #计算距离ry
        r = math.sqrt(rx*rx + ry*ry) #计算距离r
        if r < rmin:
            rmin = r #更新距离的最小值
        vx += (rx/r/r/r*qi[1])*H#计算速度vx
        vy += (ry/r/r/r*qi[1])*H#计算速度vy
    x += vx*H #计算位置x
    y += vy*H #计算位置y
    #print("{:.7f}{:.7f}{:.7f}{:.7f}{:.7f}".format(t,x,y,vx,vy))
    #追加现在位置
    xlist.append(x)
    ylist.append(y)
    if rmin < RLIMIT:
        break #无限接近电荷时终止
#图表表示
plt.plot(xlist,ylist)
plt.show()

电荷的位置和值

一个10个电荷在原点,一个5个电荷在第四象限。
Q = ( ( ( 0.0 , 0.0 ) , 10.0 ) , ( ( 5.0 , − 5.0 ) , 5.0 ) ) Q = (((0.0,0.0),10.0),((5.0,-5.0),5.0)) Q=(((0.0,0.0),10.0),((5.0,−5.0),5.0))

运动电荷的初始位置和速度

vx,vy,x,y = -2,-1,2,2

计算距离

r = math.sqrt(rxrx + ryry) #计算距离r

计算速度

vx += (rx/r/r/r*qi[1])*H#计算速度vx

这里连续除以三次r,应该是对应上面的库伦定理。

打印

打印的时刻,坐标,坐标方向的速度。截取了开头结尾的部分数据。

0.00000002.00000002.0000000-2.0000000-1.0000000

0.01000001.97991501.9899037-2.0084992-1.0096312

0.02000001.95974391.9797095-2.0171083-1.0194191

0.03000001.93948561.9694158-2.0258299-1.0293684

0.04000001.91913901.9590210-2.0346668-1.0394842

0.05000001.89870271.9485233-2.0436221-1.0497717

...

19.97000005.96697103.0464395-0.41700850.8300425

19.98000005.96278023.0547222-0.41908360.8282736

19.99000005.95856863.0629873-0.42115650.8265034

20.00000005.95433633.0712346-0.42322750.8247318

模拟

相关推荐
豌豆花下猫1 小时前
让 Python 代码飙升330倍:从入门到精通的四种性能优化实践
后端·python·ai
夏末蝉未鸣011 小时前
python transformers库笔记(BertForTokenClassification类)
python·自然语言处理·transformer
weixin_418813873 小时前
Python-可视化学习笔记
笔记·python·学习
Danceful_YJ3 小时前
4.权重衰减(weight decay)
python·深度学习·机器学习
Zonda要好好学习4 小时前
Python入门Day5
python
电商数据girl5 小时前
有哪些常用的自动化工具可以帮助处理电商API接口返回的异常数据?【知识分享】
大数据·分布式·爬虫·python·系统架构
CoooLuckly5 小时前
numpy数据分析知识总结
python·numpy
超龄超能程序猿5 小时前
(六)PS识别:源数据分析- 挖掘图像的 “元语言”技术实现
python·组合模式
amazinging6 小时前
北京-4年功能测试2年空窗-报培训班学测开-第四十四天
python·学习·appium
UrbanJazzerati6 小时前
Xlwings安装报错:Connection timed out & WinError 32?一招解决你的安装难题!
python