blob: ee1900a1bc391075c560bb305c44d9854c64f88b (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
//Caption:time_response_for_unit_step_function
//example 11_7
//page 475
s=%s;
syms t;
CL=(s+2)/(s^2+2*s+2)
disp(CL,"C(s)/R(s)=")
y=ilaplace(CL,s,t);
disp(y,"c(t)=");
b=denom(CL)
c=coeff(b)
//Wn^2=c(1,1)
Wn=sqrt(c(1,1))//natural frequency
//2*zeta*Wn=c(1,2)
zeta=c(1,2)/(2*Wn)//damping ratio
Wd=Wn*sqrt(1-zeta^2)//damped frequency
Tp=%pi/Wd//peak time
Mp=(exp(-(zeta*%pi)/sqrt(1-zeta^2)))*100//max. overshoot
disp(Tp,"peak time=");
disp(Mp,"max. peak overshoot=")
|
