I have matrix A and I need to find (e^(At)) where t is the sampling time. How to find that? Also what is the difference between exp(A) and expm(A)?
15 次查看(过去 30 天)
显示 更早的评论
How to find (e^(At)) that? Also what is the difference between exp(A) and expm(A)?
0 个评论
采纳的回答
James Tursa
2021-10-20
编辑:James Tursa
2021-10-20
exp(A) just takes the exp( ) of the scalar elements individually. I.e., exp(A(1,1)), exp(A(1,2)), etc.
expm(A) takes the matrix exponential. I.e., conceptually the equivalent of eye(n) + A + A^2/2! + A^3/3! + ...
These are two very different operations. From your description it looks like you should be using expm( ).
更多回答(2 个)
KSSV
2021-10-20
A = rand(2) ;
t = linspace(0,1) ;
iwant = zeros(2,2,100) ;
for i = 1:100
iwant(:,:,i) = exp(A*(t(i))) ;
end
exp might give less accurate values for [0,1]; in that case expm will be accurate
0 个评论
Bruno Luong
2021-10-20
编辑:Bruno Luong
2021-10-20
Math notation e^(At) is computed in MATLAB with
expm(A*t)
or equivalently
expm(A)^t
0 个评论
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!