blob: 6e2ab4348ffd4e93e63a7156b3191dd0696075e3 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
//caption:obtain_zero_state_response
//example 9.10.12_(ii)
//page 395
s=%s;
syms t
A=[0 1;-2 -1]
B=[0 1]'
x0=[1 0]'
[r c]=size(A);//size of matrix A
p=s*eye(r,c)-A;//s*I-A where I is identity matrix
q=det(p)//determinant of sI-A
r=inv(p)//inverse of sI-A
m=r*B*(1/s);
//for calculating zero state response
x=[0;0]
x(1,1)=ilaplace(m(1,1),s,t);
x(2,1)=ilaplace(m(2,1),s,t);
disp(x,"zero input response of the system,x(t)=");
|
