summaryrefslogtreecommitdiff
path: root/257/CH7/EX7.3/example_7_3.sce
diff options
context:
space:
mode:
Diffstat (limited to '257/CH7/EX7.3/example_7_3.sce')
-rw-r--r--257/CH7/EX7.3/example_7_3.sce21
1 files changed, 21 insertions, 0 deletions
diff --git a/257/CH7/EX7.3/example_7_3.sce b/257/CH7/EX7.3/example_7_3.sce
new file mode 100644
index 000000000..10eaa8a71
--- /dev/null
+++ b/257/CH7/EX7.3/example_7_3.sce
@@ -0,0 +1,21 @@
+p=3/5
+q=poly([10 7 1],'s','coeff');
+G=p/q
+H=1
+y=G*H
+F1=1/(1+y)
+
+F=1/(1+y);
+disp(F,"1/(1+G(s)H(s))=")
+syms t s;
+Ko=limit(s*F/s,s,0) //Ko=Lt s->0 (1/(1+G(s)H(S))
+d=diff(s*F/s,s);
+K1=limit(diff(s*F/s,s),s,0) //K1=Lt s->0 (dF(s)/ds)
+K2=limit(diff(d,s),s,0) //K2=Lt s->0 (d2F(s)/ds)
+
+a=(6+5*t+6*(t^2)/2)
+b=diff((6+5*t+6*(t^2)/2) ,t)
+c=diff(b,t)
+e=Ko*a+K1*b+K2*c //error by dynamic coefficient method
+disp(e,"e = ")
+