ode

0熱度

1回答

我想通過在Matlab中編寫代碼來求解微分方程組。我在這個論壇發帖，希望有人能夠以某種方式幫助我。我有一個10個耦合微分方程組。它是一種載體宿主的流行病模型，它捕捉人類與昆蟲種羣之間疾病的傳播。由於它是一個簡單的微分方程組，因此我使用求解器（ode45）來處理非僵硬問題類型。有10個微分方程，每個代表10個不同的狀態變量。有兩個函數具有10個耦合ODE的相同系統。一個叫NoEffects_de

3熱度

1回答

Matlab：CVODE接口，需要正常的Matlab ODE參數

CVODE是流行的SUNDIALS package中的兩個ODE求解器之一。有一個Matlab接口提供，但用法與標準的Matlab ode*接口不一樣。必須先初始化解算器，然後遍歷所需的每個時間點，最後釋放內存。顯然，它不是像ode15s這樣的替代品。有沒有人寫過CVODE的包裝器，它需要Matlab ODE求解器的常用參數？

0熱度

1回答

常規微分方程式表示中的誤差

UPDATE 我在試圖找到在鏈接LE中給出的Lyapunov指數。我試圖通過採取以下方案來理解我的情況。這些都是一組常微分方程（這些只是用於測試如何與COS和sin作爲ODE工作） f(1)=ALPHA*(y-x); f(2)=x*(R-z)-y; f(3) = 10*cos(x); 和x=X(1); y=X(2); cos(y)=X(3); F1裝置dx/dt; F2 dy/dt和f3在

0熱度

1回答

integrate.odeint給出了兩個非常不同的答案時，它不應該

我有以下的Python MWE（代碼如下解釋） #!/usr/bin/python from scipy import integrate from math import * import numpy import matplotlib.pyplot as plt def base_equations(y,t,center): return [5*exp(-(t-cente

0熱度

1回答

matlab ode15s在GA內調用時不起作用

我有一個襯套的數學（參數）模型。模型的輸入是位移（正弦波），以力作爲輸出（以各自的時間步長）。該模型的參數將被確定使用遺傳算法和力方程要解決使用ode15s解算器。求解器給出假設當一組參數（不使用GA）但是從GA調用時我得到以下錯誤n節數倍的輸出：警告：矩陣是奇異的，靠近以單數或嚴重縮放。結果可能不準確。 RCOND = NaN。最後警告：如果在t = 1.570796e + 000。

0熱度

1回答

更改ODE求解器中使用的函數內的參數

是否可以使用ODE求解器（例如ode45），並仍然能夠在被調用函數內「更改」參數值？例如，如果我是使用下面的功能： function y = thisode(t, Ic) % example derivative function % parameters a = .05; b = .005; c = .0005; d = .00005; % state variables

18熱度

5回答

使用自適應步長與scipy.integrate.ode

scipy.integrate.ode的（簡短）文檔說，有兩種方法（dopri5和dop853）具有步進控制和密集輸出。查看示例和代碼本身，我只能看到從集成商獲得輸出的非常簡單的方法。也就是說，看起來你只是通過一些固定的dt步進積分器，在那個時候獲得函數值，然後重複。我的問題具有相當可變的時間尺度，所以我想只需在任何時間獲取值就可以獲得所需的公差值。也就是說，早期事情變化緩慢，所以輸出時間步驟可

0熱度

1回答

通過設置tspan = [to：very_small_step：tf]，它是否會影響ode45解算器的步長？

我知道ode45求解器具有由Matlab程序本身控制的自適應步長。下面的描述是通過Matlab website給出：指定TSPAN有兩個以上的元件並不影響>的解算器使用來遍歷從TSPAN的間隔中的內部時刻的步驟（1）至TSPAN（結束）。 ODE>套件中的所有解算器都通過基本公式的連續擴展來獲取輸出值。雖然>求解器不一定準確地到達tspan中指定的時間點，但是在指定時間點產生的解>與在內部時間點

6熱度

1回答

如何使用內部閾值來解決ODE？

我已經包含了一些賦有以下功能： myfunction <- function(t, state, parameters) { with(as.list(c(state, parameters)),{ if (X>20) { # this is an internal threshold! Y <- 35000 dY <- 0 }el

0熱度

1回答

Matlab - > scipy ode（複雜）函數翻譯

我正在學習python，numpy和scipy。我不知道是否有可能在MATLAB到Python翻譯這樣的功能： function [tT, u ] = SSolve5TH(n, t, t0,tf,u_env,utop_init, utop_final,ubottom,te_idx) options = []; [tT,u] = ode23s(@SS,t,u_env,options,@B);