blob: 2c0d0c89bb127434204f232213253e8ab8f3716c (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
|
clc;
p=30*10^6; // rated power of alternator
v=11000; // rated voltage of alternator
zs=0.005+0.70*%i; // p.u synchronous reactance
Ef=1.5; // p.u. excitation EMF
ia=1; // p.u. armature current
vt=1; // p.u. terminal voltage
t1=Ef^2-(real(zs)*ia)^2-(imag(zs)*ia)^2-1;
t2=sqrt((2*ia*real(zs))^2+(2*ia*imag(zs))^2);
t3=atand((2*ia*real(zs))/(2*ia*imag(zs))); // terms needed to find out power factor
pf=cosd(asind(t1/t2)-t3);
printf('Power factor is %f lagging\n',pf);
de=acosd((ia*abs(zs)^2-Ef^2-vt^2)/(-2*Ef*vt));
printf('Load angle is %f degrees',de);
|