Getting Nan after using expm
4 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
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 Commenti
Risposta accettata
Jan
il 1 Giu 2021
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.
Più risposte (0)
Vedere anche
Categorie
Scopri di più su Mathematics and Optimization in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!