numerical-integration

0熱度

1回答

後解決與ODE45的ODE我用下面的代碼的工作，並不能找到什麼是錯的： xx = 0:1/50:1; v = 3.*exp(-xx)-0.4*xx; xq = xx; vq = @(xq) interp1(xx,v,xq); tspan = 0:1/50:1; x0 = 3; [~, y2] = ode45(@(t,x)vq(x), tspan, x0); 我得到y2 = [3;

8熱度

2回答

與分析解決方案相比，ODE45和Runge-Kutta方法的絕對誤差

如果有人能夠幫助解決以下問題，我將不勝感激。我有以下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))

1熱度

1回答

迭代標誌

我求解ODE像 Y ''（T）+ A Y（T）+ B = 0 用Matlab的ode45求解器。它迭代直到找到滿足錯誤標準的解決方案。我想在每個步驟中讀取之前的值y進行比較。這樣做的最好方法是什麼？我的主要文件是 [t,y] = ode45(@odefnc,tspan,[0 0]',options); 和odefnc是 function dx=odefnc(t,x) ... if

2熱度

1回答

Matlab中的擺動彈簧ODE系統 - 如何使位置矢量跟隨路徑？

我對Matlab很陌生。我有一個腳本，使用ode45和arrow.m來顯示一個擺動的彈簧在Matlab中的質量隨着它在三維空間中移動的運動。該方案是差不多做我想做的。現在鑽石的密度顯示出有效的彈簧速度（除了當ode45需要個人最喜歡的數字樣本時），並且速度幾乎可以精確地與功能的步長一起考慮（至少在速度我的電腦正在運行代碼）。我想要做的是，我在代碼中註釋的位置矢量只顯示在質量的瞬時位置，即曲線的終

0熱度

1回答

Matlab的ode45不發送初始條件到函數

我有一個快速的問題，我希望有一個相當簡單的解決方案。我有一個函數，如下所示： function [dx] = myFunction(x) 凡x是用含有的結構的系統的（位置和動量）四個元素的列向量。我試圖用ode45幾個粒子在這個配置空間數值積分，就像這樣： [Tout,Yout] = ode45(@(x,t)myFunction(x),[t0 tfinal],x0); 每次我嘗試這個我失敗

1熱度

2回答

數值二重積分中的R

我試圖計算二重積分數值在R值我的代碼和我想代碼公式如下： Expectedrsquare<-function(theta=.05,sigma=.0135,alpha=0.05,r0=0,Maturity=20){ library(mvtnorm) DoubleIntegration<-integrate(function(s) { sapply(s, function(s) {

0熱度

2回答

函數在數值積分中的誤差

我正在尋找一個函數，該函數使用高斯積分或辛普森積分在數值積分中產生重大誤差。

2熱度

1回答

Matlab的：ODE45輸出不正確的用力彈簧質量阻尼器

這裏是我的Matlab代碼來解決二階常微分方程的質量彈簧阻尼系統： function Spring clear all; close all; options=odeset('RelTol',1e-6); p0 = [1 0]; %initial position and velocity [t,p] = ode45(@SpringFunction, [0 20], p0, optio

6熱度

1回答

集成與陽性面積numpy的陣列只

我希望能夠計算下面的積分使用numpys trapz功能 numpy.trapz([-1, 1]) # returns 0 但我不希望允許負的區域。有沒有一種有效的方法來做到這一點，還是我必須查找最小點並手動轉換數組？ numpy.trapz(numpy.abs([-1, 1]))有意義嗎？

1熱度

1回答

求解一個ODE給出功能時作爲離散值-matlab-

我有以下ODE： x_dot = 3*x.^0.5-2*x.^1.5 % (Equation 1) 我使用ode45來解決它。我的解作爲dim（k x 1）的向量給出（通常k = 41，它由tspan給出）。另一方面，我做了一個近似於（1）模型的模型，但爲了比較第二個模型的準確性，我想通過ode45來解決它（解決第二個ODE）。我的問題是，這第二賦被給定的離散： x_dot = f(x) %