blob: 03acf26457088ab334aa20f05098f23cb940cbdb (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
q=poly([0 1 1],'s','coeff');
G=10/q //gain FACTOR=10
H=1
T=G/(1+G*H)
Stg=(1/(1+G*H))
Stg= (-1+ %i*1)/(49+%i*1) //at s= %i*w where w=1
disp(abs(Stg),"sensitivity at w=1 is ")
//sensitivity wrt H
T=-G*H/(1+G*H)
Sth=-50/(49+%i*1) //at s= %i*w where w=1
disp(abs(Sth),"sensitivity wrt H at w=1 is ")
|