diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /257/CH9 | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
Diffstat (limited to '257/CH9')
29 files changed, 196 insertions, 0 deletions
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)
|