Интерполяция наипростейшими дробями Текст научной статьи по специальности «Математика»

Таким образом, имеем следующий список бифурцирующих наборов (ЬГ-раскладов) критических точек функции Ш: (9,12,4), (5,12, 8), (8,8,1), (4,4,1), (5,8,4), (5,4,0), (1,4,4), (1, 0,0). Им соответствуют (с точностью до поворота на угол п/4) графы, приведенные на рис. 4 (комплексы Морса [11]).

Рис. 4. Комплексы Морса

Библиографический список

1. Изюмов Ю.А., Сыромятников В.И. Фазовые переходы и симметрия кристаллов. М: Наука, 1984. 247 с.

2. Широков В.Б., Юзюк Ю.И., Dkhil B., Леманов В.В. Феноменологическое описание фазовых переходов в тонких пленках В.Н. BaTiO3 // Физика твердого тела. 2008. Т. 50, вып. 5. С. 889-892.

3. Даринский Б.М., Сапронов Ю.И., Царев С.Л. Бифуркации экстремалей фредгольмовых функционалов // Современная математика. Фундаментальные направления. М.: МАИ, 2004. Т. 12. С. 3-140.

4. Красносельский М.А., Бобылев Н.А., Мухамади-ев Э.М. Об одной схеме исследования вырожденных экстремалей функционалов классического вариационного исчисления // Докл. АН СССР. 1978. Т. 240, № 3. С. 530-533.

5. Darinskii M.M., Sapronov Yu.I., Shalimov V.V. Phase transitions in crystals characterized by polarization and deformation components of the order parameter // Ferroelectrics. 2002. V. 265. P. 31-42.

6. Даринский Б.М., Дьяченко А.А., Сапронов Ю.И.,

УДК 517.538.52, 517.538.7

ИНТЕРПОЛЯЦИЯ НАИПРОСТЕЙШИМИ ДРОБЯМИ

Е.Н. Кондакова

Владимирский государственный университет, кафедра функционального анализа и его приложений E-mail: kebox@mail.ru

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

Чаплыгин М.Н. Фазовые переходы в доменных границах ферроиков // Известия РАН. Сер. физическая. 2004. Т. 768, № 7. С. 920-926.

7. Борисович Ю.Г., Звягин В.Г., Сапронов Ю.И. Нелинейные фредгольмовы отображения и теория Лере -Шаудера // Успехи мат. наук. 1977. Т. 32, вып. 4. С. 3-54.

8. Арнольд В.И., Варченко А.Н., Гусейн-Заде С.М. Особенности дифференцируемых отображений. Классификация критических точек каустик и волновых фронтов. М.: Наука, 1982. 304 с.

9. Брекер Т., Ландер Л. Дифференцируемые ростки и катастрофы. М.: Мир, 1977. 208 с.

10. Сапронов Ю.И., Хуссаин М.А. Угловые особенности гладких функционалов в задачах о прогибах упругих балок и зарождении нелинейных волн // Труды Воронеж. зимн. мат. школы. Воронеж: Изд-во Воронеж. ун-та, 2004. С. 155-167.

11. Постников М.М. Введение в теорию Морса. М.: Наука, 1971. 568 с.

Interpolation by the Simplest Fractions

E.N. Kondakova

Vladimir State University,

Chair of Functional Analysis and its Application E-mail: kebox@mail.ru

The interpolation by means of real simplest fractions is considered. There are offered a different ways of intepolating simpest fractions construction with distinct real nodes. Necessary and sufficient conditions of existence and uniqueness of interpolating simplest fractions are received. Interpolation of constants is in detail investigated; in this case the estimation of an error of interpolation on Chebyshev‘s system of nodes is received.

Ключевые слова: интерполяция, наипростейшие дроби.

Keywords: interpolation, simplest fractions.

© Е.Н. Кондакова, 2009

ВВЕДЕНИЕ

Определение. Наипростейшей дробью (н.д.) п-й степени называется рациональная функция вида

п 1

Лп (г) = -----------------------------------------, п > 1, (1)

г — гк к = 1 к

где точки гк € С — полюсы функции Лп — не обязательно геометрически различны. Таким об-

п

разом, н.д. есть логарифмическая производная некоторого многочлена фп(г) = С П (г — гк), т.е.

к=1

Лп (г) = ^п(г)/^п(г).

Как аппарат приближения н.д., по-видимому, впервые применялись в работах [1] и [2] (аппроксимация аналитических функций в интегральных пространствах Ь1 (С, р ■ шев2) относительно плоской меры шев2 с определенным весом р(г) на ограниченных областях С С С). В работе [3] предложена явная конструкция н.д., аппроксимирующих функции в равномерной метрике, и доказан аналог теоремы Мергеляна (т.е. всякую функцию, непрерывную на не разделяющем плоскость компакте К и аналитическую в его внутренних точках, можно сколь угодно точно аппроксимировать на К посредством н.д.). Дальнейшее развитие эта тематика получила в других работах [4-12]. Оказалось, что аппроксимативные свойства н.д. и многочленов во многом сходны. Так, при аппроксимации посредством н.д. справедливы аналоги классических теорем Джексона, Уолша. Вместе с тем имеются и существенные отличия, связанные в основном со своеобразной нелинейностью н.д. Например, операции вычитания н.д. или умножения их на не натуральное число выводят из класса н.д. Отметим также, что н.д. наилучшего равномерного приближения может быть неединственной, соответствующие примеры построены П.А. Бородиным (устное сообщение).

Аппроксимации нашли различные приложения в численном анализе и в теории потенциала [1012]. Последнее связано с тем, что н.д. (1) задают (с точностью до постоянных множителей и операции комплексного сопряжения) плоские поля различной природы, создаваемые равновеликими источниками, расположенными в точках гк. Задача аппроксимации здесь состоит в том, чтобы определить расположение источников гк, создающих заранее заданное поле.

Конструкция, предложенная в [3], использует кратную интерполяцию Паде. В данной заметке изучается интерполяция посредством н.д. с простыми вещественными узлами. В разделе 1 предлагается несколько способов построения интерполирующих н.д. В разделе 2 рассматривается интерполяция вещественных констант на отрезке действительной оси. Основная трудность здесь возникает при оценке погрешности интерполяции, это связано с тем, что полюсы интерполирующих н.д. могут близко подступать к отрезку интерполяции или даже попадать между узлами интерполяции. В разделе 3 получены некоторые необходимые и достаточные условия существования, единственности интерполяционных н.д. Отметим, что интерполяционные н.д. не всегда существуют, а если существуют, то необязательно единственны; и для существования, и для единственности важны также свойства интерполируемой функции. Другими словами, важен выбор не только узлов, но и таблицы интерполяции.

1. ПОСТРОЕНИЕ ИНТЕРПОЛЯЦИОННЫХ НАИПРОСТЕЙШИХ ДРОБЕЙ

1.1. Поскольку дробь (1) имеет п свободных параметров {гк}, то интерполирующую н.д. естественно строить по п узлам (для многочлена степени п таких узлов п + 1). Пусть задана таблица интерполяции (хк,Ук), к = 1,...,п, п > 2, с попарно различными вещественными узлами интерполяции хк и вещественными значениями ук. Требуется найти многочлен <^п (х) = хп + дп-1хп-1 +... + q0 с вещественными коэффициентами, такой что для н.д. Лп(г) = ^п(г)/фп(г) выполнены условия интерполяции

Лп (хк)= Ук, к = 1,..., п. (2)

Для этого перепишем условия (2) в виде ^п(хк) — ук^п(хк) = 0 и получим систему линейных уравнений относительно неизвестных qj:

fyi yixi - 1

У2 У2х2 - 1

Vyn ynxn 1

yi^n 1 - (n - ік Л / qo \

У2^П -1 - (n - 1)xn -2

Уп x

n —1

qi

^nxn 1 - y1xn ^ nxn-1 - y2xn

- (n - 1)xn V \qn—1/

\nxn 1

(3)

Уп xn/

Через Ап = Ап(ж і, уі,..., жп, уп) обозначим матрицу этой системы. Будем записывать систему (3) кратко в виде Ап ■ q = Вп. Отсюда сразу получается

Теорема 1.1. Решение Лп = ФП/Фп интерполяционной задачи (2) существует и единственно тогда и только тогда, когда det Ап = 0. При этом коэффициенты многочлена фп получаются как решение системы (3).

Приведем еще один способ построения интерполяционной н.д. Положим

Dn (x,y) :=

У1 У1 x1 - 1

y2 y2x2 - 1

У1 xn 1 - (n - 1)xn 2

У2xn-1 - (n - 1)xn —2

y1 xn - nxn 1

n П — 1

У2 xn - nx2

yn yn xn 1 •

y yx - 1 •

,71 ^y»n 1 _____ __ l'\o">n 2 У qr,n _____________ ■то nr

n n n n n n

yx

n—1 - (n - 1)xn—2

yxn - nx

n-1

n

n-1

(4)

Предположим, что главный минор ап, находящийся в левом верхнем углу определителя (4), отличен от нуля, т.е. ап = det Ап = 0.

Теорема 1.2. Если det Ап = 0, то детерминантное уравнение (ж, у) = 0 задает единственную интерполяционную н.д. у = Лп(ж), удовлетворяющую условию (2).

Доказательство. Имеем Дп(ж&,ук) = 0 при всех к = 1,...,п, поскольку в определителе £п(ж&,ук) имеется две одинаковые строки. Следовательно, все точки (ж&, ук) лежат на кривой £п(ж,у) = 0. Покажем, что функция у = Л(ж), неявно задаваемая детерминантным уравнением, является наипростейшей дробью. Разложим определитель Дп(ж, у) по последней строке и запишем уравнение (ж, у) = 0 в виде

«1 у + «2 (уж — 1) + а3(уж — 2ж) + ... + ап(ужп — пжп ) = 0,

где — алгебраические дополнения к соответствующим элементам последней строки,

ап = det Ап = 0. Отсюда получаем

y=

nx

n-1

+ (n - 1)qn—1xn 2 + • • • + 2q2x + q1

xn + qn—1 xn—1 + • • • + q2x2 + q1x + qo

где qk = ak/an. Таким образом, детерминантное уравнение Dn(x, y) = 0 действительно определяет некоторую н.д. R(x). Единственность следует из теоремы 1.1. □

1.2. Рассмотрим вопрос об интерполяции константы f (x) = с = const. В этом случае имеем yk = с, и задача (2) принимает вид

R ( ) = Qn(xk) = с

Rn (xk) Qn(xk) с

k = 1, • • • , n,

(5)

где многочлен фп(ж) = жп + дп - 1 жп 1 + ... + до требуется определить. Система (3) в данном случае принимает вид Апд = Вп, где элементы матриц Ап = (а&>то) и Вп = (6&) вычисляются по формулам

«к,т = еж^-1 — (т — 1)жт-2, && = пжП-1 — сжП, к, т = 1,..., п.

Несложно проверить, что элементарными преобразованиями столбцов определитель матрицы Ап сводится к определителю Вандермонда. При этом det Ап = еп П (ж& — жт) и, значит, det Ап = 0.

1<&<ш<п

Отсюда и из теоремы 1.1 получается

Теорема 1.3. Решение Лп (ж) интерполяционной задачи (5) при е = 0 всегда существует и единственно. При этом коэффициенты многочлена можно получить как решение системы (3).

В случае интерполяции констант для вычисления этих коэффициентов можно получить более удобные формулы. Докажем следующее утверждение.

Теорема 1.4. Коэффициенты многочлена (х) = хп + дп-1жп-1 + ... + до из задачи (5) вычис-

ляются по формуле

з

qn-j =J](-l)fccfc j j I,' afc> j = l,...,n, ao = 1, (6)

k=0 ( j )-

где ak, k = 1,..., n суть элементарные симметрические многочлены

(x1 7 x2 7 ... 1 xn) ^ ^ xjl xj2 ' ' ' xjk . (7)

Доказательство. Заметим, что из (5) получается тождество

П

(x) _ cQn(x) = _cJJ(x _ xk). (8)

k = 1

Перепишем его в виде

n n

У Л^П k ((n _ k + 1)qn-k + 1 _ cqn-k) = _c ^ ](_1)kakxn k 1 qn = 1.

k=1 k=1

Далее, приравнивая коэффициенты при одинаковых степенях переменной x, получаем рекуррентное соотношение между коэффициентами многочлена Qn:

qn = 1, qn-j = c-1(n _ j + 1)qn-j+1 + (_1)j aj, j = 1, ...,n. (9)

Отсюда последовательно находим

qn-1 = c-1 (nqn _ ca1) = c-1 n _ a, qn-2 = c-2n(n _ 1) _ c-1a1 (n _ 1) + a2,

j

и, вообще, qn-j = c-jn(n _ 1)... (n _ j + 1) + ^(_ 1)kck-j(n _ k)... (n _ j + 1)ak, j = 1,n.

k = 1

Тем самым формула (6) доказана. □

В заключение приведем еще одну формулу для вычисления Qn(x), которая легко получается как решение дифференциального уравнения (8):

Qn(x) = -с ^ J e CP (t)dt + qoj eca:, P (t) = (t - xi) ••• (t - Xn), (10)

где параметр q0 подбирается так, что выражение Qn (x) в (10) не содержит экспонент. А именно

n

qo = E(-1)k+1 ^k(n - k)!ck-n-i. k=0

Если, в частности, узлы Xk являются корнями многочлена Чебышева, то P(t) = 21-n cos(n arccost), и из (10) получаем

Л 1 л 2 2 2 1 . 3 3 2 ( 6 3 \ ( 6 3 \

Qi (x) = x + -, Q2 (x) = x + - x + - -, Q3 (x) = x + - x + -- x + - —

с с с2 2 с \ с2 4 / \ с3 4с /

и т.д.

2. ИНТЕРПОЛЯЦИЯ КОНСТАНТ

Хорошо известно, что чебышевская система узлов {xk} является в определенном смысле наилучшей в случае интерполяции многочленами [13]. Используем ее для интерполяции константы f (x) = с посредством н.д. на отрезке [-1,1]. Итак, пусть узлы xk Є [-1,1] являются нулями многочлена Чебышева Tn(x) = cos(narccosx), т.е.

X

21-пТп(ж) = (ж — жі) ■ ■ ■ (ж — жп), жк = оо^(2П(—1 + 2к)^ , к = 1,..., п.

Через Лп(ж) обозначим интерполяционную н.д. и положим

, ч о і \ УП(ж) — С^п(ж)

Рп(ж) = Я„(ж) — С = ------п ( )------,

Уп (ж)

где Уп(ж) = жп + дп-іжп-1 + ... + до. Тогда, очевидно,

п

П (х — ) гр ( Ч

р (х) = С к=1 = С Тп (ж) (11)

Рп(ж) = —С Уп(ж) = — 2п-1 опм. (11)

Теорема 2.1. Пусть с — положительная постоянная, 0 < с < 15/31. Тогда погрешность интерполяции этой постоянной посредством н.д. по чебышевской системе узлов оценивается следующим образом:

и и с 1 с 1 — с ^ ,1П.

ІІРпІІ < 2^ ■ |дп(—1)| < 22п_1п! ■ Г—2С’ п - . ( )

Предварительно докажем две леммы.

Лемма 2.1. Если Уп(—1) — а21-п с произвольной постоянной а > 1, то функция Уп(ж) монотонно возрастает на отрезке [—1,1].

Доказательство. Если выполнено предположение леммы, то, переписав равенство (11) в виде

УП(ж) = сУп (ж) — с21-пТп (ж), (13)

получим УП(—1) — сУп(—1) — с21-п — с(а — 1)21-п = 2є, где є > 0. Выберем величину 6 = 6(е) > 0,

такую что |УП(ж) — УП(у)| < є при всех ж, у Є [—1,1] с |ж — у| < 6. Пусть натуральное N удовлетворяет

условию 2N-1 < 6. Разобьем отрезок [—1,1] на N равных отрезков

Дк = [жк-1 ,жк], жо = —1, жк = —1 + 2 к^1, к = 1,...,N.

Пусть ж Є Д1. Тогда из неравенств УП(ж0) = УП(—1) — 2є и |УП(ж0) — УП(ж)| < є находим

УП(ж) — УП(жо) — |УП(жо) — УП(ж)| — 2є — є = є > 0.

Таким образом, УП(ж) > 0 на Д1 и, значит, функция Уп(ж) возрастает на Д1. Отсюда получаем, что Уп(ж1) > Уп(ж0) — а21-п и из (13) получаем оценку УП(ж1) — сУп(ж1) — с21-п — с(а — 1)21-п = 2є. Дословно повторив предыдущее рассуждение для ж Є Д2, получим, что

(ж) — (ж1) — |УП(ж1) — УП(ж)| — 2є — є = є > 0.

Так что и на Д2 производная УП(ж) положительна, и функция Уп(ж) возрастает. Повторив аналогичную процедуру N раз, придем к заключению леммы 2.1. □

Лемма 2.2. При 0 < с < 1/2 имеет место неравенство

1 — 2с 1 — 2с

Уп(—1) > ------Сп!с-п > -----Сп!2п > 0. (14)

1 — с 1 — с

Доказательство. Положим Т(ж) = 21-пТп(ж). Тогда из равенства (13) последовательно получаем УП (ж) = сУп(ж) — сТ (ж), УП (ж) = с2 Уп(ж) — с2Т (ж) — сТ ' (ж)

и, наконец, УПп)(ж) = п! = спУп(ж) — спТ(ж) — сп-1 Т'(ж) — ... — сТ(п-1)(ж). Отсюда находим

п- 1

-п

<Э„(—1) — п!с-п — £С-|Т<к)(—1)|- (15)

к=0

Для оценки суммы в (15) воспользуемся неравенством В.А. Маркова [14] для max-нормы на отрезке [a, b] производной многочлена P степени n, ограниченного на [a,b] некоторой величиной M:

.. „(k)ll M2k 2kk! 2 2 2, „ ,2, M2k 2kk! n(n + k - 1)!

IIP (k)y < Ti--ГГ Т^ГТГП2 (n2 - 12) ... (n2 - (k - 1)2) = Tl-ГГ77ГГТ7 -TTT^ •

11 11 “ (b - a)k (2k)! v (b - a)k (2k)! (n - k)!

В нашем случае имеем [a, b] = [-1,1], P = T, M = 21-n и, следовательно,

im(kw ,м „ k -^k n k! (n + k - 1)!

|T(k)(-1)| < ckak, ak :=c k2k2^Щ! V - k)! •

Пусть 0 < c < 1. Тогда

an-1 = c-n+1n!, ak = c 22k + A, < c< 1, k = 0,..., n - 2, (16)

ak+i n2 - k2

откуда an-2 < can-1, an—3 < c2an-1, и, вообще, ak < cn—k—1 an-1. Отсюда с учетом того, что an-1 = c-n+1 n!, получаем

n —1 n —1 „

Qn(-1) > n!c—n - V ak > n!c—n - an—^ cn—k—1 > n!c—n-----------------— n!c—n = n! c—n > 0,

1 - c 1 - c

k=0 k=0

что и доказывает лемму 2.2. □

Доказательство теоремы 2.1. Заметим, что при 0 < с < 15/31 и n > 2 выполнение условия леммы 2.1 следует из первого неравенства в (14). Поэтому утверждение теоремы 2.1 сразу получается из (11) и лемм 2.1 и 2.2:

|р (x)| < с 1 < с __________________1____ < ____с____ 1 - с □

П - 2n—1 |Qn(x)| 2n—1 |Qn(-1)| 22n—1n! 1 - 2c'

Теорема 2.2. Пусть 0 < c < 1/2. Тогда все полюсы интерполяционной н.д. Rn (x), n > 2, лежат вне круга |z| < 1.

Доказательство. Воспользуемся следующей оценкой корней многочлена.

Теорема (Энестрем - Какейя, [15]). Если коэффициенты многочлена Qn(x) = xn + qn_1 xn—1 + + ... + q0 строго положительны, то все его корни лежат в кольце

qk ^ I I ^ qk

mm --------- < |z| < max ----------, qn = 1.

k=°,...,n—1 qk+1 k=0,...,n—1 qk+1

Нам понадобится только первая оценка. Достаточно показать, что все коэффициенты qj многочлена Qn из (11) положительны и i5j := qn_j/qn—j+1 > 1, j = 1, ...,n, qn = 1. Воспользуемся формулой (6). Для этого заметим, что поскольку многочлены Tn(x) являются либо четными, либо нечетными функциями, то

[ n ]

T(x) = 21—nTn (x) = xn + £ a2j xn—2j,

j=1

где [■] — целая часть. Поэтому симметрические многочлены ak (7) от корней xk многочлена T с нечетными номерами равны нулю. С учетом этого формула (6) принимает вид

[ 2 ]

_ j n! | 2k_____________________j (n - 2k)! • i

qn—j = c j7------------ттг + > c ---------ГТ7-a2k, j = 1,...,n.

(n - j)! k=1 (n - j)!

Учтем, что все |хк| < 1, тогда из (7) получаются неравенства |ст&| < СП (биномиальные коэффициен ты). Поэтому при 0 < с < 1 и при всех ] = 1,... ,п имеем

qn _j > c _

(n - j)!

( [ 2 ] c2k \

2-£

k=0 (2k)!

n!

>c _j--------------—(2 - ch c) > 0. (17)

(n - j)!

n

Теперь докажем, что ^ > 1, з = 0, ...,п. Воспользуемся рекуррентными соотношениями (9). При нечетных з = 2т + 1 имеем =0 и, следовательно, при 0 < с < 1/2 из (9) получаем: ^ = (п — з + 1)с-1 > 2. Пусть теперь з = 2т. Тогда из (17) и из неравенств ^ | < СП получаем

§ = п — з + 1 + ^ > п — з + 1 СП (п — з + 1)! с^-1 =

= (п - з + 1М - - 77^—-г-т = (п - 3 + ВД(с).

?п-^+1 с п! (2 — еЬ с)

1 с^-1

с з! (2 — еЬ с)

При з =2 имеем (с) > Р2(1/2) > 1, отсюда ¿2 = (п — 1)^2(с) > (п — 1) > 1. Из того, что функция (с) возрастает с ростом з > 2, находим ^ > 1 ■ (с) > Р2 (с) > 1. □

3. УСЛОВИЯ СУЩЕСТВОВАНИЯ И ЕДИНСТВЕННОСТИ ИНТЕРПОЛЯЦИОННЫХ Н.Д.

Рассмотрим сначала пример интерполяции с таблицей (х1 ,у1), (х2, у2). В силу теоремы 1.1 соответствующая интерполяционная н.д. Л2(х) существует и единственна тогда и только тогда, когда определитель ¿(х15 уь х2, у2) = det А2 отличен от нуля, где

л л I \ /^^1 У1х1— Л

А = А (Х1 ,У1,Х2,У2)= I.

\У2 У2х2 — 1 у

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

Пример. Для таблицы (0, —2), (1, 2) существует бесконечно много различных интерполяционных н.д. второго порядка. Они имеют вид

^ ^ х — Ь

#2(Ь; х) = 2-

?2 — 2Ьх + Ь ’

где Ь — произвольный отличный от нуля параметр. Здесь не выполнено условие теоремы 1.1 поскольку й(хі,уі,х2,у2) = 0, то есть вторая точка таблицы лежит на графике у = (х — 1/2)-1, проходящем через первую точку. Кроме того, система уравнений А2q = В2 имеет бесконечно много решений. Если же взять на этом графике вместо (1, 2) любую другую точку (х2,у2), то интерполяционных н.д. вообще не существует, так как указанная система не имеет решений.

Перейдем к случаю произвольного п.

Определение. Таблица интерполяции (хк, ук), к = 1, ...,п, с попарно различными узлами хк называется допустимой, если через точки (х&, у&) Є К2 проходит график единственной интерполяционной н.д. у = Лп(х). Другими словами, допустимая таблица удовлетворяет условию существования и единственности теоремы 1.1.

Пусть задана допустимая таблица интерполяции (х^, ук), к = 1,...,п. Тогда по теореме 1.2 детерминантное уравнение ^п(х,у) = 0 определяет единственную интерполяционную н.д. у = Лп(х). Выберем какую-либо новую точку (хп+і,уп+і), не лежащую на графике этой н.д. Тогда £п(хп+і,уп+і) = 0, или, что то же самое, det Ап+і = 0, где Ап+і = Ап+і(хі,уі,... ,Хп+і,уп+і) (см. обозначение после системы (3). По теореме 1.1 новая таблица (х^,ук), к = 1,...,п + 1, также является допустимой.

Пусть теперь точка (хп+і,уп+і) лежит на графике у = Лп(х), т.е. det Ап+і = 0. Тогда система уравнений (3) (где следует заменить п на п +1) по критерию Кронекера - Капелли [16] либо имеет бесконечно много решений, если ранги расширенной матрицы и матрицы Ап+і равны, либо, в противном случае, не имеет их вовсе.

Сравним ранги. Поскольку таблица (х&, у&), к = 1, ...,п, является допустимой, то ранг матрицы Ап+і равен п. Поэтому ранг расширенной матрицы системы (3) (где следует заменить п на п + 1) равен п, если и только если равен нулю определитель матрицы

Dn+1 (xn+1,yn+1 ) —

У1 • • • У1 хП 1 — (n- 2 1 П1 Л 1

У2 • • • У2хП-1 — (n- 1)хП-2

yn • У хП 1 • • УПХП (n- 1 ) хП 2 1)xn

П+1 • • • У«.+1 xn+l (n - 1)хП-2

(n + 1)хП — y1 ХП+1 (n + 1)xn — y2 ХП+1

(n + 1)xn — УПХП+1

,.n+1

■"n+1,

Значит, для существования бесконечного числа интерполяционных н.д. степени п +1 для таблицы (хк, ук), к = 1,...,п + 1, необходимо и достаточно, чтобы новая точка (хп+1 ,уп+1) удовлетворяла уравнению det +1 (х, у) = 0. Как и в теореме 1.2 доказывается эквивалентность

det Dn+1 (x,y) — 0

У — Rn+1(x)

(18)

где ДП+1 (х) есть некоторая н.д. степени п + 1.

Из предыдущих рассуждений получается

Теорема 3.1. Пусть (хк, ук), к = 1,... , п, является допустимой таблицей интерполяции.

Новая таблица (хк, ук), к = 1,..., п + 1, является допустимой тогда и только тогда, когда точка (хп+1 ,уп+1) не лежит на графике интерполяционной н.д. у = Лп(х), построенной по первоначальной таблице (хк, ук), к = 1,..., п.

Для новой таблицы (хк , ук), к = 1,..., п + 1, существует бесконечно много интерполяционных н.д. Дп+1 тогда и только тогда, когда точка (хп+1,уп+1) лежит на пересечении двух графиков у = Дп(х) и у = ЯП+1(х) (см. (18)).

Для новой таблицы (хк, ук), к = 1,..., п + 1, не существует ни одной интерполяционной н.д. Лп+1 тогда и только тогда, когда точка (хп+1, уп+1) лежит на графике у = Лп (х), но не лежит на графике у = ЛП+1. (х).

П

Работа выполнена при финансовой поддержке АВЦП РНПВШ (рег. номер 2.1.1/5568) и гранта РФФИ (проект 08-01-00648).

Библиографический список

1. Chui C.K. On approximation in the Bers spaces // Proc. Amer. Math. Soc. 1973. T. 40. С. 438-442.

2. Chui C.K., Shen X.C. Order of approximation by electrostatic fields due to electrons // Constr. Approx. 1985. T. 1. С. 121-135.

3. Данченко В.И., Данченко Д.Я. О равномерном приближении логарифмическими производными многочленов // Теория функций, ее приложения и смежные вопросы: Материалы школы-конф., посвящ. 130-летию со дня рожд. Д.Ф.Егорова. Казань, 1999. С. 74-77.

4. Долженко Е.П. Наипростейшие дроби // Теория функций, ее приложения и смежные вопросы: Материалы V Казанск. междунар. летней школы-конф. Казань, 2001. С. 90-94.

5. Косухин О.Н. Об аппроксимативных свойствах наипростейших дробей // Вестник Моск. ун-та. Сер. 1. Математика. Механика. 2001. № 4. C. 54-58.

6. Бородин П.А. Оценки расстояний до прямых и лучей от полюсов наипростейших дробей, ограниченных по норме Lp на этих множествах // Мат. заметки. 2007. Т. 82, № 6. С. 803-810.

7. Бородин П.А., Косухин О.Н. О приближении наипростейшими дробями на действительной оси // Вест. Моск. ун-та. Сер. 1. Математика. Механика. 2005. № 1. C. 3-8.

8. Данченко В.И., Данченко Д.Я. О приближении наипростейшими дробями // Мат. заметки. 2001. Т. 70, № 4. С. 553-559.

9. Данченко В.И. Об аппроксимативных свойствах

сумм вида ^2 ХкЛ,(Хкг) // Мат. заметки. 2008. Т. 83,

к

№ 5. С. 643-649.

10. Данченко В.И. Оценки производных наипростейших дробей и другие вопросы // Мат. сб. 2006. Т. 197, № 4. С. 33-52.

11. Фрянцев А.В. О численной аппроксимации дифференциальных полиномов // Известия Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2007. Т. 7, вып 2. С. 39-43.

12. Фрянцев А.В. О полиномиальных решениях линейных дифференциальных уравнений // УМН. 2008. Т. 63, № 3(381). С. 149-150.

13. Демидович Б.П., Марон И.А. Основы вычислительной математики. М.: Физматгиз, 1963. 660 с.

14. Фихтенгольц Г.М. Основы математического анализа. СПб.; М.; Краснодар: Лань, 2005. 512 с.

15. Прасолов В.В. Многочлены. М.:Физматлит, 2002. 453 с.

16. Курош А.Г. Курс высшей алгебры. СПб.: Лань, 2004. 432 с.