問題タブ [nonlinear-optimization]

次の機能を最小化しようとしました。

x+y=1 かつ 0<=x,y<=1 の場合。Rsolnp パッケージで gosolnp を使用するには、まず cons を定義して eqfun 引数で使用します。

次に、gosolnp 関数を適用しました。

答えは x = 0 および y = 1 である必要がありますが、gosolnp を実行するたびに試すことができるように、その時点で func がほぼ 0 である新しいポイントが得られます (厳密にはそうではありません)。Mathematica と Maple はこの関数の最適化を非常に高速に行い、x = 0 および y = 1 という真の答えを返しますが、代わりに R で実行するたびに、正しくない新しい解が得られます。

alabama または DEoptim で spg() などの別の最適化関数も試しましたが、問題は解決されませんでした。

だから私の質問は：

1-Rでfuncを最小化できる解決策はありますか?

2- R と Mathematica の精度に違いはありますか? Mathematica では正確な答えが得られるのに R では得られないのはなぜですか?

前もって感謝します

問題: 非線形方程式で複数の既知の値に対して未知数を解くことができない

(I) は不明で、(v) は既知の値で、残りは定数です
。vの 1 つの値についてこれを解くことができます。
（2）同様に、定数（Rs、Io、Il、Rsh、a）が8760行91列の異なるデータフレームにある場合、（I）をどのように解決しますか？

論文からプログラミングの定式化があり、特定の問題を解決するためのツールを提供したいと考えています。著者はそれを線形計画法 (LP) インスタンスと述べていますが、私にはよくわかりません。公式は次のようになります。

関数を介してプログラムしようとしましたcplexqcp（二次制約のためですが、制約にはx_i^2変数が含まれていません）。しかし、私は受け取りCPLEX Error 5002: Q in %s is not positive semi-definite errorます。これは非凸制約を伴う非線形計画法のインスタンスですか? ツールで解決しCPLEXたり、ツールを使用したりできNLPますか? 私はスタッフの初心者ですLP/NLP（スタッフに関するコースは受講していません）ので、質問の回答の詳細を説明するのを本当に歓迎します.

本当にありがとう。

私は、既知の直径と種のマッピングされた樹木のプロットで、.5m^2 の「ごみトラップ」のネットワークで収集されたごみの量を記述するモデルの適合を最適化しようとしています。選択したモデルには、ごみ生産の相対成長スケーリングと、ごみ移動距離の指数関数的減衰という 2 つの要因があります。

ただし、トラップ データには複数のツリーからの入力が含まれています (これは、タイトルで言及されている「欠落レベル」です)。

これまでのところ、シミュレートされたデータでも非常に複雑な結果が得られました。直径と距離の組み合わせが十分にあれば、関数はある程度十分に制限されるはずです。この分析は、私が模倣している記事で実行されました。ログ（Obs.Litter）の分析も試しましたが、これが道だと思います。しかし、ログ バージョンをコーディングした方法で、パフォーマンスが向上すると期待される結果が得られたかどうかはわかりません。

この時点で、このタイプの「隠れたプロセス」を使用した非線形回帰のフィッティングまたはモデルのフィッティング問題の経験が豊富な人から、あらゆる種類のアドバイス (コードベースまたは概念的) を探しているだけだと思います。データ シミュレーションのコードとさまざまな可能性を以下に示します。OpenBUGS のベイジアン階層モデルを使用してこれらのパラメーターを推定することで、有益な事前情報のみを使用して、もう少し成功しました。

関数を最大化するベクトルp（つまり）の値を見つけたいと思います。 p(1),p(2),p(3),...A(p)

私はそれを行うためにMATLABを使用していますが、fsolveそれが役立つと思われるものを見つけました。だから私は機能を作りましたA：

そして、次のような連立方程式を解く必要があります。

結果のpベクトルが私の問題の解決策になります。

fsolveこの連立方程式 ( ) をどのように解くことができp0=[1 1 1]ますか?

拘束された K3 (K3 は速度定数) を使用した化学動力学モデル (2 組織コンパートメント モデル) があります。

出力特性をプロットするために、化学動力学モデルとともにプラズマ機能をモデル化しました

以下のコードから速度定数 k3 を推定したいと思います

すべてのパラメータの初期推定値を以下に示します

上記の非線形運動モデルコードからK3パラメーターを推定する方法を教えてください

上記の関数では、a1、a2、a3、b1、b2、b3、td、tmax、k1、k2 のパラメータ値は一定のままです。

時間 t 値の変化に伴って K3 値がどのように変化するかを知りたいです。このため、最初に t 間隔 t=0:0.1:60 で K3 を推定したいと思います。

どんな助けでも大歓迎です