differential-equations

0熱度

1回答

我想正確地在Matlab中輸入微分方程，所以我可以ode45解決它們，但我不能讓我的代碼工作。這是我進入代碼： function ydot=Untitledrt(z,y) ydot = zeros(6,1); %y(1)=A %y(2)=B %y(3)=C %y(4)=D %y(5)=P %y(6)=T m1 = 6; m2 = 9; m3 = 5; k1 = 6; k2

1熱度

1回答

隨機熱方程 - Fortran

我有一段代碼，它解決了一維隨機熱方程的週期性邊界條件。隨機項是高斯白噪聲。我的問題是，我正確實施噪音？高斯噪聲被定義爲具有均值爲零，並且第二時刻告訴我們在任何時間對的值是相同分佈的且在統計上是獨立的。 !I first define arrays for the uniform and gaussian random numbers real, dimension (-1:n) :: u,

2熱度

2回答

C++程序已停止工作 - 求解常微分方程

我正在編寫一個C++程序來爲一個大學任務找到一階微分方程的解。該程序啓動，然後一旦我輸入迭代次數，我會得到錯誤消息「歐拉的method.exe已停止工作」。這是我的代碼： #include <functional> #include <vector> using namespace std; double f_r(double x, double r)

3熱度

1回答

比較Boost.Odeint vs Scipy.integrate.odeint？

我最近關注boost.odeint庫，我對可能性和可配置性的數量感到驚訝。然而，廣泛使用scipy.integrate.odeint（這實質上是fortran中的ODEPACK的包裝），我不知道他們的表現如何比較。我知道boost.odeint也帶有並行化，這在scipy中是不可能的（據我所知），這會提高性能，但我要求單核心案例。但是，因爲我必須在這種情況下將boost.odeint（使用cyt

0熱度

2回答

從MATLAB ODE求解如何比較結果

我解微分方程系統用Matlab這樣 tspan = [0 10]; res = ode15s(@(t,x) func,tspan,x0); ，我有另外一個系統，它比@func指定的slighlty不同，解決像 res2 = ode15s(@(t,x) func2,tspan,x0); 我想通過計算每個時間步的差異來比較結果。但由於Matlab解算器是可變時間步長，因此res.y（解）矩

2熱度

1回答

Matlab：微分方程（ode45）：我能逆轉tspan以獲得更好的初始條件嗎？

我使用ode45在Matlab中求解/繪製一個二階微分方程。我的tspan是從0到0.25。但是接近於零的初始條件不明確（斜率變爲無窮大，複數值）。 0.25附近的條件明確定義（斜率和值都爲零）。問題：我可以扭轉tspan，並使用了「最終狀態」的初始條件？我知道我可以做到這一點（看下面的代碼），我得到一個看起來像我期望的情節，但這是一個有效的一般做法嗎？我在這種情況下幸運嗎？ ode45提

1熱度

1回答

解決wolfram mathematica中的微分方程組的非線性系統

如何求解非線性微分方程組並得到該解的繪圖？該系統沒有初始條件。例如， X '=（X + Y）^ 2 - 1 Y'= -y^2 - X + 1

1熱度

1回答

R包deSolve：錯誤checkFunc

我對R很新，我試圖用de deSolve包來解決一個微分方程組。我總共有10000次的時間步數，每100次的步數我必須改變變量的值。我試過代碼給我下面的錯誤： Error in checkFunc(Func2, times, y, rho) : The number of derivatives returned by func() (2) must equal the length

0熱度

1回答

Python多個初始條件

我對編程相當陌生。我在課堂上選擇了一些R和Matlab，但我試圖提高我對Python的熟悉程度。我試圖編寫一個模型，該數學模型將給出多個初始條件下的微分方程組的數值近似。我目前擁有的代碼如下： import matplotlib.pyplot as plt # parameters sigma=10 gamma=1 alpha=50 delta=2.1 # initial condi

-1熱度

1回答

在R中使用deSolve的差分Eq

由於之前不清楚，我對此表示歉意。我現在更瞭解這個功能，但可以在幾個方面使用一些幫助。我想找回轉換（X）與音量（V）的關係，或者反過來也可以。在我看來，傳統的「時代」術語是我想用0-1的X序列取代的，X是轉換，記得如此界定0和1.0 下面，rw是反應速率，並且是任何給定時刻的分壓的函數，它們被描述爲Pw，Px，Py和Pz，它們本身是初始條件（P.w0，v.0）的函數，轉換也是X. 預先感謝您 rm