問題タブ [nonlinear-functions]

SAS-proc nlin について質問があります。

10000回のシミュレーションの手順を実行しています。それらの多くは収束せず、間違った結果をもたらします。

この反復が収束しなかったことを示すバイナリ変数を出力テーブルに追加したいと考えています。

誰もそれを行う方法を知っていますか?

どうもありがとう、

ペリー

方程式を満たす正方行列 A を見つける必要があります

ALA = -17/18A -2(ALL + LAL + (LL).A) + 3(AL + LA) -4L.LL + 8L.L - 44/9L + 8/9*(ID)

、ここで、L は対角行列 L = {{2/3,0,0,0},{0,5/12,0,0},{0,0,11/12,0},{0,0 ,0,2/3}}。

A が 2 次元と 3 次元の場合は答えを見つけることができますが、4 次元以上に問題があります。

実際には、行列 A も方程式 AA = A を満たす必要がありますが、適切な行列 L を使用すると、上記の方程式だけで十分です。

これは私のコードです。

A = テーブル[a[i,j],{i,1,4},{j,1,4}]

B = ALA

ID = IdentityMatrix[4]

M = -17/18A -2(ALL + LAL + (LL).A) + 3(AL + LA) -4L.LL + 8L.L - 44/9L + 8/9*(ID)

diff = (B - M)//ExpandAll//Flatten (ここで 16 の非線形連立方程式を取得します)

A1 = A/.Solve[差分 == 0][[1]]

このコードをしばらく実行した後、計算するのに十分なメモリがありませんというエラーが発生します。

この場合、16 個の方程式と 16 個の変数があります。エントリのいくつかはパラメータですが、結果が得られるまでどのエントリかわかりません。

とにかくこの問題を解決する方法があるかどうかはわかりません。理論的に可能な合理的な（おそらく整数）答えが必要です。

この問題は、行列方程式または他の方法で解決できますか? これに関する 1 つの問題は、方程式と変数が多すぎることです。

R で ggplot を使用し、関数を使用してプロットを作成しましnlsた。これは非線形回帰であったためです。グラフからデータを抽出したいと思います。具体的には、グラフの 2 つの曲線の値を知りたいです。値は 0 で、R を使用して実行する必要があります。

出来ますか？

データポイントペアのセットがあります(y,x)。フォームを使用して関数を適合させたい

c の値を求めたい。Matlablsqcurvefitはこれに対応していません。ローカルオプティマに引っかかっているようです。

それを行う方法に関する提案はありますか？

ありがとう！

pythonを使用してコーディングするには、次の問題があります。

x、y、z、t、HF、M1F、M2F の 7 つのパラメータがあります。ユーザーはこれら 3 つのいずれかを入力する必要があり、プログラムは残りを計算する必要があります。

私が持っている関係は次のとおりです。

問題の解決を試みます。

7 つのパラメーターがあり、ユーザーは 3 を入力する必要があります => 4 つのパラメーターが残ります。つまり、4 つの未知数を持つ 4 つの非線形方程式系を解くことがすべてです。

scipy.optimize非線形方程式系を解くために使用できるオンラインを読みました。しかし、私は最初の推測が必要です。

問題の物理に戻ると、次の初期条件があります。

を使用して問題を解決できるように、これらの初期条件はどのように初期推定を得るのに役立ちますscipy.optimizeか?