matlab警告imaginary,MATLAB上机作业西安电子科技大学

第三次上机作业

班级：07127020 学号：07127020 姓名：方正阳 符号运算：

1.求符号矩阵

111213

212223

313233

a a a

A a a a

a a a

??

??

=??

??

??

的行列式值和逆，所得结果应采用

“子表达式置换”简洁化。

>> clear

>> A=sym('[a11,a12,a13;a21,a22,a23;a31,a32,a33]')

A =

[ a11, a12, a13]

[ a21, a22, a23]

[ a31, a32, a33]

>> DA=det(A);

>> IA=inv(A);

>> x_DA=subexpr(DA,'D')

x_DA =

a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31

>> x_IA=subexpr(IA,'D')

D =

1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)