28.1 Вычисление значений многочленов ¶
Значение многочлена, представленного вектором c, может быть вычислено в точке x очень просто, как показывает следующий пример:
N = length (c) - 1; val = dot (x.^(N:-1:0), c);
Хотя приведенный выше пример демонстрирует, насколько легко вычислить значение многочлена, это не самый устойчивый алгоритм. Для многочленов большей степени следует использовать более элегантные алгоритмы, такие как метод Горнера, что и делает функция Octave polyval.
В случае, когда x является квадратной матрицей, многочлен, заданный c, всё ещё определён. Как и когда x является скаляром, очевидная реализация легко выражается в Octave, но и в этом случае более элегантные алгоритмы работают лучше. Функция polyvalm предоставляет такой алгоритм.
-
:
y =polyval(p, x)¶ -
:
y =polyval(p, x, [], mu)¶ -
:
[y, dy] =polyval(p, x, s)¶ -
:
[y, dy] =polyval(p, x, s, mu)¶ -
Вычисляет значение многочлена p в заданных значениях x.
Если x является вектором или матрицей, многочлен вычисляется для каждого элемента x.
Если присутствует mu, вычисляет многочлен для (x - mu(1)) / mu(2).
Помимо вычисления значения многочлена, второй результат представляет интервал прогнозирования y +/- dy, который содержит как минимум 50% будущих прогнозов. Для расчёта интервала прогнозирования необходима структурированная переменная s, полученная из
polyfit.См. также: polyvalm, polyaffine, polyfit, roots, poly.
© 1996–2023 The Octave Project Developers
Permission is granted to make and distribute verbatim copies of this manual provided the copyright notice and this permission notice are preserved on all copies.
Permission is granted to copy and distribute modified versions of this manual under the conditions for verbatim copying, provided that the entire resulting derived work is distributed under the terms of a permission notice identical to this one.Permission is granted to copy and distribute translations of this manual into another language, under the above conditions for modified versions.
https://docs.octave.org/v9.2.0/Evaluating-Polynomials.html