Expm(x):求矩陣x的指數,這個矩陣x只能是方陣;在matlab中,采用pade近似算法進行計算。
舉壹個matlab中的例子,看看:
x =[1 2 3;4 5 6];%計算壹個2行3列的矩陣的exp,得到每個元素的exp。
& gt& gtexp(x)
ans =
2.7183 7.3891 ?20.0855
54.5982 ?148.4132 ?403.4288
& gt& gtexpm(x)?%如果使用這個函數,會得到壹個錯誤,因為此時X不是方陣。
使用== >時出錯;時光
內部矩陣尺寸必須壹致。
== >中的錯誤;expm & gtPadeApproximantOfDegree at 118
A2 = A * AA4 = A2 * A2A6 = A2 * A4
== >中的錯誤;在39時爆炸
f = PadeApproximantOfDegree(m _ vals(end));
& gt& gtexpm([1 2 3;4 5 6 ;7 8 9])%此時expm()的內部是壹個3行3列的方陣,說明結果是正確的。
ans =
1.0e+006 *
1.1189 1.3748 1.6307
2.5339 3.1134 3.6929
3.9489 4.8520 5.7552