Файл исходных данных MAIN. 
DAT

nwK=60 !Кол. витков в Кат-ке

rKatvn=15d-3 !радиус Кат. внутренний, м

tWKr=0.5d-3 !толщина Витка Катушки в направлении радиуса, м tWKz=10d-3 !толщина (высота)Витка Катушки в направлении оси Z, м tlzwwr=ld-4 !толщина Изоляции Виток-Виток в направлении радиуса, м shrzbK=ld-3 !шаг разбиен.(желаемый) для витков Катушки, м udSoprK=1.7d-8 !удельное Сопротивление (материала) Катушки, Омм s0=0.2d-3, [начальный зазор между катушкой и диском (по металлам).

hDisk=6d-3, [толщина Диска (в направлении оси Z), м rDvn=ld-2 [радиус Диска внутренний.

tDr=35d-3 [радиальная толщина диска (tDr=rDvn-rDvn).

! rDnar=45d-3 [радиус Диска наружный.

shrzbD=ld-3 !шаг разбиен. Исходный (желаемый) для Диска, м.

udSoprD=2.4d-8 [удельное Сопротивление (материала) Диска, Ом-м YdPID=2.8d3 [Удельная плотность материала диска, кГ/мЗ.

Amas=0.5d0 [Масса, присоединённая к диску (масса мех., приводимого в дв.).

roProv=1.7d-8 [удельное сопротивление (материала) Провода, Ом-м DiamPr=3d-3 [Диаметр Провода, м.

TizPr=ld-3 [Толщина изоляции Провода, м.

DIProv=0.5d0 [Длина Провода (расстояние от конд-ра до кат-ки), м ht=2d-6 [шаг по времени расчёта хар-к ИДМ, сек.

tkon=5d-3 [время окончания циклов при расчёте хар-к ИДМ, сек.

Uc0=600d0 [начальное напряжение на конденсаторе, В.

emk=500d-6 [ёмкость конденсатора, Ф.

LTiristor=1.21d-7 [Индуктивность Тиристора, Гн RTiristor=2.33d-3 [Активное сопротивление Тиристора, Ом xkon=5d-3 [ход конечный (плоный) диска и присоединённой массы.

nlntervX=3 !кол. разбиений Хода конечного на участки для постр.сплайнов.

Fm0=400d0 [Начальная сила механической хар-ки.

zl=3d4 [жесткость на первом участке мех. хар-ки.

xl=ld-3 [длина первого участка (по аргументу) мех. хар-ки.

Fml=600d0 [Сила в начале второго участка мех. хар-ки z2=6d4 [жесткость на втором участке мех. хар-ки.

kratnPrt=15 [кратность вывода на Печать таблицы результатов счёта.

poligrf=.true. [указывает, выводить графики каждой величины на «своём» поле [(ecnn.true.), или на одном поле всех величин (eoin.false.) meth=2 !PARAM (12)=meth-MeTOfl интегрирования 1-Адамса-Мултона 2-Гира miter=l !PARAM (13)=miterМожет иметь значения 0,1, 2, 3, где.

!0-Метод функциональных итераций. При этом д.б.РАЯАМ (19)=0 !1-Метод хорд с Якобианом, предоставляемым пользователем !2-Метод хорд с Якобианом разделённых разностей !3-Метод хорд с Якобианом, замещенным диагональной матрицей.

[При этом PARAM (19) д.б.=0, т. е.матрица АА не используется iatype=l !PARAM (19)=iatypeтип матрицы АА; Может быть 0,1, 2, где !0-матрица АА не используется !1-матрица AA=const (nMeeT место на этапе трогания).

[2-матрица АА с элементами, ЗАВИСЯЩИМИ от времени.

{ВНИМАНИЕ: В процессе счёта ПРИ ДВИЖЕНИИ диска программа {устанавливает его в 2.

! 1-nwK 2-twKr 3-twKz 4-rKatvn 5-tDr 6-rDvn 7-hDisk.

ZNIZ= 50, 3d-4, 5d-3, ld-2, 20d-3, 5d-3, 6d-3 !Мин. зн-я переменных.

ZVRX=95, ld-3, 10d-3, 2d-2, 8d-2, 5d-2, ld-2 !Макс. зн-я переменных.

OPTPAR=t, t, f, t, t, t, f.

kDopViz=400 {Задаваемое число расчётов целевой функции при отпимизации {После достижения счётчиком этого значения {программа предлагает завершить или продолжить оптимизацию.

Поиск оптимального ИДМ выполняется из точки с начальными данными, указанными в файле данных. Заметим, что в файле отсутствует запись iStart=2. указывающая на автоматический поиск оптимального варианта из второй стартовой точки. Это объясняется тем, что:

• 1. эффективность поиска оптимального варианта методом деформируемого многогранника выше, чем градиентным методом наискорейшего спуска;
• 2. случайный характер формирования начального комплекса определяется значениями границ (верхней и нижней) разрешённого диапазона изменения переменных;
• 3. влияние начальной точки исходных данных при поиске оптимального варианта на результат (на найденную оптимальную точку) меньше, чем при поиске оптимума градиентным методом.