initial commit / add all books

1 files changed, 63 insertions, 0 deletions

diff --git a/3793/CH14/EX14.4/exp_14_4.sce b/3793/CH14/EX14.4/exp_14_4.sce
new file mode 100644
index 000000000..e4de23a44
--- /dev/null
+++ b/3793/CH14/EX14.4/exp_14_4.sce
@@ -0,0 +1,63 @@
+clear;

+clc;

+z1=2.00;

+z2=(-.250);

+z3=1.200;

+z4=11.00;

+z5=(-1.00);

+//z1,z2,z3 are currents and z4,z5 is voltage

+//if v2 is short circuited

+z11=(1/(2+(24/10)));

+z22=(6/10)*z11;

+z33=(4/10)*z11;

+z44=z33*6;

+z55=z22*4;

+//when v1 is short ckted;

+z2a=(1/(4+(12/8)));

+z1a=(6/8)*z2a;

+z3a=(2/8)*z2a;

+z4v=z3a*6;

+z5v=z2a*4;

+H=[z11 -z22; -z22 z2a; z33 z3a; z55 z4v; -z55 z5v];

+

+H1=[H]';

+W=[100 0 0 0 0; 0 100 0 0 0; 0 0 75 0 0; 0 0 0 75 0; 0 0 0 0 75;]

+G=H1*W*H;

+

+z=[z1;z2;z3;z4;z5];

+x=(inv(G))*H1*W*z;

+zh=H*x;

+z=[z1;z2;z3;z4;z5];

+eh=z-zh;

+fh=100*eh(1)^2+100*eh(2)^2+75*eh(3)^2+75*eh(4)^2+75*eh(5)^2;

+k=3;

+alphha=.01;

+R=[.0100 0 0 0 0; 0 .0100 0 0 0; 0 0 .0133 0 0; 0 0 0 .0133 0; 0 0 0 0 .0133];

+F=H*inv(G)*(H')*inv(R);

+Rtemp=[diag(F)];

+

+er1=eh(1)/sqrt((1-F(1,1))*100);

+er2=eh(2)/sqrt((1-F(2,2))*100);

+er3=eh(3)/sqrt((1-F(3,3))*75);

+er4=eh(4)/sqrt((1-F(4,4))*75);

+er5=eh(5)/sqrt((1-F(5,5))*75);

+mprintf("for k=3 and alpha=.001 estimated value of function f is %f, which is greater than critical value of 11.35. therefore it contains some bad data\n",fh);

+mprintf("standardized error are er1=%f , er2=%f, er3=%f, er4=%f, er5=%f and we can see that er4 is greater so it should be rejected\n ",er1,er2,er3,er4,er5);

+H(4,:)=[];

+R(4,:)=[];;

+R(:,4)=[];

+W(4,:)=[];

+W(:,4)=[];

+G1=H'*W*H;

+z(4)=[];

+xh1=inv(G1)*H'*W*z;

+zh1=H*xh1;

+eh1=z-zh1;

+fh1=100*eh1(1)^2+100*eh1(2)^2+75*eh1(3)^2+75*eh1(4)^2;

+mprintf("for k=2 and alpha=.001 estimated value of function f is %f, which is greater than critical value of 9.21. therefore the measurement data is inaccurate and need to be more precise\n",fh1);

+

+

+

+

+

+