diff options
Diffstat (limited to '1523/CH7/EX7.19')
-rwxr-xr-x | 1523/CH7/EX7.19/ex7_19.sce | 31 |
1 files changed, 31 insertions, 0 deletions
diff --git a/1523/CH7/EX7.19/ex7_19.sce b/1523/CH7/EX7.19/ex7_19.sce new file mode 100755 index 000000000..584719d1f --- /dev/null +++ b/1523/CH7/EX7.19/ex7_19.sce @@ -0,0 +1,31 @@ +//Graph Theory : example 7.15 :(pg 7.34 & 7.35)
+Q=[1 -1 0 0;0 -1 1 1];
+printf("\nQ=");
+disp(Q);
+printf("\nThe KCL equation in matrix form is given by");
+printf("\nQ.Yb.(Q^T).Vl=Q.Is-Q.Yb.Vs");
+printf("\nQ.Yb.(Q^T).Vl=Q.Is");//Vs=0
+Yb=diag([5,5,5,10]);
+Is=[-10;0;0;0];
+printf("\nYb=");
+disp(Yb);
+printf("\n(Q^T)=");
+disp(Q');
+printf("\nIs=");
+disp(Is);//current entering into nodes is taken as negative
+x=(Q*Yb);
+printf("\nQ.Yb=");
+disp(x);
+y=(x*Q');
+printf("\nQ.Yb.(Q^T)=");
+disp(y);
+z=(Q*Is);
+printf("\nQ.Is=");
+disp(z);
+printf("\nLoad voltages:");
+M=[10 5;5 20];
+P=inv(M);
+N=[-10;0];
+X=(P*N);
+disp(X);
+printf("\nvl1=-1.14 V \nvl2=0.28 V");
|