如何在 scipy 集成模塊“solve_ivp”中正確設置“rtol”和“atol”以求解具有未知解析解的 ODE 系統？

[英]How to correctly set the 'rtol' and 'atol' in scipy integration module 'solve_ivp' for solving a system of ODE with unknown analytic solution?

我試圖使用 solve_ivp 在 Python 中重現 ode45 求解器的一些結果。 盡管所有參數、初始條件、步長以及“atol”和“rtol”（即 1e-6 和 1e-3）都相同，但我得到了不同的解決方案。 這兩種解決方案都收斂於周期性解決方案，但類型不同。 由於 solve_ivp 使用與 ...

如何在 Matlab 中將變量的矩陣和列向量相乘？

[英]How do I multiply a matrix and a column vector of variables in Matlab?

我有一個矩陣“F”（3 * 20 雙）和一個單元格類型“ans”（20*1），由 (1,x,y,z,xy,yx,...) 形式的變量組成。 我想將 F 和 ans 相乘以獲得方程組並使用它來獲得函數 稍后我將使用它來解決 ODE。 這給了我錯誤 現在，我無法將 celltype 轉換為 do ...

如何在這個 matlab 中用 runge-kutta 替換 ode45 方法？

[英]how to replace the ode45 method with the runge-kutta in this matlab?

我嘗試了一切，到處尋找，但找不到解決我問題的方法。 這是我嘗試用 runge kutta 方法替換 ode45 的代碼，但它給了我錯誤。 我的 runge kutta function： ...

求解一階方程的 7 個方程 ode45 Matlab

[英]Solving 7 equations of first order equations ode45 Matlab

我正在嘗試求解具有七個微分方程的系統。 而且我很難理解 ode45 求解器。 這些是方程式： ω2_dot = -0.75 ω1 ω3 ω1_dot = 0.75 ω2 ω3 + 0.2 ω3_dot = 0 q1_dot = 1/2(ω1q4 + ω2q3 - ω3q2) q2_dot ...

MATLAB ODE45：為什么解析解和數值解之間存在很大差異？

[英]MATLAB ODE45: Why is there a big divergence between analytical and numerical solution?

目前我試圖理解 ode45。 所以我想解決一個練習。 微分方程為：y=y'，y(0)=10。 我寫了這段代碼： 我知道解析解是指數函數。 所以我把它插入plot來驗證解決方案。 （exp function 需要向上移動 9。） 但是有一個分歧我無法解釋。 我理解錯了什么？ 紅色：解析解，藍色：數 ...

將矩陣作為輸入傳遞給 scipy.integrate.solve_ivp python

[英]Passing matrices as input in scipy.integrate.solve_ivp python

我正在解決下面給出的 2DOF 彈簧質量阻尼器系統： 這些是 2 個控制方程 我已經通過以下方式解決了它： 現在，這是通過將 2 個控制方程轉換為 4 個方程來完成的，如下所示： 這個問題是我必須分別寫每個方程（如 function F） Matlab 有一種使用 Ode45 functio ...

解決方案結構 sol=ode45 與 [t,y]=ode45 相比具有不同的時間步數

[英]Solution Structure sol=ode45 has different number of time steps compared to [t,y]=ode45

作為 Matlab 的初學者，我正在使用 test.m 探索 ode45 function，其中包含要求解的 ode 方程。 在編輯器中，我通過 2 個選項將其分別稱為 function。 在一種情況下，我有定義的選項 1 在選項 2 中，我要求一個結構 output 但是，與選項 1 相比，結構選 ...

從 Function 調用術語時遇到問題

[英]Trouble Calling Terms from Function

我有一個定義 function 的代碼，然后我嘗試在另一個表達式中使用我在 function 中定義的變量。 當我這樣做時，我收到一個錯誤，即： 我不確定為什么 phi 是未定義的，因為我在 if/else 語句中有它。 用我的（縮短的）代碼解釋可能會更好： 我計算 phi1、phi2 和 if/e ...

不確定如何使用“ODE45”命令修復錯誤

[英]Not Sure How to Fix Errors Using 'ODE45' Command

我無法弄清楚如何修復我的腳本，特別是使用 ODE45 命令。 這是我到目前為止所擁有的： 我最近雖然問題出在我的變量上，這就是為什么我將它們都定義為全局的。 當我嘗試在命令 window 中運行以下命令時： 我收到這些錯誤： 我曾嘗試在 Mathworks 和其他網站上對其進行研究，但無法弄清楚問題 ...

SciPy RK45 (solve_ivp) 是否准確計算 function 評估的數量？

[英]Does SciPy RK45 (solve_ivp) count the number of function evaluations accurately?

我想針對 SciPy RK45測量我自己的 ODE 積分器的性能。 因此，我需要確切知道 RK45 使用的右側 function 評估的數量。 有誰知道 sol.nfev 的數字是否准確，即沒有重復？ 例如，如果 RK45 拒絕一個步長並重復該步，那么該步的初始評估 f(t,x) 是否會計算多次？ ...

對於不同的初始條件值，使用 matlab 中的 ODE45 function 求解非線性方程

[英]Solve a nonlinear equation using ODE45 function in matlab for different values of initial conditions

我已經編寫了一個腳本來計算和求解一個簡單的倒立擺系統。現在假設我想用不同初始條件值的 ODE45 function 求解系統的非線性動力學方程。如何使用 for 循環求解 state不同初始條件值的 X 向量？我寫了一個 for 循環來做到這一點，但我無法得到我想要的答案。請幫助我。這是我的 fun ...

將 ode45 代碼從 MATLAB 轉換為 python 代碼

[英]Convert ode45 code from MATLAB to a python code

如何使用 python 解決此 MATLAB ode 問題這是帶有 BC 的 IVP：F'''+FF''-F'^2+1=0 F(0)=F'(0)=0, F'(oo)=1 當前的 matlab 代碼將生成以下 plot 它與教科書的解決方案相同： 問題是我需要使用 python 重新編碼相同的問題 ...

使用 ode45 在 MATLAB 中求解具有離散值的微分方程

[英]Solving differential equations with discrete values in MATLAB using ode45

我有一個微分方程- 使用 ode45，我得到了 L(x) 的解。 我有一個由 L_val 表示的 L(x) 值數組。 使用這個解決方案，我打算解決另一個微分方程。 我怎樣才能解出 w(x)？ 特別是因為 L(x) 不是 x 的函數，而是一個離散值數組。 ...

MATLAB ode45 OutputFcn 監控 if 循環的變化值

[英]MATLAB ode45 OutputFcn to monitor changing value of if loop

我嘗試保存/查看在 ode45 微分方程求解過程中由 if 循環更改的變量 m。 if 循環在所有其他方程之前在方程文件中： 我已經查看了 ode45 的 OutputFcn Option 的 matlab 描述，還閱讀了https://de.mathworks.com/help/deep ...

如何使用 ODE45 求解具有困難結束條件的微分方程組？

[英]How to use ODE45 to solve for a system of differential equations with a difficult end condition?

我的方程組可以寫成- y1' = F_1(x,y1,y2) y2' = F_2(x,y1,y2) 其中 F1 和 F2 是 x、y1、y2 的一些函數。 y1 和 y2 是 x 和 y1'=dy1_dx 和 y2'=dy2_dx 的函數。 初始條件為 x=0，y1=y2=0 ...

使用帶有 Octave 的 Matlab 工具箱“PDToolbox”：錯誤消息

[英]Using Matlab toolbox “PDToolbox” with Octave: Error messages

我想在八度音階下使用 Matlab 中的“PDToolbox”。 我正在使用https://github.com/carlobar/PDToolbox_matlab中的版本，但我也嘗試過發布 v1。 我在 ubuntu 18.04 下使用倍頻程 4.2.2。 安裝工具箱后，我首先嘗試在測試目錄中運行 ...

ode45, 參數錯誤, 生長動力學

[英]ode45, error in odeargument, growth kinetics

我正在嘗試解決這個與微藻生長動力學相關的 ODE 系統，我附上了我的 function 和腳本。 我收到一些與 arguments 相關的錯誤。 請讓我知道我的錯誤在哪里謝謝:) 該系統由 7 個 ODE 和 7 個對應的 state 變量和 25 個參數組成 Functionfunction dx ...

odefun 的工作原理（matlab）

[英]How odefun works (matlab)

我試圖弄清楚如何解決 matlab 中的微分方程 [t,y] = ode45(odefun,tspan,y0)。 我有一個odefun的例子，但我不確定這段代碼的哪一部分是做什么的，或者我只是不懂語法。 eq(1,1)=, eq(2,1)=, eq(3,1) = - 是向量嗎？ 是第一次辯論時間 ...

如何使用 ODE45 的輸入 u 進行模擬 - MATLAB

[英]How can I simulate with a input u with ODE45 - MATLAB

我有兩個 function 手柄： 如何使用控制變量 u 在 MATLAB 中模擬它們？ 控制變量 u 隨時間變化。 ...

odeint：無法根據規則“安全”將數組數據從 dtype('complex128') 轉換為 dtype('float64')

[英]odeint: Cannot cast array data from dtype('complex128') to dtype('float64') according to the rule 'safe'

以下代碼給出了錯誤：Cannot cast array data from dtype('complex128') to dtype('float64') based on the rule 'safe' 但是，如果我從func返回真實值（而不是復雜值），例如： 然后odeint工作沒有任何錯誤。 ...