当前位置:网站首页>Solutions of ordinary differential equations (2) examples
Solutions of ordinary differential equations (2) examples
2022-07-05 08:44:00 【Royal fox God double blazing】
use SciPyのscipy.integrate.odeint Solve simultaneous differential equations , In time t(tの The scope is 0~ Just a few seconds , example :t=0~2 second ) Reconciliation x(t), y(t), z(t) Make a picture .
As t=0 Initial conditions of ,x(t=0) = -10, y (t=0) = 0, z (t=0) = 35.0. also , About the coefficient a,b,c, You can try a = 40, b = 5, c = 35 And a = 40, b = 10, c = 35 The situation of . And time division Δt Take a smaller value appropriately . in addition ,Δt If it is too small , The amount of calculation will become very large . Adjust according to the errors after the actual implementation .
Here are the solutions python Program :
import numpy as np
import scipy.integrate as sciin
import matplotlib.pyplot as plt
# Compare with the independent variable t The function name of the related derivative function is placed in F in
def f(F, t,params):
x,y,z = F
f_values = [a*(y-x),(c-a)*x-x*z+c*y,x*y-b*z] # Write separately x,y,z The derivative of is equal to the formula on the right of
return f_values
# Amplification coefficient
a = 40
b= 5
c = 35
# Put the above three coefficients into parameters in
parameters = [a,b,c]
# Set up x,y,z The initial value of the
x0 = -10
y0 = 0.0
z0 = 35
# Put the initial value into Y0 in
Y0 = [x0,y0,z0]
# Starting point , The end point , Interval setting
tStart = 0.0
tStop = 2
tInc = 0.01 # interval
# Summarize the above to t in
t = np.arange(tStart, tStop, tInc)
# sciin.odeint Explain ODE
solution = sciin.odeint(f, Y0, t, args=(parameters,))
# Make a picture
plt.figure(figsize=(9.5, 6.5))
plt.plot(t, solution[:, 0], color='black')
plt.plot(t, solution[:, 1], color='green')
plt.plot(t, solution[:, 2], color='red')
plt.xlabel('time, t' , fontsize=14)
plt.ylabel('theta(t)', fontsize=14)
plt.show()
I'm a genius .
边栏推荐
- 319. Bulb switch
- Mathematical modeling: factor analysis
- My university
- [noi simulation] juice tree (tree DP)
- Daily question - input a date and output the day of the year
- Bluebridge cup internet of things competition basic graphic tutorial - clock selection
- Business modeling of software model | stakeholders
- leetcode - 445. Add two numbers II
- Tips 1: Web video playback code
- Basic number theory -- Euler function
猜你喜欢
UE pixel stream, come to a "diet pill"!
leetcode - 445. Add two numbers II
Halcon color recognition_ fuses. hdev:classify fuses by color
Example 003: a complete square is an integer. It is a complete square after adding 100, and it is a complete square after adding 168. What is the number?
319. 灯泡开关
Example 004: for the day of the day, enter a day of a month of a year to judge the day of the year?
猜谜语啦(4)
Guess riddles (6)
Example 002: the bonus paid by the "individual income tax calculation" enterprise is based on the profit commission. When the profit (I) is less than or equal to 100000 yuan, the bonus can be increase
RT-Thread内核快速入门,内核实现与应用开发学习随笔记
随机推荐
Esphone Feixun DC1 soft change access homeassstant
Lori remote control commissioning record
One dimensional vector transpose point multiplication np dot
Daily question - input a date and output the day of the year
猜谜语啦(10)
Business modeling of software model | overview
696. Count binary substring
Example 007: copy data from one list to another list.
Array,Date,String 对象方法
287. 寻找重复数-快慢指针
猜谜语啦(142)
Lori remote control LEGO motor
One question per day - replace spaces
319. Bulb switch
Classification of plastic surgery: short in long long long
猜谜语啦(6)
多元线性回归(梯度下降法)
猜谜语啦(3)
[noi simulation] juice tree (tree DP)
Is the security account given by Yixue school safe? Where can I open an account