Getting Nan after using expm
3 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
raha ahmadi
el 1 de Jun. de 2021
Comentada: raha ahmadi
el 1 de Jun. de 2021
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 comentarios
Respuesta aceptada
Jan
el 1 de Jun. de 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.
Más respuestas (0)
Ver también
Categorías
Más información sobre Mathematics and Optimization en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!