python中sympy庫(kù)求常微分方程的用法
更新時(shí)間：2020年04月28日 11:46:15 作者：t4ngw
這篇文章主要介紹了python中sympy庫(kù)求常微分方程的用法，文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來(lái)一起學(xué)習(xí)學(xué)習(xí)吧
問(wèn)題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庫(kù)求解微分方程的解
程序，如下
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庫(kù)求解微分方程的解
程序，如下
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庫(kù)求常微分方程的用法的文章就介紹到這了,更多相關(guān)python sympy常微分方程內(nèi)容請(qǐng)搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家！
相關(guān)文章
這篇文章主要介紹了python 如何讀取列表中字典的value值，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。如有錯(cuò)誤或未考慮完全的地方，望不吝賜教2022-02-02
這篇文章主要介紹了Python類和實(shí)例的屬性機(jī)制原理詳解,文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2020-03-03
這篇文章主要介紹了python操作ie登陸土豆網(wǎng)的方法,涉及Python操作頁(yè)面元素的技巧,具有一定參考借鑒價(jià)值,需要的朋友可以參考下2015-05-05
Python的異常問(wèn)題(Python Exception)
這篇文章主要介紹了Python的異常問(wèn)題(Python Exception),具有很好的參考價(jià)值,希望對(duì)大家有所幫助,如有錯(cuò)誤或未考慮完全的地方,望不吝賜教2023-08-08
這篇文章主要為大家詳細(xì)介紹了Python中圖像通用操作的實(shí)現(xiàn)，例如：圖像旋轉(zhuǎn)、圖像縮放等，文中的示例代碼講解詳細(xì)，需要的可以參考一下2023-07-07
Python利用openpyxl類實(shí)現(xiàn)在Excel中繪制樂(lè)高圖案
在商場(chǎng)看到一個(gè)超級(jí)瑪麗的樂(lè)高圖，感覺(jué)使用excel的顏色填充也能畫出來(lái)。所以本文將借助openpyxl類實(shí)現(xiàn)在Excel中繪制樂(lè)高圖案，需要的可以參考一下2022-12-12
Python實(shí)現(xiàn)連接MySql數(shù)據(jù)庫(kù)及增刪改查操作詳解
這篇文章主要介紹了Python實(shí)現(xiàn)連接MySql數(shù)據(jù)庫(kù)及增刪改查操作,結(jié)合實(shí)例形式詳細(xì)分析了Python安裝PyMySQL模塊及使用PyMySQL模塊進(jìn)行mysql數(shù)據(jù)庫(kù)的連接、增刪改查等相關(guān)操作技巧,需要的朋友可以參考下2019-04-04
Python數(shù)據(jù)可視化：箱線圖多種庫(kù)畫法
這篇文章主要介紹了Python數(shù)據(jù)可視化箱線圖多種庫(kù)畫法,本文給大家介紹的非常詳細(xì)，具有一定的參考借鑒價(jià)值，需要的朋友可以參考下2019-11-11
python使用PyFetion來(lái)發(fā)送短信的例子
這篇文章主要介紹了python使用PyFetion來(lái)發(fā)送短信的例子,需要的朋友可以參考下2014-04-04
python實(shí)現(xiàn)猜數(shù)字游戲(無(wú)重復(fù)數(shù)字)示例分享
這篇文章主要介紹了python實(shí)現(xiàn)猜數(shù)字游戲(無(wú)重復(fù)數(shù)字)示例,需要的朋友可以參考下2014-03-03