<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />皆さんこんにちは 、
JMP で微分方程式を解くことができるかどうか知りたいですか? はいの場合 - その方法を示す例を投稿するか、これを説明するブログ投稿を提供してください。 いいえの場合、誰かが理由を説明できますか?
注: これを実現するために Matlab を起動したくありません。
一番
ウダイ
JSLには派生機能があります。微積分をすべて忘れてしまったので、最後の 3 つの用語を説明することはできませんが、小さなテスト ハーネスを作成したところ、うまくいきました。
// 本当に醜い式を作る
fx = Expr(
A * x ^ W + B * x ^ Z + C * x + D * Sin( E * x ) + F * x * Log( G * x ) + Tan( x ) * Exp( x )
);
// JMP に x に関する微分を求めます。素数の p。
fxp = Derivative( Name Expr( fx ), x );
W * x ^ (W - 1) * A
+ Z * x ^ (Z - 1) * B
+ C
+ E * Cosine( E * x ) * D
+ F *ログ ( G * x ) + G / (G * x) * F * x
+ 1 / コサイン ( x ) ^ 2 * Exp ( x ) + Exp( x ) *日焼け(×)
テスト ハーネス:
A = 3;
B = 4;
C = 5;
D = 6;
E = 7;
F = 8;
G = 9;
W = 10;
Z = 11;
// ある時点で小さなテストを行います, "center"
センター = 7;
// 元の式を評価して小さな間隔で傾きを計算する
間隔 = .0000001;
半間隔 = 間隔 / 2;
x = センター - ハーフインターバル;
y1 = fx; // 中央の下側
x = センター + ハーフインターバル;
y2 = fx; // 中央のハイサイド
傾き = (y2 - y1) / 間隔; // 導関数の近似
Write( "\!nslope: " || Char( 勾配 ) );
x = センター; // 中点
yd = fxp; // 正確な導関数
Write( "\!nderivative: " || Char( yd ) );
Write( "\!n% エラー: " || Char( 100 * (勾配 - ヤード) / ((勾配 + ヤード) / 2) ) );
スロープ: 13639522056.5796
派生: 13639522109.8688
パーセント エラー: -0.000000390696846545
かなり急です！ゼロにはるかに近い中心は傾斜が緩くなっています。
導関数を計算する方法、または少なくとも一次導関数の式を保存する方法は少なくとも 2 つあります。
最初の方法として、数式を保存するか、数式列を作成する必要があります。導関数を取得するには、目的の列に移動し、列ヘッダーを右クリックして [数式] をクリックします。式エディタが開くと、式が表示されます。導関数を取得する変数を選択する必要があります。
2 番目の方法は、非線形プラットフォームです。再びモデルを作成する必要がありますが、予測式を保存する代わりに、一次導関数の式を保存します。
非線形プラットフォーム モデルのダイアログで、2 次導関数や数値のみの導関数を指定することもできます。
HTH
明細書
JSL Scripting Guide にセクションがあります。
ヘルプ==>書籍==>スクリプト ガイド
デリバティブの計算を扱います。JMP 12 Scripting Guideでは、265ページから始まります。
