%-- 10/04/19 12:35 --% s = tf('s') T1 = 1/(s+1) bode(T1) penddemo %-- 10/04/19 13:06 --% s = tf('s') T1 = 1/(s+1) bode(T1) S1 = ss(-1, 1, 1, 0) bode(T1, 'b', S1, 'r--') ctrlpref bode(T1, 'b', S1, 'r--') ctrlpref T3 = 1/(s^2 + 0.1*s + 1) M = 1 K = 1 D = 0.1 A = [0, 1; -K/M, -D/M] B = [0; 1/M] C = [1, 0] S3 = ss(A, B, C, 0) help bode tmp = ss(-1, 0, 0, 1) tmp = ss([], 0, 0, 1) K1 = 1 K2 = 1 M D T4 = 1/(M*s^2 + D*s + K1+K2) %-- 10/04/26 13:12 --% A = [1, 2; 0, 3] eig(A) [V, D] = eig(A) help eig inv(V)*A*V %-- 10/05/10 14:11 --% ls pwd ls ex23 G x t %-- 10/05/17 12:57 --% A = [-1, 10, 3; 1, 5, 3; 4, 13, 7] eig(A) B = [1; 1; 1] C = eye(3) P = ss(A, B, C, 0) x0 = B initial(P, x0) eig(A) help initlal help initial [T, D] = eig(A) T\A*T A = [5, -4; 6, -5] B = [1; 1] p = [-3; -4] F = place(A, B, p) A = [1, 2; -3, -4] B = [2; 3] p = [-2; -3] F = place(A, B, p) Uc = ctrb(A, B) det(Uc)