python裡怎麼樣求解微分方程

2021-09-10 19:55:05 字數 1425 閱讀 6053

有很多大學生問我,學習python有什麼用呢?我說:你至少可以用來解微分方程,如下面的例子,就是解決微分方程:

y"+a*y'+b*y=0 

**如下:

#y"+a*y'+b*y=0 

from scipy.integrate import odeintfrom pylab import *def

deriv

(y,t):

# 返回值是y和y的導數組成的陣列    a = -2.0    b = -0.1

return array([ y[1], a*y[0]+b*y[1] ])time = linspace(0.0,50.0,1000)yinit = array([0.0005,0.2])     # 初值y = odeint(deriv,yinit,time)figure()plot(time,y[:,0],label='y')    #y[:,0]即返回值的第一列,是y的值。label是為了顯示legend用的。plot(time,y[:,1],label="y'")     #y[:,1]即返回值的第二列,是y』的值xlabel('t')ylabel('y')legend()show()

輸出結果如下:

再分享一下我老師大神的人工智慧教程吧。零基礎!通俗易懂!風趣幽默!還帶黃段子!希望你也加入到我們人工智慧的隊伍中來!

有很多大學生問我,學習python有什麼用呢?我說:你至少可以用來解微分方程,如下面的例子,就是解決微分方程:

y"+a*y'+b*y=0 

**如下:

#y"+a*y'+b*y=0 

from scipy.integrate import odeintfrom pylab import *def

deriv

(y,t):

# 返回值是y和y的導數組成的陣列    a = -2.0    b = -0.1

return array([ y[1], a*y[0]+b*y[1] ])time = linspace(0.0,50.0,1000)yinit = array([0.0005,0.2])     # 初值y = odeint(deriv,yinit,time)figure()plot(time,y[:,0],label='y')    #y[:,0]即返回值的第一列,是y的值。label是為了顯示legend用的。plot(time,y[:,1],label="y'")     #y[:,1]即返回值的第二列,是y』的值xlabel('t')ylabel('y')legend()show()

輸出結果如下:

python裡怎麼樣通過函式名稱來獲取函式位址

如果你想通過函式的名稱來獲取函式的執行位址,可以像下面這樣實現 file builtin import example 2.py def getfunctionbyname module name,function name module import module name return geta...

python效率怎麼樣 如何優化python的效率

優化python的效率的方法 1 優化演算法時間複雜度 2 減少冗餘資料 3 合理使用copy與deepcopy 4 使用dict或set查詢元素 5 合理使用生成器 generator 和yield。優化方法 1 優化演算法時間複雜度 演算法的時間複雜度對程式的執行效率影響最大,在python中可...

學python做資料分析怎麼樣?

用python做資料分析是一件再合適不過的事情,我們舉個簡單的例子。例如收集資料的方法是使用問卷,實驗軟體 例如,psychopy,opensesame 和觀察。當問卷和實驗軟體使用數字應用程式時,我們當然也會以數字檔案格式 例如,excel電子 和逗號分隔,csv,檔案 獲取資料。如果資料集非常小...