這個用戶定義的函數如何與scipy集成odeint？

我正在研究一個電氣工程項目，我們將RLC電路中的電容上的電壓繪製爲時間的函數。我在以下教程中使用我的代碼：Python-ODE for RLC Circuit。我的代碼功能齊全，但我只想更好地理解用戶定義的函數和SciPy集成函數的工作原理。我已閱讀文檔，但我仍然有一些具體問題。這裏是我的代碼：這個用戶定義的函數如何與scipy集成odeint？

from scipy.integrate import odeint # for integrating ode 
import matplotlib.pyplot as plt # for plotting commands 
import numpy as np #For loadtxt and linspace modules 

R,L,C,tl,numplot = np.loadtxt('UserInput_RLC_Values.txt',unpack=True) 

def RLC(A,t): 
    Vc,m=A 
    V = 2.0 #voltageSource 
    res=[m,(V-Vc-(m*R*C))/(L*C)] 
    return res 

time = np.linspace(0.0,tl,numplot) 
vc,m = odeint(RLC,[0.0,0.0],time).T

這裏是我的問題：

什麼樣的作用是參數值A和T的演奏在功能RLC？爲什麼將A分配給Vc（電容器電壓）和m（v'的虛擬變量）？甚至用過嗎？
當我在odeint模塊裏面調用RLC函數時，爲什麼沒有任何參數？
odeint模塊末尾的.T是什麼意思？我在文檔中找不到任何關於它的信息。

感謝您對此的任何見解。我希望這不是一個懶惰的問題，但我認爲完全理解這將幫助我很多。

來源

2016-11-30 Griggs Hutaff

傳遞到odeint的ODE函數需要格式爲f(state, time)，其中time是自變量，state是狀態向量。相應的（矢量值）微分方程是y'(t)=f(y(t),t)

）對於需要高階方程轉換爲一階系統的標準ODE求解器。這裏使用m作爲動量（？）變量，m = dVc/dt。然後dm/dt = d²Vc/dt²是二階導數。

如果輸入是state = [ Vc, m ]輸出需要是[ dVc/dt, dm/dt ] = [ m, Eqn(t,Vc,m) ]其中d²Vc/dt² = Eqn(t,Vc,dVc/dt)是原始方程的一種形式。

t只是需要出現在參數列表中，即使ODE是時間不變的。
）您正在傳遞函數作爲參考/地址，而不是函數值
），意思是「轉置」。 odeint的結果是狀態向量列表，它們是數字對。轉置它將它轉換成一對可以分配給一對變量的列表。

來源

2016-11-30 16:56:20 LutzL

這個用戶定義的函數如何與scipy集成odeint？

回答

相關問題