python中sympy庫求常微分方程的用法
更新時間：2020年04月28日 11:46:15 作者：t4ngw
這篇文章主要介紹了python中sympy庫求常微分方程的用法，文中通過示例代碼介紹的非常詳細，對大家的學(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庫求解微分方程的解
程序，如下
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庫求解微分方程的解
程序，如下
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)文章
python中k-means和k-means++原理及實現(xiàn)
本文主要介紹了python中k-means和k-means++原理及實現(xiàn)，文中通過示例代碼介紹的非常詳細，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2022-05-05
這篇文章主要給大家介紹了關(guān)于python時間time模塊處理的相關(guān)資料，文中通過示例代碼介紹的非常詳細，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2020-10-10
這篇文章主要介紹了python中顯存回收問題解決方法，文章提供一個實例問題和處理的思路，具有一定的參考價值，需要的小伙伴可以參考一下2021-12-12
今天教各位小伙伴怎么用Python連接oracle,文中附帶非常詳細的圖文示例,對正在學(xué)習(xí)的小伙伴們很有幫助喲,需要的朋友可以參考下2021-05-05
Python使用ChainMap實現(xiàn)組合數(shù)據(jù)魔法實例探究
這篇文章主要為大家介紹了Python使用ChainMap實現(xiàn)組合數(shù)據(jù)魔法實例探究,有需要的朋友可以借鑒參考下,希望能夠有所幫助,祝大家多多進步,早日升職加薪2024-01-01
今天小編就為大家分享一篇python自動化UI工具發(fā)送QQ消息的實例，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2019-08-08
基于python?win32setpixel?api?實現(xiàn)計算機圖形學(xué)相關(guān)操作(推薦)
這篇文章主要介紹了基于python?win32setpixel?api?實現(xiàn)計算機圖形學(xué)相關(guān)操作,這次的主要分為2個主要模塊，一個是實現(xiàn)畫線，畫圓的算法，還有填充的算法，以及裁剪的算法，需要的朋友可以參考下2021-12-12
這篇文章主要介紹了python3.7 的新特性詳解，文中通過示例代碼介紹的非常詳細，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2019-07-07
Python中requests、aiohttp、httpx性能比拼
本文主要介紹了Python中requests、aiohttp、httpx性能比拼，文中通過示例代碼介紹的非常詳細，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2023-06-06
以制作一個網(wǎng)站為例，聊一聊小白如何入門Python，具有一定的參考價值，感興趣的小伙伴們可以參考一下2018-03-03