write(1,8)
8 format(////,10x)
write(1,*) ' расчет действующих и амплитудных значений токов '
write(1,*) ' графики изменения токов '
write(1,*) ' '
write(1,9) idad,id1,idbc,id2,id3,idcd
9 format(////10x,' действующее значение тока в цепи ad idad=',f13.5,' (a)',/
*10x,' действующее значение тока в ветви 1 id1=',f13.5,' (a)'/
*10x,' действующее значение тока в ветви bc idbc=',f13.5,' (a)'/
*10x,' действующее значение тока в ветви 2 id2=',f13.5,' (a)'/
*10x,' действующее значение тока в цепи 3 id3=',f13.5,' (a)'/
*10x,' действующее значение тока в цепи cd idcd='
,f13.5,' (a)'/)
write(1,55)imad,im1,imbc,im2,im3,imcd
55 format(////10x,'амплитудное значение тока в цепи ad imad=',f13.5,' (a)'/
*10x,' амплитудное значение тока в ветви 1 im1=',f13.5,' (a)'/
*10x,' амплитудное значение тока в ветви bc imbc=',f13.5,' (a)'/
*10x,' амплитудное значение тока в ветви 2 im2=',f13.5,' (a)'/
*10x,' амплитудное значение тока в ветви 3 im3=',f13.5,' (a)'/
*10x,' амплитудное значение тока в ветви сd imcd=',f13.5,' (a)'/)
write(1,10)
10 format(////,2x,65('-'))
write(1,*) ' iad=@ , i1=# , ibc=$ , i2=%
*, i3=^ , icd=& '
call grafik(miad,mi1,mibc,mi2,mi3,micd,100,70,1)
close(unit=1)
! расчет действующих значений на зажимах цепи и на реактивных элементах
12 tt=0
l=0
do 13 l=1,60
tt=tt+dt
muad(l)=umad*sin(w*tt+fiad)
muc2(l)=umc2*sin(w*tt+fi2-pi/2.)
13 mul3(l)=uml3*sin(w*tt+ficd+pi/2.)
call integral(muad,t,60,td)
udad=td
call integral(muc2,t,60,td)
udc2=td
call integral(mul3,t,60,td)
udl3=td
if(k.ne.0) go to 18
open(unit=2,file='sas2.res')
write(2,14)
14 format(////,10x)
write(2,*) 'расчет действующих значений напряжения на c2 и l3 '
write(2,*) 'графики изменения напряжения на элементах c2 è l3 '
write(2,*) ' '
write(2,15) udad,udc2,udl3
15 format(10x,’действующее значение напряжения в цепи ad,b udad=',f12.3,/
*10x,'действующее значение напряжения на конденсаторе c2,b udc2=',f12.3,/
*10x,'действующее значение напряжения на катушке l3,b udl3=',f12.3,/)
write(2,17)
17 format(///,2x,64('-'))
write(2,*) ' '
write(2,*) ' udad=# , udc2=% , udl3=& '
call grafik(muad,muad,muc2,muc2,mul3,mul3,60,70,2)
close(unit=2)
! расчет действующих значений напряжений на элементах c2 и l3
18 k=k+1
if(f.eq.200) go to 26
if(k.ne.1) go to 23
open(unit=3,file='sas3.res')
write(3,19)
19 format(////,10x,' таблица действующих значений параметров c2,l3,i2d,
* idad ' )
write(3,*) ' '
write(3,20)
20 format(2x,72('-'))
write(3,21)
21 format(2x,'i k i umad,b i udc2,b i udl3,b i idc2
*,а i idl3,а i')
write(3,22)
22 format(2x,72('-'))
23 write(3,24) k, umad,udc2,udl3,id2,idcd
24 format(2x,'i ',i2,' i ',f10.3,' i ',f10.3,' i ',f10.3,' i ',f10.3,
*' i ',f11.3,' i', /)
write(3,25)
25 format(2x,72('-'))
26 if(k.eq.26) go to 29
b=b+1
if(k.gt.26) go to 27
umad=umad+(q2-q1)/25.
go to 28
27 umad=umad+(q2-q1)/25.
28 uudc2(b)=udc2
uudad(b)=udad
uudl3(b)=udl3
iidc2(b)=id2
iidl3(b)=idcd
if(k.eq.50) go to 31
go to 1
29 b=0
umad=640.
write(3,30)
30 format(2x,72('-'))
write(3,*) ' '
write(3,*)' udc2=#, udl3=%, , idl3=^, idc2=& '
call grafik(uudc2,uudc2,uudl3,uudl3,iidl3,iidc2,25,70,3)
close(unit=3)
! скачок частоты со 111 гц до 200 гц
f=200
go to 1
31 open(unit=4,file='sas4.res')
write(4,32)
32 format(///,10x)
write(4,*)' действующие значения параметров на с2 и l3 '
write(4,*) ' udl3=#, udc2=%, idc2=^, idl3=& '
call grafik(uudl3,uudl3,uudc2,uudc2,iidc2,iidl3,25,70,4)
close(unit=4)
stop
end
! подпрограмма integral
subroutine integral(m,t,n,td)
real m
dimension m(n)
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.