diff options
Diffstat (limited to '257/CH10')
| -rw-r--r-- | 257/CH10/EX10.1/example_10_1.sce | 32 | ||||
| -rw-r--r-- | 257/CH10/EX10.2/example_10_2.sce | 14 | ||||
| -rw-r--r-- | 257/CH10/EX10.3/example_10_3.sce | 10 | ||||
| -rw-r--r-- | 257/CH10/EX10.4/example_10_4.sce | 12 | ||||
| -rw-r--r-- | 257/CH10/EX10.5/example_10_5.sce | 14 | ||||
| -rw-r--r-- | 257/CH10/EX10.6/example_10_6.sce | 14 | ||||
| -rw-r--r-- | 257/CH10/EX10.7/example_10_7.sce | 15 | ||||
| -rw-r--r-- | 257/CH10/EX10.8/example_10_8.sce | 11 | ||||
| -rw-r--r-- | 257/CH10/EX10.9/example_10_9.sce | 18 |
9 files changed, 140 insertions, 0 deletions
diff --git a/257/CH10/EX10.1/example_10_1.sce b/257/CH10/EX10.1/example_10_1.sce new file mode 100644 index 000000000..6f69e6e58 --- /dev/null +++ b/257/CH10/EX10.1/example_10_1.sce @@ -0,0 +1,32 @@ +s=%s
+Mr=2 //given
+omegaR=3 //given
+zeta=roots(16*s^4 - 16*s^2 + 1) //Mr=1/(2*zeta*sqrt(1-zeta^2))
+zeta(3,1)=0.933
+zeta(2,1)=0.0669
+disp(zeta)
+
+omegaN=omegaR/(sqrt(1-2*(0.2588)^2))
+disp(omegaN,"omegaN = ")
+
+TF= (omegaN^2)/poly([omegaN^2 2*0.2588*omegaN 1],'s',"coeff")
+disp(TF," transfer function = ")
+
+omegaD=omegaN*sqrt(1-(0.2588))
+Tr=(%pi-(atan(sqrt(1-(0.2588)^2)/0.2588)))/(omegaD)
+disp(Tr,"Tr = ")
+
+Tp=%pi/omegaD;
+disp(Tp,"Tp= ")
+
+Ts=4/(0.2588*omegaN)
+disp(Ts,"Ts = ")
+
+Tosc=2*%pi/omegaD
+disp(Tosc,"Tosc = ")
+
+N=Ts/Tosc;
+disp(N,"number of oscillations = ")
+
+Mp=%e^(-0.2588*%pi/(sqrt(1-(0.2588)^2)))
+disp(Mp,"Mp = ")
\ No newline at end of file diff --git a/257/CH10/EX10.2/example_10_2.sce b/257/CH10/EX10.2/example_10_2.sce new file mode 100644 index 000000000..ee36e843b --- /dev/null +++ b/257/CH10/EX10.2/example_10_2.sce @@ -0,0 +1,14 @@ +Mp=0.12 //given from table
+zeta=0.5594
+Tp=0.2
+
+omegaN=%pi/(Tp*sqrt(1-zeta^2));
+disp(omegaN,"omegaN = ")
+
+Mr=1/(2*zeta*sqrt(1-zeta^2))
+disp(Mr,"Mr = ")
+
+
+omegaR=omegaN*(sqrt(1-2*(zeta)^2))
+disp(omegaR," omegaR= ")
+
diff --git a/257/CH10/EX10.3/example_10_3.sce b/257/CH10/EX10.3/example_10_3.sce new file mode 100644 index 000000000..a5a181fec --- /dev/null +++ b/257/CH10/EX10.3/example_10_3.sce @@ -0,0 +1,10 @@ +
+s=%s;
+G=10/(s*(s+10))
+T=G/(1+G)
+disp(T,"T = ")
+
+//compare A*sin(w*t) and 10*sin(8*t)
+A=10;
+w=8;
+disp("c(t) = A*10/(sqrt((10-w^2)^2 + 100*w))*(sin(8*t-atan(10*w/(10-w^2))))")
diff --git a/257/CH10/EX10.4/example_10_4.sce b/257/CH10/EX10.4/example_10_4.sce new file mode 100644 index 000000000..1864171a2 --- /dev/null +++ b/257/CH10/EX10.4/example_10_4.sce @@ -0,0 +1,12 @@ +Mp=16.2
+Tp=%pi/(5*sqrt(3))
+zeta=0.5
+
+omegaD=%pi/Tp
+omegaN=omegaD/(sqrt(1-zeta^2))
+disp(omegaD,"omegaD = ")
+disp(omegaN,"omegaN = ")
+
+disp(omegaN*(sqrt(1-2*zeta^2)), " omegaR = ")
+disp(1/(2*zeta*sqrt(1-zeta^2))," Mr = ")
+
diff --git a/257/CH10/EX10.5/example_10_5.sce b/257/CH10/EX10.5/example_10_5.sce new file mode 100644 index 000000000..e2f1bae77 --- /dev/null +++ b/257/CH10/EX10.5/example_10_5.sce @@ -0,0 +1,14 @@ +s=poly(0,'s');
+omegaR=7
+
+zeta=abs(roots(poly([0.043766 0 -1 0 1],'s',"coeff")))
+disp(zeta)
+disp("but for zeta>0.7 , Mr does not exist, so neglect higher value")
+zeta=0.2141
+disp(zeta)
+
+omegaN=omegaR/(sqrt(1-2*(zeta)^2))
+disp(omegaN,"omegaN = ")
+
+disp(omegaN*sqrt(1-2*zeta^2 + sqrt(2-4*zeta^2+4*zeta^4))," B.W. = ")
+
diff --git a/257/CH10/EX10.6/example_10_6.sce b/257/CH10/EX10.6/example_10_6.sce new file mode 100644 index 000000000..2e164637f --- /dev/null +++ b/257/CH10/EX10.6/example_10_6.sce @@ -0,0 +1,14 @@ +Mr=1.1
+omegaR=11.2
+
+zeta=abs(roots(poly([0.2066 0 -1 0 1],'s',"coeff")))
+disp(zeta)
+disp("but for zeta>0.7 , Mr does not exist, so neglect higher value")
+zeta=0.54
+disp(zeta)
+
+omegaN=omegaR/(sqrt(1-2*(zeta)^2))
+disp(omegaN,"omegaN = ")
+
+TF=omegaN^2/(poly([0 2*zeta*omegaN 1],'s',"coeff"))
+disp(TF," TF = ")
\ No newline at end of file diff --git a/257/CH10/EX10.7/example_10_7.sce b/257/CH10/EX10.7/example_10_7.sce new file mode 100644 index 000000000..5f03883dc --- /dev/null +++ b/257/CH10/EX10.7/example_10_7.sce @@ -0,0 +1,15 @@ +s=%s
+
+//G=k/(s*(s*tau+1)) and T=G/(1+G)
+omegaN=12
+
+zeta=roots(poly([0.22225 0 -1 0 1],'s',"coeff"))
+disp("but for zeta>0.7 , Mr does not exist, so neglect higher value")
+zeta=0.578
+disp(zeta)
+
+tau=1/(2*sqrt(144)*zeta)
+disp(tau,"tau = ")
+k=144*tau;disp(k,"k = ")
+
+disp(omegaN*sqrt(1-2*zeta^2 + sqrt(2-4*zeta^2+4*zeta^4))," B.W. = ")
\ No newline at end of file diff --git a/257/CH10/EX10.8/example_10_8.sce b/257/CH10/EX10.8/example_10_8.sce new file mode 100644 index 000000000..4e5c422e8 --- /dev/null +++ b/257/CH10/EX10.8/example_10_8.sce @@ -0,0 +1,11 @@ +s=%s;
+G=100/(s*(s+8))
+T=G/(1+G)
+disp(T,"T = ")
+
+//compare denominator with s^2+2*zeta*omegaN + omegaN^2
+omegaN=10
+zeta=0.4;
+disp(1/(2*zeta*sqrt((1-zeta^2)))," Mr = ")
+disp(omegaN*sqrt(1-2*zeta^2)," omegaN = ")
+
diff --git a/257/CH10/EX10.9/example_10_9.sce b/257/CH10/EX10.9/example_10_9.sce new file mode 100644 index 000000000..f99a25731 --- /dev/null +++ b/257/CH10/EX10.9/example_10_9.sce @@ -0,0 +1,18 @@ +s=%s
+//G=k/(s*(s+a)) and T=G/(1+G)
+Mr=1.04
+omegaR=11.55
+
+disp("for zeta>0.7 , Mr does not exist, so neglect higher value")
+zeta=0.6021
+disp(zeta)
+
+omegaN=omegaR/(sqrt(1-2*(zeta)^2))
+disp(omegaN,"omegaN = ")
+
+
+k=omegaN^2
+disp(k,"k=")
+disp(2*zeta*omegaN," a = ")
+disp(omegaN*sqrt(1-2*zeta^2 + sqrt(2-4*zeta^2+4*zeta^4))," B.W. = ")
+disp(4/(zeta*omegaN)," Ts = ")
\ No newline at end of file |