differential-equations

1熱度

3回答

我用熱傳導建模物理系統，並做數值計算，我需要用三對角矩陣求解線性方程組。我正在使用此算法得到結果：http://en.wikipedia.org/wiki/Tridiagonal_matrix_algorithm但我擔心我的方法很簡單，並且不是最優的。應該用什麼C++庫以最快的方式解決該系統？我還應該提到矩陣不會經常改變（只有等式的右邊部分被改變）。謝謝！

2熱度

1回答

解決ODE的

在Matlab中我知道如何直至甚至出現解決一個ODE（例如，解決ODE直到y座標爲零）： http://www.mathworks.com/help/techdoc/ref/odeset.html#f92-1017470 是否有一個比喻在八度？下面是相關的代碼： lsode_options('Events', @events); t0 = linspace (0, 20, 1000); [t

3熱度

1回答

什麼「邊界條件」可以使矩形看起來像一個圓圈？

我正在用MATHEMATICA求解具有周期或自由邊界條件的方形域上的時間和空間上的四階非線性偏微分方程(t, x)。 WITHOUT使用保角變換，什麼邊界條件在邊緣或拐角我可以使用使正方形域「似乎」等用於我的非線性偏微分方程的圓形結構域，其是笛卡爾？我會不喜歡用的選項有：保角變換改變我的公式極/圓柱座標？這是我追求純粹出於利益，以防萬一有人尖叫血腥謀殺，如果被誤解爲功課問題！ :P

5熱度

1回答

如何訪問modelica中的時間變量

我想在Modelica中對顯式時間相關的微分方程建模。說我想在Modelica語言來 Real x(start=1); der(x) = t; 建模。訪問時間變量t的正確方法是什麼？

2熱度

1回答

python求解複數變量的微分方程

夥計，是否有可能在python中用複數變量解決ODE問題？等式I具有以下形式 dx/dt = -a x -i y(t) 其中y（t）是已知函數，a是已知數，並且i是-1的根。我試圖使用odeint（），但它提供了許多錯誤消息。我猜odeint（）不適用於複雜的變量。因此，一種方法是將x的實部和虛部分開，並將原始ODE視爲兩個耦合的ODE。但我也想知道是否有一些更方便的方法來完成這項任務

5熱度

3回答

解決一個延遲微分方程（DDE）系統約束非負的解決方案

在MATLAB中，ode45有一個參數NonNegative約束的解決方案是非負。 They even wrote a paper about how this method works以及它是如何愚蠢的，只要它變成負數就將y_i設置爲0，因爲這通常不起作用。現在，MATLAB也有dde23用於求解延遲微分方程，但這個積分器沒有等效的參數NonNegative。不幸的是，我的任務是爲現有ODE

0熱度

1回答

非線性微分方程的求解

我一般不使用Mathematica，我需要它與其他程序進行比較。我想解決三個微分和非線性方程組的系統。爲此我使用Dsolve。當我放置非線性項（指數）時，一切都會出錯。這裏是我的代碼： equa = {x'[t] == z[t] - Exp[y[t]], y'[t] == z[t] - y[t], z'[t] == x[t] + y[t] - z[t], x[0] == 0, y[0]

-1熱度

1回答

在matlab中的差異

我需要找到一個對象的加速度，公式給出的文本是a = d^2（L）/ d（T）^ 2，其中L =長度和T =時間我計算了這個通過使用這個公式MATLAB a = (1/(T3-T1))*(((L3-L2)/(T3-T2))-((L2-L1)/(T2-T1))) 或 a = (v2-v1)/(T2-T1) 但IM沒有得到正確的答案，任何機構可以告訴我如何用MATLAB任何其他方法來找到（一）

9熱度

2回答

如何使用mathematica繪製斜坡場？

我試圖用mathematica繪製一些微分方程的斜場，但無法弄清楚。說我有方程 y' = y(t) y(t) = C * E^t 我該如何繪製斜坡場？我發現了一個例子，但方式複雜，我瞭解 http://demonstrations.wolfram.com/SlopeFields/

0熱度

1回答

調整幅度調製/波形整形

我有一些可以進行幅度調製並繪製它的工作代碼。但是我試圖改變調製的方式看起來像（y變量）所以它看起來像一個雞蛋形狀。我發現了一個公式/網站看起來不錯 http://www16.ocn.ne.jp/~akiko-y/Egg/index_egg_E.html 但我不知道如何將其轉換成MATLAB /倍頻代碼改變變量y %test_amplitude modultaion fs=1000; t=l