differential-equations

4熱度

1回答

我正在尋找一種在Matlab中繪製斜場的方法。這裏就是我要找：我有一個公式 dy/dx = f(x,y) 或 dx/dt = f(x,y) dy/dt = g(x,y) ，我希望把它畫在一個不錯的方式因爲關於它的唯一answer沒有回答我的問題，它花了我一些時間如何做到這一點。因爲這不是我一直在MATLAB中做的事情（很可能直到下一次我需要它，我會忘記它）我正在爲我創建一個備忘錄如

1熱度

1回答

如何繪製並求解楓葉12的微分方程？

考慮初始條件爲$y\left(0\right)=1$的微分方程$y^{\prime}=y-2$。一）使用歐拉法有4個步驟尺寸0.2的估計$y\left(0.8\right)$ 我知道如何做手工這一點;然而，我已經安裝了楓樹12，並試圖找出如何用楓樹做到這一點，然後製作一個圖形顯示該功能的每一步。有什麼建議麼。我試圖在楓樹上尋找，但它一直指向我的功能新版本的楓樹軟件，我沒有。我發佈了這個問題作爲

0熱度

1回答

使用龍格庫塔解決貝塞爾函數

我正在爲我的一個類的任務，我應該寫一個代碼使用我選擇的程序（我選擇了Matlab）來解決貝塞爾函數微分方程使用4階Runge-Kutta方法。作爲參考，貝塞爾函數DE是： x^2 *（J_n）''+ x *（J_n）'+（x^2-n^2）* J_n = 0。我能夠通過此分離成兩個耦合的第一順序的DE：（J_n） '= Z_n和（Z_n）' +（1/X）* Z_n + [（X^2-正^ 2）/

0熱度

2回答

如何在Matlab中表示和求解這個微分方程組？

的方程如下： DXI（T）/ dt的= -c * XI（T）*義（噸） DYI（T）/ dt的= A *Σ{i = 1至其中a，b和c是任何你想要的常數值，例如a = 1，其中a，b和c是任意常數值，例如a = 1， b = 2，c = 3。（t）* y1（t））表示從i = 1到n的和，例如n = 3：x1（t）* y1（t）+ x2（x1（t）（t）* y2（t）+ x3（t）* y3（t

4熱度

1回答

在Python中繪製軌道軌跡

如何在python中設置三個主體問題？如何定義函數來解決ODE？三個方程是 x'' = -mu/np.sqrt(x ** 2 + y ** 2 + z ** 2) * x， y'' = -mu/np.sqrt(x ** 2 + y ** 2 + z ** 2) * y，並 z'' = -mu/np.sqrt(x ** 2 + y ** 2 + z ** 2) * z。寫爲6階，我們有 x' =

1熱度

1回答

python執行3d圖時會凍結

從plotting orbital trajectories開始，我們有以下代碼。價值已經改變爲已知的IC工作。如果此代碼是正確的（它不能是雖然），它會產生運行這段代碼只是凍結我的電腦或輸出絕對錯誤的地塊。有人可以幫我找到解決這個問題的方法嗎？ import numpy as np from scipy.integrate import odeint import matplotlib.p

-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