Решение уравнения с одним неизвестным, страница 11

Перед выполнением этапов вычислений целесообразно перевести Mathcad  в режим ручных вычислений (Tools-Calculate-убрать-галочку-у-команды-Automatic Calculation): нажатием клавиши <F9> следует запускать вычисления как отдельных программ при отладке, так и всего комплекса при окончательном расчете.

Этапы вычисления

1.  Ввод в РДМ уравнения (5.4) 12-го порядка, векторов параметров a, b, c, вектора порядка уравнения vp (в виде входной таблицы для сокращения габаритов объекта) и вычисления полного количества корней (60) для последующего контроля правильности составления программ.

2.  Программа вычисления корней уравнения (вектора vr12) посредством решателя polyroots(v) для векторизационного сочетания значений параметров: для расчета используются значения компонент параметров с одинаковыми номерами. Программа почти совпадает с программой рис. 5.3 для одного параметра: в полиномном векторе v1 компоненты всех трех параметров имеют один и тот же индекс(номер) j (сравните с программой рис. 5.3 для всех сочетаний компонент двух параметров, где индексы компонент разные). Контроль правильности составления программы производится вычислением длины вектора vr12.

3.  Программы формирования расчетных векторов rva, rvb, rvc и столбца с номерами строк Таблицы ИТИ. Для каждого расчетного вектора следует вычислять его длину, контролируя тем самым правильность введенной программы.

4.  Формирование заголовка таблицы.

5.  Вызов таблицы ИТИ.

6.  Формирование названия таблицы в виде текстового блока, располагаемого над таблицей.

 


Рис. 5.5. Этапы квазивекторизованного вычисления зависимости                                                                             корней уравнения (5.4) от 3 параметров

5.1.4. Определение погрешности вычисления корней                                   решателем polyroots(v)       

Абсолютная погрешность dr вычисления корня r уравнения (5.0) определяется из линейного отрезка ряда Тейлора функции уравнения f(x) в окрестности точного значения корня rт :