differential-equations

4熱度

1回答

我對Julia很新，我目前正在學習如何用它來求解微分方程。我試圖通過Christopher Rackauckas運行一個簡單的預製代碼，但我得到了一個錯誤。代碼可以找到here。我在這裏也寫： using DifferentialEquations alpha = 0.5 #Setting alpha to 1/2 f(y,t) = alpha*y u0 = 1.5 prob = ODEP

0熱度

1回答

離散微分方程並在每個時間步中找到參數的值

我想通過在t上使用for循環來找到f(t)的值。管理f的公式爲： df/dt = cos(f) g(x,y,t) x和y在二維座標是。我寫了下面的代碼： dt=1; a=2; f(1,1)= a*rand(1,1); for t = 2:100 f(1,t)= cos(f(1,t-1)) g(x,y,t) dt; end 但我不確定此方法是否正確，或者如果可以寫得更好。任何

1熱度

2回答

兩個For循環同時工作？ bvp4c

我想解決使用bvp4c（邊界條件）的微分方程系統。我使用兩個for循環來解決bvp4c的功能，但問題是循環先完成，然後只有最後一個值用於第二循環的迭代，是否有辦法讓它們同時工作？這意味着第一個for循環的第一次迭代使用第二個循環的第一次迭代值（而不是最後一次）。感謝 function RTrajfoll(X,Y) clf; for i = 1:length(X)-1

3熱度

1回答

爲什麼DifferentialEquations.jl中的tableau和explicit解算器？

我正在查看DifferentialEquations.jl程序包。在 DiffEqDevTools/src/ode_tableaus.jl我可以看到的桌面 Midpoint和RK4。但我也可以看到這些方案的明確代碼 OrdinaryDiffEq/src/integrators/fixed_timestep_integrators.jl。我有點期待使用tableaux而不是是一個明確的求解器。

0熱度

1回答

desolve微分方程不工作

我寫了一些代碼，有一個邏輯增長組件（即隨着N接近'承載能力'增長速度較慢，直到它達到'承載能力'停止增長）。但是，當我在R中運行它似乎並沒有工作。一些人口最終大於承載能力。我看了數學和一切OK。所以我認爲問題在於dN/dt只是針對每個人羣計算一次。有誰知道如何解決這個問題？任何幫助將不勝感激！實施例以下代碼： library('optimbase') library('deSolve')

1熱度

2回答

在八度中，是否有辦法爲一個變量求解兩個變量中的方程

我有一個MATLAB腳本，它使用拉普拉斯變換解決了非均勻的一階線性IVP。（在這個例子中，腳本設置爲解決IVP ，） syms x(t) s X; a0 = -3; x0 = 4; rhs = t^2; lhs = diff(x,t) + a0*x; ode = lhs - rhs Lx = X; LDx = s*X - x0; LHS = LDx + a0*Lx; RH

0熱度

1回答

Swift 3中的Lotka-Volterra模型

我正在嘗試實現一個rk4函數來求解2個微分方程。我有這樣的代碼，實現了龍格庫塔4方法： //RK4 method func rk4_func(y_array: [Double], f_array: [(([Double], Double) -> Double)], t_val: Double, h_val: Double) -> [Double] { let length

2熱度

2回答

定義一個函數用於從數據文件求解方程

我對Python和實際上任何基本的編程語言都是全新的，我將Mathematica用於我所有的符號和數字計算。我正在學習使用python並發現它非常棒！這是我正在努力解決的一個問題，但沒有任何線索而陷入困境！我有一個數據文件例如 0. 1. 0.01 0.9998000066665778 0.02 0.9992001066609779 ... .. 哪些只是{噸，COS [2噸

1熱度

1回答

根據條件改變微分方程中的參數？

我想根據位置條件改變移動機器人的輸入輪速度？ if d > 0 & d < 0.4 p.WR = 0.51; p.WL = 0.5; elseif d > 0.4 & d < 0.8 p.WR = 0.5; p.WL = 0.51; elseif d == 0.4 p.WR = 0.5; p.WL = 0.

-1熱度

1回答

在desolve

二進制開關我在試圖解決在R A ODE跑起來對抗一個問題。我有一個參數Q，當它比另一個參數h更大時停止流入它。直到它到達的地方發生切換，然後停止運行，並給我的消息的時間點的ODE正常工作： DLSODA- At current T (=R1), MXSTEP (=I1) steps taken on this call before reaching TOUT In abov