Python抛体运动

最新推荐文章于 2024-07-31 17:05:31 发布

River Chandler

最新推荐文章于 2024-07-31 17:05:31 发布

阅读量907

点赞数 2

分类专栏：经典力学与分析力学文章标签： python

本文链接：https://blog.csdn.net/Chandler_river/article/details/123810204

版权

经典力学与分析力学专栏收录该内容

13 篇文章 5 订阅

订阅专栏

import numpy as np
from math import pi,sin,cos,radians
import matplotlib.pyplot as plt
def main():
    XX = []
    YY = []      
    angle=eval(input("Enter the lauch angle(in degrees):"))
    vel=eval(input("Enter the initial velocity(in meters?sec):"))
    h0=eval(input("Enter the initial height(in meters):"))
    time=eval(input("Enter the time interval:"))
    k=eval(input("Enter the k"))
    n=eval(input("Enter the n"))
    m=eval(input("Enter the mass")) 
    theta=radians(angle)
    xvel=vel*cos(theta)
    yvel=vel*sin(theta)
    for i in range(0,n+1):
        xpos=0
        ypos=h0
        print(i)
        if i != 0:
            
            while ypos>=0:
                xpos=xpos+time*xvel
                XX.append(round(xpos,3))
                xvell=xvel-time*k/m*xvel*(xvel**2+yvel**2)**((i-1)/2)
                yvel1=yvel-time*9.8-time*k/m*yvel*(xvel**2+yvel**2)**((i-1)/2)
                ypos=ypos+time*(yvel+yvel1)/2
                YY.append(round(ypos,3))
                yvel=yvel1
                xvel=xvell
            
            plt.plot(XX,YY)
            print("\nistance traveled meters.",format(xpos))
        else:
            while ypos>=0:
                xpos=xpos+time*xvel
                XX.append(round(xpos,3))
                yvel1=yvel-time*9.8
                ypos=ypos+time*(yvel+yvel1)/2
                YY.append(round(ypos,3))
                yvel=yvel1
            
            plt.plot(XX,YY)
            print("\nistance traveled meters.",format(xpos))      
    plt.show()
main()

物理受到的阻力 $f=-k*v^n$