Разработка алгоритма и программы расчета электрической цепи (массив индуктивностей - 39 мГн, напряжение цепи - 140 В, частота - 260 Гц), страница 6

pause

GOTO 125

do 136 i=1,50

136  write(*,435)

898  FORMAT(2x,'I',F8.3,'I',F8.3,'I',/,25('_'))

126  PAUSE

STOP

END

SUBROUTINE GRAFIK(V,W,X,Z,Y,N,K)

DIMENSION V(N),W(N),X(N),Z(N),Y(N),LV(500),LH(500)

DO 35 I=1,K

35  LH(I)='-'

WRITE(*,*)

XMAX=X(1)

XMIN=X(1)

DO 36 I=1,N

XMAX=AMAX1(XMAX,V(I),W(I),X(I),Z(I),Y(I))

36  XMIN=AMIN1(XMIN,V(I),W(I),X(I),Z(I),Y(I))

WRITE(*,37) (LH(I),I=1,K)

37  FORMAT(1X,127A1)

RAS=XMAX-XMIN

N0=INT(K/2)

DO 39 I=1,N

N1=INT((V(I)-XMIN)*K/RAS)+1

N2=INT((W(I)-XMIN)*K/RAS)+1

N3=INT((X(I)-XMIN)*K/RAS)+1

N4=INT((Z(I)-XMIN)*K/RAS)+1

N5=INT((Y(I)-XMIN)*K/RAS)+1

DO 38 L=1,K

38  LV(L)=' '

LV(N0)='I'

LV(N1)='@'

LV(N2)='#'

LV(N3)='$'

LV(N4)='&'

LV(N5)='*'

39  WRITE(*,40) (LV(Q),Q=1,K)

40  FORMAT(1X,127A1)

RETURN

END

subroutine arg(x,fi,ampl)

complex x,j                                       

real xr,xc,fi,ampl

PI=ACOS(-1.)

j=(0,1)

xr=x

xc=-j*(x-xr)

ampl=sqrt(xr**2+xc**2)

IF(XR.EQ.0.AND.XC.LT.0)THEN

FI=-PI/2

GOTO 1

ELSE IF(XR.EQ.0.AND.XC.GT.0) THEN

FI=PI/2

GOTO 1

ELSE IF(XR.EQ.0.AND.XC.EQ.0)THEN

FI=0

GOTO 1

ELSE

fi=atan(XC/XR)

END IF

1    return

end

SUBROUTINE GRAFIK1(V,W,N,K)

DIMENSION V(N),W(N),LV(500),LH(500)

DO 35 I=1,K

35  LH(I)='-'

WRITE(*,*)

XMAX=V(1)

XMIN=V(1)

DO 36 I=1,N

XMAX=AMAX1(XMAX,V(I),W(I))

36  XMIN=AMIN1(XMIN,V(I),W(I))

WRITE(NU,37) (LH(I),I=1,K)

37  FORMAT(1X,127A1)

RAS=XMAX-XMIN

N0=INT(K/2)

DO 39 I=1,N

N1=INT((V(I)-XMIN)*K/RAS)+1

N2=INT((W(I)-XMIN)*K/RAS)+1

DO 38 L=1,K

38  LV(L)=' '

LV(N0)='I'

LV(N1)='@'

LV(N2)='#'

39  WRITE(*,40) (LV(Q),Q=1,K)

40  FORMAT(1X,127A1)

RETURN

END

subroutine mgv(amp,fi,freq,x,n,dt,ZOEF)

real x(n)

pi=acos(-1.)

t=0

do 1 i=1,n

x(i)=ZOEF*amp*sin(2*pi*freq*t+FI)

t=t+dt

1 continue

return

end

subroutine mv(amp,fi,freq,x,n,dt)

real x(n)

pi=acos(-1.)

t=0

do 1 i=1,n

x(i)=amp*sin(2*pi*freq*t+FI)

t=t+dt

1 continue

return

end

!    Метод трапеций

subroutine trap(m,t,n,rez)

real m

dimension m(n)

h=t/n

s=(m(1)**2+m(n)**2)/2

do 11 i=2,n-1

11    s=s+m(i)**2

s=s*h

rez=sqrt(s/t)

return

end

!     Метод Симпсона

subroutine simp(m,T,n,rez)

real m

dimension m(N)

H=T/N

S1=0.

S2=0.

DO 61 I=2,N-1,2

S1=S1+m(I)**2

61  CONTINUE

DO 62 I=3,N-1,2

S2=S2+m(I)**2

62  CONTINUE

rez=(m(1)**2+4.*S1+2.*S2+m(N)**2)*H/3.

rez=SQRT(rez/T)

RETURN

ENd

SUBROUTINE GRAFIK2(V,W,X,Z,Y,N,K)

DIMENSION V(N),W(N),X(N),Z(N),Y(N),LV(500),LH(500)

DO 35 I=1,K

35  LH(I)='-'

WRITE(4,*)

XMAX=X(1)

XMIN=X(1)

DO 36 I=1,N

XMAX=AMAX1(XMAX,V(I),W(I),X(I),Z(I),Y(I))

36  XMIN=AMIN1(XMIN,V(I),W(I),X(I),Z(I),Y(I))

WRITE(4,37) (LH(I),I=1,K)

37  FORMAT(1X,127A1)

RAS=XMAX-XMIN

N0=INT(K/2)

DO 39 I=1,N

N1=INT((V(I)-XMIN)*K/RAS)+1

N2=INT((W(I)-XMIN)*K/RAS)+1

N3=INT((X(I)-XMIN)*K/RAS)+1

N4=INT((Z(I)-XMIN)*K/RAS)+1

N5=INT((Y(I)-XMIN)*K/RAS)+1

DO 38 L=1,K

38  LV(L)=' '

LV(N0)='I'

LV(N1)='@'

LV(N2)='#'

LV(N3)='$'

LV(N4)='&'

LV(N5)='*'

39  WRITE(4,40) (LV(Q),Q=1,K)

40  FORMAT(1X,127A1)

RETURN

END

SUBROUTINE GRAFIK12(V,W,N,K)

DIMENSION V(N),W(N),LV(500),LH(500)

DO 35 I=1,K

35  LH(I)='-'

WRITE(4,*)

XMAX=V(1)

XMIN=V(1)

DO 36 I=1,N

XMAX=AMAX1(XMAX,V(I),W(I))

36  XMIN=AMIN1(XMIN,V(I),W(I))

WRITE(NU,37) (LH(I),I=1,K)

37  FORMAT(1X,127A1)

RAS=XMAX-XMIN

N0=INT(K/2)

DO 39 I=1,N

N1=INT((V(I)-XMIN)*K/RAS)+1

N2=INT((W(I)-XMIN)*K/RAS)+1

DO 38 L=1,K

38  LV(L)=' '

LV(N0)='I'

LV(N1)='@'

LV(N2)='#'

39  WRITE(4,40) (LV(Q),Q=1,K)

40  FORMAT(1X,127A1)

RETURN

END

7 РЕЗУЛЬТАТЫРАСЧЕТА

ISHODNIE DANNYE

________________________________________

I   R1   I   R2   I   R3   I   C1   I   C2   I   L    I   F    I   Q1   I   Q2   I   U

________________________________________

I 48.0000I460.0000I 40.0000I .000032I   .0032I   .0390I270.0000I 40.0000I   .0044I140.0000I

tok v cepi I=      (2.552293E-01,7.909554E-03)

Uad        Il          Ic      ISUM     UBC

___________________________________________

I    .000I    .002I    .010I    .008II   4.205I

-------------------------------------------------------I  23.637I    .013I    .064I    .051II   4.001I