differential-equations

    -1熱度

    1回答

    是否有可能做這樣的事情 ---- ------- --->|+ |------->| 1/s |-----o---------> | -| | | | ---- ------- | /|\ | -------------------- 在MATLAB ODE m文件 沒有SIMULINK的使用情況如何? 如果可能,它是如何完成

    -1熱度

    1回答

    Boost.Odeint沒有解決兩點邊值問題(TPBVP)的代碼。我無法找到任何C++實現。你知道任何開源實現嗎? 謝謝!

    2熱度

    1回答

    我正在嘗試創建一個SIR模型(求解帶有兩個參數的3個微分方程),並在MATLAB上繪製解決方案,並且不斷收到錯誤'沒有足夠的輸入參數'。我真的不知道自己出錯的地方。下面是代碼: function dx = sir(t,x) dx=[0; 0; 0]; beta = .003; delta = 1; dx(1)= -beta *x(1)*x(2); dx(2)=beta*x(1)*x(2)

    1熱度

    2回答

    我正在使用MATLAB的ODE套件來解決確定神經元行爲的13個微分方程組。現在我想添加第二個神經元,它由同一組微分方程操作,但受到第一個神經元的影響。更重要的是,這第二個神經元也會影響第一個神經元。 (這兩個單元之間的前饋和反饋。) 有沒有一種方便的方法來做到這一點?我可以將微分方程分佈在兩個函數文件中,還是必須將它們複製到原始文件之下,以便在同一個文件中存在更長的方程列表?我希望能夠每個單元格有

    8熱度

    2回答

    如果有人能夠幫助解決以下問題,我將不勝感激。 我有以下ODE: dr/dt = 4*exp(0.8*t) - 0.5*r ,r(0)=2, t[0,1] (1) 我在兩種不同的方式解決(1)。 通過Runge-Kutta方法(4階),並通過ode45在Matlab中。我比較與分析解決方案,它由下式給出了兩個結果: r(t) = 4/1.3 (exp(0.8*t) - exp(-0.5*t))

    3熱度

    1回答

    這是以下形式的線性常微分方程的N維繫統的自學: dx/dt=Ax 其中A是該系統的係數矩陣。 我知道你可以通過確定A的所有特徵值的實部是否定的來檢查穩定性。如果存在任何純粹虛構的特徵值A,則可以檢查振盪。如果存在任何純粹虛構的特徵值A,則我在閱讀的書中的作者引入了用於檢測系統的穩定性和振盪的Routh-Hurwitz準則。這似乎是比計算特徵值更有效的計算捷徑。 如果您現在可以快速找到特徵值,那

    0熱度

    1回答

    我正在嘗試求解二階微分方程。我已經設法在matlab中使用ode45來解決它的一個角度(這取決於角度的函數中有一個因素)。在我的程序中,與ode45的輸入參數也有相同的角度。 我想要的是解決不同角度的方程。我測試了一個for循環,我只需通過ode45運行給定的角度(作爲輸入參數),並繪製出不同的曲線。 但是我不確定它是多麼正確,因爲我沒有改變ode45解決的函數的角度。 所以問題是我如何爲每個循環

    1熱度

    1回答

    我已經使用matlabFunction定義了一個函數。這裏是代碼: matlabFunction([a16;-((1+x16^2)/(2*x16))*a16],'vars',{x16,[a16]},'file','DE_19') 這似乎工作。但是,當我嘗試使用ode45解決由matlabFunction定義的微分方程時,出現錯誤。下面是代碼: [x,y] = ode45(@(x16,Y) DE

    1熱度

    2回答

    我有一個單元格數組,包含函數(帶有函數句柄),我想在for循環中評估它們。我要評估微分方程: X1 '= X2,X2 '= AX2-BX1 我的代碼是這樣的: init = [0,0]; F = {@(x1,x2) x2,@(x1,x2)(a*x2-b*x1)}; X0 = init; for i=1:10 X = X0 + c*F(init(1),init(2));

    1熱度

    1回答

    我打算使用Phabricator進行團隊代碼審查,是否可以只安裝和使用差異?我們需要的僅僅是代碼審查工具。謝謝!