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

---
 1151/CH4/EX4.16/example16.sce | 22 ++++++++++++++++++++++
 1 file changed, 22 insertions(+)
 create mode 100755 1151/CH4/EX4.16/example16.sce

(limited to '1151/CH4/EX4.16')

diff --git a/1151/CH4/EX4.16/example16.sce b/1151/CH4/EX4.16/example16.sce
new file mode 100755
index 000000000..29b86c783
--- /dev/null
+++ b/1151/CH4/EX4.16/example16.sce
@@ -0,0 +1,22 @@
+ // bode p l o t is given in the question we have to find transfer function
+s=%s;
+sys = syslin('c',48/(s*(s +20)*(s^2+2.52*s +16)))
+bode(sys)
+show_margins (sys ,'bode' )
+gm= g_margin (sys)
+pm= p_margin (sys)
+disp (gm ," g a in margin=")
+disp (pm ," phas e margin=")
+if (gm <=0 | pm <=0)
+printf (" system is unstable ")
+else
+printf (" system is stable ")
+end
+printf("a) At w=1 rad/sec slope changes -20 db/decade indicating a term K/s \n");
+printf("b) At w=4 rad/sec slope changes from -20 db/decade to -60 db/decade indicating a term (1+2d/w*s+s^2/w^2)in denominator(net change -40 bd/decade)\n ");
+printf("c) At w=10 rad/sec slope changes from -60 db/decade to- 40b/decade indicating a term (1+s/10)in numerator\n");
+printf("d) At w=20 rad/sec slope changes from --40 db/decade to- -60b/decade indicating a term (1+s/20)in denominator\n");
+printf("e) its also mention 20 logK=3.521")
+k=1.5;
+disp(K,"K=");
+printf("the given transfer function is G(s)=48/(s*(s +20)*(s^2+2.52*s +16))");
-- 
cgit