mnkpoly

Синтаксис:

P = mnkpoly(X, Y, S, n);

Аргументы:

X – входной массив координат х для таблично-заданной функции Y(X),

Y – входной массив координат y для таблично-заданной функции Y(X),

S – входной массив среднеквадратичных отклонений для каждой координаты таблично-заданной функции Y(X),

n – степень апроксимирующего полинома, целое число.

Описание:

mnkpoly(A, B) – функция возвращает массив коэффициентов полинома степени n, полученного в результате аппроксимации таблично-заданной функции Y(X), с заданным среднеквадратичным отклонением S для каждой точки. Вычисление производится при помощи сингулярного матричного преобразования (SVD). Если значения среднеквадратического отклонения не известны, их можно задавать равными единице.

Входные массивы X, Y, S могут задаваться:

P = mnkpoly (X, Y, S, n);

P = mnkpoly ([x₁,x₂,x₃,x₄],[x₁,x₂,x₃,x₄], [s₁,s₂,s₃,s₄], n);

P = mnkpoly ([0, 1, 6, 9],[-1, 5, 0, 4], [0.1, 0.5, 0.7, 0.4], 2);

Результат:

P – выходной массив, содержащий коэффициенты полинома, полученного в результате аппроксимации таблично-заданной функции:

p(t) = p₀+p₁t+p₂t²+…+p_ntⁿ

Выходной массив имеет размерность n+1.

Пример:

initialization
X = [0, 1 , 6, 9];
Y = [-1, 5, 0, 4];
S = [0.1, 0.5 0.7 0.4];
P = mnkpoly(X,Y,S,2);
end;
Z = polyval(P, X);