Getting error by inverting a matrix

3 次查看（过去 30 天）

显示更早的评论

John Miller 2012-7-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/43930-getting-error-by-inverting-a-matrix

在 MATLAB Online 中打开

Hello,

I changed a little bit the parameters of my matrix J, but now I dont get the inverse of J.

Error message:

Error using mupadmex
Error in MuPAD command: Error: Exponent overflow. [poly]
Evaluating: rnormal
Error in sym/privUnaryOp (line 1700)
            Csym = mupadmex(op,args{1}.s,varargin{:});
Error in sym/inv (line 1128)
        X = privUnaryOp(A, 'symobj::inv');
Error in detofJac (line 10)
Ji=inv(J);

Could it be, that my matrix J now with the paramter changes is singular or is too complicated?

Thank you very much for you help.

7 个评论
显示 5更早的评论隐藏 5更早的评论

John Miller 2012-7-18

it doesnt matter what I intent to do with its invers, I just need it.

And J is a 2x2 matrix: With disp(J); the command window gives me the following: (I just selected and copied from the command window) If I change the parameters in this matrix back to normal , matlab gives me the inverse matrix....

[ (843937847447871284485259215221*atan((7591295679393379*y1)/281474976710656 - 6029284073134799/4503599627370496))/20282409603651670423947251286016 - (843937847447871284485259215221*atan((7591295679393379*y1)/281474976710656 + 63269/100000))/20282409603651670423947251286016 - (2008016153376818114561722882034279487857531019*y1^(2039279169394285/70368744177664))/11150372599265311570767859136324180752990208 - (843937847447871284485259215221*y1*(7591295679393379/(281474976710656*(((7591295679393379*y1)/281474976710656 + 63269/100000)^2 + 1)) - 7591295679393379/(281474976710656*(((7591295679393379*y1)/281474976710656 - 6029284073134799/4503599627370496)^2 + 1))))/20282409603651670423947251286016 - (18521919337695825475900947825254439*y1^(21237/1250))/3094850098213450687247810560000, 0] [ 0, (15798256513974917917454560974767*atan((6339*y2)/250 - 13423/12500))/324518553658426726783156020576256 - (15798256513974917917454560974767*atan((6339*y2)/250 + 2842928789975017/4503599627370496))/324518553658426726783156020576256 - (14746468850402715516729304781749667*y2^(5894254897307163/562949953421312))/158456325028528675187087900672000 - (171587469704128069312579344587508042432114690929*y2^(8162746177111353/281474976710656))/91343852333181432387730302044767688728495783936 - (15798256513974917917454560974767*y2*(6339/(250*(((6339*y2)/250 + 2842928789975017/4503599627370496)^2 + 1)) - 6339/(250*(((6339*y2)/250 - 13423/12500)^2 + 1))))/324518553658426726783156020576256]

Kye Taylor 2012-7-18

在 MATLAB Online 中打开

If j is just two-by-two, why not compute the determinant by hand to see if it equals zero.

If your matrix is

A = [a,b;c,d],

then

det(A) = a*d-b*c

If det(A) ~= 0, the inverse is equal to

Ainv = [d,-b;-c,a] / det(A);

Star Strider 2012-7-18

From the symbolic matrix posted earlier, when I tried to calculate its determinant, I got the same error that calculating the inverse generated, probably for the same reason.

I earlier suggested direct left-matrix-divide (\) and right-matrix-divide (/) to avoid the need to explicitly calculate the inverse, but that seems to have gone nowhere.

请先登录，再进行评论。

请先登录，再回答此问题。