?????????????? ?????? ? ????????? ?????????? ??????? ?????

Подпрограмма расчета переходного процесса
subroutine obrabotka
use obmen
real iL, L, iL0
! Ввод параметров
!'R1= ','R2= ','L= ','C= ','E= ','dt= ','Tsum= '
R1=vhod(1); R2=vhod(2);L=vhod(3);C=vhod(4);E=vhod(5);
! Ввод расчетного шага dt и времени расчета Tsum
dt=vhod(6); Tsum=vhod(7);
! Вычиcление числа расчетных точек Nt
Nt=floor (Tsum/dt); ! вычисление наибольшего целого от деления
if (Nt>nsum) then
print *,'NT= ',nt
CALL DWGMSG ('Предупреждение:|&
&Слишком маленький шаг|устанавливается число расчетных точек|nt=nsum');
Nt=nsum-1
endif
Uc0=0.; iL0=0.; ! Задание начальных условий
t(1)=0.
iLgr(1)=iL0; Ucgr(1)=Uc0
do n=2, Nt ! Задание цикла по n
! Расчетные формулы
iL=(E-Uc0+L/dt*iL0)/(R1+dt/C+L/dt+R2)
Uc=dt/C*iL+Uc0
! присвоение новых начальных значений для напряжения и тока
iL0=iL; Uc0=Uc;
! формирование массивов тока, напряжения и времени для вывода на график
iLgr(n)=100 * iL; Ucgr(n)=Uc; t(n)=t(n-1)+dt;
end do
return
end subroutine obrabotka
154