Предыдущая страница реферата | 8  9  10  11  12  13  14  15  16  17  18 | Следующая страница реферата

Для того, чтобы иметь возможность хотя бы качественно, но быстро, оценивать правильность работы модели необходимо осуществить визуализацию решения. Поэтому в модуле вывода на каждом шаге вычислим входную и выходную переменные усилителя и заполним этими данными очередной столбец массива D.
В этот же столбец запишем текущие значения времени Т . Массив D передадим через общую область в главный модуль, а оттуда подпрограмме построения графиков KRIS. В автогенераторе некоторое время длится процесс самовозбуждения. Нас интересует процесс установившихся колебаний, поэтому запись данных в массив будем делать только начиная с некоторого момента времени TZ. Эта величина и счетчик точек также включим в общую область.

SUBROUTINE PRIN(T,Z,F,N,IER)
С Подпрограмма вывода результатов шага.

DIMENSION Z(N*4),F(N*4),D(4,15000)

COMMON K,TZ,TAU,D

IF(T.GE.TZ)THEN

K=K+1
С Вычисление значения переменной входа U1.

CALL US(2,Z(3),U1)
C Вычисление значения переменной выхода U2.

CALL US(3,Z(3),U2)
С Заполнение массива.

D(1,K)=T
С Выход усилителя будет изображаться на графиках кривой номер 1.

D(2,K)=U2
С Вход усилителя будет изображаться на графиках кривой номер 2.

D(3,K)=U1

ENDIF

RETURN

END

13 Главный модуль решения системы уравнений

В главном модуле в соответствие с требованиями подпрограммы метода
Рунге - Кутта ARK объявим массивы для решения системы третьего порядка.
Имена массивов сохраним такими же, как имена формальных параметров подпрограммы ARK. Зададим нулевые начальные условия и равные для всех интегральных переменных весовые коэффициенты погрешности. Из исходного файла будем вводить: время начала записи данных в выходной массив TZ , параметр [pic] , время начала интегрирования ТN, время конца интегрирования ТК, максимальный шаг интегрирования НМ задаваемую погрешность ЕР.

DIMENSION Z(12),RAB(9),F(12),D(4,15000)
С Главный модуль решения системы уравнений

EXTERNAL FUN,PRIN

COMMON K,TZ,TAU,D
С Задание начальных условий и весовых коэффициентов погрешности.

DO 1 K=1,3

Z(K)=0.

1 Z(K+3)=0.33333

READ*,TZ,TAU,TN,TK,HM,EP

K=0
С Решение системы.

CALL ARK(HM,TN,TK,Z,RAB,F,3,EP,FUN,PRIN,IER)
С Вывод результатов в форме графиков и таблиц.

Рекомендуем скачать другие рефераты по теме: реферат на тему работа курсовые работы, изложение 4.





Предыдущая страница реферата | 8  9  10  11  12  13  14  15  16  17  18 | Следующая страница реферата