Примечание!Что ж, график у нас теперь есть. Как же определить, какая именно точка соответствует решению нашего уравнения? С этим поможет трассировка. Щелкните по получившемуся графику правой кнопкой мыши и выберите Trace. У вас на графике появится этакий " прицел", который можно будет передвигать мышкой, и окно, в котором будут отображаться координаты точки, находящейся в середине прицела (см. скриншоты). Прицел этот можно, кстати, двигать и при помощи клавиатурных стрелочек, но все равно в любом случае передвигаться он будет исключительно и только вдоль кривой, отображающей ход нашей функции, нули которой мы и ищем. Думаю, основную идею вы уже уловили: передвигая крестик-" прицел" вдоль кривой, можно в окне трассировки увидеть точку, Y-координата которой будет равна нулю. Здесь, правда, может возникнуть небольшая проблема. В нашем уравнении, например, корень получился не целым, а трассировать функцию по X мы можем только с шагом, равным единице. Конечно, для того, чтобы локализовать корень для последующего применения функции root, в нашем случае достаточно и такой точности: достаточно найти интервал, на котором функция будет менять свой знак (то есть ее значение будет меняться с положительного на отрицательное или, напротив, с отрицательного на положительное). Но для других уравнений может случиться такое, что на единичном интервале может лежать и несколько решений — действительно, почему бы и нет, ведь даже для квадратных уравнений несложно придумать случай, когда такое может случиться. Какой из этого может быть выход? Самый простой и очевидный — уменьшить шаг трассировки по оси X, чтобы можно было более точно искать нулевые точки по оси Y. Сделать это в MathCAD очень просто. Помните, как мы задавали диапазон для нашей переменной x? Давайте кое-что поменяем и запишем его следующим образом: x: = -10, -9.9..10. Теперь уже шаг трассировки у нас будет не единица, а одна десятая (0, 1). Для того, чтобы уменьшить его еще в десять раз, нужно записать x: = -10, -9.99..10. Общий принцип здесь вполне очевиден: мы записываем следующую точку диапазона следом за первой и, таким образом, указываем шаг, с которым MathCAD будет строить (а в итоге — и трассировать) наш график. Для более наглядного отображения данных имеет смысл добавить на график координатную сетку. Сделать это можно, если кликнуть по графику дважды, а затем в появившемся окне (см. скриншот) установить флажки напротив пунктов Grid Lines в категориях X Axis и Y Axis. Первый флажок устанавливает координатные линии для оси X, а второй, соответственно, для оси Y. Единственное, что также рекомендую — так это сразу поменять и цвет линий координатной сетки, кликнув по квадратикам рядом с этими пунктами — по умолчанию он задан кислотно-зеленым, что не очень удобно для глаз. Хотя это уже, конечно, дело вкуса. 5 Задать как приближение значения точек пересечения х1, х2, х3. В примере х1=-0.9, х2=0.2, х3= 0.7. 6 Вычислить значение корней с помощью формул: root (f(x1), x1), root (f(x2), x2), root (f(x3), x3). Полученные значения корней такие: х1=-0.92, х2=0.21, х3= 0.721 (рис. 18).
Рисунок 18 – Результат нахождения корней с использованием функции root II Для уравнения 1. Создать вектор из коэффициентов уравнения, используя панель управления Matrix (Матрица) (рис.19) и задав один столбец и четыре строки для коэффициентов уравнения.
Рисунок 19 – Диалоговое окно для определения вектора из коэффициентов уравнения
Вектор из коэффициентов уравнения будет иметь следующий вид
2. С помощью встроенной функции r: = polyroots (v) найти корни уравнения и представить их в виде вектора rT, транспонированного по отношению к r, то есть преобразованного из столбца в строку.
3. Создать циклы для переменной х и количества найденных корней:
4. Построить графики для функции и определить функцию в точках корней. В точках корней значения функции равны нулю. 5. Определить значения корней на графике (рис. 20).
Рисунок 20 – Результат нахождения корней с использованием функции polyroots
III Для уравнения 1. Записать левую часть уравнения
2. Поставить логический знак «=» и в правой части записать 0.
|
найти корни на интервале [-1.1, 7.1], шаг изменения переменной х равен 0.1.
.
