040

40 Rozdział 4

if nargin < 2

Zodb=10+20i; % domyślna wartość impedancji odbioru

end

if nargin < 1

Esk=400/sqrt(3); % domyślna wartość skuteczna sem zrodla

end

j=sqrt(-l); % definicja operatora części urojonej w=2*pi*50; % pulsacja prądu 3-fazowego % jednostki podstawowe Up=Esk; Zp=abs(Zodb); Ip=Up/Zp;

%

za=l+j*l; % obciążenie 100% zb=0.8 +j*0.8; % obciążenie 80% zc=0.5; % obciążenie rezystancyjne 50% zn=Zn/Zp;

% obliczenia prądów i napiec Ea = sqrt(2);

Eb = sqrt(2)*exp(j*( -2*pi/3) )

Ec = sqrt(2)*exp(j*(    2*pi/3)    )

ya=l/za; yb=l/zb; yc=l/zc; yn=l/zn;

Un = (ya*Ea+yb*Eb+yc*Ec)/(ya+yb+yc+yn);

Ua=Ea-Un; Ub=Eb-Un; Uc=Ec-Un;

Ia=ya*Ua; Ib=yb*Ub; Ic=yc*Uc;

% wykresy wektorowe U=[Ua Ub Uc]

compass(U), title('napięcia fazowe¹); disp(' Dalej? Wciśnij ENTER¹); pause;

I=[Ia Ib Ic]

compass(I), title('prądy fazowe'); disp(' Dalej? Wciśnij ENTER'); pause; disp('prądy fazowe');

disp (	la,	A '	; Ia*Ip
disp (	Ib,	A ’	; Ib*Ip
disp (	Ic,	A '	; Ic*Ip
disp (	napięcia		fazowe'
disp (	Ua,	V '	; Ua*Up
disp (	Ub,	V '	; Ub*Up
disp (	Uc,	V '	; Uc*Up

U=U*Up; % wartości w V

1=1*Ip; % wartości w A

disp('KONIEC funkcji []=o3f()’);

return