Cтраница 3
В подпрограмме, реализующей решение указанной системы методом квадратных корней с удвоенной точностью ЕС ЭВМ, используются только модули факторизации и вычисления решения. Итерационное уточнение здесь не выполняется по следующим причинам. Во-вторых, поскольку для итерационного уточнения невязку необходимо вычислять на повышенной точности, то в данном случае пришлось бы использовать расширенное машинное слово ЕС ЭВМ ( с мантиссой в 112 бит или с 28 шестнадцатиричными разрядами), что можно реализовать лишь специальными подпрограммами. А это и слишком громоздко, и в большинстве случаев неоправданно, как указывалось выше. [31]
Из табл. 8.8 видны три интересных момента: ( 1) Величина невязки не меняется заметно, даже если достигнуто значительное улучшение точности вычисленного решения. Итерационное уточнение стоит довольно дорого; затраты на уточнение превосходят на 50 % затраты на обычное решение системы с 14х14 - матрицей, обладающей хорошими свойствами, и на 25 % - для системы с 40х40 - матрицей, обладающей неплохими свойствами. Оценка числа обусловленности почти бесполезна. Нет готового объяснения тому факту, что ошибки в решениях, полученных подпрограммой LINEQ1 без итерационного уточнения, соответственно лучше ошибок в решениях, полученных пакетом LINPACK; этот факт имел место для всех 10 просчитанных примеров, за исключением трех случаев, когда были достигнуты одинаковые результаты. [32]
Если вычисленное решение не выдерживает описанную выше проверку, то вырабатывается предупреждающее сообщение и подпрограмма пробует применить итерационное уточнение. Итерационное уточнение дорого как с точки зрения времени выполнения подпрограммы, так и с точки зрения размеров требуемой памяти. Поэтому LEQT2F применяет подпрограмму итерационного уточнения, только если она необходима. Итерационное уточнение последовательно улучшает ответ до тех пор, пока он не станет правильным с точностью выполнения машинных операций. Если матрица столь плохо обусловленная, что процесс уточнения не сходится, то вырабатывается окончательная ошибка. Если значение IDGT равно 0, то LEQT2F опускает проверку точности и выполнение итерационного уточнения. [33]
![]() |
Параметры системы САЦ-СКЦ-КГМО для t 7 ч. [34] |
Таким образом, использование полученной априорным усреднением по пространственной координате величины отношения расхода к давлению в линеаризованных уравнениях нестационарной газопередачи может приводить к большим погрешностям вычислений. Итерационное уточнение параметров линеаризованных систем, основанное на усреднении по двум координатам, позволяет иметь решения, близкие к решениям нелинейных систем, полученных методом сеток. [35]
Отметим, что малые компоненты вектора х, как правило, имеют низкую относительную точность. Первый шаг итерационного уточнения позволяет получить оценку чувствительности решения системы уравнений. [36]
Полученные в результате детального расчета значения потерь давления во всех теплообменных аппаратах и трубопроводах установки сравниваются с ранее заданными при расчете тепловой схемы величинами гидравлических потерь. При необходимости производится итерационное уточнение результатов расчетов тепловой схемы установки и теплообменных аппаратов. [37]
Это наводит на мысль использовать процедуру асе inverse для решения систем уравнений. Ведь если процедура итерационного уточнения сходится и I - АХ Цоэ для матрицы X, принятой за обратную, существенно меньше единицы, то это является строгим доказательством того, что X дает правильно округленный результат обращения. [38]
![]() |
График пяти низших частот.| График коэффициентов демпфирования. [39] |
Решение поставленной задачи па собственные значения осуществляется методом Мюллера без выделения комплексного параметра ( о в явном виде. На каждом шаге итерационного уточнения величины т методом ортогональной прогонки решается краевая задача (3.181), (3.182), собственные числа определяются из условия равенства нулю определителя системы, вычисляемого методом Гаусса. В расчетах варьируется мгновенный модуль упругости Е2, параметры ядра релаксации А, а, р, диаметр d4 стекло-пластикового кожуха при неизменной его толщине и фиксированных значениях других параметров. [40]
Решение системы АХ В найдено как результат последовательного решения систем LY В и U X Y. Матрица В сохранена для выполнения итерационного уточнения. Матрица X размера пХ г размещена на месте Y. При обращении к процедуре числовой массив х необходимо задать. Например, его можно отождествить с массивом Ь, поскольку они имеют одинаковую размерность. [41]
Решение системы АХ В найдено как результат последовательного решения систем LY В и DUX Y. Матрица В сохранена для выполнения итерационного уточнения решения X. Матрица X размера nX r размещена на месте Y. Например, его можно отождествить с массивом Ь, поскольку они имеют одинаковую размерность. [42]
Для решения системы уравнений RX V использован метод обратной подстановки. Для первого приближения X необходимо провести итерационное уточнение. Матрица невязок В - АХ должна быть вычислена с двойной точностью. Матрица коррекций Z вычисляется по столбцам. Первый шаг итераций будет выполнен всегда. [43]
В отличие от них в процедуре ortho 2 итерационное уточнение не реализовано, зато в ней предусмотрена существенная экономия памяти за счет того, что массив исходной матрицы А использован для размещения промежуточных результатов, а затем и элементов обратной матрицы. [44]
Метод: Матрица А факторизуется на нижнюю и верхнюю треугольные матрицы L и U, и затем последовательно решаются системы LZB и UXZ. Применяется накопление скалярных произведений с удвоенной точностью и итерационное уточнение, поэтому решение очень точное всякий раз, когда S на выходе равно нулю. [45]