summaryrefslogtreecommitdiff
path: root/3889/CH7/EX7.4/Ex7_4.sce
diff options
context:
space:
mode:
Diffstat (limited to '3889/CH7/EX7.4/Ex7_4.sce')
-rw-r--r--3889/CH7/EX7.4/Ex7_4.sce32
1 files changed, 32 insertions, 0 deletions
diff --git a/3889/CH7/EX7.4/Ex7_4.sce b/3889/CH7/EX7.4/Ex7_4.sce
new file mode 100644
index 000000000..260d60321
--- /dev/null
+++ b/3889/CH7/EX7.4/Ex7_4.sce
@@ -0,0 +1,32 @@
+//Example 7.4
+//page 436
+//Control Systems: Principles and Design
+//M Gopal, Second Edition, Tata McGraw-Hill
+//Chapter:Compensator Design Using Root Locus
+
+xdel(winsid())//close all graphics Windows
+clear;
+clc;
+//transfer function
+s=%s;
+P=1/((s)*(s+3)*(s^2+2*s+2));
+
+//Root locus plot using evans root locus
+f=figure()
+evans(P)
+title("Root locus of 1/((s)*(s+3)*(s^2+2*s+2))",'fontsize',5)
+h=legend('');
+h.visible="OFF"
+zoom_rect([-2 -1.5 2 1.5]*kpure(P)/3)
+a=gca();
+a.x_location="origin"
+a.y_location="origin"
+a.parent.background=8;
+legends(['root locus';'';'asymptotic directions';'open loop poles'],[2,3,1,-2],with_box=%f,opt="ur",2.8)
+sgrid();
+K=poly(0,'K')
+R= routh_t(P,K);
+disp(R,'Routh Array=')
+kval= kpure(P)
+disp(kval,'Limiting Gain Kmax=')
+