問題タブ [integral]

微積分の先生は、台形則を使って与えられた区間の定積分を計算するプログラムを教えてくれました。プログラムされた関数が算術関数のように入力を受け取り、出力を生成することは知っていますが、その逆を行う方法がわかりません。つまり、出力が与えられたときに入力を見つけます。

問題は次のように述べています。

「t=0 から t=9 までの移動距離を見積もるには、刻みの数 n を変化させる台形則を使用します。n > D の場合、台形和がこの制限 (468) の 0.01 単位内にある数 D を見つけます。 ."

電卓を使って「プラグアンドチャグ」で制限を推定しましたが、通常の代数関数を使用すると、次のように簡単に実行できることがわかります。

limit (468) = 変数 x の代数式 (x について解く)

ただし、プログラムされた関数に対してこれを行うにはどうすればよいですか? 出力が与えられた場合、プログラムされた関数の入力をどのように決定しますか?

間隔 0 と 9 の間の多項式 (x^2+11x+28)/(x+4) の定積分を計算しています。計算機の台形規則関数は、間隔 0 と 9 の間の定積分を計算します。与えられた数の台形 n を使用します。

全体として、これを行う方法を知りたい：

n について解く: 468 = trapezoidal_rule(a = 0, b = 9, n);

私の TI-83 での trapezoidal_rule(a, b, n) のコード:

私はこの構文にもコンピュータ アルゴリズムにも精通していないので、誰かがこのコードを代数方程式に変換するのを手伝ってくれたり、そうする方向を教えてくれることを期待していました。

編集: これは私の宿題の一部ではありません。単なる知的好奇心です。

私は物理学の授業に取り組んでおり、現在このセクションで立ち往生しています。試してみましたが、うまくいきませんでした。ここで本当に助けが必要です。

台形の法則について、質問: 方程式 f(x)=x4*(1-x)44/(1+x^2) の積分値は何ですか?

これは私が試したコードですが、答えが得られません

sympy.Eq(func, var)ラムダで使用して生成した方程式を使用しようとしています。「追加」オブジェクトのリストを返しているようですが、使用方法がわかりません。strを使用する前に型キャストを試みましたevalがfrom sympy import I、問題と関係があると思います。

コードとエラーメッセージは次のとおりです。

py solid.py -vy = (x-1.0) (x-3.0) *2 1.0 3.0 -軸 y --メソッド ディスク

実行例を次に示します。

トレースバック (最新の呼び出しが最後):
ファイル "solid.py"、136 行目、印刷ボリューム (args.var、args.equation、args.a、args.b、args.axis、args.m) 内 ファイル "solid. py"、75 行目、ボリューム func = eval('lambda y: ' + '(' + func[0] + ')**2') TypeError: 'str' オブジェクトと 'Add' オブジェクトを連結できません

そして、最初にそれを str() に型キャストしようとすると、次のようになります:

ファイル「solid.py」、136 行目、印刷ボリューム (args.var、args.equation、args.a、args.b、args.axis、args.m) ファイル「solid.py」、76 行目、ボリュームintegral = integrate.quad(func, a, b) ファイル "/usr/lib/python2.6/dist-packages/scipy/integrate/quadpack.py"、185 行目、quad retval = _quad(func,a,b) ,args,full_output,epsabs,epsrel,limit,points) ファイル "/usr/lib/python2.6/dist-packages/scipy/integrate/quadpack.py", 行 249, in _quad return _quadpack._qagse(func,a ,b,args,full_output,epsabs,epsrel,limit) quadpack.error: 指定された関数は有効な浮動小数点数を返しません。

また、print func[0]与えます：

2.3333333333333 + (0.296296296296298 - y/2 + (-0.0877914951989024 + (0.592592592592595 - y) 2/4) (1/2)) (1/3)*(1/2 - I*3 (1/2)/2) + 0.444444444444444/((1/2 - I*3* (1/2)/2) (0.296296296296298 - y/2 + (-0.0877914951989024 + (0.592592592592595 - y) 2/4) (1/2))**( 1/3))

以前に完全なスクリプトを提供できなかったことをお詫び申し上げます。

MATLAB で次の積分を計算する方法を知りたいと思いました。

関数ファイル を作成しましたpn。これは次のとおりです。

とスクリプトbegin.m：

積分の答えは「1」でなければなりません。

方程式のテキスト入力を取り、それを定積分として評価するのに苦労しています。に渡す呼び出し可能な関数が必要scipy.integrateです。

Fortran や C と比較して、Matlab で Feynman path Integral を実行するのは時間がかかることを私は知っています。

しかし、誰かが経路積分による調和振動子のMatlabコードを持っていますか? 私はウェブ上で（そしてMatlabフォーラムでさえ）見つけることができませんでした。

Matlab に変換する方法がわからない Fortran コードの下 (私は初心者です) ありがとう、ジョニ

私は今論文を書いています。私は積分で立ち往生しました。私の関数は次のように定義されています

myfun =(exp(t*Q)*V*x)(j);

ここでQ、 とVは行列 (n*n) で、x要素が 1 のベクトルです。計算後、そのベクトルの j_th 要素を取得し、関数を t に対して積分する必要があります。

quadをmatlabで使いたいです。ただし、ポイントは、内部マトリックスが同じサイズではないことを報告することです。ここのAは数字じゃないから?...

これどうやってするの。そうしないと、t 自体に対してループを作成することしかできず、これは非常に低速です。

ありがとう

ここでは、2つのベクトルの積分の計算を実装するためにあなたの助けを借りたいと思います。Rに関連する積分計算のページをチェックしました。しかし、数学のトレーニングがほとんどないので、それでも自分でそれを行うことはできません。

私の目的は、この文のアイデアを実装することです。「位置ごとにレート推定値をプロットすると、遺伝子地図はこのプロットの積分にすぎません。」これは、私が変数（レート、ポジション）を持っていることを意味し、各ポジションには独自のレートがあります。各ポジションのレートの積分を計算したいと思います。ここでは、位置は単調に増加しています。

このタスクは、数学の計算に精通している人にとってはそれほど複雑であってはなりません。それで、それについての指示/指示を教えていただけますか？

前もって感謝します。

次の二重積分の評価に問題があります。

明らかに、cos 関数の振動性がこれらの問題を引き起こしています。f と g をより大きな数に増やすと、matlab は特異点と統合が失敗したことを訴えます。

x の -pi から pi および y の -pi から pi の制限を超えています。

私は使用しています：

助けやアイデアをください???

quad/quadl を使用して計算された積分をプロットする方法を誰かが知っているかどうか、またはこれが可能かどうかを知りたいです。

trace パラメーターをゼロ以外に設定できることを読みました。これにより、各反復の情報が提供されますが、その情報を使用して積分をプロットする方法と方法がわかりません。

ありがとう。