differential-equations

0熱度

2回答

我訪問了gnu gsl網站，我沒有找到那裏的例子來解微分方程直觀（特別是因爲它使用二階微分方程）。 https://www.gnu.org/software/gsl/manual/html_node/ODE-Example-programs.html#ODE-Example-programs 有人可以指導哪裏可以找到一個描述性的指南，以解決一個非常簡單的第一階微分方程。例如，supoose我的函

1熱度

1回答

MATLAB從dsolve公式中找出最大點

我正在求解兩個微分方程並繪製第二個微分方程的結果。我想在圖中的最大點找到值。但是，也許因爲我使用的是sym，我無法確定如何執行此操作。我知道在數學中，你只需找到一階導數並將其設置爲零。 % Constants k1 = 1/10; k2 = 1/1.3; k3 = 1/1.3; k4 = 1/(1/48); k5 = 1/3.9; t = 1:5; % Initial concent

-1熱度

1回答

Runge Kutta 4在C++中用於Lorenz系統

有人可以在我的代碼中找到錯誤嗎？它可以在特殊的點上工作，但容差不適合該方法。我的錯誤是非常簡單的，所以我不認爲它的問題。請幫忙。 #include <math.h> #include <stdlib.h> #include <stdio.h> #include <conio.h> #include <time.h> #include <iostream> #include <strin

0熱度

2回答

在MATLAB中定義並求解嵌套的非線性ODE

我給出微分方程y'' = -g + a(t)/m和a(t) = k*y'^2，其中y是t（時間）的函數。我的初始條件是y(0) = 600;和y'(0) = 0; 在MATLAB我知道如何與 ydd = diff(y,t,2) == -g + a(t)/m; 定義y''但我的事實，這是一個「套」非線性微分方程，我迷路了不太清楚如何定義它，更不用說在MATLAB中解決它。

1熱度

1回答

Matlab Picard方法 - 將現有的symfun分配給矢量

我正在研究matlab中的picard方法。這需要在一個被代入的函數上進行多次迭代以成爲積分多項式。現在我有一個現有的多項式SYMS x，它是由一些向量定義： for i = 1:degree+1 polynomial = symfun(polynomial + a(i) *x^(i-1), x); end syms t; y = symfun(zeros(1,maxItera

0熱度

2回答

我有麻煩，在OCaml的製作差分功能（不是語法問題）

0熱度

1回答

製作一次狀態變量reacjes閾值

作爲玩具例如，考慮微分方程，DX/DT = P X ODE = function(t, X, p){ list(p*X) } 人們可以從包deSolve沒有問題解決這個使用lsoda。 times = seq(0, 20, by=.1) out = lsoda(y = 1, times, ODE, p = 2) 不過，我想lsoda停止計算閾值已經達到後，說xMax = 1e5。

0熱度

1回答

無法使odeINT接受由循環生成的方程

我想使用Python的odeINT來集成由循環生成的多組方程。方程式都是耦合的，所以必須通過一次調用odeINT同時進行集成。問題是初始條件（「y0」）必須是列表或矩陣列表（而不僅僅是列表）。 odeINT給出了這個錯誤：「初始條件y0必須是一維的」。我想知道如何解決這個問題。這是一個代碼示例;非常感謝任何想法。 class network: def __init__(self):

0熱度

2回答

如何指定求解R中的微分方程的最終值（而不是初始值）

我想解決R中的微分方程（deSolve？），我沒有初始條件，但只有最終條件的狀態變量。如何才能做到這一點？典型的代碼是：ode（times，y，parameters，function ...）其中y是初始條件，函數定義了微分方程。

0熱度

1回答

Python中用於pde的有限差分求解器的交替方向隱式方法

我正在執行Alternating direction implicit method來解決FitzHugh–Nagumo反應擴散模型。我已經爲它找到一個Python實現例如在blog，但我認爲這是在該方法的錯誤 - 在模板這裏介紹：它不應該是一半時間步反應項乘以˚F？