автореферат диссертации по электротехнике, 05.09.05, диссертация на тему:Разработка и исследование метода синтетических схем для вентильных цепей на основе системных методов

Введение.%ч. \

1. Выбор метода анализа вентильных цепей. ю

1.1. Задача анализа вентильных цепей.до

1.2. Обоснование методов реализации основных блоков алгоритма расчета вентильных цепей.

1.3. Методика интегрирования уравнений вентильных цепей. 2-й

1.4. Выводы.

2. Разработка методики использования системных методов для построения синтетических схем.

2.1. Разработка и исследование метода переходных характеристик для двухполюсных синтетических схем. 2.I.I. Формальное построение разностной схемы.4/

2.1.2. Связь метода переходных характеристик с интегралом Дюамеля.

2.2. Использование метода переходных характеристик для многополюсных подсхем.

2.3. Метод переходных характеристик и степенные разложения, системные методы.б

2.4. Выводы.8Z

3. Алгоритм решения некорректных задач.

3.1. Некорректные, задачи, возникающие при построении вентильных'макромоделей, и методы их решения.М

3.2. Разработка тестовых схем для алгоритмов анализа * электрических цепей.ЗД

3.3. Использование тестовых схем для исследования погрешности в решении от введения G-3 и R-Е ветвей.

3.4. Разработка алгоритмов стягивания смежных узлов.НО

3.5. Выводы.Н

4. Реализация разработанных алгоритмов в комплексе программ расчета вентильных схем.Н

4.1. Общая структура комплекса программ РВС.Ив

4.1.1. Общий алгоритм работы.ш

4.1.2. Обработка подсхем постоянной структуры.ш

4.1.3. Обработка подсхем переменной структуры.М .« »

4.1.4. Процедуры ввода, вывода, контроля и редактирования информации.

4.1.5. Процедуры управления состоянием вентилей.J

4.2. Использование комплекса программ для расчета схем вывода мощности из канала МГДГ.'.М

4.2.1. Работа комплекса программ при использовании прямых методов.*

4.2.2. Возможности использования матричного процессора для анализа переходных процессов в инверторе МГД-генератора.

4.3. Использование метода синтетических схем для » расчета цепей с распределенными параметрами.W

4.4. Выводы.

Основные результаты работы.Ю\

Интенсивное развитие методов теоретической электротехники идет, преаде всего, в направлении все более широкого использования математического моделирования для проектирования, исследования и анализа совместной работы сложных электротехнических устройств. Это определяет актуальность исследований, направленных на создание новых методов анализа динавлических режимов электрических цепей с нелинейными элементагли, ориентированных на использование современных Э Ш . Перспективным направлением развития методов анализа вентильных цепей является иерархическое, многоуровневое макромоделирование, при котором на первый план выступает разработка эффективных математических моделей электротехнических устройств. Базой для построения таких моделей мол:ет служить метод синтетических схем, основанный на использовании неявных методов численного интегрирования. Привленение метода синтетических схем для решения задачи анализа вентильных цепей весьма специфично и ставит целый рдц проблем, определяелшх особенностями вентильных цепей. Важнейшими из этих особенностей являются: - большая, чем у цепей постоянной структуры, длительность временного интервала исследования вентильной цепи; - необходимость совместного учета, наряду с медленными, определяеглыми длительностью всего временного интервала исследования цепи процессами, быстро протекающих процессов, способных изменить состояние вентилей и принципиальным образом исказить дальнейшее поведение всей системы в целом; - многократность изменеьшя в процессе расчета структуры анализируемой цепи, при этом ни тип образовавшейся в результате очередной коглмутации структуры, ни время, ни последовательность возникновения различных структур заранее неопределимы, - 5 В связи с этим важнейшей проблемой при численном решении уравнений вентильных цепей является проблема их жесткости. Одним из наиболее значительных результатов на пути ее решения явилась разработка Ракитским Ю.В. системных методов интегрирования, обладаЕь щих значительнывш прешлуществами перед другими методами в быстродействии, точности и устойчивости интегрирования при численном решении жестких систем дифференциальных уравнений. Трудность их использования для анализа сложных цепей связана с резким возрастанием вычислительных затрат с увеличением размерности задачи.Это обстоятельство обуславливает актуальность разработок, направленных на использование системных методов в рамках диакоптического подхода к анализу цепей, для построения на их основе синтетических схем многополюсников.Другой важной проблемой при макромоделировании вентильных цепей произвольной структуры является проблема создания методов решения некорректных задач. Основополагащшли в этом направлении являются работы Деьшрчяна К.С, Бутырина П.А., Шипулиной Н.А., Мустафы Г.М., Короткова Б.А., в которых решение некорректных задач производится на основе топологического анализа графа всей цепи. Продолжительность топологического анализа делает практически невозможным использование этих методов при анализе сложных цепей с большим числом вентилей и ставит задачу дальне1&1его развития методов решения некорректных задач.Актуальными являются также вопросы использования современных средств вычислительной техники, в частности, матричного процессора (ЕС-2345) для повышения эффективности расчетов вентильных цепей. Использование матричного процессора целесообразно лишь при выполнении логически простых операций, связанных со сложением, сравнением и углнояением элементов векторов и матриц. В связи с этим возникает проблема исследования особенностей применимости - 6 итерационных методов решения узловых уравнений синтетических схем, так как использование именно этих методов позволяет в наибольшей степени реализовать преиьщества матричного процессора.Ряд важных объектов электроэнергетики таких, как линии электропередач, кабельные линии, требуют учета распределенности их параметров. Создание синтетических схем таких объектов требует специального исследования ввиду принципиального отличия их свойств от свойств электрических цепей с сосредоточенными параметрами.Целью работы является исследование возможностей и проблем использования системных методов интегрирования для создания эффективного по вычислительным затратам програшлно-алгоритмического аппарата анализа переходных процессов в вентильных схемах произвольной структуры на базе метода синтетических схем.Указанная цель работы предполагает необходимость решения следущих задач: 1. Выбор способа моделирования вентильных элементов и соответствущего ему метода формирования уравнений вентильных цепей, дающего возможность строить на их основе эффективные алгоритмы анализа переходных процессов в вентильных цепях произвольного вида.2. Выявление тех методов численного интегрирования, применение которых целесообразно для интегрирования уравнений вентильных цепей, и выбор наиболее оптимального по вычислительным затратам.3. Исследование и выбор различных стратегий интегрирования.4. Разработка универсального подхода к построению синтетической схемы многополюсника на основе системных методов интегрирования.5. Создание эффективных макромоделей подцепей, состоящих только из вентильных элементов (вентильных макромоделей).6. Создание метода решения некорректных задач, возникащих при анализе вентильных цепей, позволящего не проводить топологический анализ всей схемы после каадой коммутации.7. Разработка класса тестовых схем для исследования и сопоставления различных алгоритмов анализа переходных процессов в электрических цепях и проведение с помощью тестовых схем анализа точности метода решения некорректных задач, 8. Разработка и создание универсальной програмгш анализа вентильных цепей, базирущейся на разработанных методах,и исследование ее эффективности при расчетах переходных процессов в цепях с большим числом вентилей.9. Исследование возможностей использования для расчета вентильных цепей нового вычислительного средства - матричного процессора, путем применения итерационных методов.Ю . Разработка вопросов использования метода синтетических схем для анализа линейных и нелинейных систем с распределенными параметрами и исследование особенностей его программной реализации.Исследования по указанным вопросам излагаются в работе в следущем порядке.Б первой главе формулируется задача анализа вентильной цепи, определяется область эффективного пртленения всей разрабатываемой методики. Исследуются вопросы выбора стратегии интегрирования, проводится критический анализ методов численного интегрирования, выделяются методы, использование которых целесообразно при анализе вентильных цепей.Во второй главе рассматриваются вопросы использования системных методов численного интегрирования для анализа вентильных цепей. Обобщены на случай произвольного многополюсника методы переходных характеристик. Показана глубокая связь метода переход- 8ных характеристик с системшлш методами, основанная на Паде-аппроксимации матрицы состояния многополюсника. Глава завершается построением универсальной синтетической схемы тлногополюсника в базисе узловых напряжений.В третьей главе рассматриваются методы нахоадения решения некорректных задач, не требущие проведения полного топологического анализа графа всей цепи. Точность этих методов исследуется с помощью специально построенных тестовых схем. Приводится также алгоритм получения точного решения для тестовых схем.В четвертой главе рассматриваются вопросы реализации методики анализа вентильных цепей в созданном автором комплексе програшл.Приводится подробное описаьше алгоритмов работы комплекса программ, Приводятся результаты расчетов режшлов пуска различных типов инверторов строящейся 1ЩЭС-500. Исследуются вопросы применения итерационных методов решения систем линейных уравнений к анализу вентильных цепей. Проводится оценка воз!ложностей использования матричного процессора. В заключении главы рассмотрены вопросы построения синтетических схем для систем с распределенными параметрами.На защиту выносятся следующие, полученные в работе; основные результаты: - предложен новый универсальный подход к определению параметров и структуры синтетической схемы многополюсника на основе систегшых методов интегрирования; - разработан новый алгоритм расчета токов и напряжений всех ветвей, включая ветви, содержащие вентили,в вентильных цепях произвольной структуры; - предложен новый класс тестовых схем, позволшщих всесторонне оценивать и сравнивать различные алгоритмы расчета электрических цепей; - 9- разработан новый программно-алгоритглический подход, позволяющий в рамках метода синтетических схем эффективно производить переформирование уравнений вентильной цепи; - исследованы возможности применения нового итерационного метода и новых вычислительных средств к анализу вентильных цепей. -iO'

Основные результаты диссертационной работы докладывались на:

- совместной сессии Научных советов АН ССОР по комплексным проблемам "Методы прямого преобразования тепловой энергии в электрическую", "Теплофизика", "Научные основы электрофизики и электроэнергетики" (Москва, 1981 г.);

- У1 Всесоюзной межвузовской конференции по теории и методам расчета нелинейных цепей и систем (Ташкент, 1982 г.);

- Ш Всесоюзной научно-технической конференции "Проблемы преобразовательной техники" (Киев, 1983 г.);

- Всесоюзном научном семинаре "Автоматизация проектирования электротехнических устройств и систем" (Москва, 1983 г.);

- научно-технической конференции в ЛПИ им.М.И.Калинина (1982 г.);

- научных семинарах кафедр "Теоретические основы электротехники" МЭИ и ЛПИ им.М.И.Калинина.

По материалам диссертации опубликованы следующие работы:

I. Коровкин Н.В., Савинский П.Л. Разработка тестовых задач для алгоритмов анализа электрических цепей. - Изв. вузов. Электромеханика, 1980, № 4, с.361-364.

2. Бутырин П.А., Карташев Е.Н., Коровкин Н.В. Применение метода табличного моделирования для расчета схем вывода мощности из МГД-канала. - Б кн.: Моделирование и расчет магнитных полей и электродинамических усилий в электрических машинах и аппаратах. Омск, 198I, с.14-21.

3. Коровкин Н.В. Средства повышения эффективности программ расчета схем вывода мощности МГД-генератора. - В кн.: Моделирование и расчет магнитных полей и электродинамических усилий в электрических машинах и аппаратах. Омск, 1981, с.22-27.

4. Демирчян К.С., Ракитский Ю.В., Бутырин П.А., Карташев Е.Н., Коровкин Н.В. Проблемы численного моделирования процессов в электрических цепях. - Изв. АН СССР. Энергетика и транспорт, 1982, № 2, с.94-114.

5. Демирчян К.С., Бутырин П.А., Карташев Е.Н., Коровкин Н.В. Математическое моделирование мостовых преобразователей. - Электронное моделирование, 1982, J3 2, с.51-58.

6. Коровкин Н.В. О выборе метода численного интегрирования уравнений электрических схем с переменной структурой. - В кн.: Исследование в области теоретических основ электротехники и инженерной электрофизики. Иваново, 1982, с.61-63.

7. Демирчян К.С., Бутырин П.А., Карташев Е.Н., Коровкин Н.В. Макромоделирование преобразовательных систем применительно к ЮТ-генератору. - В кн.: 71 Всесоюзная конференция по теории и методам расчета нелинейных цепей и систем. Ташкент, 1982, т.1, с.6-8.

8. Коровкин Н.В., Селина Е.Е. Расчет на ЭВМ электрических цепей с распределенными параметрами, зависящими от частоты.

В кн.: Сложные электромагнитные поля и электрические цепи. Уфа, 1982, с.81-89.

9. Коровкин Н.В., Селина Е.Е. Расчет переходных процессов в

-2 Daцепях с индуктивными связями на основе метода дискретных резис-тивных схем. - В кн.: Сложные электромагнитные поля и электрические цепи. Уфа, 1982, с.90-94.

10. Демирчян К.С., Бутырин П.А., Карташев Е.Н., Коровкин Н.В. Организация макромоделирования слогшых вентильных схем. - В кн.: Ш Всесоюзная научно-техническая конференция по проблемам преобразовательной техники. Киев, 1983, т.5, с.39-40.

11. Демирчян К.С., Бутырин П.А., Карташев Е.Н., Коровкин Н.В. Математическое моделирование преобразовательных систем МГД-гене-раторов. - В кн.: Ш Всесоюзная научно-техническая конференция по проблемам преобразовательной техники. Киев, 1983, т.5, с.18-19.

12. Коровкин Н.В., Королева Т.И. Метод генерации тестовых схем для исследования и сравнения алгоритмов, использующих для анализа электрических схем метод табличного моделирования. -Электронное моделирование, 1983, JS 6, с.103-105.

13. Коровкин Н.В., Селина Е.Е. Эффективный алгоритм численного моделирования волноеого процесса в электроэнергетических устройствах. - В кн.: Электромагнитные процессы в электрических машинах и аппаратах. Омск, 1983, с.103-109.

14. Коровкин Н.В., Королева Т.И. Метод тестирования численных методик, используемых в расчетах электроэнергетических устройств. - В кн.: Электромагнитные процессы в электрических машинах и аппаратах. Омск, 1983, с.109-113.

1. А.с. 479127 (СССР). Устройство для моделирования нелинейной характеристики индуктивной катушки / Виноградов С.Е., Ницен-ко Е.М. - Опубл. в Б.И., 1975, Га 28.

2. А.с. 580562 (СССР). Устройство для моделирования переменной емкости / Виноградов С.Е. Опубл. в Б.И., 1977, № 42.

3. Костенко М.В., Ефимов Б.В., Зархи Н.И., Гумерова Н.И. Анализ надежности грозозащиты. Л.: Наука, Ленингр. отд-ние, 1981. - 126 с.

4. Костенко М.П., Нейман Л.Р., Бладзевич Т.Н. Электромагнитные процессы в системах с мощными выпрямительными установками. -М.-Л.: АН СССР, 1949. 107 с.

5. Поссе А.В. Схемы и режимы электропередач постоянного тока. Л.: Энергия, Ленингр. отд-ние, 1973. - 303 с.

6. Разработка мероприятий по обеспечению заданогокачества электрической энергии с экспериментальной проверкой на У-25, № гос.per. 77073604, Отчет о НИР, 1980, с.157.

7. Пищиков С.И., Сергеенков Б.Н., Шейнкман B.C. Многоэлементная инверторная установка для фарадеевского МГД-генератора с секционированными электродами. В сб.: Магнитогидродинамический метод получения электроэнергии. М., Энергия, 1968, с.100-112.

8. Бутырин П.А. Разработка принципов макромоделирования вентильных цепей на основе методов решения некорректных задач. -Дисс. . канд.техн. наук. Л., 1980. - 236 с.

9. Демирчян К.С., Бутырин П.А., Карташев Е.Н., Коровкин Н.В. Математическое моделирование мостовых преобразователей. Электронное моделирование, 1982, № 2, с.51-58.

10. Дикур Д.П. Метод моделирования на ЦВМ вентильных преобразовательных схем. Л.: Изв. НШПТ, 1970, сб.16, с.46-53.

11. Дижур Д.п. Метод моделирования на ЦВМ вентильных преобразовательных систем. Л.: Изв. НИИПТ, 1970, сб.17, с.44-50.

12. Кучеренко Н.Е., Толстосумов А.С. Алгоритмы и программы анализа вентильных преобразователей. В сб.: Автоматизация проектирования в электронике, 1979, вып.19, с.15-24.

13. HayJ.L.j ИLncjOiCLtti N. Q. ^LncunU б1пгих.(оЛогъ of rriLitti -- conire.tteb fiydc systems dig it а-t comfLutez, . J£££ Conf.

14. Rec. Souse'c. Ing. Comp-ui:. Jlppeic., Cermet , CotozaA) , -J969.

15. Пеиг. Уоы . N. У. , 1969 . p.p. S/2-S3S

16. Акбаров Г.А. Разработка алгоритмов и программ анализа переходных процессов в вентильных преобразователях. Дисс. . канд. техн. наук. - Ташкент, 1980. - 231 с.

17. Нейман Л.Р., Демирчян К.С. Теоретические основы электротехники. Л.: Энергоиздат, т.2, 1981. - 415 с.

18. Круг К.А. Электромагнитные процессы в установках с управляемыми ртутными вентилями. М.-Л.: ОНТИ, 1953. - 115 с.

19. И. Pack-, с. Re&QU.£et, J. Sch-опек.- Simulation Jhitomatijiie de aorvLrtbtuss-eLL'Ls static^nes slvl Ozdinateui. ^fiio^zcLmme SfiCSO, Optimisation" RafbfwtZ J) <3 R S T N ?4. Я D&4 a

20. H. Fach , С. Be6ou£et, J. Schonen.- Л method, ofjinvLuPatiotv of thy zistor static coriirzztetJ.3JS JJrtnuoX meeting 12 th Los ЛпдеШ S9??

21. Con-j-ezen-ce Record, p.p. JJJ4.

22. Дижур Д.П. Вычислительная устойчивость и погрешность цифровых моделей преобразователей. Л»: Труды НИИПТ, 1976, сб.23, с.92-108.

23. Иванов В.К. 0 линейно некорректных задачах. Докл. АН СССР, 1962, т. 145, №2, с.270-277.

24. Тихонов А.Н. 0 решении некорректно поставленных задач и методе регуляризации. Докл. АН СССР, 1963, т.151, JS 3, с.501--504.

25. Бутырин П.А., Коротков Б.А., Шипулина Н.А. Особые случаи расчета цепей с ключевыми элементами. Изв. вузов. Электромеханика, 1978, 1й 12, с. 1380-1382.

26. Шарапов И.М., Мустафа Г.М., Конев Ф.Б. Алгоритм моделирования на ЦВМ одного класса вентильных схем. Изв. вузов. Электромеханика, 1977, №6, с.660-668.

27. Демирчян К.С., Бутырин П.А. Теоремы эквивалентности решений одного класса некорректных задач теории электрических цепей. Электронное моделирование, 1981, № I, с.14.

28. Чахмахсазян Е.А., Бармаков Ю.Н., Гольденберг А.Э. Машинный анализ интегральных схем. М.: Советское радио, 1974. -272 с.

29. Сигорский В.П., Петренко А.И. Алгоритмы анализа электронных схем. М.: Советское радио, 1976. - 608 с.

30. Петренко А.И., Власов А.И., Тимченко А.П. Табличные методы моделирования электронных схем на ЦВМ. Киев: Вища школа,1977. 188 с.

31. Чуа Л.О., Лин Пен-Мин. Машинный анализ электронных схем. М.: Энергия, 1980. - 640 с.

32. Демирчян К.С., Ракитский Ю.В., Бутырин П.А., Карташев Е.Н., Коровкин Н.В. Проблемы численного моделирования процессов в электрических цепях. Изв. АН СССР. Энергетика и транспорт, 1982, № 2, с.94-114.

33. Шакиров М.А. Преобразования и диакоптика электрических цепей. Л.: ЛГУ им.А.А.Жданова, 1980. - 185 с.

34. Беллерт С., Возняцки Г. Анализ и синтез электрических це-* tпей методом структурных чисел. М.: Мир, 1972. - 434 с.

35. Бутырин П.А., Коротков Б.А., Попков Е.Н. Формирование tуравнений электромагнитных процессов цепей переменной топологической структуры. Л.: Труды ЛПИ, 1977, $ 357, с.66-72.t

36. Коротков Б.А., Попков Е.Н. Цифровая модель графа и некоIторые ее применения к решению электротехнических задач. Электронное моделирование, 1982, № I, с.86-92.

37. Попков Е.Н. Методика формирования на ЭВМ уравнений пере/ходных процессов и исследование электрической системы с переменной структурой. Дисс. . канд.техн.наук. - Л., 1982. - 176 с.

38. Нейман Л.Р., Демирчян К.С. Теоретические основы электро-tтехники. Л.: Энергоиздат, т.1, 198I. - 533 с.. * *

39. Бахвалов Н.С. Численные методы. М.: Наука, т.1, 1975. -*631 с.

40. Хемминг Р.В. Численные методы. М.: Наука, 1972. - 296 с., .

41. Сеа-г С. IV. NununeiLccL£ initial. irafue fii&Her^ in otdlaciTij. clLJjeietxtlaC cyLtaXioru. ~ &ncj?e.uroocL QtLfj-s, N. Y. : Szentccd- Hate, Inc., i9 9 J

42. Современные численные методы решения обыкновенных дифференциальных уравнений / Под ред. Дж.Холла, Дж.Уатта. М.: Мир, 1979. - 312 с.

43. Ильин В.Н. Основы автоматизации схемотехнического проектирования. М.: Энергия, 1979. - 380 с.

44. Моделирование и оптимизация на ЭВМ радиоэлектронных*устройств / Под ред. Бененсона Э.М. М.: Радио и связь, 1981. - 271 с.

45. Белов Б.И., Норенков И.П. Расчет электронных схем на ЭЦВМ. М.: Машиностроение, 1971. - 143 с.

46. Анисимов Б.В., Белов Б.И., Норенков И.П. Машинный расчет элементов ЭВМ. М.: Высшая школа, 1976. - 336 с.

47. Schlckm.cuv И., 3n,U£jtolLon, •System. of a, fr-ontinaxi tletuTDzji Jfna,tysLs ramme, IEEE Зга,ги. on, CitCLLit Якету } p.p. 398 386 - flttcj. 19?0,

48. Конев Ф.Б., Ярлыкова H.E. Методы численного решения систем дифференциальных уравнений, применяемых в цифровых моделях вентильных преобразователей. М.: Информэлектро, 1978. - 49 с.

49. Карташев Е.Н.,Разработка и исследование алгоритмов анализа электрических цепей на базе системных методов. Дисс. . канд.техн. наук. - Л., 1977. - 342 с.

50. Luucjcz. IV. ft c\lie,tiori fob fl- stc^S-LUt^ of Иаеси, nu,ltistef Inte^toJlon- fotmu,ta.e Convfu-Uny ,, 3, -196S /p.p. SS0-2S550. 2)ak^LList Л ffveciae. s±ccH£Lty firofrtern /or lin-EO^t rbu,£Ustefv methods. /963. /V3. 6IT p. 27-^9

51. Ek£e fi. I. On, Shde aff-rooumatLons to the, &x.fLone.ntLa.L and. Л Sta-Ste me.lh.ocLs fa г the пилыт-е. tica.t Sotuilon, of- LtiitLaC- uaiue. ffLoblems. - Re&. Pefit. CSR.R 20i0. Vnv. of Ufatettoo

52. EncjtcuvcL R. бьгоъ елИпъаЛ^ foz Rungz- JtuMa. гг/Де Solution to Systzrm of O.D.E.'s.- Cotn/tut. CJ., v ■ 166-169.

53. Ракитский Ю.В., Устинов C.M., Чернорудский И.Г. Численные методы решения жестких систем. М.: Наука, 1979. - 208 с.

54. ZPwthbo М., RoUuon / On the st&UtLty сиъ cuttLLtcucjj of one stefv methods fox. Sotuiruj stiff system of otdLnauj dJ-ffeben-Ua^ eyLLatiotu. - fried:fи £om{x. , p. iA5- IffZ

55. SkcLrrbfbLntL L.F. , Waits И.А. Jhtocn Cm^duLt one-ste/t methods. tTlalh,. Comfi. 1969. V. £3 , 73J-?4o

56. Sha.rn.fUne L. F. , Walts H. A. fi ~ sta£U BtocM. imfvtidi One- stefv methods.- В IT , 19 72 , v 12 , гбг-Лбб

57. Lcun-6ett J. J). Comjuxlcdjlonat methods in DirtincLiy cUffezentLat ec^uMJ-Lotu . J. Wi£ey. : /V. У. , J 9 73

58. Ракитский 10.В. Методика последовательного увеличенияjniaraj-* . интегрирования систем обыкновенных дифференциальныхуравнений. ДАН СССР, 1972, т.207, В 4.

59. Ракитский Ю.В. Новые численные методы решения систем обыкновенных дифференциальных и разностных уравнений. Труды

60. ЛПИ им.М.И.Калинина, 1973, IS 332, с.87-93.* * »

61. Ракитский Ю.В., Кириллова Л.К., ЗимницкиЁ В.А. Новые методы численного расчета переходных процессов в колебательных системах. Прикладная механика, 1971, т.10, вып.7, с.32-38.

62. Дешрчян К.С., Волков В.М., Карташев Е.Н. Сравнительный анализ методов численного интегрирования при расчете переходных процессов в электрических цепях. Электричество, 1976, № 9.

63. Демирчян К.С., Ракитский Ю.В. Новые методы расчета переходных процессов в электрических цепях. Электроника и моделирование, Наукова думка, 1976, вып.2.

64. Люк Ю. Специальные функции и их аппроксимации. М.: Мир,1980. 608 с.

65. Хованский А.Н. Приложение цепных дробей и их обобщений к вопросам приближенного анализа. М.: Гостехиздат, 1956. - 203 с.

66. Коровкин Н.В. 0 выборе метода численного интегрирования уравнений электрических схем с переменной структурой. В кн.: Исследование в области теоретических основ электротехники и инженерной электрофизики, Иваново, 1982, с.61-63.

67. Коровкин Н.В., Савинский П.Л. Разработка тестовых задач для алгоритмов анализа электрических цепей. Изв. вузов СССР. Электромеханика, 1980, с.361-365.

68. Коровкин Н.В., Королева Т.И. Метод генерации тестовых схем для исследования и сравнения алгоритмов, использущих для анализа электрических схем метод табличного моделирования. -Электронное моделирование, 1983, В 6, с.103-105.

69. Коровкин Н.В., Королева Т.И. Метод тестирования численных методик, используемых в расчетах электроэнергетических устройств. В кн.: Электромагнитные процессы в электрических машинах и аппаратах, Омск, 1983, с.109-113.

70. Королева Т.И. Разработка алгоритма и исследование возможностей повышения эффективности комплексного анализа линейных электрических цепей на ЦЕМ. Дисс. . канд.техн. наук. - Л., 1977. - 263 с.

71. Уилкинсон Дж.Х. Алгебраическая проблема собственных значений. М.: Наука, 1970.

72. Зимницкий В.А. 0 моделировании на ЦВМ процессов в прецизионных гироскопических устройствах. Дисс. . канд.техн.наук. - Л., 1975. - 163 с.

73. Парлетт Б. Симметричная проблема собственных значений.

74. М.: Мир, 1983. 382. с. * «

75. Фаддеев Д.К., Кублановская В.Н., Фаддеева В.Н. Линейные алгебраические системы с прямоугольными матрицами. В кн.: Современные численные методы (Материалы мездународной летней школы по численным методам, Киев, 1966). -М.: 1968, вып.1, с. 16-75.

76. Математическое моделирование синхронных и асинхронных вращащихся машин, силовых выпрямителей и инверторов и их систем управления и регулирования. Л.: Отчет по НИР, 1975,9380 . 156 с.

77. Лунц А.Г. Приложение матрично-булевой алгебры к анализу и синтезу релейно-контактных схем. ДАН СССР, 1950, т.70, № 3.

78. Кристофидиес Н. Теория графов. М.: Мир, 1978, с.321.

79. Коровкин Н.В. Средства повышения эффективности программ расчета схем вывода мощности МГД-генератора. В кн.: Моделирование и расчет магнитных полей и электродинамических усилий в электрических машинах и аппаратах, Омск, 1981, с.22-27.

80. Демирчян К.С., Бутырин П.А., Карташев Е.Н., Коровкин Н.В. Организация макромоделирования сложных вентильных схем. В кн.: Ш Всесоюзная научно-техническая конференция по проблемам преобразовательной техники, Киев, 1983, т.5, с.39-40.

81. Демирчян К.С., Бутырин П.А., Карташев Е.Н., Коровкин Н.В. Макромоделирование преобразовательных систем применительно к

82. МГД-генератору. В кн.: У1 Всесоюзная конференция по теории и методам расчета нелинейных цепей и систем, Ташкент, 1982, т.1, с.6-8.

83. Демирчян К.С., Бутырин П.А., Карташев Е.Н., Коровкин Н.В. Математическое моделирование преобразовательных систем МГД-гене-раторов. В кн.: Ш Всесоюзная научно-техническая конференцияпо проблемам преобразовательной техники, Киев, 1983, т.5, с.18-19.

84. Бутырин П.А., Карташев Е.Н. Комплекс программ анализа »электрических цепей с вентильными элементами. П Всесоюзное научно-техническое совещание по САПР САУ и АСУ, Челябинск, 1978, с.100.

85. Апряткина И.Б., Коровкин Н.В. Повышение эффективности программ расчета трехмерных электростатических полей в электрооптических устройствах. Изв. АН СССР, 1983, й 2, с.139-142.

86. Демирчян К.С., Ракитский Ю.В. НоЕые методы оптимизации численных расчетов электрических цепей и полей. Изв.АН СССР. Энергетика и транспорт, 1981, J3 2, с.46-52.

87. Белопольская Е.Б. Реализация двухпараметрического итерационного метода расчета электромагнитных полей. Изв. АН СССР. Энергетика и транспорт, 1981, $2, с.53-56.

88. Тьюарсон Р. Разреженные матрицы. М.: Мир, 1977. - 183 с.

89. Фаддеев Д.К., Фадцеева В.Н. Вычислительные методы линейной алгебры. М.: Физматгиз, I960. - 656 с.

90. Процессор матричный ЕС 2345. Метод доступа. Руководство программиста ШКЧ.990.172 РП1. 1973. 83 с.

91. Процессор матричный ЕС 2345. Принципы работы и функциональное описание. Руководство программиста ШКЧ.990.172 РП. 1979. 144 с.

92. Гусев Ю.П. Установившиеся режимы и электромагнитные переходные процессы в электрической части МГДГ с многомостовой инверторной подстанцией. Дисс. . канд.техн. наук. - М., 1984. - 265 с.

93. Самарский А.А., Николаев Е.С. Методы решения сеточных уравнений. М.: Наука, 19 . - 592 с.

94. Коровкин Н.В., Селина Е.Е. Расчет переходных процессов в цепях с индуктивными связями на основе метода дискретных схем. -В кн.: Сложные электромагнитные поля и электрические цепи, Уфа, 1982, с.90-94.

95. Коровкин Н.В., Селина Е.Е. Расчет на ЭШ электрических цепей с параметрами, зависящими от частоты. В кн.: Сложные электромагнитные поля и электрические цепи, Уфа, 1982, с.81-63.

96. Коровкин Н.В., Селина Е.Е. Эффективный алгоритм численного моделирования волнового процесса в электроэнергетических устройствах. В кн.: Электромагнитные процессы в электрическихмашинах и аппаратах, Омск, 1983, с.103-104.• » •

97. Александров Г.Н. Коронный разряд на линиях электропередачи. М.-Л.: Энергия, 1964. - 227 с.1. АКАДЕМИЯ НАУК СССРт

98. ИНСТИТУТ ВЫСОКИХ ТЕМПЕРАТУР127412. Москва. И-412. Коровинское шоссе. ИВТАН. Телефон 485-83-45. Телетайп 113274 ИВТАН. Телекс 7659 ИВТАН SU.1. На ^ Проректорупо научной работе ЛПИим. М.И.КалининаJ1. АКТ

99. Настоящий акт не является основанием для финансовых расчетов эжду участниками договора.1чальник отдела,к.т.н.г. научный сотрудник, к. т.н.1. Честноков В.В,1. Карташев Е.Н.У

100. О внедрении алгоритмов к программ анализа электрических цепей, разработанными ассистентом кафедры ТОЭ ЛШ имени Ю, Калинина Коровкиннм Н.В.

101. DCL (U26E0ITtL26SAVE|U26KEYe,EDIT0R,DR,F,WL!B.RFRAGiCH0L,iwn6, 0f

102. CALL L26E0IT(AM CALL u26save (a) j CALL L26KEYG (a) j 0(

103. GET L!ST(LEIT,FUR,MET0D) i 0<

104. PUT SKIP EDIT { »LE IT 'iLEIT , rFllR *,fur i '"ETOD'iNETOO j (a (5 ) , 0<a(8).2(a(5),f|6j))i 01

105. GET LIST ( IPED, IPE .NRTS ,NWTS ,N J VEAU »DLIT Ir^P I IvERT , ipR INT I и lEPSZE IT } ! 01 PUT EDIT<'I0tD',lPED,'IPEMPE,'NRTSr,NRTS, 'NWTS'.NKTS, 'N I VE AU «, N I VEAU , 0(

106. DLlTIMP',DUTIf1P,'lVeRT',JVERT,4PR!NT',lPRjNT,'H',Ht 01

107. EPSZEIJ ',EPSZE I T j (SKIP ,s (A , F(2 ) , X ( ?)) , 2 (A , E I 9 , 2) , X<3 ) ) ) | 0(

108. QMEG=3, 1М59Ь25 IlFi=50i 01

109. JPE1=CEIL(lPE/lpEO)+1K0; 01lSCHET = CElL(DI.JT|Mp/(1800f5*h)); 0i

110. AVLCCATE R T s , И у s j 0!call edit0r(leit»nrtsirtslnwts|wt5|nrnax|nhhax|nl,n2|nd|nnn|hmiALLO = h;hnew = h;halbh = h/2;igd(3) = iprlnt| igd(4)=nrmax« tcd(5)cnrts5 0allocate NWi d3 1=1 TO nwtsj 01

111. V>'TsmtU(16)=^C0SD(WTS(I)tU(l6)); 01

112. NW { J, | i ) = 2 j NW(2|I)s4; N'« (3,1 ) =3J N W ( 4 , I } = 5 ; 01

113. CAUL DR ( NRTS , RTS i NW TS i WIS i NViMAX , U , TT ) { 0{

114. END; H = И N j if KOM = l THEN TU = -TT; Fl5E TL = TT i 01

115. MI = 1 I MOD(MI »IPHP)sZ | TL<0 THEN 00 ? 0i

116. F mij>ife1 THEN GO TO mgr i 0i 00 1 = 1 TO I2j GR(f1I J, I) =ABS (WTS ( I) ,T ( 15) ) i ENDS

117. TTT (П! J) =TL; и J J = M j J * l ; E M D ; GO TO BEG; И GR ! CAU F ( I PF. I , GR , TTT ) } 051. OUT: PROC(TG » Mp i 0(

118. DCL TG, tfI | J , J, IS , I IS , К » 0f1. MOD(M I i50Ja0 THEN 0f

119. PtJT SKJP EDI T ( 'ВРЕМЯ/ЭЛ, ГРАД/ Ml ',(1 1 = 1 TO 20)) 0f

120. X(1) ,A(21) ,F(1) ,19 F(5) | i 0{

121. PUT SKIP ED IT(TG,MJ) (F(12 , 4 ) , F (8) ) } 0<

122. J I S = 2 | DO 1 = 1 то NWTSi IS = 0 ; DO J = 1 TO 6| 0f

123. WTS(n,PAR(J)B4'B THEN DO! 0f1.cIS + l i 0<

124. IS = 1 THEM K = IIS i ELSE K = 1i 0<

125. PUT EDIT(NNW(J) ) (F (K) ) i END; 01

126. ENDi IF is = 0 THEN JIS=I!S+5; ELSE I IS = 6^ IS j 0{end; end out; 0 j1. FORM! PROCEDURE! 0t1. комаре THEN DO; N F I = 0 в ; 011. IALL0 = 1 THEN DO; 0(

127. METODrl THEN FREE JAK.UHlLFj 011. f1 ETОDг2 THEN FREE 5PIL; 0(

128. MET0D=3 THEN FREE 0 L D , U H I L F i 0<

129. FREE PRAW, F J КС I OA i END? 0(lALL0=i;DO JSNH6 To NN6 5 wsu( I )=4 » DO J = 2 f 3 I 4 i 0(1. ( WTS ( I ) «PAR(J? 1 ) = 'B)£ (WIS ( I ) ,PAR(J + 2) = Г0 'B» THEN 0( DO; NFIsNFJ+l; WSU ( I )=WSU ( I } F I (NFI , 1 )BKTs( I ) .NNU ( J) } END!

130. ENDJENDI DO 1=1 TO NFI; DO J=i TO NRTSJ 01

131. RTS(J) ,NNU(2)=FI (Ы) THEN DC » FI(I,2)sJ| GO TO ЦК 3 ! END; 0<end ; мкз; end; ilf=nu do ipni to nwmax; do j = i то nfi; 0<

132. I=F1(J,U THEN DOi FlUDsMi GO TO M К 2 ; END! 0<

133. ENDiFILU )-ILFi IL. F s I L F * 1 i MK2 J E^Di N2=lLF«i; nD = N2->NH15 0<

134. ALLOCATE Wj, GRAF, USM, FIKCIJAi 0<

135. CALL WL IB ( NHT5 , WTS ,N J, ?N2 I VEAU , I VERT ,ND , W I ,F IL) J . 0<1.A?. = M1^1B; 01nu=NI; s<do к с i to NWTSi 0<n 1 1 = n 1 1 *w 5 u ( к " p 5 01

136. JPER , IGD (2j=0 } INN=1» IGD( 1 )=KOM | 0(if met0d = 2 i met0d = 3 then sy( j ) ,ns = i ra2 ! 00

137. С ALU RFRAGUGO|PRAH.JAK,RTS,FIL|FIKCIJAiUiH» fj p * »TT ,sy »METOD) ! 00

138. J P E Я = J G D ( 2 )} ENDS 001. ELSE DO! 02

139. N N = 0 i l G D И ) P К Q fj i lGn(2)=lPE«i 00

140. ИЕТ0О = ? I MfcT0D = 3 THEN SY ( l ) , NS= I RftZ f 02

141. CAUL RFRAG(IGD,PRAtt,JAK,RTSlFIl.»FIKClJA,U,H,'PRWf,TT1SY,flETODM 02j PER s J GD(2) { end; 021. METODsl THEN 02

142. CALL CHOL ( JRAZ , IPRIWT , IPEP i INN1 tPRAW ,UHIUF) i ELSE 021. ncT0D=2 THEN DO S 02

143. IP£R"=J THEN IF INN""t THEN GO TO OBR 5 ELSE GO TO TRAP i 02f TT»i,5«H THEN GO TO TRAP! 021. SPIL-r0 rB i NE = 0j 021. DO 1=1 TO JRAZnlj 02

144. SPIL ( l\I ) = M'Bi HEAP = SY ( I I .«EXT! 02

145. START? IF HEAP*=0 THEN DO» 02

146. NS = 5YIHEAP),NS5 SP I L( I tNS) iSPIL<NS»I) = ' 1 rB ! 02

147. NE = NEMJ HEAP = SY (HEAP) .NEXT 5 GO TO START; EN^i 0*

148. END} SPIUIRAZiIRAZfs'l'B? N E E = 0 ! 021. DO 1 = 1 TO IRAZMf 021. DO J= J + 1 TO I RAZ } 02

149. SPIL ( Jt I ) fi M 'B THEM 02do к = I * 1 TO J RAZ } 021. 5PIL<| ,«);'! 'B THEN IF SPIL I J , К)45 ' I 'B THEN DO I 02

150. SPIL(J»K)='1'B; NEE=NEE*H END; 02

151. END; END} END; NEBNE+NEE/2i 02

152. PUT EDIT? fNEs ' i N£} ( A 1 F ( 6 ) ) | 0C1. MR ( 1 ) , К = 3. i 0fi

153. DO 1 = 1 TO IR A 7," i ) NEE = e; . 02

154. DO JrJ + 1 TO IRAZ; 0f IF SPjL(I»J) = flfB THEN D05 NEEaNEE+П UC|K)=J} K=K+JI END; END J 0£1. UR(I+1)=UR{I}*NEEj 02

155. END; UR(IRAZ)=Mf* l 5 0{ TRAPJ IPERaP,lNN?l UR,UC,SPIL HE изменились v ♦/do 1 = 1 то ne; us!j)вйе0; end; a<

156. HEAPsSY(heap) , next 5 go to spurt 5 EnD{ 01end; do 1 = 1 to i r a z 1 old( i )=praw ( i ) ; 0(

157. Do 1=1 to iraz; old(i)=prah{i); 01startu(i).praw(d=startu(i)+uhilf1i'*sy(I).y; 01end; NORM I END 1 0)if !1l<6 then put skip edit(praw) (f(«,l) ) } 010 ' 01

158. DO 1 = 1 TO NRMAXi U( I )= pRASM I ) » END} 0(1. К = 0 ; 0(

159. PO I = N 1 TO N2 *■ NF jt i 0(

160. FIL(I»>3 THEN и( I )sPRAW(F1KCI JA (FlL( I ))) { ЕцsE 0<

161. DO i К = К * i I П =F I (К , 2 M 0f

162. U (RTS ( I I ) , NNU {? ) ) =U ( R T Г, ( I I) , NNU ( 1 ) ) +RTS { I J | ,T (4 ) I

163. RTS ( I!) ,ME3 = tRTS ( I I ) ,UAB | 0(1. END » END i 0f1. WТОК=01 0(1. DO 1=1 TO NRTS; 0{

164. I = ABS(«TS(I) ,NNU(1) ) 1 0f

165. J = abs(RTS(I>,nnu(2) ) ; 0(1.u А В = R T s { I) , L A 0 i 011. IL АВ > Я T MEN 01

166. HI,RTS(I)tT(2)=(U(II)^U(IJ))#RTS(I),T(5URTS<I),T<6i; 0< ELSE DOJRTSU) tLAU = rRTS( I) ,LAB}HI ,RTS( J } ,T<2>=0| GO TO ML I { ENDi 01wtok{ i i?=wtok(i j 0<

167. WTOK(I J)=WTOK(J J)+HI ! 0(

168. ML 1 I IF I LAB5 11 I UAB = 12 THEN IF RTS(I I,N N U{4)=1 THEN 0<

169. RTS( I ).PAR<6)=RTS(П ,T(3)+H#HI/RTS ( П,PAR<3) ; ELSE 01 RTS<n»pAR(6|=RTS(IJ ,1(31*14/(2 *RTS(UiPAR(5)))*<Hl*RTSm,Tm)l 0<

170. ILAB = H THEN IF RTS(I)tNNU(4)=1 THEN 0<

171. РЕГУЛЯТОР ТОКА ПИ^РЕГУЛЯТОР 0!

172. WTs ( >U ( 16) -BETA В УСТ. РЕЖИМЕ WTS()U(16)* -COS (Е5ЕТ A MIN В ГРАД 01

173. WTs Л Т | 16}-СДВИГАЕМАЯ ГРАНИЦА ИМПУЛЬСОВ УПРдВ/|ЕНИЯ«КТ5 U К U I 01

174. WTSOTU5) * МОДУЛЬ ТОКА МОСТА » W Т S ( J Т ( 15) 01

175. WTS(!NNU(8)=ЧИСЛО ШАГОВ С ИМПУЛЬСОМ УПРАВЛЕНИЯ«ИТSUNNU(9) 0'

176. WTS()NNU(9)-НОМЕР ВЕНТиЛЯ, ЖДУЩЕГО И"П, УпрдвЛЕНИЙ^Нts()nnu(9) 0(

177. О Б 4 И Е ГЛ06, ПЕРЕМЕННЫЕ Н Т S ( ) U ( 1 5 ) „ Т ( М И Л /|ИС Е К ) 0'

178. J SCHE Т-ЧИСЛО ШАГОВ» КОТОРОЕ W Т 5 ( ) N N U ( 6 ) » Т О К УСТАВКИ 0<

179. НЕОБХ, СДЕЛАТЬ С НМЛ. УПР, WTS DNNU (7 j irKP+ 1000 0! Н я иА Г ИНТЕГРИРОВАНИЯ FUHR ( »3)гЕрS , КР*ЕРS . ( ИНТ ( ЕРS) )/Т , 01 ТТ-ТЕКУМЕЕВРЕПЯ ш/ 01

180. QS=TT*18E35 IF JP=1 THE" /* РЕГУЛЯТОР ТОКА + / 011. DO 1=1 ТО N*TS? 01

181. ALF=(WTS(I),U(16)r30MWTS(I).T(l5)/5Ke^252r3,fl*wTS(i),U(),6); 0!

182. ALFM5H then ALF = 15PJ 01 IF 1 = 1 THEN PUT EDIT(ALF)(F(B) ) }IF I>1 I I<6 THEN PUT EDtT(ALF) (F(A ) ) | 01

183. CALL LSjENDi IF IP~2 THEN /* ПИ^РЕГУЛЯТоР */ 01 DO 1 = 1 TO NWTS i И J=FLOAT ( WTS ( I ) .NNU(6} ) jAUF = FL0AT (WTS(J) . nnu( 7) )/1000«01

184. FUHR( I i t) = (HI-WTS(J) ,T (\5 ) )/HIi FUHR(J,2)=ALF*FUHR(Iit) I 01

185. FUHR( I |3) = (FUHR( I ,3) * H * F U H R ( I t 1) ) / (WJS( I ) ,U( 1 3 j * 1E » 3 ) I 01hi=fuhr(1,2)+fuhr(1,3); 0!

186. H i > 0 THEN hi ,FUHR(! »3)=0{ 0!

187. WTS ( н . IMP ( JL) = Г0 'B THEN DO j Si

188. ALFsQSn(WTS(I),T<16)+ALF)\ 01

189. WTS(J)tT(J)sWTS(J),T(J)*(WTS(I),T«J*7)»-WTS(I)tT{J|)*C5j 0 0

190. WTS(I».UtJ»"*fTS<n,U(J)*{WTS(!).U<J*7}«WTS<I)fU<0})*PI 00

191. ENDiEND5VST,y = USTMUnUST)+D; ME 7 j £N0 WMKM6I HE? END HOCTt 00i- • v v -> i ep

192. Ч1В! PROC (nwts , Wts,NJ ,N2 ,N J VEAU , I VERT f ND.WI , Fit,) REORDER» 00! DCL, <NWTS,N*,N2fNIVEAU,:VERT,ND.n,I?»!?,n,I,J,IJ,K) BJNFJXEP, 00!

193. BH.WI (*,*)) В I T(8) AL I GNED » 00'

194. Fat*) , I UF ) BIN FIXED. 0011.WTS<*> ALICNEDi 0012 ( Ц А В » NU i NNU !9) , NP ) , $0',2 (PAR' 15) , ItfPU5) ) fUTU) , 00i 2 ITI16I , U( 16) ) !1. DO I=N1 TO N2-4$ 00{

195. WI(I»I + l)fWIU»I+2),WI(I,I+3),Wl(I,U4)» 00<1.( I t I ) s 'B 5 ENP'. DO I4N2-3 TO N2! 00£

196. DO J = N 2" 3 TO N2! W J ( I , J ) ? ff? fB ; END} WJ(!«l>s'l'B} END t 00«00 К s1 TO NWTsi 00f

197. WM6: I 1 spIц(WTS « к »,NNU(1) ) ; I2SFIU(WTS(K),NNU(3)J ; 00$ DO 1 = 1 TO 6? IF WTS(K) , P A. R (I) s ' 0 ' 0 THEN

198. I< = 3 THEN ООП I=FIU (WTS (K) , NNU( 1*1 ) ) U i , Ц ) ( П , I 1 ) ^BjENDi 00«

199. EUSF DO 5 I l4PIt<WTS(K).NNU(1-2)) | W I ( J 2 , I J | , К П I I , I 2 ) = ' 1'B f END) END I 00?

200. F WTS(K) ,PAR(7r-'0'R THEN WI (11,12) ,WI II2,lUn'UB{ 0021.UP = K i 206111! DO K=M TO !2i DO !;Ц TO 121 00?

201. BHsW I(I|K)i 00 j=ll TO 12! 00?

202. PC L LEITUNG CHAR (6) , ULF I J »K , ICN , I , 1Ц INDEX , 02

203. ND.NRTS ,N1 ,N2 ,NWTS,NRMAX ,NWMA>f) BIN FIXED(15,0j, 021. А Я S В U J L T I N i 0 2

204. SYSIN FILE EXTERNAL, SYSPRINT F IU£ EXTERNAL PRjNT, 02

205. H,HALBH,FUR) DEC FLOAT(6) , 021 RTS(*)I 022 (LAB, NUf NNU(4)» NP, PAR(6), T(6)), 021 WTS{*) ALIGNED» 022 (LAB» NUf NNU(9) ( NP) , 022 (PARI 15) | IMP < 15 ) ) В IT ( 1 ) . 022 (T ( 16) , U( 16) ) , 0C

206. NNN(*i*| i II I (0"N W И A X ) CTLI 021 3 ALIGNED, г (LAB,NU , NNU(9 ) ,MP» (PAR ( 15 ) , IMP( 15) ) 8lT(l)i 02

207. T(16) ,U(16) ) , ERR CHAR{ 1 ) ,ML (9 ) LABEL I N IТ CALL Fl ( L E I T U N G ) » 0«on error begin; 02on error system; 02

208. PUT L1ST (LEI TUNG,H,NRTS,NWTS,ND,Nl,N2iNRMAx»NWHAXjI I J»RTS,WTS) ; END! 0«

209. F J I F R О С ( L E J T U N Q ) | 0 2t* Ml ВООЛ WTS И RTS СПИСКОВ 02

210. M2 СОКРАЩЕННЫЙ ВВОД WTS И RTS СПИСКОВ 0{

211. M9 РЕДАКТИРОВАНИЕ RTS И WTS СПИСКОВ 0$

212. МЗ ВЫ90П RTS СПИСКА, M4 ВуВОД WTS СПИСКА */ 0«

213. PCL L E I T liN G CHAR(S) } 0{

214. MUl)=Mli Г11(2}=П9; M L ( 3 ) = M 3 j ML И) iML(5)=M5 } 0t

215. JF LEITUNGs 'Т" = И„OUT ' THEN ML(411 02

216. JF LEITUNGs»Тяй, qyTf THEN DO; M L ! 1 !я П 2 ? ML(4|sM4; END; 0*

217. JF LEITUNG='T = 0, RTS? THEN г1ИП=М2» 0«

218. JF LE I TUNG= 'T"" = » , WTS ' THEN ML(3)sM4J 0£

219. JF L E I T U N G = ' T = (■-•, WTS' THEN DO! ГЦ. ( 1 ) sM2 t ML ( 3 ) =M 4 J END? 0(-IF LEITUNGa »T-*sBg NO? THEN M L ( 3 ) в M 5 i - . - .- 0«1. lEITUNG='T = 0, NOP THEN DO ; МЦ1)=М2; ГЦ.(?)вН5{ END; END FI| 0tfjl; ppoc; err='0", icn=i; 01

220. DO WHILE(ERR"S'*'J ; GET ed!T(ERR){A(1)){ E Np i 0f

221. END FIL v INDEX, ICN = 0! 0<

222. И 7 J INDEXSJNDEX+1 i GO TO ML( INDEX) » 01пз; put skip; do 1=1 то nrts; 01

223. PUT EDJT( '* ' ,I, '* ',RTS(I) ) (Л( 1 ) ,F(?) I A( 1 ) ,7 F(4j,SKIPil2 E(10,2))} 0(1. END; GO TO мв; 0i

224. M4j PUT SKIP; DO I=1 TO NWTS! 01

225. PuT EDI M '* ', I > '* SWTS ! I ) ) ( A ( 1 ) ,F (3 ) , A { 1 I , 12 f ( A j f 2 ( x ( 4 ) , i 5 B(l))i 01

226. X(5),2E(8il),SKJP,30E(8,l))5 E NO J 011. GO TO MbJ 011. M2J DO 1=1 TO NUTS; 011. GET L1ST{IL , ILF) J . 0j

227. IL>IO0 I IL < 0 I IL F > 9 I IL F <2 THEN GO TO MF ? 01

228. W т S( I) ,L А В = I L $ WTS { I ) , MU= ILF ; 0(

229. DO J=1 TO 9 ; q E T L J ST ( I L ) » I F IL<0 THEn GO To flF i WTS ( I J , NNU I J J = IU {END 5 011. GFT L I ST ( IL ) ; 011. !L<1 I IL > 15 THEN GO TO F } 011. WTS(l),Np=JLJ 011. DO J = 1 TO 15; 01

230. WTS(I) ,PAR(J) |WTS( I) ,IMP(J) = 'BiIF J< = IL THEN GET L I ST ( УTS( I ) . PAR ( 0)) ! 01

231. WTS(I),u(j)fWTs(j) ,T(J)=0; 01

232. End; /* СДВИГАЕМАЯ ГРАНИЦА имп, УПРАВЛЕНИЯ, 01параметры пи-регуляторл т, бетд мin */ 02

233. GFT LIST (WTS ( I / ,Т ( 16) ,ViTS ( I ) ,U ( 15) ,WTS ( I ) ,U ( 16) ) J 021. GO ТО MGUT » 02

234. ИГ! PUT EDIT( 'ERROR', 'ERROR'WTS', I ,J) 02

235. SKIP,A,S«|P,A,X(30) .A, 2 F{8) ) { 02 CALL FIL;1. MGиT; END i 021. DO 1=1 ГО NRTSi 021. GET LIST ( IL« ILF) I 021. !L>1^ I IL <0 I 11 F < 2 • J LF>4 THEN CO TO HFR 5 02

236. RTS ( I ) ,LA3 = ILJ RTS(I) ,NU=ILr; 021. DO Jsl TO i; 02

237. J > 11, F THEN R T S ( I ) • N N U ( J ) = 0 | EtSE 02 DO {GET L I ST ( I L) i I F ABS(IL)>1E3 THEN GOTO MFR | RTS ( !) , NNU(J> = IUIEND}END|02

238. GET L I S T j I L ) i IF Ц<0 ! I L >6 THEN GO TO tfFR | R TS ( I j , NP = JL } 02

239. J=1 TO 6J IF J>IL THEN RTS ( I) , PAR(J)=0J ELSE 02

240. GET LIST(RTS( I ) .PAR(J) ) ; RTS I J) ,T(J)=05 END i GO TO MFGUTl 02

241. П F R J PUT t'D IT ( 'ERROR'« 'ERROR 'RTS S',I,J) 022 (SKIP,A) ,X(3P) ,A,2 F(8| ) ? 02call fil; hfgut; end; go to mb; 02

242. Ml! DO 1=1 TO Nrtjs; GET LIST(IL,ILF)i 02

243. I L> 1( !L<Z | IL F > 9 J I L F <2 THEN GO TO f', ERROR J 0Й

244. HT5 ( I ) , LABsIL! WTS ( I ) .NU=ILF i 0Й DO. j=l TO 9 ? GET UST{IL)»IF IL<C THEN GO TO MERROR jWTS { I ) >NNU ( J)-a j 011. END; GET LI s т ( JU 5 (M

245. I L < 1 | IL, > 1 5 THEN GO TO и E RR0 R 5 W T S ( I ) , NP = 11 i DO J=1 TO 1 5 5 0« WTS ( I ) ,PAR ( J) , #( T 5 ( } ) f IMP ( J) = '0 'B i IF J< = I L THEN GET L I ST{WTS( I ) , PAR(J)) 50*

246. END i GET LIST ( (WTS ( I ) t IMP (v)) DO J = 1 TO IL))5 it

247. DO J=1 TO 16 iIF J> J L THEN WTS(I) ,T(J) iWTS{I) ,U(J)=05 Hi

248. F IL>100 I U<v ! ILF <2 | IL F > 4 THEN GO TO MRS . 0t

249. RTS ( I ) ,LAB=IL( RTS(I).NU=ILF 5 0«

250. DO jsl TO 4; IF J>ILF THEN RTS( I ) ,NnU(J)S05 ELSE 0{ DO InET L J 5 T(IL) »IF A0S{ I L)> 1 E 3 THEN GO TO MRJRTS( I j,NNU(J)=ILiEND}END\ dl

251. GET LIS T( I L ! I IF I L < 1 i I Ц > 6 THEN GO TO MR 5 RTS(I) , NP=ILl Ък

252. DO j=i TO 6; IF j>il THEN RtS(j),pAR(j)=0'} 0f

253. ELSE GET LlST(RTS(l()fPAR(J))! END; 0(

254. GET LIST(CRTS(I)fT(J) DO J=1 TO 6))j GO TO П6J 0<

255. MR| PUT EDIT! 'ERROR', 'RTS',I ,0) ( 5 К I P , A , X { 3 0 ) ,A,2 F(8)|} M6| END} 0f

256. GO TO M8} 0( M9| IF lС N s1 THEN DO} PUT ED IT( 'РЕДАКТИРОВАНИЕ HE ПРОИЗВОДИЛОСЬ, ТАК КА0«

257. К БУЛИ ОбНАРУиЕНЦ ОШИБКИ В ИСХОДНОЙ цН <? ОРМ а ц ИИ г) { S К.1 р » x ( 2 0 ) , а ) 5 0<go то ме; end; nrmax ,nwmax = 0; halbh=h/2j 0(

258. Do 1 = 1 TO NWTSj WTS(I) ,NNU(9)=2? Do J=1 TO WyS 1 ? ) , NU i 0(

259. KsWTS ( I ) »NMU (0) } IF К > n v- m A x THEN nwmAX=KJ ends END} 0< DO 1 = 1 TO NRTSipO J=1 TO RTS ! IJ ,NU5 К= RTS ( I) ,NNU(Jj jIF K<0 THEN DO! 0(

260. K?pKJ IF К > N R M A X THEN NRMAX = Ki 0(end; end; end; ni=nrmax*i; n2 = nwmax » nd=n2*n i<m ; 0<

261. I L = i » I IU = 2 THEN RTS { I ) ,PAR(2)e^RTS( I ) .PAR(2) 5 000if i l = 12 THEN RTS( I ) ,PAR(4)s^RTS( I ),PAR {4 ) ; 00?

262. IL = U | Il, = |2 THEN RTS( I ) ,T (3) =^RTS( I ) ,T I3M 0001. i l = 7 | Use THEN RTS(I) ,PAR<3>SnRTS(J).PAR<Jj | END! END! 002

263. DO 1 = 1 TO N R T S J a=RTS(l ) ,LA3! 0001. 1Ц=1 I IJHEN DO! 000

264. RTSn.T(5) = l/RTS(I),PAR(U!RTS(I).T(6)=RTS«lJ,T(5)*RTS{I),PAR(2M 00»

265. END! JF Il, = 5 | IL s 7 THEN DO} 00Я

266. RTS(IbpAR(5)BH/(RTS(I),PAR(2)*H*RTS(I),PAR(U)l 000

267. RTS(I)fPAR(6J=MALBH/(RTS(nlPAR(2UHAl.BH«RTS(I),pAR(Uit END! 0001. jl=h then DO! 000

268. END! NNN = 0; DO laj ТО N W Т $ } 0001. 4TS ( J) ,LAB = 2 THEN NNN(2 , J)=NNN(2 , 1 ) + 1 J 000

269. WTS ( I ) f L,AB = " THEN NNN (2 ,3 ) =NNN (2 ,3 )J 0001. «TSU'iUASaS THEN NNN ( 2 . 5 ) =NNN ( 2 , 5 ) + 1 . 000

270. JF #ts(j»,LAtfai THEN NNN(2,4)=NNN(2 , 4)*Ц 000

271. WTS( I ) ,LA8SJ THEN NNN(2,2)=NNN(2,2) + l! . END! 000

272. JbFcl; DO K = i TO 5! DO I = ILF TO NNN ( 1, К ) * I I.F в }f 000

273. KTS( П «ЬАЧ^зК THEN DO J=IM TO NWTSJ 000

274. WTS ( J) ,U8 = K THEN DO ! 000

275. В= f T 5( I ) » 4TS( I )BWTS(j) ! WTS(J>=Bl GO TO M92! END! 000

276. END! M92J END! UFs UF*NNN ( 1 , К ) ! END! 000 NNN(1, 1)s i ! 00 I=2 i 3 , 4 , 5 J

277. NNN( 1, I )=NNN(2, M>M! NNN(2,I)=NNN(1,I)+NNN(2 , I j„ 1 ! £ND ! 000

278. M8| GO TO M 7 , M5: END EDITOR} 000pro3><nhi ,nni .wts , ^тзк,и,н> ^eorde^; ocl розоз external;

279. OcL ( NlHl ,NN1 , I NO ) BIN FIXED,

280. N>U<#) ,WTOKU) ,UH!i JS.WS.WNEW) , I WTS(*) ALIGNED»2 (la3i nj » mnj ( 9 ) » чр) . 2 (pari 15) ,ир( 15) ) 3it( 1) , 2 T( 16 ) ,2 и (16);1. DO IND-siHi TO NN 1 5

281. WNEW.^TS( IND) ,T(2)=^TDK<WTS( IMD> , I? ) 5и И , WTS < I No ) . u ( 2 ) = j'( ITS ( JN3 j . nnu ( 1 ) ) J ( WTS ( IMQ) . MN J ( 2) ) 51. Us=rfti<IN3).U(l)}

282. Do I = 4 • 5 IF HT9 I I 3 F W I ( I DO I = 7 • 9 IF *TS I I=P W I ( I DO 1=10, IF HTS i I = P wid1. WTS < < WTS ( < WTS< <i f i f i f

283. WTS(K).NMU(П > » I? = FlL(WrS(K),MNU(5)); L ( W T S ( К ) . N MU ( 9 ) ) ; 3 !1. K) .p4r( I 3' THEN OQ 5

284. N ВЕНТИЛЕЙ <АГОаНЭЯ Г=>УПП^t #/со1. Wil:i 1sF fl(wts(к) . n n j ( 1) ) ; i2=fluwts(k) .mnjutsu) .nj) > i1. DO isi TO WTSU).Nj-i;

285. PjT sк I эi PО I т < W I ) ( s.< I (N0)3(1) ) ;1. ivert's» then do;

286. DD I = M 1 To N2; WI |I,!|SM'B! END i L ; DDI < = N1 T 3 N2 ; DO I= N1 TO N2 J1. B^t=J(iJ i l i- . . . .1. D3 J=N1 TO N2J1. кi (i, j) = '0'B thevj1. 3H = '1 * 3 rHHNif wi(k,j)*4'8 then rf'm i . j) s'l 'b ;end L;e^D ;1. ЕчО WL13«;

287. WM12J PROC(NH12,NN12,WTS,WTOK,UiH,ZEIT,STT»STUW) REORDER', DCU I WTSf*) ALIGNED,2 (LAB,NU,NNU(9) , NP t (PAR,IMP) (15) В I T ( 1 ) , ( T , U ) ( 1 6 ) ) , I I NO , I ,K,L,NH12|NN12,X(7) ) BIN FIXEDU5.0),

288. T,H,0,HILF, (WTOKiU) (*) , (ZEIT ,STT,STUW) (*,*) ) DEC F L О A T ( 6 ) i DO I N D s N P I 2 TO NN 31 2 ?

289. DO 1=1 TO 14! /* запись исходных ТОКОВ */

290. STT(I NO,I)sfcTS(I N 0 J,T(I) } /* И НАПРЯЖЕНИЙ ВЕНТИЛЕЙ */ STUW(IND,I)SWTS(I NO J ,U( I ) ! /* (СТАРЫХ) ДЛЯ БУДУЩЕГО +/1. ВЫЧИСЛЕНИЯ ZEIT12 */

291. ВЫЧИСЛЕНИЕ И ЗАПИСЬ В WTS(*),Т(*) НОВЫХ ТОКОВ ВЕНТИЛЕЙ */ WT0K(WTS(IND).NNU(5))=-HT0K(WTS(IND),NNU(1))^WT0K(WTS(IND),NNU(?))rWT0K(VTS(IND).NNU(3))-TWT0K(WTS(IND),NNU(4));1. DO L = 4 , 10; 0 = 0 ?

292. DO I=L,L + 1,L + 2| wTs ( 1ND? , T(J)=-lE20i IF WTS(I NO J.PAR(I) = 'B THEN wтS< I NDI,T(I)=0{ ELSE IF WTS I IND),PAh(1-3)s»0»в THEN

293. L = 4 THEN hTS( I NDI♦T(I I= WTOK(WTS(J ND) ,NNU( I"2) ) |

294. ELSE WTS(IND),T(I|=WTOK(WT$<IND) ,NNU < 1-41);

295. QsD + WTS(IND) , T ( I ) i END; IF D>-5E19 THEN GO TO WH; /* НЕТ СКВОЗНЫХ ГОРЕНИЙ */

296. D<t1,5E20 THEN GO TO wMl}/#ECTb ПАРНЫЕ СКВОЗНЫЕ ГОРЕНИЯ*/ ELSE

297. L = 4 THEN IF WTS( I ND> ♦ PAR( 1 3) = ' 1 'B THEN CO TO KHlj

298. ELSE К = W T S( IND) ,NNU(5) ; ELSE IF WJS { INO) f PAR ( M) s M 'B THEN GO TO W M 1 ;

299. ELSE K = WTS( IND) .NNU (9) }т = w т о к (к);

300. WTS ( IND) ,T (L) S-1E20 THEN DO!

301. TS(IND) ,T (L)=T-WTS ( IND) ,T (LM )-WTS { IND) tT<L*2U GO TO WMJ END! ELSE IF WTS( IND).T(L + U="1E20 THEN DO!

302. WTS ( IND) ,T (L + l ) sT^WTS ( IND) , T (L j-WTS ( IND) , T <L*2 ) » GO TO ни; END', ELSE

303. WTS(IND) .7 (L + 2) =-lE20 THEN DO

304. WTS ( I ND ) . T (1*2) =T-WTS ( IND) ,T (L ).WTS ( I N D ) , T (L + U I CO TO WMi E N D I1. ELSE GO TOw m1 i два или более сквозных горения, возможно с x ( i 3 ) и1и х(н)#/ip l = 4 then if w t s ( ind) .par ( 13) = '0 '8 then x(7.)s0{1. ELSE X(7)52;

305. ELSE IF WTS( IND).PAR ( 14) = '0 '8 THEN X<7)=0J1. ELSE X ! 7)=2 ;

306. DO 1=Lг3 TO L + 2 i IF WTS I I ND1.PAR(I) = '0 '8 THEN X(I-L*4)=0j1. ELSE X ( I rL*4 ) = t\

307. END! H!LFsXll)*Xl4hM2}*XI5!tX(3)*l(|6)*X|7)lif l s 4 then н ilfs(wwtok(wts( ind) .nnu(5) )nx ( 1 )*п-0 . 5*x {4))* wtok(wts(ind),NNU(2))tX(2i*(1-0.5*X(5))+ w t о к ( w t s ( ind) | nn u i 3 ) )чх(3)*(1^0»5*х(6) ) * wtok(wts(Ind),nnu(4)i)/hilf;

308. ELSE HILF=(-WTOK(WTS(lND).NNU(9))^X{t)*(l-0,5*i<(4))* WTOK(WTS( INO) »NNU(6) )"Х(21*!1",0|5*Х{5) ) * WTOK(WTS(IND),NNU(7))'x(3)*Ur0t5*X(6))*

309. WTOK (WTS( INO),NNU(8 J) )/Н ILF i DO !=LiL*'l|L + 2{

310. L = 4 THEN WTS( IND)fT(I!=0,5*WTOK(WTS(IND).NNU(I"2))+HILF;

311. ELSE wts ( IND) ,T ( I ) =0,5«WTOK (WTS ( IND) ,nnu ( I-4 ) )+HH.F ;1. END 5

312. WTS ( IND) t T(L) "WTS (IND? . T ILH )-WTS ( IND) . T (L + 2) i

313. ELSE IF WTS(IND).PAR(14) = r0ffi THEN WTS(IHD),T( 1 4)s0; ELSE WTS(lND)fTU4)=-"WT0K(WTSUND)*NNU(9))r

314. WTS(IND),TU)-WTS<IND},T(LM!9*TSUND),Ta*2!; FND| /* END ЦИКЛА По L /* ВЫЧИСЛЕНИЕ И ЗАПИСЬ В WTS(*),U(*) НОВЫХ НАПРЯЖЕНИЙ ВЕНТ ИЛЕ-1 */ ОО 1=1 ТО 6!

315. I< = 3 THEN /* АНОДНАЯ ГРУППА */ DO;

316. WTS(INO),U(I)sU(WTS(IND),NNU<lj}-U(WTS(IND).NNU(I+l)); WTS(IND),U(I*6)=U(WTS(IN0),NNU(5))'-u(WTS(IND).bNU(l-.5));ENDi

317. ELSE /* КАТОДНАЯ ГРУППА */ DO;

318. WTS(IND).U(I)sU(WTS(IND),NNU(Is2))'-U(WTS(IN0) , NNIJ (5 ) ) { WTS(IND),U(I+6)=U(WTS(IND)lNNU(I*2))BU(WTS(IND).NNU(9));ENDi end;

319. END; /* ЦИКЛА ПО IN D */ END WM12J

320. Wvi: PR03<< mHn.^nn.WTS » «TDK ,U »H,2EITN) HEORHR! TCL PCHOs external;

321. WmE*i»wTS(IND) .Т(1«-в)=й|ТЭК(ИТ5(1 nJD) ,n|NJ{ Г) ) ; U^.4TS<INO>,U(neisJMTS(lND},N\|U(WrS<lND>.*jy)H1. UMTS<IND).N 4U( I )) ;

322. WS=*TS(IND).U( J ) 5 US = </TS ( IND) . T < I ) 5

323. TS ( I ND) , ИР ( I ) ' 1 fB & UH>t> § US«P3 THEN ZEiTN(IND,I)=H#Us/(US-JH) + iE-U.J 1-- JH = 0 & WNEW < э0303 THEN ZEiTN(IND,I)=H*W3/<WMEM~WS)-lE»10;1. Емо WNJ;wA: pro:'(mHa,nna»wts«»(tdk,u«h,zeita) reorder; tcl po^ds external;

324. DCL Z в I т A ( * , * ) DE3> -LOAT ( 5) ; DCL < MHA , NN A , I , I N 3;) BIN F I X£D i

325. H»U{*) . WTOK ( * ) » UH> ff NEW • NS f JS I i i WTS|*J alic^ej»2 (LAB* NJ1MNJ(9) , MP) , 2 ( 9 A R ( 1 5 ) • HP ( i 5 ) ) 9 I T {1) , 2 T( 16 ) , 2 U ( IS ) ;do indpnHa to nna;

326. Do 1 = 1 TO WTS ( IND) . v)U-»l 5

327. WME^»WTS(IND).T(H8)=«/TDK(WTS(IMD),vlNJ(I)); U^.^TS ( IND) .U( Нб) U ( "TS ( IND) . NMU( I) )

328. JMTS( INQ) .N^U(WrS( IND) . MM) ) Ius-#|ts ( inq) . и ( i >; WS=WTS(INO>.T(I)i1. *|TS ( IND) . HP ( I ) = f 1 'B A US<30 & U4M THEN

329. E I T A ( IND, I ) = H * U S / (US"JH) + ie-1J'5 I- JH = 0 4 WNErfOO^OS THEN ZElTA(IND,I)=H»W5/<WME^-rtS)-lE«i0; e>id W a i - .process; 0<

330. U26SAVE! PROC(A) REORDER! 0f

331. DCL A Deo FLOAT,SYSPRINT FILE EXTERNAL PRINT} 0(

332. PUT EDITf'K ПРОГРАММЕ ПОДКЛЮЧЕН РАЗДЕЛ L26SAVE, ' 0('ПРОЦЕДУРЫ? F « ВЫВОД ГРАФИКОВ,' 0i

333. DR ВЫВОД ТАБЛИЦ РЕШЕНИЙ,' 0S' CHOL г РЕШЕНИЕ СИСТЕМ ДИН, АлГ, УРАВНЕНИЙ,'' RFRAG г БИБЛИОТЕКА МАКРОМОДЕ/.Ef\ ВЕТВЕЙ НЕИЗЧЕНЯПЩ0<

334. БйСЯ СТРУКТУРА ; 'Н5 (SKIP,X(30) .А» ) i 0<1. END L26SAVEJ 051. PROCFSS! 01

335. FI PRPC<{РЕ|GRИТТ) REORDER; 0(

336. DCL (IPE,I»J» I R. I I »К) 3 I N Г I ХЕП( 1 5) * 0(

337. SYSPRINT TILE EXTERNAL PR I NT i ( ABS , FLOQR ) В U I L T J N i 0(

338. GR(*|*) .T T T(*)) DEC FLOAT; 0(

339. DCL (MAX * I 2) »MIN(12) ,STEP(12) ) DEC ( 6) FLOAT « 01

340. S(1li) ,S3( 12) ) CHAR(l) ALIGNED, (KLEI N,GROSS,HI I DEC FLOAT(6) } 0t

341. GR( I t J) >ПАХ ( JJ THEN MAX ( 0 ) =GR ( I , J ) 0{

342. DO 1=1 TO 125 02 PUT EDIT (MJN( I ) , ПОАГ ПО КРИВОЛ ',SS(U,' P ABgH f , sT£P ( I ) , M AX ( \ ) j (R ( MF ) ) \ 0i

343. PUT EDIT(S) (SKIP,X(8) , HI A ( 1 ) ) ! 02end ; 02

344. PUT S К IP(2) tPIT( 'ВРЕМЯ') (A(5)) i 021. JF IR>0 THEN DO; 02

345. DO I = 1T0 It; ПАХ ( I ) =!<LE IN ; KLE I N = KLE IN*H I ; END} 02 PUT EDIT((MAX(J) DO 1 = 1 TO U) ) (E(9»2) , E( l П»2) , 8 E (\\, 2 ) , E(8 , UI 5 021. ENDI 021. DO 1=1 TO IPEi 02

346. Ssf '! IF I П>0 THEN DO К = 1 TO 1U BY 11! S(K} = ','} END | 021. DO J:1 TO JZi 02

347. STEP(J!сИ THEN GO TO LU J 02

348. I = FlnOR ( ( G R ( I , J ) M I N ( J ) )/STEP(0j + lfL"5) + l| 021. S(in=ss(J); 02ll! end; hj = 1s000*ttt ( i ) i 02put edit(hi.s) ( sk i p • f ( a , 2 ) ,x(1) ,1u a(l))lend} end f} 02process; 02

349. PR J PROC!NRTS»RTS,NwTS,HTS,NWMAX,U»TT) REORDER; 02

350. DCL (NRTS,NKTS,NWMAX) BIN FIXED! 15) f 02

351. MFi FORMAT(SKIP , 10 E( 1 2 , 4) ) j 02

352. PUT ED Iт( (J DO J = K TO K + 9) MR(flFi) ) » 01

353. PUT E D I т ( (и ( 0 ) DO J = K TO К + 9 ) ) ( R ( M F ) > 5 01к = к * 10 ; 0(1. End; 01

354. PUT EDIT ( * J DO J = K TO K + Ml ) ) (P (MF1 ) ) ! 0<

355. PUT EDIT((U<J) DO J=K TO K*MI)»(R«MFI)» 0'flsFLOOR (NRTS/10) I Ml=MOD(NRTS i 10) 5 K = 0R5 0<1. H|:HMj 0i1. К a I 5 0 (

356. PUT F.DlTt'TOKn ВЕТВЕЙ ') (R (MF2 ) ) 5 0t1. DO J=1 TO П} 0(

357. PUT EDIT( (J DO J = K TO K + 9) ) (R (MF 1 ) ) i 0«

358. PUT EDIT ( (RTS ( J) . T { 1 ) DO J = K TO К + 9-) ) ( R I MF ) ) 5 0gк = К * 1 0 } 0{1. END} 04

359. PUT ED I T( (J DO J = K TO K + M1) MR(MF 1 ) ) i 0«

360. PUT EDIT( (RTS(J)fT( 1) DO JcK TO К+ M 1)) (R(MFIM 0<

361. PUT EDIT ('ТОКИ И НАПРЯЖЕНИЯ НА В EH Т И Л Л X ' ) ( R (fi F 2 ) ) { 0{0{

362. PUT ED I т и I DP 1 = 1 ТО 6)>(SKIP,6 f F < 10) »Х( 10) ) ) } 0£ро i=i то nwts; m=wts(I!.np; 0в

363. PUT EDIТ I I WTS \ I ) fT(J) ,WTS ( 1 ) ,u \ J) DO J=1 TO Ц))

364. SKIP, (2*H)EJ 10,3) ) S 02end;

365. PUT EDIT ( 'СОСТОЯНИЕ ВЕНТИЛЕЙ И ИМПУЛЬСЫ УПРАВЛЕНИЯ О {R(MF2)) | 021. PUT SK I Р I 02ро i=i то nwts; • 02

366. PUT EDIT((WTSuJ,PAR(J) DO J=1 To 6!)(X(9),6 B(lj)5 02

367. PUT EDJT< (ИТБ ( J ) , IMP ( J) DO J = 1 TO 6I)(X(3!,6 8 ( 1 ) ) iENDiEND OR 5 02process; 02

368. RAKITI PROC ( IRAZ ,SY ,STARTU ,PRVI ,HI ) REORDER} 02

369. PCL (I.HEAP , JRAZ. IGL) BIN FlXED(l5,0j, 02

370. SYSPR1NT FILE EXTERNAL PR 1 NT , t MOD • SUM ? BuitTiN. 0ЙfJB,A(,F,NR,AP2,BET,HI , (STARTUjPRW) ( * ) f DEC FL0aT(6), 00

371. AP,H,R,D) ( I RA Z > DEC FLOAT(6 ) С T L I 00

372. S Y (*) , 2 ((NS.NEXT) BIN FIXED(15,0), Y DEC Fi0AT(6)M 00

373. MULT? PRPC(V»w)5 /* Y*Vя>W */ 00

374. DCL (V,W)(*) DEC FL О A T(6) J W = V 5 00

375. DO 1 = 1 TO IRAZmlj HEAP = SY(J),NEXT . IF НЕАРя=0 THEN DO 5 00

376. START! W (SY (HEAPj (NS) =W (SY (HEAP) ,Ns)+SY (HEA.P) ,Y*V ( U 5 00

377. W(i)=VMp+SY(HEAP)(Y*V|SY(HEAP).NS)jHEAP = SY(HEAP),NEXTI 00

378. HEA?" = 0 THEN GO TO START} END J 00end; end hulti allocate ap,h,r,d{ 00

379. CHOLI PROC ( IRaZ,IPRINT , IPER , INN , jak ipraw»UHILP) REoRPERi 00

380. DCL { IRAZ , JPR1NT , IPER , INN) BIN FI XED{ 1 5) , SQRТ BUlLTJN, ' 00

381. PRO-ESS''FЯJAK,K,^,*PRAW!K•, 5 END;pRAWm=Hll,F*UHlLF<b JENDI EN0 QHOL j RFRAG! ^®^||^D|PRAW|jAKjRTS,FlLiFlKCIJAiU|H,STR|TTfSY|f1ET0D) REORDER!

382. DCL ( UD(»i , f i L (*) ,FIkcIJa <*> ) BIN f i xfo ( i s j .

383. ON ERROR BEGIN; PUT SKIP LIST(K,NNUi.NNU2.il TT H ,HR , I , J ,HEAP) i END I /* П E T 0 D = I ПЕТЭО ХОЛЕЦКОГО БЕЗ УПАКОВКИ МАТрИЦЫ,

384. MtTOD ХОЛЕЦКОГО С УПАКОВкОй ПАТРИЦЫ.3 mfctod ракитского с упаковкой матрицы, */ tor! proc( i , j,ylj, jijj ; ' ' pcl (i , j, i 1 ) bin fixeo(15,0) ,

385. YIJ , J I J) DEC FLOAT(6) i \F I>0 THEN PRAM I)=PRAW( n-jlj; IF THEN P4AW(J)=PRAW(J)*jrj{if sTRr'prw then go to fin; if j-j then go to f i n 5 sy(j|,y=sy(j>•г+уij; if i>0 thendo i sy( j ) ,y--sy( i ) , y + y i j;

386. I>J THEN D0 I 1 = 1 ; I=j; J=n . ^ N D } IF SY( D t.MEXT=0 THEN DOJ

387. SY(I),N£XTsHEAP5 SY(HEAP) , Y = чY|J i

388. ELSE*<HEa?'fNSsJ; HEAP=HEAP+1I GO TO FJNJ END;

389. NNU I >«<3 THEN Do; JAK(NNU 1 , NNU 1 >=JAK{NNU1 , NNUП+PARC I J А К(N N U1 | N N U 2) , JAK(NNU2,NNU1)=JAK(NNU1 , NNU2UPARGiENDi

390. JAK(NNU2.NNUZ)= JAK(NNU2,NNU2)+ PARG i 001. FIN; END PRE5 00

391. PRAW = 05 41F 5TR=fJ pr then /* ОБНУЛЕНИЕ ЯКОБИАНА 001. METOD=l THEN >JAK = 0; 00

392. ELSE DOi JRAZ = SY( 1 ),NS5SY ( *) fNEXT=0JHEAP = JRaZ+1 ; 00

393. DO 1=1 TO IRAZ; SY(I).NS=Ii SY(J).Yc0; END? 00end; 001. KOMsl THfcN 00do; 00

394. МЕТКА (2 I ) ,МЕТКА (22 ) =RE } 00

395. МЕТКА(23)=ELL; МЕТКА(24)=EC ; 00

396. МЕТКА(25)=ERL; METК A(26)=ERC ; 00

397. METKA(27)=ERLEfMETKA(28)sERlESJN5 00

398. МЕТКА(30)=ELE5 МЕТКА(3 1 )=ERLC ' МЕТКА ( 32)=ERLCES ! 00

399. МЕТКА(33)=ELES! ПЕТКА(34)=ELCSS 001. МЕТКА(29)sREOSS? 001. END? ELSE DO; 00

400. МЕТКА ( 1 ) |"ETKA (2 ! = RE \ METKA(3)=LLl METKAMIsC} METKAf5)=RU 00

401. METKA(6)5RC; МЕТКА(7)=RLE; 00

402. МЕТКА < 10 ) =IE ; MET KA (1 1 ) =RUC ! МЕТКА ( 12 ) =RL,CES ; 00

403. МЕТКА ( 8 ) fRLES IN ; flETKA (9) rREOSS i 00

404. МЕТКА(13)=LES; МЕТКА( 14)=LCS ; 00end; 00

405. HALBH = H/2I OBRHE = l/H; OBRH= 1/HALBH ; 00do k = i to nrts; 00

406. NNUlPRTS(K) .NNUUJ ; NNU2 = RTS(K) ,NNU(2) 5 001.1sA6S(NNU i i i I2=ABS(NNU2) ; 001. HIL = U(nnU(l2) ; 00

407. NNU2<0 THEN NNU2=-NNU2; 001. ELSE 00

408. DO» U=FIL(NNU2); IF I I > 0 THEN NNU2 = FlKCI JA (FlU (NNU2 ) ) } 00else go to rs; 00end; 001. NNU1P0 THEN GO TO RS 5 00

409. NNUK0 THFN NNUl = r.NNUli 00else nnui=fikcjja(fii.(nnui) ); 00rsi pare = b; if kom = i then do; igd<2> =1;go to метка(RTS(Ki ,lab+20) $ 00

410. END! ELSE GO то МЕТКА(RTS(K) .LAB) ; 00

411. RF! PARGrRTSjK) ,T(5) J 001. PARt =RTS(К) ,PAR(2) i 00

412. PARJ = RTS ( К) , T ( 6) ;GO TO R? 00

413. I PARGsRTS { <) tPAR(6! 5 00

414. PARJ=RTS(K)fT<l>+PARC*HIL;C0T0Ri 00

415. C« PARG = RTS(K)tPAR( 1)/HAL8H{ 00

416. PAR J = -RTS (K ) , T ( I )-PA.RG*H IL 5 GO TO R; 00

417. RL!PARG=RTS(K),PAR(6)i 00 PARJ=(HlL40BRH*RTS(K),PAR{2)"RTS(K),PAR(ni*RTS(K!,TU))*PARG$ 00go то r; 00

418. RCJ HILF = RTS{K) t P AR ( 2 ) i 00

419. HI=HJLF*RTS(K) > P A R ( 1 ) i 001. HF=HI+HAL8Hi 001. PARG = H ILF/HFJ 00

420. PARJs ПН1-НАЦВН) *RTS (K) ,T ( 1) «HILF*H JLJ/HFIGO TO RJ 00

421. RLE|HILF=RTS|K|,PAR(2)i 001. PARE=RTS(K)tPAR(3)} 00

422. HF=HALbH*RTS(K) ,PAR(1 | J 001. H I SH I LF' + HF j 001. PARG = HALBH/H I i 00

423. PARJ^( (HILF-HF)*RTS(К) ,T ( 1 )*HIL*HALBH*H*PARE)/HI {CO TO R; 001. 5 PAREsRTS(K),PAR(2) 5 00

424. PARG = HALBH/RTS(К ) .PAR( 1 ) 5 00

425. PARJsRTS(Kj,T(1)+PARG*(2*PAPE*HIL)iGO TO RJ 00

426. HI l=h ilf + rlF + h J { 0<1. PARGsHALUM/hiii 0<1. PARFsRTSUI » p AR(5) i 0f

427. SI pARAsRTS(K).PAR(2) i 0<

428. PARF = RJS(«) . P AR(3 ) J 0£

429. PARGsRTS(ft ) ,par(6 ) » 0£

430. PARErPAKA*COS(ONEG*(TT+H)+PARF) i 0£ parj = rt5 (К) ,T ( П+pa!»G» (HIL + pare*para*cos (omec*tt*parf) ) igo TO r| 0£

431. S! PARG = RTS(ft| , PAR(6 ) ! RTS(К),N^U(4)= 05 0£

432. PARFpRTS(ft),PAR (3) I PARA = RTS(КJ . PAR ( 2 ) \ 02

433. HI=RTS(K).PAR(l) i PARU = "RTS(K). t ( 3 ) j 02

434. PAR£sPARA*COS(OMEG*fTT + H)+PARF)+PARU| 0 £

435. PA^J=(HILF*RTS(K) fT(t)+H*PARE)/HF5 GO TO Rj 02

436. ELEJ PARE=RTS(tf),par(2 j j 00

437. PARG=H/RTS(K) jPAR ( 1 ) » 00

438. PAR J=RTS (f< ) . T (\) +PARG*PARE i GO TO Ri 00

439. ERLC! HILF = RTS(K) .PAR (2) ; RTS ( К) . NNU(4)= \\ 00

440. PARt = r-RTS (К ) ,T (3) } 00

441. HF = H(LF + H*(RTS(k).PAR(1)+h/PTS{K),PAR(3)M 001. PARGsH/HF? 00erlces; hiuf=rts(kj,par<2); 02

442. UCHOL ! proc ( IRAZ |n'e .0 ,u |praw ,uc ,u" , ipr int , iper » inn) REORDER» 02

443. DO i = l то н» и 11)=u(i j/нi; end; 00do m l = l to н; i = u с(m l) i li=ur(i); h lsUR( i * i )Г i ; 00

444. Материалы о внедрении.4 отр.