From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001 From: prashantsinalkar Date: Tue, 10 Oct 2017 12:27:19 +0530 Subject: initial commit / add all books --- 3793/CH14/EX14.2/exp_14_2.sce | 38 ++++++++++++++++++++++++++++++++++++++ 1 file changed, 38 insertions(+) create mode 100644 3793/CH14/EX14.2/exp_14_2.sce (limited to '3793/CH14/EX14.2') diff --git a/3793/CH14/EX14.2/exp_14_2.sce b/3793/CH14/EX14.2/exp_14_2.sce new file mode 100644 index 000000000..1e5287b65 --- /dev/null +++ b/3793/CH14/EX14.2/exp_14_2.sce @@ -0,0 +1,38 @@ +clear; +clc; +z12=complex(.05,.20); +z23=complex(.075,.25); +c1=.025; +c2=.005; +w1= (.1568*10^(-4)); +w2= (.1679*10^(-4)); +w3= (.0668*10^(-4)); +w4= (.0702*10^(-4)); +W=[w1 0 0 0; 0 w2 0 0; 0 0 w3 0; 0 0 0 w4]; +v1=1.05; +v2=1.05; +v3=(1.05); +h1=(v1/z12); +h2=(v2/z12); +h3=(v2/z23); +h4=(v3/z23); +H=[h1 0 0 0; 0 h2 0 0; 0 0 h3 0; 0 0 0 h4]; +H1=conj(H); +D=H1*W*H; +D1=real(D); +A=[1 -1 0; -1 1 0; 0 1 -1; 0 -1 1]; +B=[-1 0; 1 0; 1 -1; -1 1]; +b=[1;-1;0;0]; +E=(B')*D; +f=E*B; +s1=complex(.50,-.12); +s2=complex(-.48,.10); +s3=complex(.80,-.40); +s4=complex(-.78,.38); +S=[s1;s2;s3;s4]; +vm=(inv(H))*(conj(S)); +vb=inv(f)*E*(vm-(b*v1)); +V=[v1;vb]; +printf("V = ") +disp(V); + -- cgit