Если набор экспериментальных данных получен со значительной погрешностью, то не имеет смысла использовать интерполяцию полиномами Лагранжа и сплайнами для обработки результатов. В этом случае необходимо провести аппроксимирующую кривую, которая не проходит через экспериментальные точки, но в то же время отражает исследуемую зависимость, сглаживает возможные выбросы за счет погрешности эксперимента.
Обозначим узлы исходной таблицы данных через xj, где 0 ≤ j ≤ n – номер узла. Считаем известными значения экспериментальных данных в узловых точках ƒ(xj) = ƒj. Введем непрерывную функцию φ (х) для аппроксимации дискретной зависимости f(xj). В узлах функции φ(х)и f(x)будут отличаться на величину
ε j= φ(xj) – f(xj). Отклонения ε jмогут принимать положительные и отрицательные значения. Чтобы не учитывать знаки, возведем каждое отклонение в квадрат и просуммируем квадраты отклонений по всем узлам:
n n
Q = ∑ ε j ² = ∑[φ(xj) - ƒ(xj)]² (3.1)
j=0 j=0
Метод построения аппроксимирующей функции φ (х)из условия минимума величины Q получил название метода наименьших квадратов (МНК).
Наиболее распространен способ выбора функции φ(х)в виде линейной комбинации
φ(х) = c0 φ 0 (x)+ c1 φ 1 (x) + ... + cm φm(x), (3.2)
где φ 0 (х), φ 1 (x),...,φ m(х) — базисные функции; m ≤ n;
c0(x), c1 (x),…,cm (x)- коэффициенты, определяемые при минимизации величины Q.
Математически условия минимума квадратов отклонений Q запишем, приравнивая к нулю частные производные от Q по коэффициентам 0≤k≤ m:
n
∂ Q⁄∂ c0 = 2∑ [c0 φ 0 (xj) + c1 φ 1(xj) +…+
j=0
+ cm φ m(xj) - ƒj] φ 0 (xj) = 0,
n
(3.3) |
j=0
+ cm φ m(xj) - ƒj] φ 1 (xj) = 0,
…………………………………….
n
∂ Q⁄∂ cm = 2∑ [c0 φ 0 (xj) + c1 φ 1(xj) +…+
j=0
+ cm φ m(xj) - ƒj] φ m(xj) = 0.
Из системы линейных алгебраических уравнений (3.3) определяются все коэффициенты сk. Система (3.3) называется системой нормальных уравнений. Матрица этой системы имеет следующий вид:
(3.4)
и называется матрицей Грамма. Элементы матрицы Грамма являются скалярными произведениями базисных функций:
n
(φ j, φ k) = ∑ φ j(xj) φ k(xj). (3.5)
j=0
Расширенная матрица системы уравнений (3.5) получится добавлением справа к матрице Грамма столбца свободных членов
, (3.6)
где скалярные произведения, являющиеся элементами столбца, определяются аналогично (3.5):
n
(φ j, ƒ) = ∑ φ j(xj) ƒ j, (3.7)
j=0
Отметим основные особенности матрицы Грамма, полезные при программной реализации алгоритмов МНК:
1) Матрица симметрична, т.е. a i j = a ji, что позволяет сократить объем вычислений при заполнении матрицы.
2) Матрица является положительно определенной, следовательно, при решении системы нормальных уравнений методом исключения Гаусса можно отказаться от процедуры выбора главного элемента.
3) Определитель матрицы будет отличен от нуля, если в качестве базиса выбраны линейно независимые функции φk(x), при этом система (3.3) имеет единственное решение.
При обработке экспериментальных данных, определенных с погрешностью εв каждой узловой точке обычно начинают с аппроксимации функцией φk(x), представимой одной-двумя базисными функциями. После определения коэффициентов ckвычисляют величину Q по формуле(3.1). Если получится, что , то необходимо расширить базис добавлением новых функций φk(x).Расширение базиса необходимо осуществлять до тех пор, пока не выполнится условие .
Выбор конкретных базисных функций зависит от свойств аппроксимируемой функции ƒ(x), таких как периодичность, экспоненциальный или логарифмический характер, свойства симметрии, наличие ассимптотики и т.д.[8]