7 files changed, 98 insertions, 0 deletions
diff --git a/28/CH3/EX3.2/ex3_2.sce b/28/CH3/EX3.2/ex3_2.sce
new file mode 100755
index 000000000..0e04d7c71
--- /dev/null
+++ b/28/CH3/EX3.2/ex3_2.sce
@@ -0,0 +1,12 @@
+syms K;
+s=%s;
+G=syslin('c',25(s+1)/(s+5));
+p=K;
+q=s^2+s;
+J=p/q;
+F=G*J;
+T=F/(1+F); // Closed loop transfer function
+disp(T,"C(s)/R(s)")
+// sensitivity w.r.t K = dT/dK*K/T
+S=(diff(T,K))*(K/T)
+disp(S,"Sensitivity")
+\ No newline at end of file
diff --git a/28/CH3/EX3.3.a/ex3_3_1.sce b/28/CH3/EX3.3.a/ex3_3_1.sce
new file mode 100755
index 000000000..403637f6f
--- /dev/null
+++ b/28/CH3/EX3.3.a/ex3_3_1.sce
@@ -0,0 +1,15 @@
+syms K1 K t;
+s=%s;
+p=K1*K;
+q=t*s+1+(K1*K);
+T=p/q;
+disp(T,"V(s)/R(s)")
+// sensitivity w.r.t K is dT/dK*K/T
+S=(diff(T,K))*(K/T)
+// given K1=50 K=1.5
+s=0
+S=horner(S,s)
+K1=50;
+K=1.5;
+S=1/(1+K1*K)
+disp(S,"sensitivity=")
+\ No newline at end of file
diff --git a/28/CH3/EX3.3.b/ex3_3_2.sce b/28/CH3/EX3.3.b/ex3_3_2.sce
new file mode 100755
index 000000000..79f3ce513
--- /dev/null
+++ b/28/CH3/EX3.3.b/ex3_3_2.sce
@@ -0,0 +1,14 @@
+syms A K K1 t
+s=%s;
+p=K1*K*A;
+q=s*(1+(t*s)+(K1*K));
+K=1.5;
+K1=50;
+V=p/q
+v=limit(s*V,s,0)
+// given steady state speed = 60km/hr
+A=60*(1+(K1*K))/(K1*K)
+// steady error e(ss)=A-v
+v=60;
+e=A-v;
+disp(e,"e(ss)=")
diff --git a/28/CH3/EX3.3.c/ex3_3_3.sce b/28/CH3/EX3.3.c/ex3_3_3.sce
new file mode 100755
index 000000000..eae27d0c6
--- /dev/null
+++ b/28/CH3/EX3.3.c/ex3_3_3.sce
@@ -0,0 +1,9 @@
+// under stalled conditions
+syms Kg K1 D;
+A=60.8;
+A*K1=Kg*D;
+// given Kg=100 
+Kg=100;
+K1=50;
+D=(A*K1)/Kg;
+disp(D,"upslope=")
+\ No newline at end of file
diff --git a/28/CH3/EX3.3.d/ex3_3_4.sce b/28/CH3/EX3.3.d/ex3_3_4.sce
new file mode 100755
index 000000000..c6ededab5
--- /dev/null
+++ b/28/CH3/EX3.3.d/ex3_3_4.sce
@@ -0,0 +1,8 @@
+// steady speed=10km/hr
+syms K Kg D
+(((A-10)*K1)-(-D*Kg))K=100;
+A=(60.8*10)/60;
+K=1.5;
+Kg=100;
+D=((100/K)-((A-10)*K))/Kg;
+disp(D,"Down slope=")
+\ No newline at end of file
diff --git a/28/CH3/EX3.3.f/ex3_3_5.sce b/28/CH3/EX3.3.f/ex3_3_5.sce
new file mode 100755
index 000000000..d33114f39
--- /dev/null
+++ b/28/CH3/EX3.3.f/ex3_3_5.sce
@@ -0,0 +1,11 @@
+// for open loop system
+// given speed=60km/hr
+syms R K1 K;
+(R*K1*K)=60
+K1=50;
+K=1.5;
+R=60/(K1*K)
+disp(R,"Input open=")
+// for closed loop
+R=60(1+(K1*K))/(K1*K)
+disp(R,"Input closed=")
+\ No newline at end of file
diff --git a/28/CH3/EX3.3.g/ex3_3_6.sce b/28/CH3/EX3.3.g/ex3_3_6.sce
new file mode 100755
index 000000000..0e780d76d
--- /dev/null
+++ b/28/CH3/EX3.3.g/ex3_3_6.sce
@@ -0,0 +1,29 @@
+// for open loop
+syms t g s;
+s=%s;
+K1=50;
+K=1.5;
+g=20;
+V=syslin('c',((K1*K)*0.8)/(s*((g*s)+1)))
+// taking inverse laplace
+v=ilaplace(V,s,t)
+v=60(11-%e^(-t/20))
+// given v=90%
+v=0.9;
+t=-20*log(1-v);
+disp(t,"time open=")
+// for closed loop
+syms K' g'
+s=%s;
+V=syslin('c',(60.8*K')/(s*((g'*s)+1)))
+// taking inverse laplace
+v=ilaplace(V,s,t)
+// given 
+K'=75/76; 
+g'=.263;
+v=60(1-%e^(-t/.263))
+// at v=90%
+v=0.9;
+t=-.263*log(1-(v/60));
+disp(t,"time closed=")
+