Breaking down a numerical integration into two parts
1 回表示 (過去 30 日間)
古いコメントを表示
Sorry for this stupic question.
Suppose we have this integral
Of couse this can be solved analytically. But I still want to solve it numerically as:
dx = 0.01;
x = 1:dx:4;
fx = x.^2 + 5*x;
I = trapz(x,fx)
Now, we break down this integral into two part. Analytically, we write it as:
I am confused about the lower limit of the second integral when solving it numerically. Should it be 2 or . If we use 2, will not it mean that we are counting one point two times, once in and once in ?
0 件のコメント
採用された回答
Torsten
2024 年 3 月 6 日
I1 = dx/2 * f(x0) + dx * sum_{i=1}^{n-1} f(xi) + dx/2 * f(xn)
I2 = dx/2 * f(xn) + dx * sum_{i=n+1}^{n+m-1} f(xi) + dx/2 * f(x_{n+m})
->
I1 + I2 = dx/2 * f(x0) + dx * sum_{i=1}^{n+m-1} f(xi) + dx/2 * f(x_{n+m})
Thus since the endpoints in the trapezoidal rule are only counted half, all is fine.
その他の回答 (0 件)
参考
カテゴリ
Help Center および File Exchange で Numerical Integration and Differential Equations についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!