問題タブ [newtons-method]

MATLAB で Newton-Raphson 法を使用して、3 変数の 3 つの非線形システムを解こうとしています。3 つの非線形方程式を次に示します。

ニュートン・ラフソン法を使用して、、、cおよびsの値を見つける必要があります。q

これが私がこれまでに持っているものです：

誰かが私のコードをチェックしてください。いくつかのエラーがあるため、動作していません。前もって感謝します。

Matlab で Python に記述したバック ステップ コードを使用してニュートン法をコーディングしようとしていますが、Python 構文に問題があります。Matlab では約 5 回の反復が必要ですが、私の Python コードは最大反復回数 1000 までループし、バック ステップ メカニズムが機能しない (負の対数を計算しようとする) ため、ドメイン エラーが発生します。私はしばらく Python を使用していないので、何らかの構文を混乱させている可能性が高いです。

正しく動作する Matlab コードは次のとおりです。

Python コードは次のとおりです。

scipy.newton メソッドを使用して pandas データフレームを最適化しようとしています。

まず、私のデータフレームの作成は以下です。次に、関数 Px を作成します。3 番目に、別の関数 YieldCalc を作成します。ここでは、scipy.newton を使用して最適化し、Px = 0 のような Rate の値を見つけます。次に、その値を新しい列 'Yield' に追加しようとしていますが、次のエラーが発生します。どんな助けでも大歓迎です。前もって感謝します。

エラー/出力:

Newton-Raphson アルゴリズムを使用して、単精度ハードウェアを使用して IEEE-754 単精度浮動小数点値を除算しています。

次の 2 つのリンクで説明されている方法を使用しています。

1. ウィキペディア ニュートン・ラフソン課
2. 私が使っているニュートン・ラフソン法

ただし、Xi から X_3 を計算する (つまり、3 回の反復を使用する) にもかかわらず、私の答えはまだ少しずれています。なぜそうなのか不思議です。MATLAB を使用して結果を比較しています。

間違った結果の例を示す出力を次に示します。

MATLAB コードをここに添付しました。

テスト入力を含むファイルを次に示します (これは generatedfloats_for_fpdiv.txt という名前のファイルに保存する必要があります)。

mpmath.findroot の多次元ニュートン法を使用して、連立方程式を数値的に高精度に解こうとしています。システムの例を次に示します。

残念ながら、解決策は得られませんが、エラーがスローされます

ZeroDivisionError: 行列は数値的に特異です

これは、findroot がヤコビアンを計算しようとしているために発生していますか? これは、システムが過小決定されていることを意味しますか?

scipy の fmin_cg を使用して出発点を見つけましたが、解をより高い精度に磨きたいと考えています。fmin_cg の関数として、関数 f の 27 エントリの二乗和を最小化しました。

mpmath.findroot の問題を回避できない場合、このシステムを高精度で解決するためのより良い方法はありますか?

私はmatlabを初めて使用し、近似x = aを開始してニュートンラフソン法をn回反復する関数を作成する必要があります。この開始近似は反復としてカウントされず、別の要件として for ループが必要です。投稿された他の同様の質問を見てきましたが、私の場合は while ループを使用したくありません。

これは私の入力が想定されているものです：

これまでのコードは次のとおりです。

関数を呼び出そうとすると、次のエラーが表示されます。

問題は、この VPA 機能をどのように使用するかです。私のコードの残りの部分もおそらく100％正しいとは限りませんが、vpaの問題に対処したり、コードの他の部分を修正したりする助けがあれば大歓迎です.

ありがとう！