Разработка алгоритма и программы расчёта электрической цепи, страница 9

13

24

25

12

zefd=real(zef)          zadd=real(zad)

zefm=aimag(zef)    zadm=aimag(zad)

z3d=real(z3)         z4d=real(z4)

z3m=aimag(z3)    z4m=aimag(z4)

z1d=real(z1)         z2d=Real(z2)   

z1m=aimag(z1)    z2m=аimag(z2)

Xp=XL3-XC4      z4=cmplx(r1,xp)

zef   =z1*z2*z3/(z2*z3+z3*z1+z1*z2)

z1=cmplx(R2,-XC1); z3=cmplx(R5,0.)    

z2=cmplx(R4,XL2-XC2) 

Xc1=1./(w*C1*1.E-6);    XC2=1./(w*

*C2*1.E-6);     Xc4=1./(w*c4*1.E-6)     

XL2=w*L2*1.E-3    

XL3=W*L3*1.E-3

Xc1=1./(w*c1*1.E-6)

w=2.*pi*F2

T=1./f2

L=0

pri=0

I4m=umef/z3mod

FAZA

UmEF=zefmod*I1m

zadmod=sqrt(zadd**2+zadm**2)

Fi3=0

UmL3=I1m*XL3

UmC1=I2m*XC1

I2m=Umef/z1mod

I3m=UmEf/z2mod

I1m=Umad/zadmod

zefmod=sqrt(zefd**2+zefm**2)

z3mod=sqrt(z3d**2+z3m**2)

z2mod=sqrt(z2d**2+z2m**2)

z1mod=sqrt(z1d**2+z1m**2)

Fi2=ff

Fi1=ff

      FAZA

PR

MguL3(i)=UmL3*sin(w*t2-fiad-pi/2.)

MguC1(i)=UmC1*sin(w*t2+fi1+pi/2.)

mgi1(i)=i1m*sin(w*t2+Fiad)

mguad(i)=Umad*sin(w*t2-fiad)

mgi4(i)=I4m*sin(w*t2+Fi3)

mgi2(i)=I2m*sin(w*t2+Fi1)

mgi3(i)=I3m*sin(w*t2+Fi2)

t2=0

dt=t/n

fiad=ff

FAZA