最新国产好看的视频,伊人天堂AV在线,国产Aaaaaa视频,蜜臀视频在线观看一区,人妻av色图,密臀久久久精品影片,青青视频免费观看毛片,久草在线观看视,国产三级精品色情在线

python中sympy庫求常微分方程的用法

 更新時間:2020年04月28日 11:46:15   作者:t4ngw  
這篇文章主要介紹了python中sympy庫求常微分方程的用法,文中通過示例代碼介紹的非常詳細(xì),對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值,需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧

問題1:

程序,如下

from sympy import *
f = symbols('f', cls=Function)
x = symbols('x')
eq = Eq(f(x).diff(x, x) - 2*f(x).diff(x) + f(x), sin(x))
print(dsolve(eq, f(x)))

結(jié)果

Eq(f(x), (C1 + C2*x)*exp(x) + cos(x)/2)

附:布置考試中兩題

1.利用python的Sympy庫求解微分方程的解 y=f(x),并嘗試?yán)胢atplotlib繪制函數(shù)圖像

程序,如下

from sympy import *
f = symbols('f', cls=Function)
x = symbols('x')
eq = Eq(f(x).diff(x,1)+f(x)+f(x)**2, 0)
print(dsolve(eq, f(x)))
C1 = symbols('C1')
eqr = -C1/(C1 - exp(x))
eqr1 = eqr.subs(x, 0)
print(solveset(eqr1 - 1, C1))
eqr2 = eqr.subs(C1, 1/2)
# 畫圖
import matplotlib.pyplot as plt
import numpy as np
x_1 = np.arange(-5, 5, 0.1)
y_1 = [-0.5/(0.5 - exp(x)) for x in x_1]
plt.plot(x_1, y_1)
plt.axis([-6,6,-10,10])
plt.grid()
plt.show()

結(jié)果

Eq(f(x), -C1/(C1 - exp(x)))
FiniteSet(1/2)


2.利用python的Sympy庫求解微分方程的解 y=y(x),并嘗試?yán)胢atplotlib繪制函數(shù)圖像

程序,如下

from sympy import *
y = symbols('y', cls=Function)
x = symbols('x')
eq = Eq(y(x).diff(x,1), y(x))
print(dsolve(eq, y(x)))
C1 = symbols('C1')
eqr = C1*exp(x)
eqr1 = eqr.subs(x, 0)
print(solveset(eqr1 - 1, C1))
eqr2 = eqr.subs(C1, 1)
# 畫圖
import matplotlib.pyplot as plt
import numpy as np
x_1 = np.arange(-5, 5, 0.01)
y_1 = [exp(x) for x in x_1]
plt.plot(x_1, y_1, color='orange')
plt.grid()
plt.show()

結(jié)果

Eq(y(x), C1*exp(x))
FiniteSet(1)

到此這篇關(guān)于python中sympy庫求常微分方程的用法的文章就介紹到這了,更多相關(guān)python sympy常微分方程內(nèi)容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家!

相關(guān)文章

最新評論

宜兰市| 剑阁县| 德令哈市| 阳高县| 阿荣旗| 霍州市| 安塞县| 郁南县| 安福县| 绥中县| 开远市| 天长市| 星座| 梁河县| 林口县| 宝鸡市| 舟山市| 固原市| 景德镇市| 长沙县| 赤水市| 东光县| 普兰县| 泽库县| 溆浦县| 赤水市| 浦县| 呼图壁县| 乐安县| 车致| 米泉市| 马鞍山市| 建平县| 闽侯县| 阜南县| 广饶县| 和田市| 合阳县| 易门县| 东城区| 桃园市|