ch3_34.m

来自「matlab7.0从入门到精通(光盘内容)」· M 代码 · 共 16 行

A = [0    -6    -1;
     6     2   -16;
    -5    20   -10];                 %矩阵A
x0=[1;1;1];                           %初值x(0)
t=0:0.01:1;                           %时间t
xt=[];                               %计算结果x(t)
for i=1:length(t),                      %计算每个时间点上的x(t)
    xt(i,:)=expm(t(i)*A)*x0;
end;
plot3(xt(:,1),xt(:,2),xt(:,3),'-o')           %画三维图显示x(t)
grid on;                             %给图添加辅助网格
xlabel('xt(:,1)');                       %给图的坐标轴x命名
ylabel('xt(:,2)');                       %给图的坐标轴y命名
zlabel('xt(:,3)');                       %给图的坐标轴z命名
title('矩阵指数函数结果');             %给图添加标题