How to calculate expm of a badly scaled matrix?

Question

alex corno 2013-6-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/77906-how-to-calculate-expm-of-a-badly-scaled-matrix

I'm using expm for a matrix with big and small coefficients and I get NaN in the result. Is there a method to get the good matrix exponential even if the matrix is badly scaled? Normally there is a exponential for every matrix.

ALEX

2 个评论
显示无隐藏无

Matt J 2013-6-4

编辑：Matt J 2013-6-4

在 MATLAB Online 中打开

Suppose you had a simple badly scaled matrix like

     >> A=diag( [1e6,1] )
    A =
         1000000           0
               0           1

What would you like the result of expm(A) to be? The following has no NaNs. Would it be an okay answer?

    >> expmA=diag( exp([1e6,1])  )
    expmA =
           Inf         0
             0    2.7183

alex corno 2013-6-4

I think I have a different problem. My matrix isn't even invertible, there are two zero lines, moreover there are some coefficients that are of order of 1e6. In my opinion at some point Matlab is doing some inversion that's why expm doesn't work properly with badly scaled matrices. In my problem I changed the units of physical parameters (distance -> µm Stresses -> kPa) and now it works. But I'm looking for some existent function (transformation) which would enable to calculate expm for any matrix (apart for huge numbers that shouldn't appear in a physical problem there's no chance that e^100 will be a coefficient in a transfer matrix)

Alexandru

请先登录，再进行评论。

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

Answer 1

Richard Brown 2013-6-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/77906-how-to-calculate-expm-of-a-badly-scaled-matrix#answer_87726

You might this Cleve's corner instructive http://blogs.mathworks.com/cleve/2012/07/23/a-balancing-act-for-the-matrix-exponential/

In particular, have you tried the alternative versions expmdemo1 and expmdemo3 ?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Matt J 2013-6-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/77906-how-to-calculate-expm-of-a-badly-scaled-matrix#answer_87662

编辑：Matt J 2013-6-4

在 MATLAB Online 中打开

I think I have a different problem. My matrix isn't even invertible, there are two zero lines, moreover there are some coefficients that are of order of 1e6. In my opinion at some point Matlab is doing some inversion that's why expm doesn't work properly with badly scaled matrices.

No, it's not an invertibility issue. The reason for the NaNs is the big 1e6 order numbers. Below is an example of a non-invertible matrix on which EXPM works fine.

     >> expm(diag([1,0]))
    ans =
        2.7183         0
             0    1.0000

3 个评论
显示 1更早的评论隐藏 1更早的评论

alex corno 2013-6-4

Yes, OK but is there a way in Matlab(a transformation or something) to get expm even if the matrix is badly scaled? In my case I just changed the units of my unknowns but is there a better method?

Matt J 2013-6-4

编辑：Matt J 2013-6-4

在 MATLAB Online 中打开

Even if the matrix is well-scaled, you could have the same problem, for example:

    >> A=diag([1e6,1e6]); %a well-scaled matrix
    >> expm(A)
    ans =
       NaN   NaN
       NaN   NaN

The bottom line is, you must avoid matrices with large eigenvalues, and then things should be fine. Otherwise, they will not be fine.

There are any number of transformations you could consider to reduce the magnitude of the eigenvalues. Dividing the matrix by a large number is one way, but I can't know if that suits your application.

请先登录，再进行评论。