問題タブ [mupad]

forループでMatlabの2つの非線形方程式のシステムを解こうとしています。for ループが原因でプログラムの実行が遅くなる可能性があることは認識していますが、現在、1 つのインスタンスを解決するのに約 1 時間かかります。コードをデバッグすると、次の行がプログラムをハングさせることに気付きました。

プログラムを終了すると、次のようになります。

ここで、97 行目は同じコード行を参照しています。これは私がソルブ関数を呼び出す私のコードのセクションです:

私は同次座標で作業しており、muPad で記号計算を行うのが好きです。私の変換行列Tが 4x4 で、左上がRと呼ばれる 3x3 回転行列、右上がbと呼ばれる 3x1 ベクトル、右下の要素が1であるとしましょう。

したがって、次のようになります。

通常の Matlab 表記で。

ここで、muPad を展開して要素を操作するのではなく、シンボルRとbを操作する必要があります。

したがって、2 番目の変換T2を次のように定義すると、

T と T2 の積は次のように表示されます。

私はここで初めてなので、フォーマットが悪いことを許してください..（ラテックスを含める方法はありますか？）

私はコードを持っています：

Mupad が出力 e(s) を与えるのはなぜですか? Mupad 出力:

G(s)*exp(1)(s) (または G(s)*e(s) を青色で表記)

E(s) が Mupad によって自動的に変更されないようにするにはどうすればよいですか ...正直なところ、これは遅れています... ボーナス :(

....うわあ

E(s):=R(s)-B(s) エラー: 識別子 'E' は保護されています。[_割当]

私はu0ステップ関数として持っています（区分定義）：

今私は計算します：

を与える

教科書で私は答えを見つけました：

これは MuPAD のバグですか?

MuPAD では次のように単純化します。

与える：

そんなこと知ってる

so f(t)=2*t, so は Soをdirac(t+1)=dirac(t-(-1))意味し、MuPAD が返すものとは異なります。k=-12 * t * dirac(t+1) =2*(-1)*dirac(t) =-2*dirac(t)-2*dirac(t + 1)

私の間違いはどこですか？

x1 x2 m1 と m2 を解くために mupad で numeric::fsolve を使用しています。使用されるコマンドは次のとおりです。

Scase2:=numeric::fsolve([m1,g2,eqn1,eqn2],m1=0,m2=0..infinity,x1=-3.14..3.14,x2=-3.14..3.14],RestrictedSearch,MultiSolutions)

返される答えは次のとおりです。

[m1 = 0.0, m2 = -8.678897929, x1 = 3.38742203*10^(-17), x2 = 2.77006362], [m1 = 0.0, m2 = 19705.59085, x1 = -1486.62367, x2 = 157.7356704], [m1 = 0.0, m2 = -6151.321562, x1 = -6.594711892*10^(-18), x2 = 78.32275445], [m1 = 0.0, m2 = -2334.614422, x1 = -1.047226205*10^(-13), x2 = -48.81668048], [m1 = 0.0, m2 = 38.00363756, x1 = -70.90099452, x2 = 11.28424375], [m1 = 0.0, m2 = -3017.254857, x1 = -5.285351078*10^(-17), x2 = 55.38965896], [m1 = 0.0, m2 = -14353.49312, x1 = -3.299575199*10^(-16), x2 = 119.665553], [m1 = 0.0, m2 = 4746.862416, x1 = 216.8709045, x2 = 69.03215292], [m1 = 0.0, m2 = -337.6953647, x1 = 4.231726504*10^(-19), x2 = 18.53286231], [m1 = 0.0, m2 = 59888.45924, x1 = 769.621396, x2 = -244.978099], [m1 = 0.0, m2 = -754.4068652, x1 = 7.46136118*10^(-19), x2 = 26.99119714], [m1 = 0.0, m2 = -72400.20717, x1 = -1.43883627*10^(-12), x2 = -269.417393], [m1 = 0.0, m2 = -2953.978962, x1 = -3.936403445*10^(-22), x2 = -54.84860971], [m1 = 0.0, m2 = -0.6449920855, x1 = -2.627245131*10^(-15), x2 = -0.5757632436], [m1 = 0.0, m2 = -1646.776257, x1 = 3.55576955*10^(-15), x2 = -40.34209161], [m1 = 0.0, m2 = -52.36064061, x1 = -5.846045816*10^(-13), x2 = -6.937973456], [m1 = 0.0, m2 = 31302.17624, x1 = 707.7689947, x2 = 225.2898681], [m1 = 0.0, m2 = -6184.407888, x1 = 400.9347135, x2 = 127.621483], [m1 = 0.0, m2 = -30523.06702, x1 = 6.438047089*10^(-17), x2 = -175.083226], [m1 = 0.0, m2 = -212.5775725, x1 = 8.213227391*10^(-19), x2 = -14.08916721]

ただし、ソルバーによって返される回答には、指定された範囲外のソリューションが含まれています。まあ、手動で許容できる解決策を見つけることができます。これを行うためのより良い方法はありますか? ソルバーが許容範囲外の解を返すのはなぜですか?