問題タブ [integral]

quad2dネストされたtrapzループで呼び出しをバイパスする方法があるかどうか知りたいです。私の問題についてもう少し詳しく説明します。現在、この方法で二重積分の計算を実行しています。

は次のようにPHI11定義されます

次の方法でE_integral.mat取得できます。

ここでの問題は、ネストされた関数を使用して、見落としquad2dてhandle functionより実用的なアプローチを支持することは可能でしょうか?trapz

これまでのところ、次のコードを試しましたが、期待どおりの結果が得られませんでした。

画像マトリックスを指定して、合計面積テーブルまたは積分画像を作成しようとしています。それが何であるかわからない人のために、ウィキペディアから：

合計面積テーブル (インテグラル イメージとも呼ばれます) は、グリッドの長方形サブセットの値の合計を迅速かつ効率的に生成するためのデータ構造およびアルゴリズムです。

言い換えれば、イメージ/マトリックス内の任意の長方形領域の値を一定時間で合計するために使用されます。

これを R に実装しようとしていますが、コードの実行に時間がかかりすぎているようです。

このリンクからの疑似コードは次のとおりです。inは入力行列または画像で、intImgは返されるものです

そして、ここに私の実装があります

入力行列と出力行列の例:

ただし、480x640マトリックスでは、これには約 4 秒かかります。論文では、それらの次元についてミリ秒単位でかかると説明しています。

ループやインデックス作成で非効率的なことをしていませんか?

C++ で書いて R でラップすることも考えましたが、C++ にはあまり詳しくありません。

ありがとうございました

第 2 種フレドホルム積分方程式に似た方程式があります。

この方程式を解くために、特定の方程式に収束することが保証されている反復解が与えられています。ここでの唯一の問題は、この反復法を MATLAB で実装することです。

今のところ、コードの問題のある部分は次のようになります。

delta_aの関数でxあり、反復の初期値を表します。とmsの関数です。xxi

お気づきかもしれませんが、反復では (積分の前) と(積分の内側) のdelt両方に依存したいと考えています。残念ながら、このコードの書き方 (関数ハンドルを使用) では、希望どおりの数値が得られません。は定義されていないため (定義するまで)、のとの の 2 つの異なる関数として記述することもできません。では、それが積分の内部に依存していることを確認しながら、反復から数値を取得するにはどうすればよいでしょうか?xxideltxxixiintegraldeltxi

これを解決する方法について何か提案はありますか?

数値積分の使用

入力パラメータの説明: x は数値のベクトルで、残りはすべて定数です。私のコードの問題は、入力パラメーター x が使用されていないことです (これは、x がシンボルとして扱われていることを意味すると思います)。

4つの異なる変数を持つ関数があり、変数が1つだけの場合の数値積分f(x, y, z, t)を見つけたい:quad

出来ますか？象徴的な結果が必要だと思います。他に何か考えはありますか？

統合する巨大な機能があります：

シムズXY

f=(228155022448185.*(cos((2.*pi).*y)./exp(131738205584307./(35184372088832*x)) - 1)*(cos((8.*pi.*y)/ exp(131738205584307./(8796093022208*x)) - 1)*(cos((8.*pi.*y)/exp(131738205584307./(8796093022208.*x)) + cos((18.*pi) .*y)/exp(1185643850258763./(35184372088832.*x)) - 2))/((18014398509481984. *(x.^2)).*exp(x. * ((1981232555272083.*(y.^)) 2))/2251799813685248 - y./16 + 1./16)))

(x:[0,inf) と y:[0,1]) を統合する必要がありますが、quad2d と dblquad でエラーが発生します。

これらのエラーが表示される理由を説明していただけますか? どうすれば修正できますか？