問題タブ [integral]

この積分を数値的に評価することに固執しており、単純な直交スキームを実行した経験しかないので、これが素人っぽく見える場合はご容赦ください。入れ子になった積分は、Gauss-Hermite(-inf,+inf) と Gauss Laguerre スキーム (0,+inf) の変形を使用して、ガウス密度とガンマ (v/2,2) 密度で以下に示します。積分がほぼ終了し、中間ステップは良さそうに見えますが、全体の積分を評価するために重みを組み合わせる方法に行き詰まっています。問題を解決するより良い直交スキームをコード化するために、コード/その他のアイデアを変更する際の提案に本当に感謝しています。\begin{equation} \int^{\infty} {-\infty}\int^{\infty} {0} \prod_{i=1}^n\Phi \left(\frac{\sqrt{w/v }\,C{i}-a{i}Z}{\sqrt{1a {i^2}}\right)f{z}(Z)f{w}(W)dwdz \end{式}

% スクリプトは、ノード、ウェイト、パラメーターを定義し、1 つのメイン関数と 1 つのサブ関数を呼び出します

ロー=0.3; nfirms=10; h=repmat(0.1,[1,nfirms]); T=1; R=0.4; v=8; alpha=v/2;GaussPts=15;

積分式で与えられる 2 つの変数の縮退超幾何関数を計算する必要があります。

そして、数値積分を取るためにMatlabを使用しました：

結果は 54.9085 で、この値が間違っていることはわかっています。ですから、特異点が 0 の上記の積分の真の値を計算するのを手伝ってください。

Matlabで次のタイプの積分を計算したい。

これは関数 e^-(u)*u の積分であり、境界はゼロと無限大です。この積分は 1 を返す必要があります。

Matlabでこれを行うにはどうすればよいですか?

Erlang が一般的にこの種の用途に使用されないことは知っていますが、ひょっとして、積分/微分が既に記述されている数値関数 (GNU Octave の quad のように、関数と 2 つの数値を取る) はありますか?

誰かがmatlabで以下を計算する方法を教えてもらえますか:

どこ

一部の整数型では、浮動小数点値が整数の表現可能な範囲をはるかに超えている場合でも、浮動小数点型の値に最も近い値を見つけるにはどうすればよいですか。

またはより正確には：

浮動F小数点型にします (おそらくfloat、double、またはlong double)。I整数型にしましょう。

と の両方が の有効な特殊化を持っているFと仮定します。Istd::numeric_limits<>

の表現可能な値が与えられ、FC++03 のみを使用して、表現可能な最も近い値を見つけるにはどうすればよいIですか?

私は、純粋で効率的でスレッドセーフなソリューション、および C++03 によって保証されていること以外はプラットフォームについて何も想定していないソリューションを求めています。

そのような解決策が存在しない場合、C99/C++11 の新機能を使用して解決策を見つけることは可能ですか?

C99 の使用lround()は、ドメイン エラーが報告される方法が自明ではないため、問題があるようです。これらのドメイン エラーは、移植可能でスレッド セーフな方法でキャッチできますか?

注: Boost はおそらくそのboost::numerics::converter<>テンプレートを介してソリューションを提供することを認識していますが、その複雑さと冗長性が高いため、そこから本質的な要素を抽出できていないため、それらのソリューションが有効かどうかを確認できませんでした。 C++03 を超える仮定。

I(f)次の素朴なアプローチは、 の整数部分がfの表現可能な値ではない場合、の結果が C++03 によって定義されないという事実のために失敗しますI。

次に、次のアプローチを検討してください。

F(std::numeric_limits<I>::min())と の不可欠な部分がでF(std::numeric_limits<I>::max())表現できない可能性があるため、これも失敗しますI。

最後に、これも失敗するこの 3 番目のアプローチを検討してください。

今回I(f)は常に明確な結果が得られますが、F(std::numeric_limits<I>::max())よりもはるかに小さい可能性があるため、 未満の複数の整数値である浮動小数点値をstd::numeric_limits<I>::max()返す可能性があります。std::numeric_limits<I>::max()std::numeric_limits<I>::max()

F(i)変換が最も近い表現可能な浮動小数点値に切り上げるか切り下げるかは未定義であるため、すべての問題が発生することに注意してください。

C++03 (4.9 Floating-integral conversions) の関連セクションは次のとおりです。

整数型または列挙型の右辺値は、浮動小数点型の右辺値に変換できます。可能であれば、結果は正確です。それ以外の場合は、表現可能な次の低い値または高い値の実装定義の選択です。

多くのデータポイントがあり、スライディング ウィンドウの曲線の下の領域を計算したいと考えています。しかし、それはかなり速いはずです。少しググって、Java で NewtonCotes の実装を見つけましたが、より高速な方法があるかどうかはわかりません。

何か案は？

特定のステップ関数の積分である関数を定義したばかりで、取得した三角形関数を積分したいと考えています。 jerk:1,-1,1-1 一連のジャークの特定のリストに対応するステップ関数。三角形関数を取得するには、

私は次の式を使用しましたa(t)=kt+ac。k はジャーク、ac は各ステップの連続性と変化を保証する定数です。しかし、式を使用するとv(t)=(1/2)kt^2+ac*t+vc、速度に負の値が得られますが、これは正しくありません。

積分の数値を返すいくつかの scipy 関数を見つけましたが、曲線が必要です。これを統合する方法について何か考えはありますか? ありがとうございます

みんな、私は誰かがこの問題を解決する必要がありますか? そのコードをコンパイルすると、次のエラーが発生します。

私はこの部分に問題があります:

どうすればそれを修正できますか

コンデンサの電圧と電流を計算しています。次に、エネルギーも決定したいと思います。エネルギーは電力の積分ですが、電力関数を統合することはできません:

エラーが表示されます: The integrand function must return a output vector of the same length as the input vector.

誰でも助けることができますか？