From 7bc77cb1ed33745c720952c92b3b2747c5cbf2df Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Sat, 3 Feb 2018 11:01:52 +0530
Subject: Added new code

---
 3885/CH6/EX6.2/Ex6_2.sce | 45 +++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 45 insertions(+)
 create mode 100644 3885/CH6/EX6.2/Ex6_2.sce

(limited to '3885/CH6/EX6.2/Ex6_2.sce')

diff --git a/3885/CH6/EX6.2/Ex6_2.sce b/3885/CH6/EX6.2/Ex6_2.sce
new file mode 100644
index 000000000..a413259da
--- /dev/null
+++ b/3885/CH6/EX6.2/Ex6_2.sce
@@ -0,0 +1,45 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 6.2
+
+clc;
+clear; 
+s=poly(0,'s')
+//calculation of gain K
+kv=30//velocity error constant is 30sec^-1
+// open loop transfer function G(s)=K/s*(s+4)*(s+80)
+//by definition of velocity error constant applying limit s=0 in G(s) 
+disp('the value of K is 9600;')
+h=syslin('c',9600/(s*(s+4)*(s+80)))
+bode(h)
+show_margins(h) 
+xtitle("uncompensated system")
+//from the plot the phase margin of uncompensated system is 12
+//but the system requires phase margin of 33 so lag compensation required
+pm=38//choose PM of compensated system is 38 degree
+phigcn=38-180// phase of G(s) at new gain cross over frequency
+wgcn=4.7//the frequency corrosponding to  phase of -142 is 4.7 rad/sec
+agcn=16//db magnitude at wgcn
+//20log betaa=16
+betaa=10^(agcn/20)
+disp(betaa,'the value for betaa is')
+zc=wgcn/10//zero of lag compensator
+t=10/wgcn 
+disp(t,'the value for t is')
+pc=1/(betaa*t)
+disp(pc,'pole of lag compensator is')
+//transfer function of lag compensation is (s+1/t)/(s+1/betaa*t))
+hc=syslin('c', (6.3*(1+2.13*s))/(1+13.419*s))
+disp(hc,'the transfer function of lag compensator is;') 
+//open loop transfer function of compensated system is h*hc
+hcmp=syslin('c',h*hc)
+disp(hcmp,'open loop transfer function of compensated system is') 
+figure()
+bode(hcmp)
+show_margins(hcmp)
+xtitle("compensated system")
+
+
-- 
cgit