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

---
 773/CH9/EX9.07/9_07.sci | 32 ++++++++++++++++++++++++++++++++
 1 file changed, 32 insertions(+)
 create mode 100755 773/CH9/EX9.07/9_07.sci

(limited to '773/CH9/EX9.07/9_07.sci')

diff --git a/773/CH9/EX9.07/9_07.sci b/773/CH9/EX9.07/9_07.sci
new file mode 100755
index 000000000..bedd58bd9
--- /dev/null
+++ b/773/CH9/EX9.07/9_07.sci
@@ -0,0 +1,32 @@
+//calculates//
+s=%s;
+sys1=syslin('c',20/(s*(s+2)));
+syms Kt;
+sys2=Kt*s;
+sys3=sys1/.sys2;
+p=simple(sys3);
+disp(p,"G(s)=")
+H=1;
+sys=sys3/.H;
+sys=simple(sys);
+disp(sys,"C(s)/R(s)=")
+[num,den]=numden(sys)
+cof_a_0 = coeffs(den,'s',0) // coeff of den of symbolic function CL
+cof_a_1 = coeffs(den,'s',1)
+//Wn^2= cof_a_0,comparing the coefficients
+Wn=sqrt(cof_a_0) 
+Wn=dbl(Wn);
+disp(Wn,"natural frequency Wn=")           // Wn=natural frequency
+//cof_a_1=2*zeta*Wn
+zeta=cof_a_1/(2*Wn)
+zeta=0.6;   
+Kt=((2*zeta*Wn)-2)/20;
+disp(Kt,"Kt=")
+Wd=Wn*sqrt(1-zeta^2);
+disp(Wd,"Wd=")
+Tp=%pi/Wd;
+disp(Tp,"Tp=")
+Mp=100*exp((-%pi*zeta)/sqrt(1-zeta^2));
+disp(Mp,"Mp=")
+Ts=4/(zeta*Wn);
+disp(Ts,"Ts=")
-- 
cgit