diff options
Diffstat (limited to '3872/CH8/EX8.2/EX8_2.sce')
| -rw-r--r-- | 3872/CH8/EX8.2/EX8_2.sce | 17 |
1 files changed, 17 insertions, 0 deletions
diff --git a/3872/CH8/EX8.2/EX8_2.sce b/3872/CH8/EX8.2/EX8_2.sce new file mode 100644 index 000000000..8b3d83108 --- /dev/null +++ b/3872/CH8/EX8.2/EX8_2.sce @@ -0,0 +1,17 @@ +//Book - Power System: Analysis & Design 5th Edition
+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye
+//Chapter - 8 ; Example 8.2
+//Scilab Version - 6.0.0 ; OS - Windows
+clc;
+clear;
+Ip = [10; 10*(cos(120*%pi/180)+%i*sin(120*%pi/180)); 10*(cos(-120*%pi/180)+%i*sin(-120*%pi/180))]; //given column vector of phase current in A
+function [Ip1]=phaseshift(x1,x2) //Function for shifting the phase
+ [r theta]=polar(x1);
+ Ip1=r*(cos(theta+x2*%pi/180)+%i*sin(theta+x2*%pi/180));
+endfunction
+I0 = 1*(Ip(1,1)+Ip(2,1)+Ip(3,1))/3; //zero sequence current in A
+I1 = 1*(Ip(1,1)+phaseshift(Ip(2,1),120)+phaseshift(Ip(3,1),240))/3; //positive sequence current in A
+I2 = (Ip(1,1)+phaseshift(Ip(2,1),240)+phaseshift(Ip(3,1),120))/3; //negative sequence current in A
+printf('\nThe zero sequence current V0 = %f A',I0);
+printf('\nThe positive sequence current V1 = %f A',I1);
+printf('\nThe negative sequence current V2 = %f A',I2);
|