summaryrefslogtreecommitdiff
path: root/29/CH6/EX6.10.23/exa6_10_23.sce
diff options
context:
space:
mode:
Diffstat (limited to '29/CH6/EX6.10.23/exa6_10_23.sce')
-rwxr-xr-x29/CH6/EX6.10.23/exa6_10_23.sce24
1 files changed, 24 insertions, 0 deletions
diff --git a/29/CH6/EX6.10.23/exa6_10_23.sce b/29/CH6/EX6.10.23/exa6_10_23.sce
new file mode 100755
index 000000000..2abaf3686
--- /dev/null
+++ b/29/CH6/EX6.10.23/exa6_10_23.sce
@@ -0,0 +1,24 @@
+//caption:determine_sensitivity_w.r.t_K_and_H
+//example 6.10.23
+//page 192
+syms Wn zeta K H;
+s=%s;
+A=sym('Wn^2/(s^2+2*zeta*Wn*s)');
+B=(A*K);
+CL=B/(1+B*H);
+CL=simple(CL);
+disp(CL,"overall_transfer_function:M(s)=");
+disp("sensitivity_w.r.t_K:");
+a=diff(CL,K);
+b=K/CL;
+b=simple(b);
+Sk=a*b;
+Sk=simple(Sk);
+disp(Sk,"sensitivity_w.r.t_K,Sk=");
+disp("sensitivity_w.r.t_H:");
+c=diff(CL,H);
+d=H/CL;
+d=simple(d);
+Sh=c*d;
+Sh=simple(Sh);
+disp(Sh,"sensitivity_w.r.t_H,Sh="); \ No newline at end of file