問題タブ [integral]

MATLABでの二重積分の実装について質問があります。

それは知られています

を利用する

k1 = 1E-04：0.001：1E + 04; k2 = 1E-04：0.001：1E + 04; k3 = 1E-04：0.001：1E + 04;

上記の手順（F11、F22、およびF33の式を呼び出す）により、以下に示す結果が得られます。

ここで質問があります：

phi11、phi22、phi33のみを含む二重積分（またはk2とk3の入れ子になった単一積分）を使用して同じ結果を達成したいので、F11、F22、F33の式を直接呼び出す必要はありません。

今まで私はこのコードを使用しています（F11計算に制限されています）：

最後の図に示されているのと同じベンチマーク結果にはつながりません。

この問題にどのように取り組むことを提案しますか？よろしくお願いします。

よろしく、FPE

私は持っています

その後、F_33を次のように計算する必要があります

さて、私の質問に行きましょう。私は論文から次のことを知っています：

したがって、phi_33の式は次のようになります。

問題は、現在（forループ内で）使用している長い式を挿入したこの最終式をどのように利用できるかということです。

phi_33の最後の式は、処理が簡単で（特に、前者の記述における無謀な間違いのため）、独立変数である「参照渡し」（k2、k3）になります。

どんなヒントも大歓迎です。

よろしく、fpe

次のディストリビューションがあります。

左の分布の下の 0.95 の領域と右の領域の 0.05 を分離する値を見つけたいと思いますx(信頼性の片側 95% 区間)。

経験曲線を関数に当てはめ、関数を統合して目的の値を取得する必要があると思いますが、どこから始めればよいかわかりません。

これはRでどのように行うことができますか?

写真に表示されているこれら2つの方程式に問題があります。

1. 2 つの式で C(m) と S(m) を表す 2 つのベクトルがあります。これらの方程式を Matlab で実装しようとしています。連続積分演算を行う代わりに、総和を行うべきだと思います。たとえば、最初の式

   A1 = sqrt(合計(C.^2));

私は正しいですか？また、||dM|| を含む方程式 2 を実装する方法がわかりません。助けてください。

2. これらの 2 つの方程式の数学的な意味は何ですか? 最初のものは「二乗和」に関連している可能性があると思います.C(m)がベクトルの場合、この式はベクトルCの確率変数の総分散またはベクトルCの何らかの平均を測定しますか? 2つ目はどうですか？

どうもありがとうございました! A.

離散ベクトル フィールドを連続ベクトル フィールドに補間できる Mathematica スクリプトを作成しようとしています。次に、連続ベクトル場から積分曲線 (磁力線) を抽出します。

ListStreamPlotMathematica の関数が積分曲線を視覚的に表示できることは知っています．ただし、それらを構成する実際のポイントが必要です。

私は、次の効果を持つ力ずくの方法を書くことを考えました: 積分曲線は再帰的に決定できるようです。リストでポイントが選択された場合、積分曲線の「次の」ポイントは、最も近い勾配を持つポイントになります。(各「ポイント」は実際にはベクトルであるため、それぞれに独自の勾配があります。)「次の」ポイントの候補は、現在のポイントの上、下、左、または右の 4 つだけです。「次の」点が新しい開始点になり、グリッドの端に到達するまでプロセスが繰り返されます。

しかし、Mathematica の上位関数を使ったもっと簡単な方法はないのでしょうか。（私の方法が概念的に正しいかどうかもわかりません。）

このタスクに関するご意見やご提案をいただければ幸いです。

ありがとう！

誰かが私を助けて、何が問題なのか教えてもらえますか？いくつかの積分を計算する必要があり、このエラーが発生し続けます。

例：

SymPy では、不定積分に制限を適用して評価することは可能ですか?

Y_def.evalf()数を生成します。

Y_indef.evalf((theta, 0, sympy.pi/2))同じ答えが得られるようなものを探しています。

私は、ユーザーが長方形の量、開始、および停止を指定する関数の積分を見つけるためのプログラムに取り組んでいました。

注：長方形の左端の点を使用しています。

私は機能が完璧に機能しています（少なくとも、完璧に見えます）。しかし、私はそれのためにワンライナーを書くことができるかどうかを見たかったのですが、私が使用しているのでどのようにすればいいのかわかりませんeval()。これが私の元のコードです：

これは正常に機能します。

実際の答えは1000000/3ですので、私の関数はかなり良い見積もりを出します（1000個の長方形のみ）。

ワンライナーでの私の試み：

ただし、セミコロンを使用しているため、これは真のワンライナーではありません。また、それは少し遅く（数秒長くかかります、これはかなり重要です）、間違った答えを与えます：

それで、この関数にワンライナーソリューションを使用することは可能ですか？eval()実装方法とfloat(x)/n同じリスト内包表記がわかりませんでした。float(x)/n関数内で仮想の「ステップ」を実現しrangeます。

ありがとう！

特定の関数の積分を計算することを目的としたプログラムをMatlabで作成しました。積分には台形法を使うことになっています。現時点では、コードは次のようになっています。

問題は、私の関数がyの単一の値を返すのは、それがもう少しあるはずなのにだけだということだと思います。複数の値を返すように関数を書き直すにはどうすればよいですか？それとも、ここで間違っているのは何か他のものですか？

前もって感謝します！終わり