diff options
author | prashantsinalkar | 2018-02-03 11:01:52 +0530 |
---|---|---|
committer | prashantsinalkar | 2018-02-03 11:01:52 +0530 |
commit | 7bc77cb1ed33745c720952c92b3b2747c5cbf2df (patch) | |
tree | 449d555969bfd7befe906877abab098c6e63a0e8 /3885/CH6/EX6.14/Ex6_14.sci | |
parent | d1e070fe2d77c8e7f6ba4b0c57b1b42e26349059 (diff) | |
download | Scilab-TBC-Uploads-7bc77cb1ed33745c720952c92b3b2747c5cbf2df.tar.gz Scilab-TBC-Uploads-7bc77cb1ed33745c720952c92b3b2747c5cbf2df.tar.bz2 Scilab-TBC-Uploads-7bc77cb1ed33745c720952c92b3b2747c5cbf2df.zip |
Diffstat (limited to '3885/CH6/EX6.14/Ex6_14.sci')
-rw-r--r-- | 3885/CH6/EX6.14/Ex6_14.sci | 34 |
1 files changed, 34 insertions, 0 deletions
diff --git a/3885/CH6/EX6.14/Ex6_14.sci b/3885/CH6/EX6.14/Ex6_14.sci new file mode 100644 index 000000000..e1257acc2 --- /dev/null +++ b/3885/CH6/EX6.14/Ex6_14.sci @@ -0,0 +1,34 @@ +//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 6.14
+
+clc;
+clear;
+s=poly(0,'s')
+//dominent pole sd=-zeta*w=%i*w*sqrt(1-zeta^2)
+zeta=0.8//damping ratio
+w=2//natural frequency of osciilation in rad/sec
+sd=(-zeta*w)+((%i*w)*sqrt(1-zeta^2))
+disp(sd,'the dominennt pole is')
+d=abs(sd)
+disp(d,'the value of d is ')
+betaa=phasemag(sd)
+disp(betaa,'the value of betaa is;')
+h=syslin('c',20/s*(s+2)*(s+4))//given tranfer function G(s)
+//find magnitude and phase of G(s) at s=sd
+a=20/(sd*(2+sd)*(4+sd))
+ad=abs(a)
+disp(ad,'the value of ad is')
+phid=phasemag(a)
+disp(phid,'the value of phid is')
+kd=sind(phid)/(d*ad*sind(betaa))//derivative constant
+disp(kd,'the derivative constant is')
+kp=(-sind(betaa+phid))/(ad*sind(betaa))//proportional constant
+disp(kp,'the integral constant is')
+hc=syslin('c', s*((0.243/s)+0.557))//transfer function of PD controller is kpof +kd*s
+disp(hc,'transfer function of PD controller is')
+hcmp=syslin('c',h*hc)//transfer function compensated system
+disp(hcmp,'transfer function compensated system ')
|