exa2_31.m

来自「有关matlab r2007的实用教程」· M 代码 · 共 10 行

clear all
a=[-1,0,0;0,-2,0;0,0,-3]
x0=[1;1;1];
t=0:0.1:1.0;
for i=1:length(t)
    x(:,i)=expm(a*t(i))*x0;
end
plot3(x(1,:),x(2,:),x(3,:));
grid on