From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 29/CH11/EX11.8/exa11_8.sce | 28 ++++++++++++++++++++++++++++
 1 file changed, 28 insertions(+)
 create mode 100755 29/CH11/EX11.8/exa11_8.sce

(limited to '29/CH11/EX11.8/exa11_8.sce')

diff --git a/29/CH11/EX11.8/exa11_8.sce b/29/CH11/EX11.8/exa11_8.sce
new file mode 100755
index 000000000..2bd7e7fef
--- /dev/null
+++ b/29/CH11/EX11.8/exa11_8.sce
@@ -0,0 +1,28 @@
+//Caption:time_response_for_unit_step_function
+//example 11_8
+//page 476
+s=%s;
+syms t;
+G=(20)/((s+1)*(s+5))
+H=1;
+CL=G/.H
+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
+Mp=(exp(-(zeta*%pi)/sqrt(1-zeta^2)))*100//max. overshoot
+Tp=%pi/Wd//peak time
+t=(2*%pi)/(Wn*sqrt(1-zeta^2))//period of oscillation
+ts=4/(zeta*Wn)//settling time
+N=Wd/(2*%pi)*ts//no. of oscillations completed before reaching steady state
+disp(Tp,"peak time=");
+disp(Mp,"max. peak overshoot=");
+disp(t,"period of oscillation");
+disp(N,"no. of oscillations completed before reaching steady state=");
+disp(ts,"settling time=")
-- 
cgit