問題タブ [nonlinear-functions]

次の連立方程式を記号的に解こうとしています。

関数 solve(Sys,[x,y,z,w,q,r,t]) を使用して得た: 警告: 明示的な解が見つかりませんでした。

169で解決中

しかし、Mathematica で同じ連立方程式を解こうとすると、2 つの解が見つかります。私は何か間違ったことをしていますか？

ありがとう！

コスタキスの紙の解を再現しようとしています。この論文では、de Heligman-Pollard モデルを使用して、要約された死亡表を完全な生命表に拡張します。モデルには、適合する必要がある 8 つのパラメーターがあります。著者は修正された Gauss-Newton アルゴリズムを使用しました。このアルゴリズム (E04FDF) は、コンピュータ プログラムの NAG ライブラリの一部です。Levenberg Marquardt は、同じ一連のパラメーターを生成するべきではありませんか? LM アルゴリズムのコードまたはアプリケーションの何が問題になっていますか?

変数 x[1:3] のベクトルとその共分散行列 V[1:3,1:3] の推定値が与えられた場合、次のような 2 つ以上の変数に線形関数または非線形関数を適用するユーティリティを探します。

x[4] = x[1] + x[2]

x[5] = x[4] / x[3]。

x[4] と x[5] の推定値には単純な代数が必要です。

線形変換 x[4] を含む共分散行列は単純に (H * V * H') です。

H =

x[5] を追加した共分散行列は、項 a と b の 1 次テイラー級数近似で推定できます。

H =

概念的には、アルゴリズムがどのように機能するかを知っています。しかし、ユーザー インターフェイスである種の汎用方程式パーサーを試みる場合は特に、多くのコーディングが必要になります。

この問題に対応する既存の R ライブラリはありますか?

Cox 比例ハザード回帰を使用して、何らかの連続危険因子と私の好みの疾患との推定 U 字型関連を調査したいと思います。これまでのところ、私はこれを行ってきましたがR、survivalパッケージは次のようになっています

ただし、ドキュメントは残念ながらあまり冗長ではありませんが、これらのスプラインは制限されていない (つまり、テールで線形にする必要がある) ため、テールで大きな分散を示しているように思えます。必要な近似の複雑さを (df引数によって) 提供します。フィットの複雑さを自動的に調整できるパッケージのようにしたいのmgcvですが、AFAIKはcoxモデルをサポートしていません。

でそれを行うことはできますRか?

三角関数を使用した 2 次の非線形 ODE があるため、2 次導関数によっては定式化できません。例えば：

arctan(y'') のような用語が存在しなければ、次のように自分の頌歌関数を書くことができます

残念ながら、二次の非線形項 (=> b*arctan(y'') ) により、 y'' に依存して頌歌を書くことができなくなります。

このような三角法をMatlabで数値的に解く方法はありますか?

この男がMATLAB で行ったことを Python で実行したいと考えています。

anaconda をインストールしたので、numpy と sympy のライブラリがあります。これまでのところ、numpy nsolve を試してみましたが、うまくいきません。私はPythonを初めて使用し、MATLABでそれを行う方法を知っていると言わなければなりません:P.

式:

-2*log(( 2.51/(331428*sqrt(x)) ) + ( 0.0002 /(3.71*0.26)) ) = 1/sqrt(x)

通常、私はこれを反復的に解決します。左側の x を推測し、右側の x を解くだけです。解を左側に置き、もう一度解いてください。左の x が右に近づくまで繰り返します。私はどのような解決策が必要かという考えを持っています。

だから私はそれを行うことができましたが、それはあまりクールではありません。数字でやりたい。カシオの15ユーロの電卓でそのまま解けるから、そんなに複雑じゃないかな？

ご協力ありがとうございました、

編集：だから私は次のことを試しました：

そして、私はこの結果を得ます：

定数も解いているからですか？

edit2: args=() キーワードを介して定数を割り当てる必要があると思うので、次のように変更しました:

しかし今、私はこれを取得します:

とにかく、別の方程式を入れると、それが機能するとしましょう2*x*Re+(k*d)/(x+5)。そのため、方程式を変換する必要があると思います。

ここで死ぬ：log（x、10）..

edit4：正しい構文はlog10（x）です...今は動作しますが、結果としてゼロになります

Equationstomatrix 関数を使用して線形方程式を行列に変換する方法を知っています

残念ながら、equationsToMatrix非線形方程式には使用できません。複数の非線形方程式を行列に変換したい場合、これを行う方法はありますか? たとえば、次の 3 つの方程式があります。

次の結果を取得したい

X、Y、およびZを解きたい3 つの未知数を持つ 4 つの非線形方程式があります。方程式は次の形式です。

...ここでa、bとは、4 つの方程式の のc各値に依存する定数です。F

これを解決する最善の方法は何ですか？