differential-equations

0熱度

1回答

我的問題是：爲什麼rhsexpr有效，但rhsfunc沒有。謝謝。 :) m = l = 1.; \[Mu] = 100; beltv[t_] = .1; spring[x_] = 1000. (l - x); viscous[v_] := -30. (v - beltv[t]); coulomb[v_] := -25. Sign[v - beltv[t]

0熱度

2回答

擬合微分方程：如何適應一組數據到差分方程中的R

與數據集： conc <- data.frame(time = c(0.16, 0.5, 1.0, 1.5, 2.0, 2.5, 3), concentration = c(170, 122, 74, 45, 28, 17, 10)) 和我想此數據代入下面的微分等式： dC/dt= -kC 其中C將是數據集中的濃度和時間t。這也會給出k的結果。有誰能給我一個線索如何在R中做到這一點？謝謝。

4熱度

1回答

如何在matlab中繪製斜場

我正在尋找一種在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，並試圖找出如何用楓樹做到這一點，然後製作一個圖形顯示該功能的每一步。有什麼建議麼。我試圖在楓樹上尋找，但它一直指向我的功能新版本的楓樹軟件，我沒有。我發佈了這個問題作爲

1熱度

2回答

實現龍格庫塔系統的DE

我試圖在MATLAB中實現Runge-Kutta Method for Systems of DEs。我沒有得到correct answers，我不確定代碼中是否有錯誤或者我用來運行它的命令。這裏是我的代碼： function RKSystems(a, b, m, N, alpha, f) h = (b - a)/N; t = a; w = zeros(1, m);

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

0熱度

1回答

線性雙方程系統，兩個變量二階導數都在

原諒我，如果這被認爲是重新發布，但我已被建議我可能已經給出了一個不好的格式。我試圖解決兩個線性二階微分方程。我想將它們分解爲單階方程，但我無法看到兩個變量在兩個問題中都有二階導數。 (m*a)u」 + (I + m*a^2)θ」 + (d*a^2)θ’ + (K - m*g*a)θ = 0 (M + m)u」 + (m*a)θ」 = -F 我已經使用了一些既dsolve和ODE45 這是

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