skrypt

xdel;

Rl=input("rezystancja w omach: ");

XLl=input("reaktancja indukcyjna:    ");

XCl=input("reaktancja pojemnościowa:    ");

Z1=R1+(XL1-XC1)*%i

R2=input("rezystancja w omach: ^rr);

XL2=input(”XL2 w omach: ");

XC2=input("XC2 w omach:

Z2=R2+(XL2-XC2)*%i R3=input("R3 w omach: ");

XL3=input("XL3 w omach: ");

XC3=input("XC3 w omach: ");

Z3=R3+(XL3-XC3)*%i Z23=Z2+Z3 Zab=Z23/(Z2+Z3)

Zglowny=(Zl+Zab)

Zabs=abs(Zglowny)

fi_z=atan(imag(Zglowny)/real(Zgłośmy))

Esk=input("definiowanie wartości skutecznej E: "); fiu=input("definiowanie kata początkowego E: "); fO=input("Podaj częstotliwość sieci: ");

Etl=Esk*cos(fiu)+%i*Esk*sin (fiu)

iglowny=(Esk/Zabs)*cos(fiu-fi_z)+(Esk/Zabs)*sin(fiu-fi_z)*%i fii=atan(imag(iglowny)/real(igłówny))

Uab=Zab - iglowny

i2=Uab/Z2

i3=Uab/Z3

url=Rl*igłówny

ul1=XL1* %i*iglowny

ucł=t-l)*XCł*%i*iglowny

ur2=R2*i2

ul2=XL2*%i*i2

t=linspace(0,0.1,1000);

it=sqrt(2)*abs(igłówny)*sin((2*%pi*f0*t)+fii);

Et=Esk*sqrt(2)*sin((2*%pi*f0*t)+fii) fi=fiu-fii pt=Et. *it;

S=abs(iglowny)*abs(Esk)

Psk=S*cos(fi)

Qsk=S*sin(fi)

xg=[0; real(i2)] yg=[0; imag(i2)] xd=[0; real(i3)] yd=[0; imag(i3)] xl=[0; real(iglowny)] yl=[0; imag(iglowny)]

plot2d4(xd, yd, axesflag=0) plot2d4(xg, yg, axesflag=0) plot2d4(xl, yl, axesflag=0)

ur=[0; real(uli);    real(url)+real(uli);    real(urli+real(ull)+real(ucl); real(urli+real(ull)+real(ucl)+real(Uab)]

ui=[0; imag(ull);    imag(url)+imag(ull);    imag(url)+imag(uli)+imag(ucl); imag(url)+imag(uli)+imag(ucl)+imag(Uab)]