問題タブ [newtons-method]

連立方程式とヤコビ行列を含む m ファイルにニュートン法を実装するプログラムを作成しています

これが私が始めた計算ファイルです：

私は最終的に収束を確認し、エラーノルムを保存し、収束せずに最大反復回数に達したかどうかも確認したいと考えています。ここまでは大丈夫そうですか？

演習として tanh(x) の根を計算しようとしています。

アルゴリズムが最初の推測を求めるニュートン・ラフソン法を使用しています。

アルゴリズムは、約 1 より大きい初期推測では収束しないと想定されています。しかし、それに到達する前に数学範囲エラーが発生します。

これは私が使用しているコードです

Q1_6Raphson(5)たとえば、実行しようとすると、次のようになります。

私が読んだことから、数値が大きすぎると数学範囲エラーが発生します。しかし、私が理解していないのは、私のコードで呼び出されたすべての関数が入力として 5 で問題ないということです:

だから問題は何ですか？数学範囲エラーの原因は何ですか?

まず、教授が間違った質問をしたかどうかはわかりません。とにかく、F(x)~U(0,1)CDFを生成しようとしましたF(x)=1-(1+x)exp(-x)(この CDF についてはx=g(F(x))、手計算では計算できませんでした)。そしてF(x)、質問が望むものを達成するために根を計算します。

0からまでINFのルート範囲uniroot()は論外だからです。したがって、ニュートン法を使用して記述します。

次に、私のコードは次のようになります。

ただし、F(x)が小さすぎてニュートン法の 1 ステップで結果が 0 未満になると、エラーが発生します。さらに、次のようにコードを修正しました。

ゼロに近い結果を拒否したため、明らかにコードが間違っています。

したがって、私の質問は、私のコードを修正して正しく計算できるかどうか、そうでない場合は、それを行う別の方法があるかどうかです。どんな援助でも大歓迎です。

このコードをhttp://blog.shay.co/newtons-method/から取得しました。

ある場合、数値の反復回数の適切な経験則は何ですか? (たとえば、「n/2 反復を使用する」。)

Matlab で Newton アルゴリズムの実装を作成しようとしています。

式を使用して関数を呼び出すと:

次のエラー メッセージが表示されます。

何が悪いのかわかりません。たぶん、Matlabに詳しい人が私を助けてくれるでしょう。

以下は私のコードです：

Matlab との最初の接触です。

アップデート：

今、私はエラーを持っています:

私が入力した：

私のプログラム全体は正しいです（さまざまな段階でチェックしました）。ただし、このモジュールで強調表示されている行は、次のエラーを返します。

なぜこれが起こっているのかわかりません。funcPsatfloat 値を返します。有益なアドバイスをいただければ幸いです。