基于势能的平面运动模拟

运动模拟

接前面

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

里面说要讲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

模拟

相关推荐
右恩4 分钟前
Docker 实践与应用举例
python·docker
凤枭香8 分钟前
Python Scikit-learn简介(二)
开发语言·python·机器学习·scikit-learn
湫ccc8 小时前
《Python基础》之字符串格式化输出
开发语言·python
mqiqe8 小时前
Python MySQL通过Binlog 获取变更记录 恢复数据
开发语言·python·mysql
AttackingLin8 小时前
2024强网杯--babyheap house of apple2解法
linux·开发语言·python
哭泣的眼泪4089 小时前
解析粗糙度仪在工业制造及材料科学和建筑工程领域的重要性
python·算法·django·virtualenv·pygame
湫ccc9 小时前
《Python基础》之基本数据类型
开发语言·python
drebander10 小时前
使用 Java Stream 优雅实现List 转化为Map<key,Map<key,value>>
java·python·list
威威猫的栗子11 小时前
Python Turtle召唤童年:喜羊羊与灰太狼之懒羊羊绘画
开发语言·python
墨染风华不染尘11 小时前
python之开发笔记
开发语言·笔记·python