問題タブ [quad]

ネストされたscipy.integrate.quad呼び出しを使用して、2次元被積分関数を統合しています。被積分関数はnumpy関数で構成されているため、入力をループして入力ごとに1回呼び出すよりも、入力の配列を渡す方がはるかに効率的です。numpyの配列により、約2桁高速になります。

ただし、被積分関数を1つの次元のみに統合したいが、他の次元に入力の配列がある場合は、問題が発生します。「scipy」クアッドパックパッケージでは、それが何でもできないようです。 numpyは配列された入力を処理します。他の誰かがこれを見たことがありますか-そして/またはそれを修正する方法を見つけました-または私はそれを誤解していますか？クワッドから得られるエラーは次のとおりです。

私がやろうとしていることの漫画版を下に置きました-私が実際にやっていることはより複雑な被積分関数を持っていますが、これはジャイストです。

肉は上にあります-下は私のポイントを示すためにベンチマークを行っています。

これは私がやろうとしていることです:

数値的に積分された配列 (データ範囲全体にわたって) を、特定の積分範囲にわたって "quad" 積分関数で除算します。

エラーが発生します：「オペランドは形状（n）（n、2）と一緒にブロードキャストできませんでした」およびnはtループの長さです

別の言い方をすれば、ループ内の「t」の値ごとに、(1) trapz を使用して一連の数値データを統合し、それを配列に変換します。（2）特定の制限の間でクワッドを使用して関数を統合し、これも配列に変換します。次に、(3) 1 つの配列を別の配列で除算して、最終的な配列を取得します。

明らかに何かが欠けています。どうすればいいですか？

私が決めた範囲で信号を統合する必要があります。

私はこの方法を使用しました：

y は 500 マイクロ秒で 2500 ポイントの配列です。たとえば、t0 は 700 で、t1 は 1500 です。0.2E-6 は、信号の 2 つの値の間のステップです。

それが正しいか？

ラムダ関数を使用して関数を python 2.7 に、たとえば 2 つの変数の 2 番目に統合する以外の方法はありますか? たとえば、この例の関数を x に統合したいと思います。

これは、次の方法で簡単に実現できます。

これを行う別の方法があるかどうか疑問に思っています。たとえば、似たようなものを使用するなど、はるかに直感的だと思います。

誰かが代替ソリューションを手に入れましたか？(ええ、関数を次のように定義して代わりにsqrfn(x,a)使用することもできますが、それは重要ではありません)。quad(sqrfn,0,10,args(2,))

計算でさらに使用される整数値の配列を作成しようとしています。問題は、integrate.quad が (answer, error) を返すことです。フロートではないため、他の計算では使用できません。2 つの数字のセットです。

私はQUADPACKについて学んでおり、scipy.integrate.quadを使用して0から無限までの積分を計算しました。非常に良い結果が得られましたが、ソフトウェアが適用した統合方法 (QAGI、QAWF など) を知りたいです。テクニックを印刷する方法はありますか？ソフトウェアは決定木を適用しますか? お時間をいただきありがとうございます。よろしく。

エラーメッセージを生成する小さなコードがあります（数式のためだと思います）。誰でも理由を知っていますか？

エラーログ：

ありがとう！

scipy で 1 つの変数だけで複数の変数を持つ関数の定積分を計算しようとしています。これは、私のコードがどのように見えるかのようなものです-

そして、このタイプのエラーを返します:

ただし、引数に数値を入れると機能します。しかし、y を数値ではなく y のままにしておきたいので、そうしたくありません。これを行う方法を知っている人はいますか？

編集：申し訳ありませんが、私が明確だったとは思わないでください。最終結果を y の関数にしたいのですが、y はまだシンボルのままです。