程序读不懂谁帮读读

yhq966 · 发表于 2008-11-4 23:23

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

clear
close
A=[[0.2 0 1 0];[0 0.1 0 1];[1 0 0 0];[0 1 0 0]];
B=[[2 -1 0 0];[-1 2 0 0];[0 0 -1 0];[0 0 0 -1]];
[V,D]=eig(-B,A);
init=[[1];[1];[0];[0]];
c=inv(V)*init;

for i=1:1:600
t=(i-1)*0.01;x1=0;x2=0;
for j=1:1:4
x1=x1+c(j)*V(1,j)*exp(D(j,j)*t);
x2=x2+c(j)*V(2,j)*exp(D(j,j)*t);
end
xt(i,1)=real(x1);xt(i,2)=real(x2);xt(i,3)=0;
end
plot(xt)

账号		自动登录	找回密码
密码			我要加入

[编程技巧] 程序读不懂 谁帮读读

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

[编程技巧] 程序读不懂谁帮读读