Файл исходных данных 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 !радиус Диска внутренний.

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

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

udSoprD=2.4d-8 Сдельное Сопротивление (материала) Диска, Ом-м.

YdPID=2.8d3 !Удельная плотность материала диска, кГ/мЗ.

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

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 [Начальная сила механической хар-ки (см. рис.1.7).

zl=3d4 [жесткость на первом участке мех. хар-ки (см. рис.1.7) xl=ld-3 !длина первого участка (по аргументу) мех. хар-ки (см. рис.1.7) Fml=600d0 !Сила в начале второго участка мех. хар-ки (см. рис.1.7) z2=6d4 !жесткость на втором участке мех. хар-ки (см. рис.1.7).

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

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

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

!1-матрица AA=const (имеет место на этапе трогания).

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

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

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

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

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

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

kDopViz=400 (Задаваемое число расчётов целевой функции при отпимизации (После достижения счётчиком этого значения (программа предлагает завершить или продолжить оптимизацию iStart=2 (Кол-во стартовых точек при оптимизации. Может быть 1 или 2 /.

В конце файла данных имеется запись iStart=2. iStart — количество стартовых точек при оптимизации. От стартовой точки может зависеть результат оптимизации. Если пользователь задал iStart=l, то в качестве стартовой точки программа использует начальные переменные, указанные в файле данных. Эти переменные помещаются в массив OPTXAR. Затем, в процессе счёта, начальные переменные в массиве OPTXAR замещаются оптимальными переменными. Если iStart=2, го первый поиск выполняется с начальными данными, а второй — из точки, которая определяется как полусумма минимальных значений допустимого диапазона и оптимальных, полученных при первоначальной оптимизации:

где i — номер оптимизируемой переменной.

После этого выполняется новый поиск оптимальных переменных. Понятно, что старт даже из двух разных точек не гарантирует нахождение глобального минимума функции цели.