1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
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);
|
