問題タブ [derivative]

Rを使用してパラメーターをデータ セットに適合させようとしていますoptim()。目的関数では、方程式の反復根解法が必要であるGため、予測値pがG(目的関数内にネストされた) の値を 0 (または 0 に近い値) にします。可能な; 安定性のために二分法を 50 回繰り返します)。

ここに問題があります: に解析的勾配を含めたいのoptim()ですが、反復関数にはそれができないと思います。しかし、分析的勾配をあきらめる前に、ここにいる全員でこの問題を実行して、見落としている解決策があるかどうかを確認したいと思いました. 何かご意見は？

注: 二分法に落ち着く前に、他のルート解決法を試しましたが、ブラケットを使用しない方法 (ニュートンなど) はすべて不安定なようです。

以下は、問題の再現可能な例です。提供されたデータ セットと の開始値を使用するoptim()と、アルゴリズムは解析的勾配なしで問題なく収束しますが、他のデータ セットと開始値ではうまく機能しません。

ありがとう！！

私はPythonで心電図信号に取り組んでおり、セマンティックマイニングから抽出された自然周波数、減衰係数、入力信号などのパラメータから心電図信号の振動動作をキャプチャしたいと考えています。このアルゴリズムは、x'、x''、x'''、および x'''' として示される単一点での 4 次導関数を使用します。ここで、x はフーリエ変換された心電図信号です。私の質問は、周波数領域で信号の時間導関数をどのように取得するかです。これを行うには、scipy.misc.derivative や numpy.diff のように Python でどの関数を使用すればよいですか?

私は Mathematica を初めて使用し、現在クォータニオンの作業に使用しています。ある特定のケースでは、次のようなクォータニオンのノルムを含む導関数を取得したいと考えています。

Mathematica は実際にクォータニオンのノルムを計算することができますが、上記の入力が展開された導関数を生成するようにノルム関数を展開できるかどうかを誰かが知っていますか (//FunctionExpand は出力で何も変更しませんでした...)。

ヒントをいただければ幸いです。

R では、相互作用項、二乗項などを含む式の右辺の記号導関数を取得する方法が必要です。

たとえば、次の 2 つの [edit:three] 式の右辺を x について微分できるようにしたいと考えています。

y~x+I(x^2)

y~x:z

編集:y~x*z

上記の 3 つの式を入力すると、それぞれ 1+2x、z、1+z を返す関数が欲しいです。

私は次のことを試しました：

R に I(x^2) (または同様に I(x*z) など) を x^2 (それぞれ x*z)、x:z を x*z (導関数を計算する目的で、x*z (式の意味で) x+z+x*z (数学の意味で) として?

第二に、deriv() から出力を取得し、それを数式の右辺のように整形する方法はありますか? 特に、D() がこの問題を軽減し、希望する形式で出力を生成することはわかっています (ただし、D() は数式を入力として処理できません)。導関数を取得したい変数ごとに D() を何度も適用することでこれを回避できますが、そのようなすべての変数の文字列を単純に入力し、配置に適した出力を受け取るとよいでしょう。式の右辺。

ありがとうございました！

Rで生存関数の導関数を取り、それを新しい関数に保存しようとしています。

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

ポイントx = 1.0でMatlabの次の関数の導関数を計算したいと思います:

私はこれを試しましたが、py = []を返します

こんにちは、C++ で動きの小さなシミュレーションを作成しました。オイラー法、ルンゲクッタ法、ミッドポイント法との違いを学生に示したかったのです。

しかし、Rungy-Kuttaモードに切り替えると、どこかでミスをして、素材ポイントがすべて消えてしまいました。私が犯した間違いは、solveRK4 メソッドにあります。添付ファイルの MinGW Developer Studio プロジェクト。mingw コンパイラ ディレクトリに配置するライブラリ フォルダもあります。: http://speedy.sh/CvDHj/LABO3.zip

'R' ボタンを押すと、'E' が Euler に、'M' が MidPoint に、RK4 に切り替わります。

問題は、solveRK4 関数のどこで間違いを犯したかです。

結果 (リリース) は次のようになります: http://speedy.sh/h28VP/zad3.exe R を押すと、ポイントが画面から消えます。

Punkt - Point Wektor - ベクトル

メイン：

構造:

数学を扱うのは久しぶりで、少しさびています。ばかげた質問をしてもいいです。

問題: 2D 点の配列としてメモリに保存されている n 対の線があるため、明示的な関数はありません。対の線が平行かどうかを確認する必要があります。導関数が同じかどうかを確認するだけで十分なので、これは非常に簡単な作業です。

アルゴリズムでこれを行うには、(私が持っている) 関数の 2 点間の直線の傾きを確認する必要があります。極端な精度は必要ないため、次の簡単な式を使用できます。

m = (y2-y1)/(x2-x1)

しかし明らかに、これは x2 = x1 という大きな問題につながります。この場合のデフォルト値を指定できません...どうすれば回避できますか?