• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    公众号

数学建模pythonmatlab编程(喷泉模拟)

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

在无风情况下的喷泉模拟

我的python代码

import numpy as np
import random
import matplotlib
matplotlib.rcParams['font.sans-serif']=[u'simHei']
matplotlib.rcParams['axes.unicode_minus']=False
import matplotlib.pyplot as plt

tt = np.linspace(0,4,1000)
#y = -tt**2+4*tt
#plt.plot(tt,y,linestyle='--')


# 没有风
vy0=10
g=10
vx0=5
x=[]
y=[]
x0=0
y0=0
dt = 4/1000
x.append(x0)
y.append(y0)
for t in tt:
    y0 += vy0*dt
    x0 += vx0*dt
    vy0 -= g*dt
    x.append(x0)
    y.append(y0)
    if y0<0:
        break
plt.plot(x,y,linestyle='--')

 

 

在右边风情况下的喷泉模拟

我的python代码

# 有向右的风
vy0=10
g=10
vx0=5
a = 1.5
x=[]
y=[]
x0=0
y0=0
dt = 4/1000
x.append(x0)
y.append(y0)
for t in tt:
    y0 += vy0*dt
    x0 += vx0*dt
    vy0 -= g*dt
    vx0 -= a*dt
    x.append(x0)
    y.append(y0)
    if y0<0:
        break
plt.plot(x,y)

 

 

在左边风情况下的喷泉模拟

我的python代码

# 有向左的风
vy0=10
g=10
vx0=5
a = -1.5
x=[]
y=[]
x0=0
y0=0
dt = 4/1000
x.append(x0)
y.append(y0)
for t in tt:
    y0 += vy0*dt
    x0 += vx0*dt
    vy0 -= g*dt
    vx0 -= a*dt
    x.append(x0)
    y.append(y0)
    if y0<0:
        break
plt.plot(x,y)

 

 

所有代码

# -*- coding: utf-8 -*-

import numpy as np
import random
import matplotlib
matplotlib.rcParams['font.sans-serif']=[u'simHei']
matplotlib.rcParams['axes.unicode_minus']=False
import matplotlib.pyplot as plt

tt = np.linspace(0,4,1000)
#y = -tt**2+4*tt
#plt.plot(tt,y,linestyle='--')


# 没有风
vy0=10
g=10
vx0=5
x=[]
y=[]
x0=0
y0=0
dt = 4/1000
x.append(x0)
y.append(y0)
for t in tt:
    y0 += vy0*dt
    x0 += vx0*dt
    vy0 -= g*dt
    x.append(x0)
    y.append(y0)
    if y0<0:
        break
plt.plot(x,y,linestyle='--')

# 有向右的风
vy0=10
g=10
vx0=5
a = 1.5
x=[]
y=[]
x0=0
y0=0
dt = 4/1000
x.append(x0)
y.append(y0)
for t in tt:
    y0 += vy0*dt
    x0 += vx0*dt
    vy0 -= g*dt
    vx0 -= a*dt
    x.append(x0)
    y.append(y0)
    if y0<0:
        break
plt.plot(x,y)

# 有向左的风
vy0=10
g=10
vx0=5
a = -1.5
x=[]
y=[]
x0=0
y0=0
dt = 4/1000
x.append(x0)
y.append(y0)
for t in tt:
    y0 += vy0*dt
    x0 += vx0*dt
    vy0 -= g*dt
    vx0 -= a*dt
    x.append(x0)
    y.append(y0)
    if y0<0:
        break
plt.plot(x,y)
View Code

 

 喷泉散落在地上模拟

 

import numpy as np
import random
import matplotlib
matplotlib.rcParams['font.sans-serif']=[u'simHei']
matplotlib.rcParams['axes.unicode_minus']=False
import matplotlib.pyplot as plt
from scipy.stats import norm
a1_list = random.uniform(0,0.5) # 生成一组正太分布的加速度,均值是0,方差是0.5
a1_list = norm.rvs(a1_list, size=1000)
a2_list = random.uniform(0,0.5)
a2_list = norm.rvs(a2_list, size=1000)
X0=[]
X1=[]
# 对每一组加速度
for i in range(1000):
    vy0=10
    g=10
    vx0=5
    vx1=5
    a1 = a1_list[i]
    a2 = a2_list[i]
    x0=0
    x1=0
    y0=0
    dt = 4/1000
    for t in tt:
        y0 += vy0*dt
        x0 += vx0*dt
        x1 += vx1*dt
        
        vy0 -= g*dt
        vx0 -= a1*dt
        vx1 -= a2*dt
    
        if y0<0:
            break
    X0.append(x0)
    X1.append(x1)
    
plt.scatter(X0,X1)

 

炫酷喷泉三维图

from mpl_toolkits.mplot3d import axes3d
from scipy.stats import norm
a1_list = random.uniform(0,0.2) # 生成一组正太分布的加速度,均值是0,方差是0.5
a1_list = norm.rvs(a1_list, size=100)
a2_list = random.uniform(0,0.2)
a2_list = norm.rvs(a2_list, size=100)
X0=[]
X1=[]
# 对每一组
ax=plt.subplot(111,projection='3d')
for i in range(100):
    vy0=10
    g=10
    vx0=5
    vx1=5
    a1 = a1_list[i]
    a2 = a2_list[i]
    xx0=[]
    xx1=[]
    y=[]
    x0=0
    x1=0
    y0=0
    dt = 4/1000
    for t in tt:
        y0 += vy0*dt
        x0 += vx0*dt
        x1 += vx1*dt
        vy0 -= g*dt
        vx0 -= a1*dt
        vx1 -= a2*dt
        xx0.append(x0)
        xx1.append(x1)
        y.append(y0)
        if y0<0:
            break
#    X0.append(x0)
#    X1.append(x1)
    
    ax.plot(xx0,xx1,y,'--')
    
plt.show()

 


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
matlab字符串操作总结发布时间:2022-07-22
下一篇:
matlab画棋盘格程序发布时间:2022-07-22
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap