Вторая пачка часть 184


Процедура TRBAK. Входные: N (тип integer) - по¬ря¬док дей¬ствительной симметрической матрицы А; М1 и М2 (тип in¬¬teger) - граничная пара обозначения массива но¬р¬ми¬ро¬ван¬¬ных собственных векторов трех¬ди¬аго¬наль¬ной матрицы Аn-1, определенной в процедуре TRED1; AZ (тип real) - мас¬¬сив размером N?N, в котором су¬щес¬твен¬ны лишь под¬д謬агональные элементы (со¬от¬вет¬ст¬ву-ет выходному мас¬си¬ву А процедуры TRED1); Е (тип re¬al) - массив размером N?1, N - 1 элементов которого от Е[2] до Е[N] определяют вн嬬¬диагональные элементы трех¬диагональной матрицы Аn-1. Первый элемент мас¬си¬ва Е[1] равен нулю. Выходные: ZA (тип real) - массив размером N ? (М2 - М1 + 1), со¬дер¬жа¬щий но¬рми¬ро¬ван¬ные соб¬ственные векторы матрицы А, прну¬ме¬ро¬ван¬ные от М1 до М2. Остальные процедуры не имеют прямого отношения к рассматриваемой задаче, выполняют вспо¬мо¬га¬тель¬ную роль, указанную в комментариях к каж¬дой про¬це¬ду¬ре. В заключение следует отметить, что МГК - сложный с вы¬числительной точки зрения метод и требует мак¬си¬му¬ма вни¬мательности. Программа, ре¬а¬ли¬зующая МГК, довольно сложная и трудоемкая. Один расчет матрицы 25?60 за¬ни¬мает до 20 минут на ЭВМ типа IВМ-286 с частотой 12 МГц. Ввиду слож¬нос¬ти текст про¬граммы приводится пол¬нос¬¬тью. Для на¬строй¬ки про¬грам¬мы перед запуском следует укବзать пол¬ный путь, где находятся графический драйвер, файл с исходными дан¬ными, и определить значения сле¬ду¬ю¬¬щих констант: PR - если PR = 0, то результаты будут записываться в файл на магитный диск С: с именем BURKOV; если PR = 1, то результаты будут выводится на печатающее устройство; если PR = 2, то на экран дисплея; START - строковая константа, указывающая полное имя файла, в котором хранятся исходные данные; N - количество параметров, которые измеряют в ходе эксперимента, определяет количество столбцов в матрице данных; М - количество измерений параметров, определяет количество строк в матрице данных. В матрице исходных данных первые три строки со¬дер¬жат служебную информацию: название массива, где и когда получены данные; фамилию и имя исследователя; условные названия измеряемых параметров (лучше, если название будет состоять из одного или двух сим¬во¬лов
Индекс
Элементарные функции    Линейные уравнения    Нелинейные уравнения    Случайные числа


Hosted by uCoz