differential-equations

23熱度

2回答

我有一組定義3D輪廓的3D點。我想要做的是獲得對應於這個輪廓的最小表面表示(see Minimal Surfaces in Wikipedia)。基本上這需要求解一個非線性偏微分方程。在Matlab中，使用pdenonlin功能(see Matlab's documentation)幾乎是直接的。可以在這裏找到其用於解決最小表面問題的示例：Minimal Surface Problem on t

-1熱度

2回答

ODE函數的連續繪圖

我有一個ODE求解器，工程很好，很流暢，但我需要在一個圖中繪製所有圖。連接圖（1）+（3）和圖（2）+（4），我必須設置啓動和停止條件，但它不適用於我，我是在死衚衕。我試圖用x_m設置結束條件而沒有結果。 options = odeset('Events',@events); [t,y] = ode45(@ph1,[0,w_max],[0,0], options); figure(1),pl

0熱度

1回答

的Scilab Xcos，插入在Xcos圖中的差分方程

我怎樣才能插入Xcos圖中的差分方程，如： Y（K + 1）= Y（k）的[一個 SQRT（Y（K -1））] + b * y（k-1）; ？感謝，問候編輯查找互聯網上我認爲要解決我的問題，最好的辦法就是使用： scifunck block。

0熱度

1回答

解決微分方程絕對值

我想解決這種形式的方程： X」 = -A.x + B. | SIN（100 * PI * T）| 我用ODE45這樣的： function find_x t = 0:0.001:10; x0 = 0; R1 = 90000; R2 = 1000; C = 0.001; [t,x]=ode45(@rhs, t , x0); plot(t,x); function dxdt

-1熱度

1回答

整合軌道軌跡2

原來的二階微分方程是 x'' - 2 * omega * y' - omega ** 2 * x = - mue * (x + pi2 * r12)/np.sqrt((x + pi2 * r12) ** 2 + y ** 2) ** 3 - mum * (x - pi1 * r12)/np.sqrt((x - pi1 * r12) ** 2 + y ** 2) y'' + 2 * omega *

0熱度

2回答

如何數值測試振盪曲線的穩定性？

我現在有一個振盪曲線，它是一組非線性常微分方程解的一部分。隨着時間的流逝，我需要測試這條曲線的穩定性/收斂性。如何用Matlab做到這一點？這個數字看起來是這樣的：

4熱度

1回答

如何在Python中使用dorpi5或dop853

我已經瀏覽了scipy.integrate.ode，但我無法找到如何實際使用這些集成方法，dorpi5和dop853。我想嘗試整合ode integration python versus mathematica我的Python代碼與這兩種方法，看看它是如何影響結果，但不知道如何。

1熱度

1回答

在Matlab中求解微分方程，ode45

我想用Matlab中的函數ode45來求解一個帶有三個微分方程的系統。我真的不明白我得到的錯誤，我可以用一些幫助來理解我做錯了什麼。微分方程如下： F1 = -k1y1+k2(y2-y1) F2 = -k2(y2-y1)+k3(y3-y2) F3 = -k3(y3-y2) 而且我在Matlab代碼是這樣的： function dz = kopplad(t, z) global m1 m

-5熱度

1回答

計算ODE與MATLAB

微分方程： α '（T）= S（β-βα+α-Qα^ 2） β'（T）=（S^-1 ）（ - β-αβ+γ） γ'（T）= W（α-γ）初使值 α（0）= 30.00 β（0）= 1.000 γ（0）= 30.00 計算我要解決從T_0 = 0的問題至t = 10，而使用數值• = 1, q = 1和w = 0.1610 我不知道如何編寫ODE的函數，並非常感謝幫助！

2熱度

2回答

使用ODE45查找函數的最大值

我試圖在MATLAB中的微分方程系統中找到其中一個方程的位置。我試圖使用odeset的事件屬性。我如何挑選我的函數中的特定方程？ options = odeset('Events',@event); [t x tm xm ie] = ode45(@Lorenz,[0 200],I,options); function X = Lorenz(t,x) r = 15; sigma = 10