summaryrefslogtreecommitdiff
path: root/257/CH4
diff options
context:
space:
mode:
Diffstat (limited to '257/CH4')
-rw-r--r--257/CH4/EX4.1/example_4_1.sce12
-rw-r--r--257/CH4/EX4.14/example_4_14.sce17
-rw-r--r--257/CH4/EX4.16/example_4_16.sce22
-rw-r--r--257/CH4/EX4.17/example_4_17.sce7
-rw-r--r--257/CH4/EX4.18/example_4_18.sce23
-rw-r--r--257/CH4/EX4.19/example_4_19.sce6
-rw-r--r--257/CH4/EX4.2/example_4_2.sce24
-rw-r--r--257/CH4/EX4.20/example_4_20.sce17
-rw-r--r--257/CH4/EX4.21/example_4_21.sce8
-rw-r--r--257/CH4/EX4.3/example_4_3.sce23
-rw-r--r--257/CH4/EX4.4/example_4_4.sce17
-rw-r--r--257/CH4/EX4.5/example_4_5.sce22
-rw-r--r--257/CH4/EX4.6/example_4_6.sce7
-rw-r--r--257/CH4/EX4.7/example_4_7.sce8
14 files changed, 213 insertions, 0 deletions
diff --git a/257/CH4/EX4.1/example_4_1.sce b/257/CH4/EX4.1/example_4_1.sce
new file mode 100644
index 000000000..11438e9ab
--- /dev/null
+++ b/257/CH4/EX4.1/example_4_1.sce
@@ -0,0 +1,12 @@
+// F = M*s^2 +K*X + B*X*s
+syms s t V Q L C R I;
+//force-voltage method
+F=V;
+X=Q;
+M=L;
+K=1/C;
+B=R;
+
+V=I*(s*L + 1/(s*C) + R);
+
+disp("v = L*diff(i) + 1/C*int(i) + i*R") \ No newline at end of file
diff --git a/257/CH4/EX4.14/example_4_14.sce b/257/CH4/EX4.14/example_4_14.sce
new file mode 100644
index 000000000..051e3f1c1
--- /dev/null
+++ b/257/CH4/EX4.14/example_4_14.sce
@@ -0,0 +1,17 @@
+syms M1 X1 B1 K X2 M2 B2 K3 K1 K2 C1 C2 C3 R1 R2 M3 X3 B3 I1 I2 I3 L3 L1 L2 R3
+
+F=M1*X1*s^2 + B1*s*X1 + K*X1 + B2*(X1-X2)*s
+zero=M2*X2*s^2 + B2*s*(X2-X1)+K3*X2+K2*(X2-X3)
+zro=K2*(X3-X2)+M3*s^2*X3+B3*s*X3
+disp(F)
+disp(zero," 0 =")
+disp(zro," 0 = ")
+
+disp("F-V equations are")
+
+V=L1*s*I1 + R1*I1 + R2*(I1-I2)+I1/(s*C1)
+zero=L2*s*I2 + I2/(s*C3) + R2*(I2-I1)+(I2-I3)/(s*C2)
+zro=(I3-I2)/(s*C2)+L3*s*I3+R3*I3
+disp(V)
+disp(zero)
+disp(zro) \ No newline at end of file
diff --git a/257/CH4/EX4.16/example_4_16.sce b/257/CH4/EX4.16/example_4_16.sce
new file mode 100644
index 000000000..7767fceca
--- /dev/null
+++ b/257/CH4/EX4.16/example_4_16.sce
@@ -0,0 +1,22 @@
+syms K1 K2 x1 x2 x3 B1 B2 M1 M2 V Q L C1 R i1 i2 i3 C2 L1 L2
+disp("equilibrium equations are")
+F=K1*(x1-x2)
+disp(F)
+zero=M1*s^2*x2+K1*(x2-x1)+B1*s*(x2-x3)
+disp(zero)
+zro=M2*s^2*x3+B2*s*x3+K2*x3+B1*s*(x3-x2)
+disp(zro)
+
+//force-voltage method
+F=V;
+X=Q;
+M=L;
+K=1/C;
+B=R;
+disp("F-V equations are")
+V=s*(i1-i2)/C1
+disp(V,"V = ")
+zero=L1*s*i2+(i2-i1)/(s*C)+R1*(i2-i3)
+disp(zero)
+zro=L2*s*i3+R2*i3+i3/(s*C2)+R1*(i3-i2)
+disp(zro) \ No newline at end of file
diff --git a/257/CH4/EX4.17/example_4_17.sce b/257/CH4/EX4.17/example_4_17.sce
new file mode 100644
index 000000000..6c7ef98fb
--- /dev/null
+++ b/257/CH4/EX4.17/example_4_17.sce
@@ -0,0 +1,7 @@
+syms xo D1 K1 D2 K2 x1 s C1 C2 R1 R2
+zero=xo*(s*D1+K1+D2*s+K2)-x1*(D1*s+K1)
+disp((D1*s+K1)/(s*D1+K1+D2*s+K2),"xo/x1 = ")
+
+E1=i1*(R1+1/(s*C1)+R2+1/(s*C2))
+Eo=i1*(R2+1/(s*C2))
+disp(Eo/E1,"Eo/E1 = ") \ No newline at end of file
diff --git a/257/CH4/EX4.18/example_4_18.sce b/257/CH4/EX4.18/example_4_18.sce
new file mode 100644
index 000000000..2d4c44da9
--- /dev/null
+++ b/257/CH4/EX4.18/example_4_18.sce
@@ -0,0 +1,23 @@
+syms s t V q L C1 R K1 x1 x2 phi1 phi2 L1 R1 R2 L1 L2
+
+//F=K1*(x1-x2)
+//0=K1*(X2-X1) + M2*s^2*X2 + K2*X2 + B2*s*X2
+
+//F-V anolagy
+
+x1=q1;
+x2=q2;
+K1=1/C1;
+B=R;
+disp("V = 1/C1 * (q1-q2)")
+disp("0=1/C1*(q2-q1) + L2*s^2*q2 + q2/C2 + R2*s*q2")
+
+//F-I ANOLOGY
+
+M=C;
+B=1/R;
+K=1/L;
+
+disp("I=1/L1*(phi1-phi2)")
+disp("0=1/L1*(phi2-phi1) + C2*s^2*phi2 + 1/R2*s*phi2 + 1/L*phi2")
+
diff --git a/257/CH4/EX4.19/example_4_19.sce b/257/CH4/EX4.19/example_4_19.sce
new file mode 100644
index 000000000..21d4919cf
--- /dev/null
+++ b/257/CH4/EX4.19/example_4_19.sce
@@ -0,0 +1,6 @@
+syms B1 x1 x2 B2 K M s
+disp("differential equations are")
+F=B1*s*(x1-x2)
+disp(F)
+zero=B1*s*(x2-x1)+M*s^2*x2+K*x2+B2*s*x2
+disp(zero) \ No newline at end of file
diff --git a/257/CH4/EX4.2/example_4_2.sce b/257/CH4/EX4.2/example_4_2.sce
new file mode 100644
index 000000000..74a256d4d
--- /dev/null
+++ b/257/CH4/EX4.2/example_4_2.sce
@@ -0,0 +1,24 @@
+syms s t V q L C R L1 L2 X1 X2 q1 R1 R2 q2 X
+
+//F-V anolagy
+F=V;
+x=q;
+M=L;
+K=1/C;
+B=R;
+
+V=L1*q1*s^2 + R1*s*q1 + R2*s*(X1-X2)
+//0=L2*s^2*q2 + (1/C)*q2 + R2*s*(q2-q1)
+//REPLACING I/s=Q
+disp("V=L1*s*I1 + R1*I1 + R2*(I1-I2)") //LOOP 1
+disp("0=L2*s*I2 + 1/(s*C) + R2*(I2-I1)") //LOOP 2
+
+//F-I ANOLOGY
+
+phi=X;
+F=I;
+
+I=phi*(C*s^2 + 1/(R*s) + 1/L)
+//REPLACING phi=V/s
+I=V*(s*C + 1/R + 1/(s*L))
+disp("i(t)=c*diff(v) + v/R + 1/L*int(v)") //taking laplace inverse
diff --git a/257/CH4/EX4.20/example_4_20.sce b/257/CH4/EX4.20/example_4_20.sce
new file mode 100644
index 000000000..e8619fdda
--- /dev/null
+++ b/257/CH4/EX4.20/example_4_20.sce
@@ -0,0 +1,17 @@
+syms i1 i2 R1 R2 C s L x1 x2 B1 B2 M K R
+disp("loop equations are")
+V=i1*(R1+R2)+1/C*('i1-i2')/s-i2*R2
+disp(V)
+zero=L*s*i2+i2*R2-i1*R2+(i1-i2)/(C*s)
+disp(zero)
+
+//force-voltage method
+F=V;
+M=L;
+K=1/C;
+B=R;
+disp("F-V equtions are")
+F=s*x1*(B1+B2)+K*(x1-x2)-s*x2*B2
+disp(F)
+zero=M*s^2*x2+B2*(x2-x1)*s+K*(x2-x1)
+disp(zero) \ No newline at end of file
diff --git a/257/CH4/EX4.21/example_4_21.sce b/257/CH4/EX4.21/example_4_21.sce
new file mode 100644
index 000000000..ef7ed077a
--- /dev/null
+++ b/257/CH4/EX4.21/example_4_21.sce
@@ -0,0 +1,8 @@
+syms M1 X1 B1 K X2 M2 B2
+
+F=M1*X1*s^2 + B1*s*X1 + K*X1 + K1*(X1-X2)
+//M2*X2*s^2 + K1*(X2-X1)=0
+X2=X1*K1/(s^2*M2+K1)
+
+disp(X2/F,"X2/F = ")
+
diff --git a/257/CH4/EX4.3/example_4_3.sce b/257/CH4/EX4.3/example_4_3.sce
new file mode 100644
index 000000000..2d4c44da9
--- /dev/null
+++ b/257/CH4/EX4.3/example_4_3.sce
@@ -0,0 +1,23 @@
+syms s t V q L C1 R K1 x1 x2 phi1 phi2 L1 R1 R2 L1 L2
+
+//F=K1*(x1-x2)
+//0=K1*(X2-X1) + M2*s^2*X2 + K2*X2 + B2*s*X2
+
+//F-V anolagy
+
+x1=q1;
+x2=q2;
+K1=1/C1;
+B=R;
+disp("V = 1/C1 * (q1-q2)")
+disp("0=1/C1*(q2-q1) + L2*s^2*q2 + q2/C2 + R2*s*q2")
+
+//F-I ANOLOGY
+
+M=C;
+B=1/R;
+K=1/L;
+
+disp("I=1/L1*(phi1-phi2)")
+disp("0=1/L1*(phi2-phi1) + C2*s^2*phi2 + 1/R2*s*phi2 + 1/L*phi2")
+
diff --git a/257/CH4/EX4.4/example_4_4.sce b/257/CH4/EX4.4/example_4_4.sce
new file mode 100644
index 000000000..c5f4942a6
--- /dev/null
+++ b/257/CH4/EX4.4/example_4_4.sce
@@ -0,0 +1,17 @@
+syms V q L C R I phi
+
+//T=J diff(diff(theta,t),t) + K*theta + B*diff(theta,t)
+
+//F-V anolagy
+T=V;
+theta=q;
+J=L;
+K=1/C;
+B=R;
+disp("V= L*s*I + R*I + I/(s*C)")
+
+//F-C anology
+
+T=I;
+theta=phi
+disp("I= C*V*s + V/R + V/(s*L)")
diff --git a/257/CH4/EX4.5/example_4_5.sce b/257/CH4/EX4.5/example_4_5.sce
new file mode 100644
index 000000000..54a05271e
--- /dev/null
+++ b/257/CH4/EX4.5/example_4_5.sce
@@ -0,0 +1,22 @@
+syms V C1 q1 q2 L1 R1
+
+disp("equivalent systems equations")
+disp("T = K1 *(theta1-theta2)")
+disp("0=K1*(theta3-theta2) + J1*s^2*theta2 + B1*s*theta2 + K*(theta2-theta3) + B*s*(theta2-theta3) ")
+disp("0=K*(theta3-theta2) + B*s*(theta3-theta2) + J2*s^2*theta3 + B2*s*theta3 + K2*theta3")
+
+
+
+//F-V anology
+
+T=V;
+K1=1/C1;
+theta1=q1;
+theta2=q2;
+J1=L1;
+B1=R1;
+disp("FV analogy")
+disp(" V = 1/C1 *(q1-q2)")
+disp("0=1/C1*(q3-q2) + L1*s^22*q2 + R1*s*q2 + 1/C*(q2-q3) + R*s*(q2-q3) ")
+disp("0=1/C*(q3-q2) + R*s*(q3-q2) + L2*s^2*q3 + R2*s*q3 + 1/C2*q3")
+
diff --git a/257/CH4/EX4.6/example_4_6.sce b/257/CH4/EX4.6/example_4_6.sce
new file mode 100644
index 000000000..79a65b3d3
--- /dev/null
+++ b/257/CH4/EX4.6/example_4_6.sce
@@ -0,0 +1,7 @@
+syms I R C s
+
+V1=I*(R+1/(s*C))
+V2=I*(1/(s*C))
+
+disp(V2/V1,"V2/V1 = ")
+
diff --git a/257/CH4/EX4.7/example_4_7.sce b/257/CH4/EX4.7/example_4_7.sce
new file mode 100644
index 000000000..ef7ed077a
--- /dev/null
+++ b/257/CH4/EX4.7/example_4_7.sce
@@ -0,0 +1,8 @@
+syms M1 X1 B1 K X2 M2 B2
+
+F=M1*X1*s^2 + B1*s*X1 + K*X1 + K1*(X1-X2)
+//M2*X2*s^2 + K1*(X2-X1)=0
+X2=X1*K1/(s^2*M2+K1)
+
+disp(X2/F,"X2/F = ")
+
generated by cgit (git 2.34.1) at 2024-12-19 18:28:32 +0000