A = [0    -6    -1;
     6     2   -16;
    -5    20   -10];                 %矩陣A
x0=[1;1;1];                           %初值x(0)
t=0:0.01:1;                           %時(shí)間t
xt=[];                               %計(jì)算結(jié)果x(t)
for i=1:length(t),                      %計(jì)算每個(gè)時(shí)間點(diǎn)上的x(t)
    xt(i,:)=expm(t(i)*A)*x0;
end;
plot3(xt(:,1),xt(:,2),xt(:,3),'-o')           %畫三維圖顯示x(t)
grid on;                             %給圖添加輔助網(wǎng)格
xlabel('xt(:,1)');                       %給圖的坐標(biāo)軸x命名
ylabel('xt(:,2)');                       %給圖的坐標(biāo)軸y命名
zlabel('xt(:,3)');                       %給圖的坐標(biāo)軸z命名
title('矩陣指數(shù)函數(shù)結(jié)果');             %給圖添加標(biāo)題