numerical-integration

0熱

2答えて

最大値の[x = value1、x = value2、...]という形式の要素のリストの並べ替え

解決方法またはすべての方法を使用すると、ソートされていない要素のリストが取得されます。問題は、私はそのリストをソートしたいが、それは[x = value1、x = value2、x = value3]の形式であり、ソート関数はこれらのx =式でうまく動作しないためできない。これを解決する方法はありますか？ルジャンドル多項式の根を得ながら背景私の問題は、ガウスの求積法の文脈において生じます。数式

1熱

2答えて

MATLABでIntegral2関数を使用したときのエラー

私が知る限り、誰もこれを尋ねていません。私は、関数の二重積分と同じ二重積分を計算するよう求められましたが、積分の順番は入れ替わっています（つまり、最初にdydx、次にdxdyに積分します）。ここに私のコードは次のとおりです。 %Define function to be integrated f = @(x,y) y^2*cos(x); %First case. Integrati

0熱

1答えて

時間依存入力のODE、補間を使用せずにスピードアップするには？

私はODEのシステムを解こうとしています。入力励起は時間の関数です。私はinterp1をインテグレーション機能の中で使っていますが、これは非常に効率的な方法ではありません。関数内でinterp1呼び出しを必要としないsin関数への入力励振を変更すると、はるかに高速な結果が得られるため、これはそうではないことが分かります。しかし、1ステップごとに補間を行うと、収束するのに10〜20倍の時間がかかり

0熱

2答えて

mathematica：数値積分の評価順序

こんにちは、私はmathematicaの数値積分に関して問題があります。ここに私のテスト機能 Table[NIntegrate[ Boole[Sqrt[1 - cosk^2]*Sqrt[1 - cosk2^2] > Abs[a - cosk*cosk2]]/ Sqrt[(1 - cosk^2)*(1 - cosk2^2) - (a - cosk*cosk2)^2], {cosk, -1,1}

6熱

1答えて

統合は

私はtrapz機能 numpy.trapz([-1, 1]) # returns 0 不可欠使っnumpys次計算することができるようにしたいしかし、私は、負の領域を許可する必要はありません。これを行う効率的な方法はありますか、または最小点を探して手で配列を変換する必要がありますか？ numpy.trapz(numpy.abs([-1, 1]))は意味がありますか？あなたは私たちが単にnp.t

0熱

1答えて

FFTによる振動機能のフーリエ積分/フーリエ変換

f(x) = cos(x^2)およびg(k) = pi^0.5 cos((pi*k)^2 - pi/4)は、フーリエ対である。私はそれがすべてで同意すれば、結果は非常に小さなk範囲にg(k)と一致するしかしすなわちがSum[ fn * exp(2 pi * ik x_n), {n, 0, N-1} ] * Delta_x でIntegrate[ f(x) * exp(2 pi * ikx),

1熱

1答えて

MATLAB積分関数の評価数

MATLABで "積分"サブルーチンを使用して積分を評価しています。「積分」は適応型直交ルールに基づいています。 "integral"で関数評価の数を返すことはできますか？

6熱

1答えて

パフォーマンスが

を統合し、私は積分はベクトル値ある一次元、で数値積分を実行したいと思います。 integrate()はスカラー積分を許すだけなので、私はそれを何度か呼び出す必要があります。 cubatureパッケージが適しているように見えますが、1Dインテグレーションでは非常にうまく機能していないようです。 library(cubature) integrand <- function(x, a=0.01) e

1熱

1答えて

Pyletでの粒子の統合による粒子の逃げ

2つの物体の運動方程式を解くための私の休暇中のコードがありましたが、その結果は粒子が走っていてエラーがどこで見つけられなかったのですか import numpy as np import matplotlib.pyplot as plt plt.style.use('ggplot') DIM = 2 N = 2 ITER = 1000 def acc(r, v, a):

0熱

1答えて

既知のパスで実行されるアニメーションシェイプ

私は、既知のパスを追跡する車両をシミュレートしたい、私が取り組んでいるのは、パスで移動するアニメーションシェイプです。私はmatlabの下に自分のコードを書きました。車両）はパス上を移動しますが、私は最後の四角形（記憶）を取り除きたいと思います。誰でも私がそれを達成するためにできることを私に助言してください。私のコードは、私が先に説明したものを示しています r=0.25; t=0:0.1:1