blob: c290815acc9d08eabc0c014487d926b9fa449ef2 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
|
clc;
A=[-2 1;1 -2];
q0=[.5;1];
//Q=e^(A*t)*q0+integ(e^A*(t-T)*b*x,T,t0,t)
X=0;
//therefore integ term is zero
e=spec(A);//eigen values
b0=e(1)*%e^(e(2)*t)-e(2)*%e^(e(1)*t);
b1=%e^(e(1)*t)-%e^(e(2)*t);
eAt=b0*eye(A)+b1*A;
Q=eAt*q0;
disp(Q,"thus [vc1(t) vc2(t)]=");
|
