О построении интерполяционно-аппроксимационного многочлена Текст научной статьи по специальности «Математика»

НАНОСИСТЕМЫ: ФИЗИКА, ХИМИЯ, МАТЕМАТИКА, 2012, 3 (6), С. 5-15 УДК 519.65

О ПОСТРОЕНИИ ИНТЕРПОЛЯЦИОННО-АППРОКСИМАЦИОННОГО МНОГОЧЛЕНА

В. Д. Лукьянов

ОАО «Авангард», Санкт-Петербур, Россия [email protected]

Построен интерполяционно-аппроксимационный многочлен, который решает задачу аппроксимации данных методом наименьших квадратов и значения которого совпадают с заданными значениями в некоторых точках. С помощью интерполяционно-аппроксимационного многочлена решена задача аппроксимации данных пучком кривых, с подбором наилучшим образом общих точек. Приведены примеры нахождения интерполяционно-аппроксимационного многочлена.

Ключевые слова: многочлен, интерполяция, аппроксимация, метод наименьших квадратов, интерполяционно-аппроксимационный многочлен, пучок кривых.

1. Введение

При обработке экспериментальных данных в физических исследованиях возникает задача аналитического восстановления функции по измеренным значениям. Различают два метода решения этой задачи: метод интерполирования и метод аппроксимация данных [1-3].

При интерполировании функции подбирается такой интерполяционный многочлен значения, которого совпадают с заданными значениями функции при некоторых значениях аргумента, которые здесь называют узлами интерполирования. Степень интерполяционного многочлена всегда по крайней мере на единицу меньше количества узлов интерполирования. На практике стараются ограничиться интерполяционными многочленами небольшой степени.

Аппроксимацию функции многочленами используют, если значения функции в узлах аппроксимации известны приближённо или количество известных её значений велико. Здесь, в частности, при одном и том же значении аргумента может быть задано несколько значений функций, что недопустимо в задаче интерполирования.

При аппроксимации функции находят такой многочлен, который имеет значения, с некоторой погрешностью приближающий заданные значения функции. Одним из методов решения задачи аппроксимации является метод наименьших квадратов (МНК) [1-3].

Пусть некоторые значения восстанавливаемой функции известны точно, например, из физических соображений. К восстановлению аналитического представления функции при обработке экспериментальных данных, часть которых известна точно, приходим при рассмотрении, например, краевой задачи. Здесь значения функции при некоторых значениях аргумента заданы до начала решения задачи, или эти значения измеряются со значительно большей точностью, чем при прочих значениях аргумента, где погрешность значительна.

Использование традиционного МНК не решает задачу аппроксимации функции, так, чтобы её значения совпадали с заданными значениями. Эта задача в статье будет решена объединением метода интерполяции и МНК.

2. Смешанная задача аппроксимации и интерполирования функции

2.1. Постановка смешанной задачи аппроксимации и интерполирования функции

Пусть заданы значения функции У0, У2, ..., Ум в N + 1 узлах интерполирования а0, а\, ..., ам и значения функции уь .. . , ум, в М узлах аппроксимации Ъ\, ... , Ъм. Найдём многочлен Рз(х) степени Б ^ N + М, который имеет во всех узлах интерполирования заданные значения. Здесь должны быть выполнены равенства Рз(ап) = Уп при п = 0, 1, 2, ..., N. Все остальные заданные значения в узлах аппроксимации многочлена должен аппроксимировать МНК.

Многочлен Р = Рз(х), который часть данных интерполирует, а часть - аппроксимирует, назовём интерполяционно-апроксимационным многочленом.

2.2. Построение интерполяционно-апроксимационного многочлена

2.2.1. Построим интерполяционный многочлен Тм (х) степени не выше N по значениям функции У0, У2, ..., Ум в узлах интерполирования а0, аь ..., ам. Эту задачу решает интерполяционный многочлен Лагранжа, который имеет вид [1-3]

Тм (х) = ± Уп ^, (1)

п=0

где Пп(х) представляет собой произведение многочленов первой степени вида шк(х) = х — ак при к = 0, 1, ..., N, кроме к = п,

N

Пп(х) = Ши (х) = Шо(х)шг(х) . ..Шп-\(х)Шп+\(х) ...шм (х) .

к=0, к=п

Многочлен Тм (х) является искомым интерполяционным многочленом степени не выше чем N. Выполнены равенства Тм (аг) = Уг при всех значениях г от 0 до N. Равенство проверяется непосредственным вычислением с учётом того, что Пк (аг) = 0 при всех значениях г от 0 до N и к = г, так как Пк(х) при к = г обязательно содержит сомножитель

шг(х), а шг(аг) = 0.

2.2.2. Ищем искомый интерполяционно-апроксимационным многочлен Рз (х) в виде суммы построенного интерполяционного многочлена Тм (х) и некоторого многочлена Qs(х) степени Б ^ N + М

Рз (х) = Тм (х) + Qs (х). (2)

При выполнении равенств Рз(ап) = Уп имеем Qs(ап) = 0 при всех значениях п от 0 до N, так как

QS (ап) = Р3 (ап) — ТМ (ап) = Уп — Уп = 0

Таким образом, все узлы интерполирования а0, аь ..., ам являются корнями многочлена Qs(х), поэтому для этого многочлена имеем следующее представление

к

Qs(х) = П(х)^ ^кхк-\ (3)

к=1

N

где первый сомножитель — это многочлен П(х) = Л шк (х), имеющий степень N +1 и

к=0

корни во всех узлах интерполирования, а второй сомножитель — это многочлен степени не выше чем К — 1, где К = Б — N, что обеспечивает многочлену Qs(х) степень Б ^ N + М.

Коэффициенты второго сомножителя Л1} Л2, ..., Лк в формуле (3) являются искомыми величинами и вычисляются согласно МНК. Согласно МНК ищем коэффициенты многочлена из требования обеспечения этими величинами минимума функции

м

Г (Ли Л2,...,Лк) = ^ & Ль Л2,..., Лк),

т=1

где 5т(Л1, Л2, ..., Лк) —это отклонение в узле аппроксимации Ьп значения многочлена Ря (Ьт) от известного значения функции ут,

к

ёт(Лх, Л2, ..., Лк) = Ут - Ря (Ьт) = Ут - ТМ (х) - П(х) ^ Ли Хк—.

к=1

Необходимое условие экстремума функции Г — это равенство нулю всех частных производных по переменным Л1, Л2, ..., Лк.

Имеем при всех натуральных значениях г от единицы до К

дГ м к

— = -2^(Ут - Тм (Ьт) - П(Ьт)^ Л^ ) П(Ьт) Ь^1.

т=1 к=1

Приравнивая все частные производные к нулю, после преобразований приходим к системе линейных алгебраических уравнений (СЛАУ) для вычисления коэффициентов

Л1,Л2,...,Лк.

апЛ1 + а.12 Л2 + ... + а1к Лк = в1, &21Л1 + а22 Л2 + ... + а2к Лк = в2,

(4)

ак 1Л1 + ак 2 Л2 + ... + акк Лк = в к,

где

м м N у м

Оц = £ П2(Ьт)^-2, вг = £ УтЩ(Ьт)Ьт - £ оТИГ^ Ьт) Щт) Ь]п.

т=1 т=1 п=0 п( п) т=1

Перепишем СЛАУ (4) в матричном виде АХ = В, где матрицы А = [аг^],

Т т

X = [Лк, Лк-1,..., Л1] , В = [в1, в2,..., вк] , знак Т здесь означает транспонирование

матрицы, т.е. матрицы X и В — это матрицы-столбцы.

Имеем X = А-1В, если матрица А — неособенная , у которой главный определитель А = А = 0. Все коэффициенты в представлении (3) найдены, таким образом задача нахождения интерполяционно-аппроксимационного многочлена, заданного в формуле (2), полностью решена.

Отметим, что второй сомножитель в представлении (3) может быть тождественно равным нулю, если дополнительно выполнены равенства ^ (Ьт) = ут при всех т от 1 до М. В этом случае найден многочлен TN (х), значения которого совпадают в точности с заданными значениями как в узлах интерполирования, так и в узлах аппроксимации функции. Так как все отклонения 5т в узлах аппроксимации равны нулю, то это обеспечивает требуемый минимум суммы квадратов отклонений. Следовательно, искомая функция, значения которой заданы, совпадает с многочленом ^ (х).

2.1.3.Докажем невырожденность матрицы А. Непосредственное вычисление главного определителя СЛАУ А = А даёт выражение

м м м к

А ^ Е Е ... Е ™ (Ьто , Ьт-1,..., Ьтк) Д П"(Ьтг) Ь— , (5)

т-1 = 1 т2 = 1 тк=1 г=1

где W(Ьто, Ьт1, ..., Ьтк) — определитель Вандермонда,

м

W (Ьто ,ЬШ1 ,...,Ьгпм )= П (Ьт1 - Ьт3 ).

Для доказательства того, что определитель А отличен от нуля, заметим, что его величина не измениться от перенумерации узлов аппроксимации однако это изменит порядок аргументов в определителе Вандермонда и сомножителей в произведении в представлении (5) для этого определителя. Количество всех таких перенумераций равно количеству всех перестановок из К элементов и равно К!. В каждом определителе добъёмся одинакового порядка точек и преобразуем этот определитель к первоначальному определителю Вандермонда W с учётом изменения знака и сомножителей в произведении. Если теперь сложить все эти К! равных по величине определителей, то с одной стороны получим величину

К!А. А с другой стороны сумму К! слагаемых. После объединения всех кратных сумм

к

в одну кратную, вынесения общего множителя W(Ьто, Ьт1, ..., Ьтк) П +\(Ьт1), полу-

г=1

к

чится новый множитель, состоящий из К! слагаемых вида Л Ь]-1 с различным порядком

г=1 1

сомножителей и соответствующим знаком. Этот множитель по определению определителя К-го порядка равен определителю Вандермонда. Из полученного равенства найдём величину определителя А

1 м м м к

А = К ^ £... £ W 2( Ьт1 ,...,Ьтм )\{^2(Ьтг)

т.1 = 1 т2 = 1 тк = 1 ¿=1

Определитель А представляет сумму квадратов чисел среди которых есть ненулевые, поэтому его величина отлична от нуля.

2.3. Частные случаи построения интерполяционно-аппроксимационных многочленов

2.3.1. Имеем один узел интерполяции а0, в котором известно точное значение функции У0,иМ узлов аппроксимации Ь1,..., Ьм с известными значениями функции уь..., ум. Найдем линейный интерполяционно-аппроксимационный многочлен Р1 (х). Интерполяционный многочлен имеет вид Т0(х) = У0. Искомый многочлен ищем в виде

Р1(х) = Уо + А1 (х - ао). (6)

Формула (6) для многочлена Р1(х) содержит только одну неизвестную величину А1, которую определим из условия минимума функции

м

Р (^1) = ^(Ут - ¥ - А (Ьт - ао))2.

т=1

Необходимое условие минимума функции Г(А1)приводит к уравнению для определения коэффициента А1

м

Г'(А1) = -2^2 (Ут - ¥0 - А1(Ьт - ас)) (Ьт - а0) = 0,

т=1

м

или Л ^ (Ьт - ао)2

т=1

м

Е (Ут - Уо)(Ьт - ао), откуда найдем величину

т=1

А\

м

Е (Ут - У>) (Ьт - ао)

т=1

м '

Е (Ьт - ао)2

т=1

(7)

Пример 1. На рис. 1 показаны экспериментальные результаты, взятые из оригинальной статьи Хаббла [4] и использованные для вычисления постоянной Хаббла. Точками показаны экспериментальные данные: по оси абсцисс отложено расстояние от Земли до звезды, а по оси ординат — радиальная скорость звезды. Штриховая линия —это результат использования стандартного МНК при аппроксимации функцией V = аО + в, где V — величина скорости звезды относительно Земли, О — расстояние от Земли до звезды. Имеем при этом величину в = 0. Полученный результат противоречит физическому смыслу: при нулевом расстоянии получается ненулевая относительная скорость Земли. Сплошная линия показывает результат использования формулы (7), когда при обработки тех же результатов используется функция V = аО, удовлетворяющая условию V(0) = 0.

1200

800

400

У(км > •

I /л у 1

• /

4 У

■•у ф

/ • 0(ш фсек)

106

2 • 10е

Рис. 1

2.3.2. Имеем два узла интерполяции ао и а1 со значениями Уо и У1, узлы аппроксимации Ьо, Ь1, .Ьм с известными значениями функции уо, у1, ..., ум. Найдем многочлен Рз(х).

Здесь интерполяционный многочлен — линейная функция

гр ( \ Лг Х - а1 , Х - ао

Т1(х) = Уо--+ Уг

ао - а1 а1 - ао Искомый многочлен ищем в виде

Рз(х) = Т1 (х) + (Ах + А1) (х - ао) (х - а1).

Для искомых коэффициентов получим

Л1 А = Д2

X' А2 = ~Л'

А1

(8)

Рис. 2

где

м м

д = 2 52 12 (Ьт1 - ав)2(Ьт1 - ах)2(Ьт2 - ао}2(Ьт2 - ах)2(Ьт1 - Ьт2)2,

т\ = 1 т2 = 1

мм

А1 = 52 52 Ьш1 (ут1 - Т1(Ьт1)) (Ьт1 - ао) (Ьт1 - а1) (Ьт2 - ао)2(Ьт2 - а1) (Ь^1 - Ь^),

т1 = 1 т2 = 1

мм

Д2 = 52 52 (Ут1 - Т1(Ьт1)) (Ьт1 - ао)(Ьт1 - а1)(Ьт2 - ао)2(Ьт2 - а1)(Ьт1 - Ьт2 )•

т1 =1 т2 =1

Пример 2. Даны два узла интерполирования ао = 1, а1 = 2,5 со значениями У = 1, 5, У1 = 1, 0. Точки на рис. 2 показывают значения функции в узлах аппроксимации из табл. 1, значения функции в узлах интерполирования показаны квадратиками.

ТАБЛИЦА 1. Данные для аппроксимации функции

Ьт 1,1 1,2 1,3 1,4 1,6 1,8 2,0 2,2 2,3 2,4

ут 1 0,45 0,4 0,25 0,2 0,45 0,9 1,2 1,25 1,2

Вычисления по формулам (8) дают значения А1 = -3.582, А2 = 7.692. На рис. 2 показан график функции Р3(х).

3. Аппроксимация пучка кривых с общими точками

3.1. Постановка задачи об аппроксимации пучка кривых с общими точками

Применим интерполяционно-аппроксимационные многочлены для решения задачи наилучшего приближения нескольких функциональных зависимостей по методу МНК, при условии, что графики этих многочленов представляют собой пучок кривых, проходящих через несколько общих точек, абсциссы которых известны, а ординаты следует найти из условия наилучшего приближения. Требование совпадения значений всех функциональных зависимостей может вытекать, например, из физического смысла самих функциональных зависимостей.

Имеем задачу: найти Я зависимостей, каждая из которых приближает заданное множества данных и проходит через некоторые точки с известными абсциссами. Для г-й зависимости даны узлы аппроксимации Ь^, •••, ьМ^. . Количество и расположение узлов аппроксимации для разных зависимостей могут совпадать, а могут быть и различными. Для г-й функциональной зависимости в узлах аппроксимации известны значения у(г), •••, уМ}г.

Заданы узлы совпадения а0, аь ..., ам, в которых значения всех Я искомых функциональных зависимостей совпадают. Заметим, что некоторые узлы аппроксимации могут совпадать с узлами совпадения. Ординаты точек совпадения неизвестны и подбираются МНК наилучшим образом.

Отметим, что использование МНК для аппроксимации каждой функциональной зависимости по отдельности не обеспечит совпадения значений функций в узлах совпадения.

3.2. Построение интерполяционно-аппроксимационного пучка кривых с нахождением наилучшим образом общих точек

Будем искать для аппроксимации г-й зависимости интерполяционно-аппроксима-ый многочлен Р^^х). Обозначив н

узлах совпадения а0, аь ..., ам чер{ представление (2), с учётом формулы (3)

ционный многочлен Р^^(х). Обозначив неизвестные значения функциональных зависимостей в узлах совпадения а0, аь ..., ам через У0, У1,..., Ум, запишем для многочлена Р^ (х)

Кг

Р%>(х) = Тм (х) + ВД^ А^хк-\

к=1

где многочлен Тм (х) задан формулой (1), но только значения Уп здесь неизвестны и будут найдены далее.

Потребуем, согласно МНК, чтобы искомые многочлены обеспечивали минимум следующей функции Р, аргументами которой являются коэффициенты многочлена Тм (х) и все коэффициенты А^:

1 Я 1 Мг _

р = Я Е М Е (у$ - р?®)2 = ((у(г) - р^)))2}. (9)

г=1 г т=1

При введении функции Р в формуле (9) использованы два вида усреднения. Первый вид усреднения, обозначенный чертой сверху над выражением, означает усреднение в данном случае квадрата отклонений заданных значений от одной функциональной зависимости по всем узлам аппроксимации. Второй вид усреднения, обозначенный скобками (...) — это усреднение величины, полученной при усреднении первого вида усреднения по всем функциональным зависимостям. Смысл функции Р — это среднее значение квадрата отклонения заданных значений функциональной зависимости от значений искомого интерполяционно-аппроксимационного многочлена в узле аппроксимации по всем узлам аппроксимации и по всем функциональным зависимостям. Имеем общее правило восстановления формул вычисления средних значении: для каждой г-й функциональной зависимости

1 Мг

/(в(г)) = м Е f(©,

г т=1

усреднение по всем функциональным зависимостям

1 я 1 Мг

I (вм>) = я^мХ ?

г=1 г т=1

Аргументами функции Р являются искомые значения ординат общих точек функциональных зависимостей У1; У2, ..., Ум и коэффициенты всех многочленов А^ при г =1, 2, ... Я, а для каждого значения г имеем к = 1, 2, ..., Кг. Всего этих аргумен-

Я

тов N +52 Кг.

г=1

Необходимым условием экстремума функции Г является равенство нулю всех частных производных по всем её аргументам.

Для каждой г-й функциональной зависимости при ] = 1, 2, ..., Кг имеем

Мг / Кг

дГ 2 1

---£ г1 у(т) - Тм(Ь$) - ЩЬ^)^ (№)к-1АР

дА] ' т=1 \ к=1

Приравнивая к нулю каждую эту производную и опуская сомножитель 2/Я, получим для каждой г-й зависимости СЛАУ, которая с использованием матричных обозначений принимает вид

с (г)х( г) + Н = W (г). (10)

В матричном уравнении (10) введены: матрицы-столбцы

Х(г) = (А1г), А2Г), ..., АКТ) )Т; V = (Ус, У1,...,У„ )Т;

W(r) = (^(г), W2(r), ..., ШК]Г)Т, где Ш(г) = V(r)(B(r))j-1П(B(r)); матрицы

С(г) = (С^), размер которой Кг х Кг, где С^ = (В(г)У+к-2(П(В(г)))2;

Н(г) = (И{т)), размер которой Кг х N, где ИЩ = (ВМ^^^^П(В(г)). СЛАУ (10), при условии неособенности матрицы С(г), позволяет выразить матрицу-столбец Х(г) через матрицу-столбец V

Х(г) = (G(r))-1(W(r) - H(r)V). (11)

Вычислим производные функции Г по аргументам У.г при г = 0, 1, ..., N

дг 2 * 1 М вдт^ / (Г) пп(ь(тт)) К

т - V Уп - ^ У: (ьm¡))k-lAkr)

ду Я ^ Мг т=1 пг(аг) \Ут Пп(ап) ' К"т

г=1 т=1 4 ' \ п=0 4 ' к=1

Приравнивая к нулю каждую такую производную, запишем получаемую при этом СЛАУ в матричном виде

В(г)Х(г) + Р(г) V = У(г) , (12)

где введены: _

матрица-столбец У(г) = (Ус, Уь ..., Ум)Т, V = ;

матрицы Ю(г) = (Д(к)) и Р(г) = (Р^), размер которых соответственно N +1) х КГ и N +1) х (N + 1), с элементами

;(г) = / (В(г))к-1 П(В(г))Пг(В(Г)Л р(г) = /^Ш . ^п(В(г)) гк ^ (В ) ^(В ) П(аг) ГРгп = \ Пг(аг) Пп(ап) / "

Объединение двух матричных СЛАУ (10) и (12) в одну СЛАУ позволяет найти искомые матрицы-столбцы Х(г) и V. Для решения этой СЛАУ воспользуемся соотношением (11) и исключим матрицу-столбец Х(г) из СЛАУ (12), тогда получим линейное матричное уравнение для матрицы-столбца V

Б(г)^1г))-1^(г) - Н^) + Р(г^ = У(г).

После преобразований найдём матрицу V:

V = (Р(г) — В(г)^(г))-1Н(г))-1(У(г) — Б(г) (G(r))-1W(r)), а матрицы-столбцы Х(г) вычисляем по формуле (11).

4. Частные случаи аппроксимации пучка кривых с общими точками

4.1. Аппроксимация данных пучком линейных функций, исходящих из одной точки с заданной абсциссой

Пусть пучок из Я прямых исходит из точки, абсцисса которой а0. Для прямой с

(г)

номером г, где г принимает все значения от 1 до Я, даны узлы аппроксимации Ь\ , ..., ьМд , в узлах аппроксимации заданы значения у(г), ..., уММг. Уравнения прямых ищем в виде

у = р(г)(х)= Уо + А(1Г) (х - ао).

Здесь искомыми величинами являются величина У0 и все коэффициенты А^, которые ищем из условия минимума функции

1 Я 1 Мг

Р = Я Е МтХ Ут - Уо - А^ - ао)2.

г=1 т=1

Это условие приводит к СЛАУ вида (4) и (6) и здесь имеет вид

А<г) (В(г) - ао)2 + Уо (В(г) - ао) = У(г)(В(г) - ао), г =1, 2, ..., Я; Я = А[г)(В(г) - ао) + Уо = (У*) ,

(13)

где

1 Мг 1 Мг

(В(г) - ао)2 = № - ао )2, (В(г) - ао) = — £ (Ь^ - ао),

т=1 т=1

1 Мг 1 Я 1 Мг

у^да«-^ = ^ ут<ьт - ао)2. (ум) = я ^ —. £ ^.

г т=1 г=1 г т=1

Разрешая первые Я уравнений в СЛАУ (13) относительно всех величин А^ будем

иметь

(г) _ У(г) (В(г) - ао) (В(г) - ао)

А1) = 4 - Уо = при г = 1, 2,..., Я. (14)

1 (В(г) - ао)2 (В(г) - ао)2

Исключая величины А^ из последнего уравнения в СЛАУ (7), получим линейное уравнение с неизвестной величиной Уо, решая которое получим

Ао

Уо = А, (15)

где

А = 1 - ( У(В(Г) - ао)) \ , Ао = -IУ(Г)(В(Г) - ао) ■ (В(Г) - ао)'

(В(г) - ао)2 / \ ' \ (В(г) - ао)2

(В(г) - аоП _ 1 Я 1 и=11

М

(Ь(т) - ао)

<В(Г) - ао)7 Я £ М' £ (ьт> - ао)2

т=1

2

и,

Мг

У(г)(Б(г) - ао) • (В(г) - ао) \ = 1 Л _1

Е (ьт - ао) • Ё (&т - ао)

т=1

т=1

(В(г) - ао)2

Г=1

После нахождения величины Уо значение А^ для всех г вычисляем по формуле (14)

Мг

Мг

Е № - ао)2

т=1

Рис. 3

Пример 3. На рис. 3 показано практическое использование формул (14) и (15). Здесь помечены точки аппроксимации для нахождения трёх линейных функций. Задана также абсцисса точки пересечения этих прямых ао = 0.

Штриховыми линиями показаны прямые, найденные с использованием стандартного МНК для каждой прямой поотдельности. Из рис. 3 видно, что штриховые линии в одной точке не пересекаются, тем более они не пересекаются в точке с заданной абсциссой. Пучок линий, показанный сплошными линиями, уравнения которых рассчитаны по формулам (8) и (9), аппроксимирует точки с учётом, что все линии пересекаются в одной точке с заданной абсциссой ао = 0.

4.2. Аппроксимация пучком парабол, имеющих две общие точки

Имеем две точки с заданными абсциссами ао и а1, через которые проходят Я искомых парабол, аппроксимирующие заданные множества точек. Уравнение г-й параболы имеет вид

рМ(х) = Уо+ У1 ^^ + А(1г)(х - ао)(х - а1).

ао - а1 а1 - ао

Для нахождения искомых величин Уо, У1 и А11), А12), ..., А^ имеем соотношение (11), которое здесь имеет вид

А(г) = 1 (у(г)(В(г) - ао)(В(г) - а1) -

1 (В(г) - ао)2(ВМ - а1)2 ^

-Уо ^^(В^ - ао)(В(г) - а1) - У ■^^(ВМ - ао)(В(г) - ) (16) ^о(ао) ^1(а1)

при г =1, 2, ..., Я.

СЛАУ (12) принимает вид

где

Л.

ЛссУс + Л01 Yi = Фо,

ЛюУо + Л11 Y = Ф1, 1

(17)

ij

Qi(ai)Qj (aj)

(B(r))) -

a(BM)(BM - ao)(B(r) - ai) • fy(B^BM - ao)(Br - ai)

Фг

(B(r) - ao)2(BM - ai)2 1 /7YQi(B(r))) -

Qi(ai)

Y(B(r) - ao)(B(r) - ai) • Q3(B(r))(B(r) - ao)(B(r) - ai)

(ВМ - ао)2(В(г) - а1)2 Решение СЛАУ (16) даёт следующие значения для искомых величин

У = Ао У = А1

Уо = "д , У1 = Х,

где А = ЛооЛц - Ло1Лю, Ао = ФоЛп - ФоЛо1, А1 = Ф^оо - ФоЛю.

(18)

(r)

Использование формулы (16) позволяет вычислить значения коэффициентов А Пример 4. На рис. 4 показаны результаты применения формул (16) и (18) для нахождения трех парабол, проходящих через две точки с заданными абсциссами ао = 0 и а1 = 2. Каждая из парабол аппроксимирует множество точек, имеющих на рис. 4 одинаковые обозначения.

РИС. 4

Литература

[1] Демидович Б.П., Марон И.А. Основы вычислительной математики. — М.: Наука, 1970,— 664 с.

[2] Гончаров В.Л. Теория интерполирования и приближения функций. — М.-Л.:, ГТТЛ, 1954. —325 с.

[3] Кобзарь А.И. Прикладная математическая статистика. Для инженеров и научных работников. — М.:, Физматлит, 2006. — 816 с.

[4] Hubble E. A relation between distance and radial velocity among extra-galactic nebulae// Proc. Nath. Acad. Sci. USA,—1929.—15.—P. 168-173.

ON THE ISSUE OF INTERPOLATION APPROXIMATION POLYNOMIAL CONSTRUCTION

V. D. Lukyanov

Interpolation approximation polynomial which solves the problem of data approximation by least squares method was built. Its values coincide with the given values at certain points. Beam curves data approximation problem was solved with interpolation approximation polynomial by the selection of the best common points. The examples of finding the interpolation approximation polynomial were given.

Key words: polynomial, interpolation, approximation, least squares method, interpolation approximation polynomial, beam curves.

V.D. Lukyanov - JSC Avangard, Saint Petersburg, Russia, professor, D.Sc., [email protected].