Решение систем алгебраических и трансцендентных уравнений, страница 32

6.3.1.4. Корректировка структуры Суперрешателя для системы                                                  nнелинейных уравнений

Исходная структура Суперрешателя – для системы 2 уравнений.

1. Решатель системы уравнений (рис. 6.7, А).

Автоматически учитывает размерность N решаемой системы, увеличивая размерность гнездовой матрицы до N в строках выходной таблицы. Вывод координат корней по N–2 новым неизвестным производится вручную путем вставки имени выходной таблицы с индексом, соответствующим новой неизвестной, ниже решателя. Число вновь вставляемых таблиц равно числу  N–2 новых неизвестных.

2. Матрица сеансовых корней MSK (рис. 6.8, А).

Следует увеличить количество входных таблиц kzir до  N–2(zi – новая неизвестная, i = 3…N), вписав их имена в присваивание, формирующее таблицу MSK (подэтап 3).

3. Фильтр корней заданной области FKZO (рис. 6.8, Б).

В подэтапе 1 следует ввести 2(N–2) новых граничных значений заданной области по новым неизвестным. В подэтапе 2 в программе FKZO:

·  увеличить число присваиваний «MKi ¬ 0», i = 3…N;

·   увеличить выражение условия в if  на N–2 пар логических произведений (zjn < MSK_i,j) Ù (MSK_i,j < zjv),               j = 3…N;

·  увеличить столбец инструкций в if  на N–2 присваивания «MKj_chk¬ MSK_i,j», j = 3…N;

·  увеличить число аргументов MKj в функции augment(…) на N–2, j = 3…N.

4.  Фильтр неповторных корней FNPK (рис. 6.8, В).

Увеличить число присваиваний «KZO_KNK,m¬ FKZO_i,m» в столбце инструкций 3–го if на  N–2, m = 3…N.

5.  Фильтр новых корней FNK (рис. 6.8, Г).

Увеличить число присваиваний «MNK_chnk,j ¬ 0» в начале программы FNK  на N–2, j = 3…N. Увеличить число присваиваний «MNK_chnk,j ¬ FNPK_i,j » в столбце инструкций otherwise  на N–2, j = 3…N.