3.5 Алгоритм подпрограммы вычисления действующих значений методом трапеций. Подпрограмма TRAP
1 Задаются параметры: амплитуду, частоту, начальную фазу.
2 Вычисляется значение шага интегрирования h.
3 Обнуляются ячейки, куда будет заноситься значение интеграла.
4 Задаётся цикл для расчета интеграла.
5 Вычисляется определенный интеграл S
где
a, b – начальная и конечная точки отрезка интегрирования.
n – количество интервалов.
F(a) – значение ф-ции в начальной точке.
F(b) – значение ф-ции в конечной точке.
F(xi) – значение ф-ции в i-той точке.
6 Вычисляется действующее значение 
.
7 Возврат в головной модуль.
3.6 Алгоритм подпрограммы вычисления мгновенных значений функции. Подпрограмма MV
1 Создаётся массив мгновенных значений функции согласно амплитуде, углу сдвига и синусоидальному закону.
2 Конец (возврат в точку вызова).
4 ТАБЛИЦА ИДЕНТИФИКАТОРОВ
Все обозначения параметров сведем в таблицу идентификаторов (таблица 2).
Таблица 2 ― Таблица идентификаторов
|
Обозначение параметров |
Смысл параметров |
|
1 |
2 |
|
R C L U F T Q1 Q2 X ZSUM Z UM il1M IC1M am FI ZOF ud |
Массив сопротивлений Массив емкостей Массив индуктивностей Напряжение цепи АД Частота Период Начальное значение C1 Конечное значение C1 Массив реактивных сопротивлений Сопротивление участка АД цепи (полное сопротив-ление) Массив комплексных сопротивлений Массив мгновенных значений напряжения АД Массив мгновенных значений тока в L2 Массив мгновенных значений тока в емкости C1 Массив амплитудных значений Массив углов сдвига по фазе Массив масштабных коэффициентов Действующее значение напряжения АД |
Окончание таблицы 2
|
1 |
2 |
|
icd ild ilq icq UBC IL1 IC1Z |
Действующее значение тока в индуктивности L2 Действующее значение тока в емкости C1 Массив действующих значений тока в L2 Массив действующих значений тока в емкости C1 Падение напряжения на участке BC Ток в L2 Ток в С1 |
5 БЛОК СХЕМА ПРОГРАММЫ
5.1 Блок схема головной программы
5.2 Блок схема подпрограммы ARG
5.3 Блок схема подпрограммы GRAFIC
5.4 Блок схема подпрограммы TRAP
5.5 Блок схема подпрограмм MV,MGV
6 ТЕКСТ ПРОГРАММЫ
INTEGER KOL
REAL
R(3),C(2),L,FI(4),am(5),UM(100000),il1M(100000),IC1M(100000),ZOF,ud,icd,ild,ilq(100),icq(100),isumm(100),ISUMD,uBcm(100000),UBCD
COMPLEX U,J,XC(2),XL,Z(4),ISUM,ZSUM,UBC,IL1,IC1
open(unit=3,file='in.txt',status='old')
READ(3,*)r(1),r(2),R(3),C(1),C(2),L,U,F,T,DT,q1,q2
close(3)
open(4,file='out.txt',status='new')
J=(0,1) PI=ACOS(-1.)
DO 1 I=1,2
XC(I)=-J/(2*PI*F*C(i))
1 CONTINUE
UR=int(U)
XL=J*2*PI*F*L
Z(1)=R(1)
Z(2)=1/(1/XC(1)+1/XL)
Z(3)=R(2)
Z(4)=R(3)+XC(2)
DO 2 I=1,4
ZSUM=ZSUM+Z(I)
2 CONTINUE
! OPREDELENEI TOKA V SXEME
ISUM=U/ZSUM
UBC=ISUM*Z(2)
WRITE(4,*)' ISHODNIE DANNYE '
WRITE(4,687)R(1),R(2),R(3),C(1),C(2),L,F,Q1,Q2,U
687 FORMAT(2x,40('_'),/,2X,'I',3X,'R1',3X,'I',3X,'R2',3X,'I',3X,'R3', 3X,'I',3X,'C1',3X,'I',3X,'C2',3X,'I',3X,'L ',3X,'I',3X,'F ',3X,'I',3X,'Q1',3X,'I',3X,'Q2',3X,'I',3X,'U',3X,/,2X,40('_'),/,2X,'I',F8.4,'I',F8.4,'I',F8.4,'I',F8.6,'I',F8.4,'I',F8.4,'I',F8.4,'I',F8.4,'I',F8.4,'I',F8.4,'I')
write(4,*)'tok v cepi I=',isum
write(4,*)'Uad Il Ic ISUM UBC'
write(4,*)' ___________________________________________'
IL1=UBC/XL
IC1=UBC/XC(1)
call arg(IL1,fi(2),am(2))
call arg(IC1,fi(3),am(3))
call arg(isum,fi(4),am(4))
call arg(ubc,fi(5),am(5))
am(1)=real(u)
fi(1)=0
PAUSE
KOL=T/DT
125 do 132 i=1,50
132 write(*,435)
435 format(50x)
WRITE(*,*)'NU CHE??'
write(*,*)'1--vyvod grafikov 2--nahojdenie deystvuyschih snzchenyy'
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.