1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
|
//to calculate sec voltage magnitude and ph
clc;
a=6000/100; //turn ratio
R1=780;
R2=907;
X1=975;
X2=1075;
Xm=443*1000;
disp('sec open');
//Zb=inf;
V1=6500;
V22=complex(0,Xm)*V1/complex(R1,Xm);
V2=V22/a;
disp(abs(V2),'voltage magnitude(V)');
disp(atand(imag(V2)/real(V2)),'phase(deg)');
disp('when Zb=Rb');
Rb=1;
Rbb=Rb*a^2;
Zm=complex(0,Xm/1000)*Rbb/complex(0,Xm/1000)+Rbb;
R=complex(R1/1000,X1/1000)+Zm;
Vm=Zm*V1/R;
V2=Vm/a;
disp(abs(V2),'voltage magnitude(V)');
disp(atand(imag(V2)/real(V2)),'phase(deg)');
disp('when Zb=jXb');
Rb=complex(0,1);
Rbb=Rb*a^2;
Zm=complex(0,Xm/1000)*Rbb/complex(0,Xm/1000)+Rbb;
R=complex(R1/1000,X1/1000)+Zm;
Vm=Zm*V1/R;
V2=Vm/a;
disp(abs(V2),'voltage magnitude(V)');
disp(atand(imag(V2)/real(V2)),'phase(deg)');
|
