From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001 From: prashantsinalkar Date: Tue, 10 Oct 2017 12:27:19 +0530 Subject: initial commit / add all books --- 257/CH9/EX9.1/example_9_1.sce | 4 ++++ 257/CH9/EX9.10/example_9_10.sce | 3 +++ 257/CH9/EX9.11/example_9_11.sce | 3 +++ 257/CH9/EX9.12/example_9_12.sce | 8 ++++++++ 257/CH9/EX9.13/example_9_13.sce | 6 ++++++ 257/CH9/EX9.14/example_9_14.sce | 9 +++++++++ 257/CH9/EX9.15/example_9_15.sce | 3 +++ 257/CH9/EX9.16/example_9_16.sce | 3 +++ 257/CH9/EX9.17/example_9_17.sce | 3 +++ 257/CH9/EX9.18/example_9_18.sce | 3 +++ 257/CH9/EX9.2/eg_9_2.sce | 14 ++++++++++++++ 257/CH9/EX9.20/example_9_20.sce | 8 ++++++++ 257/CH9/EX9.21/example_9_21.sce | 4 ++++ 257/CH9/EX9.23/example_9_23.sce | 3 +++ 257/CH9/EX9.24/example_9_24.sce | 3 +++ 257/CH9/EX9.25/example_9_25.sce | 25 +++++++++++++++++++++++++ 257/CH9/EX9.26/example_9_26.sce | 16 ++++++++++++++++ 257/CH9/EX9.27/example_9_27.sce | 4 ++++ 257/CH9/EX9.28/example_9_28.sce | 9 +++++++++ 257/CH9/EX9.3/example_9_3.sce | 4 ++++ 257/CH9/EX9.30/example_9_30.sce | 11 +++++++++++ 257/CH9/EX9.31/example_9_31.sce | 4 ++++ 257/CH9/EX9.32/example_9_32.sce | 22 ++++++++++++++++++++++ 257/CH9/EX9.4/example_9_4.sce | 4 ++++ 257/CH9/EX9.5/example_9_5.sce | 3 +++ 257/CH9/EX9.6/example_9_6.sce | 3 +++ 257/CH9/EX9.7/example_9_7.sce | 3 +++ 257/CH9/EX9.8/example_9_8.sce | 8 ++++++++ 257/CH9/EX9.9/example_9_9.sce | 3 +++ 29 files changed, 196 insertions(+) create mode 100644 257/CH9/EX9.1/example_9_1.sce create mode 100644 257/CH9/EX9.10/example_9_10.sce create mode 100644 257/CH9/EX9.11/example_9_11.sce create mode 100644 257/CH9/EX9.12/example_9_12.sce create mode 100644 257/CH9/EX9.13/example_9_13.sce create mode 100644 257/CH9/EX9.14/example_9_14.sce create mode 100644 257/CH9/EX9.15/example_9_15.sce create mode 100644 257/CH9/EX9.16/example_9_16.sce create mode 100644 257/CH9/EX9.17/example_9_17.sce create mode 100644 257/CH9/EX9.18/example_9_18.sce create mode 100644 257/CH9/EX9.2/eg_9_2.sce create mode 100644 257/CH9/EX9.20/example_9_20.sce create mode 100644 257/CH9/EX9.21/example_9_21.sce create mode 100644 257/CH9/EX9.23/example_9_23.sce create mode 100644 257/CH9/EX9.24/example_9_24.sce create mode 100644 257/CH9/EX9.25/example_9_25.sce create mode 100644 257/CH9/EX9.26/example_9_26.sce create mode 100644 257/CH9/EX9.27/example_9_27.sce create mode 100644 257/CH9/EX9.28/example_9_28.sce create mode 100644 257/CH9/EX9.3/example_9_3.sce create mode 100644 257/CH9/EX9.30/example_9_30.sce create mode 100644 257/CH9/EX9.31/example_9_31.sce create mode 100644 257/CH9/EX9.32/example_9_32.sce create mode 100644 257/CH9/EX9.4/example_9_4.sce create mode 100644 257/CH9/EX9.5/example_9_5.sce create mode 100644 257/CH9/EX9.6/example_9_6.sce create mode 100644 257/CH9/EX9.7/example_9_7.sce create mode 100644 257/CH9/EX9.8/example_9_8.sce create mode 100644 257/CH9/EX9.9/example_9_9.sce (limited to '257/CH9') diff --git a/257/CH9/EX9.1/example_9_1.sce b/257/CH9/EX9.1/example_9_1.sce new file mode 100644 index 000000000..83b111455 --- /dev/null +++ b/257/CH9/EX9.1/example_9_1.sce @@ -0,0 +1,4 @@ +//root locus +s=%s +sys=syslin('c',(k)/(s)) +evans(sys) diff --git a/257/CH9/EX9.10/example_9_10.sce b/257/CH9/EX9.10/example_9_10.sce new file mode 100644 index 000000000..eb5db031e --- /dev/null +++ b/257/CH9/EX9.10/example_9_10.sce @@ -0,0 +1,3 @@ +s=%s +sys=syslin('c',(k)*(s+2)*(s+4)/(s^2*(s+6))) +evans(sys) diff --git a/257/CH9/EX9.11/example_9_11.sce b/257/CH9/EX9.11/example_9_11.sce new file mode 100644 index 000000000..13e358b6e --- /dev/null +++ b/257/CH9/EX9.11/example_9_11.sce @@ -0,0 +1,3 @@ +s=%s +sys=syslin('c',(k)*(s+4)/((s+2)*(s^2+2*s+2))) +evans(sys) diff --git a/257/CH9/EX9.12/example_9_12.sce b/257/CH9/EX9.12/example_9_12.sce new file mode 100644 index 000000000..66d216f0c --- /dev/null +++ b/257/CH9/EX9.12/example_9_12.sce @@ -0,0 +1,8 @@ +s=%s +sys=syslin('c',(k)/(s*(s+1)*(s+4))) +evans(sys) + +//breakaway points + +disp("break away points are") +disp(roots(numer(derivat(sys)))) \ No newline at end of file diff --git a/257/CH9/EX9.13/example_9_13.sce b/257/CH9/EX9.13/example_9_13.sce new file mode 100644 index 000000000..7bc97da9d --- /dev/null +++ b/257/CH9/EX9.13/example_9_13.sce @@ -0,0 +1,6 @@ +//angles of departure +sys=syslin('c',k*(s+2)/(s*(s+4)*(s^2+2*s+2))) +evans(sys) + +theta=180-(135+90+18.43-45) +disp(theta) \ No newline at end of file diff --git a/257/CH9/EX9.14/example_9_14.sce b/257/CH9/EX9.14/example_9_14.sce new file mode 100644 index 000000000..7ee66f365 --- /dev/null +++ b/257/CH9/EX9.14/example_9_14.sce @@ -0,0 +1,9 @@ +s=%s +sys=syslin('c',(k)/(s*(s+5)*(s+10))) +clf +evans(sys) + +//stability + +[Ki,s]=kpure(sys) +disp(Ki) \ No newline at end of file diff --git a/257/CH9/EX9.15/example_9_15.sce b/257/CH9/EX9.15/example_9_15.sce new file mode 100644 index 000000000..8bf3da37f --- /dev/null +++ b/257/CH9/EX9.15/example_9_15.sce @@ -0,0 +1,3 @@ +s=%s +sys=syslin('c',(k)/(s*(s+1)*(s+2)*(s+3))) +evans(sys) \ No newline at end of file diff --git a/257/CH9/EX9.16/example_9_16.sce b/257/CH9/EX9.16/example_9_16.sce new file mode 100644 index 000000000..8d7e381c8 --- /dev/null +++ b/257/CH9/EX9.16/example_9_16.sce @@ -0,0 +1,3 @@ +s=%s +sys=syslin('c',(k)/(s*(s+3)*(s^2+3*s+4.5))) +evans(sys) \ No newline at end of file diff --git a/257/CH9/EX9.17/example_9_17.sce b/257/CH9/EX9.17/example_9_17.sce new file mode 100644 index 000000000..c60c035d0 --- /dev/null +++ b/257/CH9/EX9.17/example_9_17.sce @@ -0,0 +1,3 @@ +s=%s +sys=syslin('c',(k)/(s*(s+3)*(s^2+3*s+11.25))) +evans(sys) \ No newline at end of file diff --git a/257/CH9/EX9.18/example_9_18.sce b/257/CH9/EX9.18/example_9_18.sce new file mode 100644 index 000000000..cf32c764d --- /dev/null +++ b/257/CH9/EX9.18/example_9_18.sce @@ -0,0 +1,3 @@ +s=%s +sys=syslin('c',(k)/(s*(s+3)*(s^2+3*s+3))) +evans(sys) \ No newline at end of file diff --git a/257/CH9/EX9.2/eg_9_2.sce b/257/CH9/EX9.2/eg_9_2.sce new file mode 100644 index 000000000..b921c35b4 --- /dev/null +++ b/257/CH9/EX9.2/eg_9_2.sce @@ -0,0 +1,14 @@ +//angle condition +s=%s +sys=syslin('c',(k)/(s*(s+2)*(s+4))) + +f=-0.75; +//disp(-atan(s,0)) +disp(-atan(f,0)-atan(f,2)-atan(f,4)) +if(-atan(f,0)-atan(f,2)-atan(f,4)==(-3.14)) + printf("yes") + +else + printf("no") +end + diff --git a/257/CH9/EX9.20/example_9_20.sce b/257/CH9/EX9.20/example_9_20.sce new file mode 100644 index 000000000..6d7cff76e --- /dev/null +++ b/257/CH9/EX9.20/example_9_20.sce @@ -0,0 +1,8 @@ +s=%s +sys=syslin('c',(s+1)*k/(s*(s-1)*(s^2+5*s+20))) +evans(sys) + +//stability + +[Ki,s]=kpure(sys) +disp(Ki) \ No newline at end of file diff --git a/257/CH9/EX9.21/example_9_21.sce b/257/CH9/EX9.21/example_9_21.sce new file mode 100644 index 000000000..49be78017 --- /dev/null +++ b/257/CH9/EX9.21/example_9_21.sce @@ -0,0 +1,4 @@ +//given characteristic equation s^3 + 9*s^2 +k*s + k +s=%s +sys=syslin('c',(s+1)*k/((s+9)*(s^2))) +evans(sys) diff --git a/257/CH9/EX9.23/example_9_23.sce b/257/CH9/EX9.23/example_9_23.sce new file mode 100644 index 000000000..7f1cc10c5 --- /dev/null +++ b/257/CH9/EX9.23/example_9_23.sce @@ -0,0 +1,3 @@ +s=%s +sys=syslin('c',k/(s*(s+4)*(s^2+4*s+20))) +evans(sys) diff --git a/257/CH9/EX9.24/example_9_24.sce b/257/CH9/EX9.24/example_9_24.sce new file mode 100644 index 000000000..94d7a1535 --- /dev/null +++ b/257/CH9/EX9.24/example_9_24.sce @@ -0,0 +1,3 @@ +s=%s +sys=syslin('c',k/(s*(s+2)*(s^2+6*s+25))) +evans(sys) diff --git a/257/CH9/EX9.25/example_9_25.sce b/257/CH9/EX9.25/example_9_25.sce new file mode 100644 index 000000000..17b53b06b --- /dev/null +++ b/257/CH9/EX9.25/example_9_25.sce @@ -0,0 +1,25 @@ +s=%s +sys=syslin('c',k/(s*(s+4)*(s+2))) +evans(sys) + + +//values of k + +[Ki,s]=kpure(sys) +disp("k should be less than") +disp(Ki) + +//frequency of oscillations +s=%s +P=s^3+6*s^2+8*s+Ki +routh=routh_t(P) +disp(routh) +disp("frequency of oscillations is") + disp(sqrt((roots(routh(2,:))))) + + //damping ratio is 0.5 given +// cos inverse of 0.5 is 60 degrees. from the root locus, the 60 degree line crosses the locus at (-0.75+j*1.25) + +f=(-0.75+%i*1.25) +disp("k for damping ratio 0.5 is") +disp(abs(f*(f+4)*(f+2))) \ No newline at end of file diff --git a/257/CH9/EX9.26/example_9_26.sce b/257/CH9/EX9.26/example_9_26.sce new file mode 100644 index 000000000..82154dd9d --- /dev/null +++ b/257/CH9/EX9.26/example_9_26.sce @@ -0,0 +1,16 @@ +s=%s +sys=syslin('c',((k)*(s^2-2*s+5))/((s+2)*(s-0.5))) +clf +evans(sys) + +//stability + +[Ki,s]=kpure(sys) +disp(Ki) + +//damping ratio 0.5 + +f=(-0.3+%i*0.55) +disp("k for damping ratio 0.5 is") +disp(abs(((f+2)*(f-0.5)))/(f^2-2*f+5)) + diff --git a/257/CH9/EX9.27/example_9_27.sce b/257/CH9/EX9.27/example_9_27.sce new file mode 100644 index 000000000..7eebb552f --- /dev/null +++ b/257/CH9/EX9.27/example_9_27.sce @@ -0,0 +1,4 @@ +//when |a|<|b| then system is stable. hence let a=2 and b=3 +s=%s +sys=syslin('c',k/((s-2)*(s+3))) +evans(sys) \ No newline at end of file diff --git a/257/CH9/EX9.28/example_9_28.sce b/257/CH9/EX9.28/example_9_28.sce new file mode 100644 index 000000000..5403989b2 --- /dev/null +++ b/257/CH9/EX9.28/example_9_28.sce @@ -0,0 +1,9 @@ +s=%s +sys=syslin('c',k/((s+16)*(s^2+2*s+2))) +evans(sys) + +//stability + +[Ki,s]=kpure(sys) +disp("k marginal is") +disp(Ki) diff --git a/257/CH9/EX9.3/example_9_3.sce b/257/CH9/EX9.3/example_9_3.sce new file mode 100644 index 000000000..e0e179131 --- /dev/null +++ b/257/CH9/EX9.3/example_9_3.sce @@ -0,0 +1,4 @@ +//magnitude condition for GH=K/s*(s=2)*(s+4) +s=-0.75 +k=-s*(4+s)*(2+s) +disp(k) \ No newline at end of file diff --git a/257/CH9/EX9.30/example_9_30.sce b/257/CH9/EX9.30/example_9_30.sce new file mode 100644 index 000000000..2d991d48a --- /dev/null +++ b/257/CH9/EX9.30/example_9_30.sce @@ -0,0 +1,11 @@ +s=%s +sys=syslin('c',(k)/((s)*(s^2+8*s+17))) +clf +evans(sys) + +//damping ratio 0.5 + +f=(-1.15+%i*2) +disp("k for damping ratio 0.5 is") +disp(abs(((f)*(f^2+8*f+17)))) + diff --git a/257/CH9/EX9.31/example_9_31.sce b/257/CH9/EX9.31/example_9_31.sce new file mode 100644 index 000000000..5deb38b5e --- /dev/null +++ b/257/CH9/EX9.31/example_9_31.sce @@ -0,0 +1,4 @@ +s=%s +sys=syslin('c',((k)*(s+2)*(s+3))/((s+1)*(s))) +clf +evans(sys) \ No newline at end of file diff --git a/257/CH9/EX9.32/example_9_32.sce b/257/CH9/EX9.32/example_9_32.sce new file mode 100644 index 000000000..979ba13bd --- /dev/null +++ b/257/CH9/EX9.32/example_9_32.sce @@ -0,0 +1,22 @@ +//given characteristic equation we get GH +s=%s +sys=syslin('c',k/(s*(s^2+8*s+20))) +clf +evans(sys) + +//stability + +[Ki,s]=kpure(sys) +disp(Ki) + + +//damping ratio is 0.95 + +f=(-1.8+%i*0.55) +disp("k for damping ratio 0.95 is") +disp(abs((f*(f^2+8*f+20)))) + + +f=(-03.6+%i*1.1) +disp("k for damping ratio 0.95 is") +disp(abs((f*(f^2+8*f+20)))) diff --git a/257/CH9/EX9.4/example_9_4.sce b/257/CH9/EX9.4/example_9_4.sce new file mode 100644 index 000000000..4092e3bf6 --- /dev/null +++ b/257/CH9/EX9.4/example_9_4.sce @@ -0,0 +1,4 @@ +s=%s +sys=syslin('c',(k)/(s*(s+2))) +evans(sys) +printf("there are 2 branches approaching infinity") \ No newline at end of file diff --git a/257/CH9/EX9.5/example_9_5.sce b/257/CH9/EX9.5/example_9_5.sce new file mode 100644 index 000000000..96b345c1c --- /dev/null +++ b/257/CH9/EX9.5/example_9_5.sce @@ -0,0 +1,3 @@ +s=%s +sys=syslin('c',k*(s+1)/(s*(s+5))) +evans(sys) diff --git a/257/CH9/EX9.6/example_9_6.sce b/257/CH9/EX9.6/example_9_6.sce new file mode 100644 index 000000000..0de0dbe92 --- /dev/null +++ b/257/CH9/EX9.6/example_9_6.sce @@ -0,0 +1,3 @@ +s=%s +sys=syslin('c',(k)*(s+1)*(s+4)/(s*(s+3)*(s+5))) +evans(sys) diff --git a/257/CH9/EX9.7/example_9_7.sce b/257/CH9/EX9.7/example_9_7.sce new file mode 100644 index 000000000..0de0dbe92 --- /dev/null +++ b/257/CH9/EX9.7/example_9_7.sce @@ -0,0 +1,3 @@ +s=%s +sys=syslin('c',(k)*(s+1)*(s+4)/(s*(s+3)*(s+5))) +evans(sys) diff --git a/257/CH9/EX9.8/example_9_8.sce b/257/CH9/EX9.8/example_9_8.sce new file mode 100644 index 000000000..d214f1c9b --- /dev/null +++ b/257/CH9/EX9.8/example_9_8.sce @@ -0,0 +1,8 @@ +//Gh=k/((s+1)*(s+2+2j)*(s+2-2j)) + +for n=0:2 + theta=(2*n+1)*180/3 //poles=3 , zeroes=0 + disp(theta) +end + +disp(cntrd=(-1-2-2-0)/3) // real(poles-zeroes)/ number of poles-zeroes \ No newline at end of file diff --git a/257/CH9/EX9.9/example_9_9.sce b/257/CH9/EX9.9/example_9_9.sce new file mode 100644 index 000000000..46438edf3 --- /dev/null +++ b/257/CH9/EX9.9/example_9_9.sce @@ -0,0 +1,3 @@ +s=%s +sys=syslin('c',(k)*(s+6)/(s*(s+2)*(s+4))) +evans(sys) -- cgit