Исследование RC-генератора синусоидальных колебаний
| Категория реферата: Рефераты по математике
| Теги реферата: класс, контрольные работы по алгебре
| Добавил(а) на сайт: Nizamutdinov.
Предыдущая страница реферата | 2 3 4 5 6 7 8 9 10 11 12 | Следующая страница реферата
Формальными параметрами подпрограммы вывода должны быть параметры
(T,X,F,N,IER), где IER- код ошибки, определяемый в подпрограмме ARK:
IER=0,- ошибки нет;
IER=1,- знак заданного начального шага не соответствует движению от начала
интервала интегрирования к его концу;
IER=2,- начальный шаг или/и длина интервала интегрирования ошибочно заданы
равными нулю;
IER=3,- шаг в процессе счёта стал более чем в 1000 раз меньше начального.
Массивы X и F в подпрограммах правых частей и вывода можно объявлять как одномерные, с регулируемым размером X(N),F(N).
В главном модуле для подпрограммы ARK должны задаваться максимальный
(он же и начальный) шаг интегрирования HM, начало TN и конец TK интервала
интегрирования, а также значение требуемой абсолютной погрешности решения
E.
Подпрограмма ARK вычисляет решение системы и в каждой точке, удовлетворяющей условиям точности, обращается к подпрограмме вывода, передавая ей значения параметров T,X,F,IER. Пользователь может запрограммировать здесь печать необходимых переменных или накопление их в дополнительных массивах для последующей обработки. (В последнем случае дополнительные массивы следует передавать в главный модуль через общую область памяти с помощью оператора COMMON). После возврата из подпрограммы вывода, ARK продолжает вычисление следующей точки решения.
SUBROUTINE ARK(HM,TN,TK,X,R,F,N,E,PRAV,OUT,IER)
C Подпрограмма автоматического выбора шага.
C HM -Задаваемый максимальный шаг.
C TN,TK -Начало и конец отрезка интегрирования.
C N -Порядок системы.
C E -Задаваемое значение абсолютной погрешности.
EXTERNAL PRAV,OUT
C PRAV и OUT имена составляемых Пользователем подпрограмм правых
частей и вывода.
C IER -Выходной код ошибки.
DIMENSION X(N,4),R(N,3),F(N,4)
C Первый столбец массива X должен при входе содержать начальные
условия,
С на выходе в нем содержится решение.
C Второй столбец массива X должен при входе содержать весовые
коэффициенты погрешности.
C Первый столбец массива F должен заполняться вычисляемыми
C в подпрограмме PRAV значениями правых частей системы уравнений.
C Остальные элементы массивов X,R,F -рабочие.
DO 3 K=1,N
3 F(K,4)=X(K,2)
T=TN
HB=2*HM
IER=0
KP=0
KLP=1
CALL PRAV(T,X,F,N)
C Вызов составленной Пользователем подпрограммы правых частей
системы уравнений.
C T -Независимая переменная системы.
IF((TK-TN)*HM)4,5,60
4 IER=1
GO TO 60
5 IER=2
60 CALL OUT(T,X,F,N,IER)
C Вызов составленной Пользователем подпрограммы вывода результатов
шага
IF(IER.NE.0)RETURN
6 H=HB/2
CALL SH(T,H,X,X(1,2),F(1,2),PRAV,N,R)
8 T1=T+H
CALL SH(T1,H,X(1,2),X(1,3),F(1,3),PRAV,N,R)
T2=T+HB
Рекомендуем скачать другие рефераты по теме: реферат на тему работа курсовые работы, изложение 4.
Категории:
Предыдущая страница реферата | 2 3 4 5 6 7 8 9 10 11 12 | Следующая страница реферата