問題タブ [integral]

http://smsoftdev-solutions.blogspot.com/2009/08/integral-histogram-for-fast-calculation.htmlから抽出した、このコードを使用しています。

現在、そのページのコードをopencv2 c ++実装にインポートしています。これは、これまでのところです。

主要

働き

プログラムはコンパイルされますが、実行すると停止し、デバッガーはSIGSEVエラーを表示し、61行目と38行目を表示します。

38-extractIntegralHist（frame、9）;

61-bins [9] = Mat :: zeros（in.size（）、CV_32F）;

すべてがうまくいくように思われるため、問題を特定できませんでした。コードを確認してください。

編集。アップデート。関数extractIntegralHistが実行されると、プログラムが停止しているように見えます。これは、（関数の）行61にあるとされています。これがメモリ処理の例外であることはわかっていますが、何が問題なのかがわかりません。さらに、デバッガーが何を言っているのかわかりません。

ありがとう。

答えるのは簡単な質問があると思います。

私はこのコードを持っています：

次に、私は定義します

今、私はこの積分を実行したいと思います

phi_iiの代わりにphi_11を使用します。

すべてのベクトル化された変数があるので、Matlabでどのように実行できますか？

よろしくお願いします。

WKR、フランチェスコ

私は次の積分を評価しようとしています：

次の多項式の面積は次のように見つけることができます。

最初にシンプソンの法則による積分を使用する

結果はarea evaluated by Simpsons rule : 11.483072 次のようになります。次のコードは、ガンマ関数を使用して上記の式の合計を評価します。

そしてこれは11.4831を返します。関数とほぼ同じ結果quad。ここで私の質問は、累積分布関数を作成して、逆CDF変換を使用してこの分布からサンプルを取得できるようにするため、このネストされたループを取り除くことができるかどうかです。（累積分布関数を作成するためにgammainc、代わりに不完全ガンマ関数を使用しますgamma）

多項式係数が異なる可能性のある密度からサンプリングする必要があり、速度が問題になります。私はすでにモンテカルロ法を使用してそのような密度からサンプリングすることができますが、スピードアップするために密度からの正確なサンプリングを使用できるかどうかを確認したいと思います。事前にどうもありがとうございました。

私は、このソリューションを適用して、データ実験から定義された任意の分布で p 値を見つけようとしています。R の密度関数を使用してこの分布を推定しました。この関数を統合して、@mpiktas によって提案されたソリューションを適用したいと思います。ただし、関数を統合するには、密度が提供する関数を定義する値を持つ2 つのベクトルxとyではなく、入力として関数が必要です。

Rのx - y値に基づいてこの数値積分を処理する方法についてのアイデアはありますか?

私はこのコードを持っています：

今私は計算する必要があります

どこ

積分を計算する手順は次のとおりです。

現在、このループを使用しています

k>400近似を無視します。私のループは正しくないと思います。

積分をどのように計算することを提案しますか？

よろしくお願いします。

WKR、フランチェスコ

PIDコントローラーの不可欠な部分を理解していません。ウィキペディアからのこの擬似コードを想定しましょう:

Integral は、最初はゼロに設定されています。そして、ループ内で、時間の経過とともにエラーを統合しています。測定値または設定値に (正の) 変更を加えると、誤差は正になり、積分は時間の経過とともに (最初から) 値を「食い尽くし」ます。しかし、私が理解していないのは、エラーが安定してゼロに戻ると、積分部分にはまだ何らかの値（時間の経過とともに積分されたエラー）があり、コントローラーの出力値に寄与することですが、そうではありません。

誰か説明してくれませんか？

私は次のことを実行したいと思います：

どこ

後者の図に示されているパラメーターは、次のように取得できます。

たとえば、

Matlab は次のエラー メッセージを取得します。

この問題を克服する方法の手がかりはありますか?

ご連絡をお待ちしております。

よろしく、 FPE

MATLAB を使用して、次のような方程式を解こうとしています。

B = alpha*Y0*sqrt(epsilon)/(pi*ln(b/a)*sqrt(epsilon_t))*0 から pi までの (2*sinint(k0*sqrt(epsilon*(a^2+b) ^2-2abcos(シータ))-sinint(2*k0*sqrt(イプシロン)*a*sin(シータ/2))-sinint(2*k0*sqrt(イプシロン)*b*sin(シータ/2)) ) シータに関して

イプシロンは未知数です。

と を使用して積分に未知数が埋め込まれた方程式をシンボリックに解く方法を知っていますが、シンボリック積分器を使用するint()とsolve()、int()この複雑な方程式には時間がかかりすぎます。quad()、quadl()およびを使用しようとするとquadgk()、未知数が積分にどのように埋め込まれているかを処理するのに苦労します。

非統合テンプレート定数を初期化しようとしています。

コードの下を見つけてください：

問題は線にあります

非積分型であるため、定数をインラインで初期化できません。

私がオンラインで読んだことから、.cppファイルでconst初期化を移動した場合、問題はなくなったようです。ただし、テンプレート化されたクラスを扱っているため、これを行うことはできません。以下に詳述するようなエラーが発生します。

誰かが私のためにそれを見てくれませんか？よろしくお願いします。

私はC++でintをInfinityに設定するを読んでいました。真の無限大が必要なときは、numeric_limits<float>::infinity();を使用することになっていることを理解しています。その背後にある理論的根拠は、通常、整数型には、IEEE 754フロートのようにNaN、Infなどの特別な状態を表すために指定された値がないことだと思います（ここでも、C ++はどちらも義務付けていません- int＆float使用は実装に任されています）。しかし、それでもmax > infinity、特定のタイプについては誤解を招く可能性があります。私は、標準でこの呼び出しの背後にある理論的根拠を理解しようとしています。infinityあるタイプに意味がない場合は、フラグの有効性をチェックする代わりに、それを禁止するべきではありませんか？