numerical-integration

2熱度

1回答

我正在使用RK4解決一個ODE系統。我正在生成一條直線圖，這似乎是由於k3_1的上限爲-3.1445e + 24。我不明白爲什麼它被封頂。 function RK4system_MNModel() parsec = 3.08*10^18; r_1 = 8.5*1000.0*parsec; % in cm z_1 = 0.0; % in cm also theta_1 = 0.0

2熱度

1回答

R繪製積分

我在R中的積分函數有一些問題。我試圖繪製積分vo，但它似乎我沒有正確地做。 t <- seq(0, 0.04, 0.0001) vi <- function(x) {5 * sin(2 * pi * 50 * x)} vo <- function(x) {integrate(vi, lower=0, upper=x)$value} test_vect = Vectorize(vo, vec

1熱度

1回答

Python中的Verlet集成導致粒子跑掉

我的下游代碼（應該是）解決兩個物體的運動方程，但結果是粒子運行方式，我無法找到錯誤的位置 import numpy as np import matplotlib.pyplot as plt plt.style.use('ggplot') DIM = 2 N = 2 ITER = 1000 def acc(r, v, a): for i in range(N - 1):

0熱度

1回答

使用Gauss-Legendre正交的雙重積分

我想使用Gauss-Legendre積分來評估雙積分。我想出了以下代碼 m=5000; %generate weights and abscissas [wx,xx]=leg(-1,1,m); [wy,xy]=leg(-1,1,m); %define function [email protected](x,y) hypergeom(-1./4,3./2,x.^2.*y.^6);

1熱度

1回答

具有更多參數和集成的函數

我有我簡單的問題，但我無法在任何地方找到解決方案。我必須整合一個函數（例如使用辛普森的規則子例程），但我不得不將一個以上的參數傳遞給我的函數：其中一個是我稍後要整合的變量，另一個是來自我不能在函數內部執行不同的計算。問題是Simpson子程序只接受f（x）來執行積分而不是f（x，y）。弗拉基米爾建議後，我修改了代碼。下面的例子： Program main2 !--------------

2熱度

1回答

ode45 Langevin方程

我有一個關於使用Matlab來計算隨機微分方程解的問題。這些方程式是this paper（PDF）中的2.2a，b，第3頁。我的教授建議使用ode45一小段時間步，但結果與文章中的不匹配。特別是時間系列和pdf。我對函數中白噪聲的定義也有疑問。這裏爲一體的功能代碼： function dVdt = R_Lang(t,V) global sigma lambda alpha W1=sig

0熱度

1回答

多個嵌套積分

我試圖整合以下形式的嵌套積分： s_min = 0; s_max = lam; t_min = 0; t_max = pi; pmax = @(x,y) y; pmin = 0; ymax = 1; ymin = @(x) x; xmax = 1; xmin = 0; eqn = integral(@(x) x.*integral(@(y) (1/y).*integral(@

-2熱度

1回答

Matlab圓周積分

這可能是一個簡單的問題。我正在嘗試理解Matlab中的數值積分。對於一個單位圓， fun = @(x) sqrt(cosd(x).^2+sind(x).^2); q = integral(fun,0,360); 答案是360.0000，而我希望它是2*pi。請注意，角度保持爲度。我很困惑。我究竟做錯了什麼？

0熱度

1回答

數值積分python

我需要在python中減少quad（）的運行時間（我正在集成幾千個積分）。我在這裏發現了一個類似的問題，他們建議做幾個集成並添加部分值。但是，這不會提高性能。有什麼想法嗎？這裏是一個簡單的例子： import numpy as np from scipy.integrate import quad from scipy.stats import norm import time

1熱度

2回答

計算梯形積分從原點到每個時間點中的R

我有一個數據幀x： head(x) # time Qfr #1 1 0.004751271 #2 2 0.005405618 #3 3 0.005785781 #4 4 0.006028213 #5 5 0.006179973 #6 6 0.006263814 我試圖計算從time = 0數值積分高達每一個時間點，即積分： \integral_{u=0}^t Qfr du