Getting Nan after using expm
7 次查看(过去 30 天)
显示 更早的评论
Hi H dont know why I get Nan in situatin A but not in B.
Thanks in advance
% situation A
r=[1e5,0,0;
0,-2e5,0;
1,1,-5e5]
R=expm(r);
%% situation B
g=[-1e5,0,0;
0,-2e5,0;
1,1,-5e5]
G=expm(g);
0 个评论
采纳的回答
Jan
2021-6-1
See the help text of expm:
Although it is not computed this way, if A has a full set of eigenvectors V with corresponding eigenvalues D then
% [V,D] = EIG(A) and expm(A) = V*diag(exp(diag(D)))/V
r=[1e5,0,0; 0,-2e5,0; 1,1,-5e5];
exp(diag(D))
% [0; Inf; 0]
This Inf is the cause of the NaN output.
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Mathematics and Optimization 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!