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/CH11/EX11.1/example1.sce | 23 +++++++++++++++++++++++ 1151/CH11/EX11.2/example2.sce | 20 ++++++++++++++++++++ 1151/CH11/EX11.3/example3.sce | 27 +++++++++++++++++++++++++++ 1151/CH11/EX11.4/example4.sce | 15 +++++++++++++++ 1151/CH11/EX11.5/example5.sce | 15 +++++++++++++++ 1151/CH11/EX11.6/example6.sce | 16 ++++++++++++++++ 1151/CH11/EX11.7/example7.sce | 16 ++++++++++++++++ 7 files changed, 132 insertions(+) create mode 100755 1151/CH11/EX11.1/example1.sce create mode 100755 1151/CH11/EX11.2/example2.sce create mode 100755 1151/CH11/EX11.3/example3.sce create mode 100755 1151/CH11/EX11.4/example4.sce create mode 100755 1151/CH11/EX11.5/example5.sce create mode 100755 1151/CH11/EX11.6/example6.sce create mode 100755 1151/CH11/EX11.7/example7.sce (limited to '1151/CH11') diff --git a/1151/CH11/EX11.1/example1.sce b/1151/CH11/EX11.1/example1.sce new file mode 100755 index 000000000..7a23cec31 --- /dev/null +++ b/1151/CH11/EX11.1/example1.sce @@ -0,0 +1,23 @@ +printf(" Given C(s)/R(s)=14/(s^2+1.4*s+14)") +printf("characterstic equation of the given system is s^2+1.4*s+14=0"); +printf("compare it with the standard second order characterstic equation s^2+2*d*w*s+w^2=0"); +w=sqrt(14); +d=1.4/(2*w); +d1=.7; +t=2*(d1-d)/w; +pt1=%pi/(w*sqrt(1-d^2)); +mo1=exp((-%pi*d)/sqrt(1-d^2))*100; +rt=(%pi-atan(sqrt((1-d)/d)))/(w*sqrt(1-d^2)); +disp(t,"Td=") +disp(pt1,"peak time(in sec)for sytem without derivative control is") +disp(rt,"rise time(in sec)for sytem without derivative control is") +disp(mo1,"maximum overshoot(in %)for sytem without derivative control is") +printf("overal transfer function with derivative control is C(s)/R(s)=14*(1+0.274*s)/(s^2+5.236*s+14)") +printf("c(t)=1-e^(-2.618*t)*cos(2.673*t)+0.455*e^(-2.618t)*sin(2.673t)") +tp=(%pi-atan(2.58))/2.673; +tr=(1/2.673)*atan(1/0.455); +mo2=1-%e^(-2.618*tp)*cos(2.673*tp)+0.455*%e^(-2.618*tp)*sin(2.673*tp); +mp=(mo2-1)*100; +disp(tp,"peak time(in sec)for sytem with derivative control is") +disp(tr,"rise time(in sec)for sytem with derivative control is") +disp(mp,"maximum overshoot(in %)for sytem with derivative control is") diff --git a/1151/CH11/EX11.2/example2.sce b/1151/CH11/EX11.2/example2.sce new file mode 100755 index 000000000..c398f90f6 --- /dev/null +++ b/1151/CH11/EX11.2/example2.sce @@ -0,0 +1,20 @@ +printf("characterstic equation of the given system is s^2+1.4*s+14=0"); +printf("compare it with the standard second order characterstic equation s^2+2*d*w*s+w^2=0"); +w=sqrt(14); +d=1.4/(2*w); +d1=.7; +Kt=2*(d1-d)/w; +ess1=2*d/w; +disp(d,"damping ratio") +disp(w,"natural frrequency(in rad/sec)=") +disp(ess1,"steady state error without derivative control = ") +printf("overal transfer function with derivative control is C(s)/R(s)=w^2/(s^2+(2dw+w^2Kt)*s+w^2)=14/(s^2+5.23*s+14") +disp(Kt,"Kt=") +rt=(%pi-atan(sqrt((1-d1)/d1)))/(w*sqrt(1-d1^2)); +pt1=%pi/(w*sqrt(1-d1^2)); +mo1=exp((-%pi*d1)/sqrt(1-d1^2))*100; +ess2=2*d1/w+Kt; +disp(pt1,"peak time(in sec)for sytem with derivative control is") +disp(rt,"rise time(in sec)for sytem with derivative control is") +disp(mo1,"maximum overshoot(in %)for sytem with derivative control is") +disp(ess2,"steady state error with derivative control = ") diff --git a/1151/CH11/EX11.3/example3.sce b/1151/CH11/EX11.3/example3.sce new file mode 100755 index 000000000..60ef2b195 --- /dev/null +++ b/1151/CH11/EX11.3/example3.sce @@ -0,0 +1,27 @@ +printf(" Given C(s)/R(s)=10/(s^2+2*s+10)") +printf("characterstic equation of the given system is s^2+2*s+10=0"); +printf("compare it with the standard second order characterstic equation s^2+2*d*w*s+w^2=0"); +w=sqrt(10); +d=1/(2*w); +mo1=exp((-%pi*d)/sqrt(1-d^2))*100; +ts1=4/(d*w); +ess1=2*d/w; +disp(d,"damping ratio") +disp(w,"natural frrequency(in rad/sec)=") +disp(ess1,"steady state error without error control = ") +disp(ts1,"settling time(in sec)for sytem without derivative control is") +disp(mo1,"maximum overshoot(in %)for sytem without derivative control is") +printf("overal transfer function with derivative control is C(s)/R(s)=10(1+s*Ke)/(s*(s+2)+10*(1+s*Ke))") +printf("characterstic equation of the error control system is s^2+s*(2+10*Ke)+10=0"); +Ke=(2*d*w-2)/10; +disp(Ke,"Ke="); +printf("characterstic equation of the error control system with Ke = 0.16 is s^2+s*3.16+10)=0"); +d2=3.16/(2*w); +mo2=exp((-%pi*d2)/sqrt(1-d2^2))*100; +ts2=4/(d2*w); +ess2=2*d2/w; +disp(d2,"damping ratio for error control system=") +disp(w,"natural frrequency(in rad/sec) of error control system=") +disp(ess2,"steady state error with error control = ") +disp(ts2,"settling time(in sec)for sytem with derivative control is") +disp(mo2,"maximum overshoot(in %)for sytem with derivative control is") diff --git a/1151/CH11/EX11.4/example4.sce b/1151/CH11/EX11.4/example4.sce new file mode 100755 index 000000000..f3b65377a --- /dev/null +++ b/1151/CH11/EX11.4/example4.sce @@ -0,0 +1,15 @@ +printf(" unity feedback system with transfer function G(s)=K/(s*(s+10))\n Determine K so that the syetm will have damping ratio0.5.For this value of K find settling time ,peak overshootand peak time for unit step input"); +printf("characterstic equation of the given system is 1+G(s)H(s)=s^2+10*s+K=0"); +printf("compare it with the standard second order characterstic equation s^2+2*d*w*s+w^2=0"); +d=0.5; +K=(10/(2*d))^2; +w=sqrt(K); +ts1=4/(d*w); +mo1=exp((-%pi*d)/sqrt(1-d^2))*100; +pt1=%pi/(w*sqrt(1-d^2)); +disp(K,"K=") +disp(d,"damping ratio") +disp(w,"natural frrequency(in rad/sec)=") +disp(pt1,"peak time(in sec)for sytem is") +disp(mo1,"maximum overshoot(in %)for sytem is") +disp(ts1,"settling time(in sec)is") diff --git a/1151/CH11/EX11.5/example5.sce b/1151/CH11/EX11.5/example5.sce new file mode 100755 index 000000000..b843949ca --- /dev/null +++ b/1151/CH11/EX11.5/example5.sce @@ -0,0 +1,15 @@ +printf(" Given G(s)=10/(s*(s+2))"); +printf("characterstic equation of the given system is 1+G(s)H(s)= s^2+2*s+10=0"); +printf("compare it with the standard second order characterstic equation s^2+2*d*w*s+w^2=0"); +w=sqrt(10); +d=1/(2*w); +ess1=2*d/w; +disp(d,"damping ratio") +disp(w,"natural frrequency(in rad/sec)=") +disp(ess1,"steady state error without error control = ") +printf("overal transfer function with derivative control is C(s)/R(s)=10/(s^2+(2+K)*s+10))") +printf("characterstic equation of the error control system is s^2+s*(2+K)+10=0"); +K=(2*d*w-2)/10; +disp(K,"K="); +ess=0.38 +disp(ess,"steady state error=") diff --git a/1151/CH11/EX11.6/example6.sce b/1151/CH11/EX11.6/example6.sce new file mode 100755 index 000000000..97e9a310d --- /dev/null +++ b/1151/CH11/EX11.6/example6.sce @@ -0,0 +1,16 @@ +printf(" Given G(s)=10/(s*(s+2))"); +printf("characterstic equation of the given system is 1+G(s)H(s)= s^2+2*s+10=0"); +printf("compare it with the standard second order characterstic equation s^2+2*d*w*s+w^2=0"); +w=sqrt(10); +d=1/(2*w); +ess1=2*d/w; +disp(d,"damping ratio") +disp(w,"natural frrequency(in rad/sec)=") +disp(ess1,"steady state error without error control = ") +printf("overal transfer function with derivative control is C(s)/R(s)=2*Ka/(s^2+s*(2+2*Ka)+2*Ka))"); +printf("characterstic equation of the error control system is s^2+s*(2+2*Ka)+2*Ka=0"); +printf("1+Kt=0.2*Ka"); +Ka=(1.98/0.4)^2; +Kt=0.2*Ka-1; +disp(Ka,"Ka=") +disp(Kt,"Kt=") diff --git a/1151/CH11/EX11.7/example7.sce b/1151/CH11/EX11.7/example7.sce new file mode 100755 index 000000000..97e9a310d --- /dev/null +++ b/1151/CH11/EX11.7/example7.sce @@ -0,0 +1,16 @@ +printf(" Given G(s)=10/(s*(s+2))"); +printf("characterstic equation of the given system is 1+G(s)H(s)= s^2+2*s+10=0"); +printf("compare it with the standard second order characterstic equation s^2+2*d*w*s+w^2=0"); +w=sqrt(10); +d=1/(2*w); +ess1=2*d/w; +disp(d,"damping ratio") +disp(w,"natural frrequency(in rad/sec)=") +disp(ess1,"steady state error without error control = ") +printf("overal transfer function with derivative control is C(s)/R(s)=2*Ka/(s^2+s*(2+2*Ka)+2*Ka))"); +printf("characterstic equation of the error control system is s^2+s*(2+2*Ka)+2*Ka=0"); +printf("1+Kt=0.2*Ka"); +Ka=(1.98/0.4)^2; +Kt=0.2*Ka-1; +disp(Ka,"Ka=") +disp(Kt,"Kt=") -- cgit