Герасимович(математический анализ)
.pdfПрименяя формулу Л агран ж а (5.9) к функции f(x) на [дс0; дс], по лучаем
y— Y = f'(l)(x — х 0) — f'(x о) (х — х0)=>
=> y - Y = ( f ' ( t ) - r ( x o ) ) ( x - X o ) ,
где £б]дсо; |
х [с :]а; |
Ь[. |
Применяя формулу Л агран ж а к функции |
f'(l) — f'(x о) |
на [дс0; |
5 ], |
находим |
y - Y = f " { b ) { l - X o ) { x ^ X o ) , Т.
|
Р и с . |
6.14 |
|
где £| в ]лг0; £[сг]а; |
Ь[. В последнем равенстве f"(%i) > 0 , а I — х 0 > 0, |
||
если |
х — х0 > 0 , или £ — х0 < 0 , если х — х0 < |
0 . Следовательно, |
|
y > Y , |
т. е. ординаты точек кривой больше ординат точек касатель |
||
ной при одной и той же абсциссе. Точки кривой у = |
f(x) на ]а; Ь[ леж ат |
||
выше точек касательной к кривой. |
|
График функции y = f(x) на ]а\ Ь[ — вогнутый (рис. 6.15). Доказательство выпуклости графика функции на ]а; Ь[ проводится
аналогично. <1
Пример 6.7. Найти интервалы выпуклости и вогнутости графика фугкции fix) =
= >/(1 +*)■ |
Д аи и ая функция определена, непрерывна и дифференцируема |
иа |
|||||
Р е ш е н и е . |
|||||||
R, причем |
|
|
|
|
|
|
|
|
ч |
|
— 2х |
г м = |
6 х2 — 2 |
|
|
|
Г(*)= „ |
, |
(1 + * 2)3 |
|
|
||
|
|
(1 |
+ х) |
|
|
|
|
Д л я определения |
интервалов |
выпуклости и вогнутости графика функции найдем |
|||||
решения неравенств f " ( x ) < . 0 |
и f"(x) > |
0 : |
|
|
|
||
|
ЗХ ~ * < 0 = > 3 * 2 — 1 < 0 = > * е ] — 1 / V ^ |
|
|
||||
|
(1 + * 2)3 |
|
|
|
|
|
|
Зл:2 — 1 > о= > зх2— 1 > о = > * е ] —■»; |
— i /V 3 [ U ] i/V 3 ; |
°о[. |
|
||||
О + * 2)3 |
|
|
|
|
|
|
|
Кривая, изображ аю щ ая |
график функции, иа |
] — оо; — \ / - \ / з [ — вогиутаи, |
на |
||||
] — 1/д/З; 1/-\/3[— выпуклая, |
иа ]1/д/3; |
оо[ — вогнутая. В точках х = |
— 1 / - ^ 3 t х = |
=1 /V 3 она меняет направление выпуклости, т. е. по определению это ее точки пере
гиба, в которых вторая производная f"(x) = 0, а в их окрестностях /" ( х) меняет знак. График функции изображен иа рис. 6.1.
Сформулируем достаточные условия существования точек пере гиба.
148
Теорема 6.7. Если для функции f(x) вто |
|||||||
рая |
производная |
f"(x) |
в некоторой точке х 0 |
||||
обращается в нуль или не существует и при |
|||||||
переходе через нее меняет свой знак, то |
|||||||
точка М{х0\ f(x о)) является точкой перегиба |
|||||||
графика функции. |
|
|
|
|
|||
р> Пусть |
f"(x0) = 0 |
или не |
существует. |
||||
Если |
И * ) < |
О |
в 0 6{хо - 0 ) |
и |
f"(x) > О |
||
в О6(*о + 0 ), |
то |
точка |
кривой |
с |
абсциссой |
||
х 0 отделяет интервал выпуклости |
от |
интер |
|||||
вала |
вогнутости. |
Если |
/ " ( * ) > О |
в |
О(,(х0_ |
— 0 ) и f"(x) < 0 в О6 (х0 + |
0 ), то эта точка от |
|
||||
деляет интервал вогнутости от интервала |
|
|||||
выпуклости кривой. В обоих случаях точка |
|
|||||
(х0; f(xо)) является точкой перегиба графика |
|
|||||
функции. <] |
|
|
|
|
|
|
Пример 6 .8 . Найти точки перегиба графика функ |
|
|||||
ции f(x).= JC4 — 4л:3. |
|
|
|
|
Р н с . 6.16 |
|
Р е ш е н и е . Д анная функция |
определена, непре |
|||||
|
||||||
рывна и дифференцируема |
на |
R, |
причем |
|
|
|
f'(x) |
= |
4л:3 - |
12л:2, f"(x) = |
12л:2 — 24*; |
|
|
f"(x) = |
0 =>л2 — 2x = 0=>xi = |
0 , x 2 = 2. |
|
Исследуем знак f"(x) в Oe(0) и Oe(2). В Оe(0 — 0) вторая производная f"(x) > 0, а в 0е(0 + 0) }"(х) < 0. Следовательно, точка 0 (0; 0) является точкой перегиба кривой (см. рис. 6.16). В Ое(2 — 0) f"(х) < 0, в Ое(2 -(- 0) f"(x\ > 0. Значит, точка М ( 2; — 16) является точкой перегиба графика функции f(x) = x* — 4х3 (рис. 6.16).
6.5. АСИМ ПТОТЫ ГРАФИКА ФУНКЦИИ
При исследовании поведения функции на бесконечности, т. е. при х -+ + о о и при x - i — оо, или вблизи точек разрыва второго рода часто оказывается, что расстояния между точками графика функции и точками некоторой прямой с теми же абсциссами сколь угодно малы. Такую прямую называют асимптотой графика.
Различают асимптоты вертикальные (т. е. параллельные оси орди нат) и наклонные. Частным случаем наклонной асимптоты является горизонтальная асимптота.
П рямая х = х 0 называется вертикальной асимптотой графика функции у = f(x), если хотя бы один из односторонних пределов в
точке х 0 равен бесконечности, т. е. lim f ( x ) = ± 0 0 |
или lim f(x) = |
х— —0 |
+0 |
=-+ - ОО .
Очевидно, что непрерывные на R функции вертикальных асимптот не имеют; такие асимптоты существуют только в точках разрыва второго рода функции у = f(x). Следовательно, для отыскания верти кальных асимптот графика функции надо определить те значения х, при которых хотя бы один из односторонних пределов функции бесконечен.
149
|
Пример 6.9. Найти вертикальные асимптоты графика |
функции |
</ = |
—я--— |
■ |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
х — |
1 |
|
Р е ш е н и е . Д ан н ая |
функция |
при |
* = ± 1 |
имеет |
|
разрыв |
второго |
рода, |
следо |
|||||||
вательно, прямые х = — 1 , х = |
1 являю тся |
вертикальными |
асимптотами |
ее графика |
|||||||||||||
(рис. 6.17). Вычисляя односторонние пределы в этих точках, имеем: |
|
|
|||||||||||||||
|
lim |
f ( x ) — |
lim |
— |
—р |
= |
+ о о , |
lim |
— —— = |
— oo, |
|
||||||
|
X -+ -— i — о |
|
|
|
|
|
X -*— i — о x |
— 1— |
i + o AT |
|
|
||||||
|
|
lifn . - |
r |
j |
= — ° ° . |
lim |
—j-J- |
|
|
|
|
|
|
||||
|
|
*“►1—0 |
X |
— |
1 |
|
|
„ . I |
I n |
x z — |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x - H + 0 |
X * |
|
|
|
|
|
|
|
|
|
П рямая |
у = k x + |
b |
называется |
наклонной |
(если |
k = |
0 — гори |
|||||||||
зонтальной) |
асимптотой |
графика |
функции |
|
y = f(x) |
при |
дс-»-+оо |
||||||||||
(дс->— оо), если функцию |
f(x) |
можно |
представить в виде f(x) = |
||||||||||||||
= |
kx + b + |
a(x), где |
а(дс) - > - 0 при дс-»- + |
оо |
(дс->— |
оо). |
|
|
|||||||||
|
Теорема 6.8. Д л я |
того чтобы график |
функции у = f(x) имел на |
||||||||||||||
клонную асимптоту у = |
kx + |
b, необходимо и достаточно, чтобы су |
|||||||||||||||
ществовали конечные пределы: |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
lim |
-i^L = |
k, |
lim (f(x) — kx) = b. |
|
|
|
(6 . 1 ) |
||||||||
|
|
x-*~± oo |
X |
|
|
x-*-dh oo |
|
|
|
|
|
|
|
|
|
||
> |
Необходимость. |
Предполагаем, |
что |
у = |
k x - \ - b — наклонная |
асимптота графика функции y = f(x). Тогда справедливо представ
ление f(x) = kx + b + «(■*)> |
гДе |
а (х ) — |
бесконечно малая при х - * - о о . |
|
Следовательно, |
|
|
|
|
l i m - M = Hm ( к + ± + Л Щ = к, |
||||
х-*-оо % |
х-*~оо \ |
X |
X / |
|
lim (f(дс) — kx) |
- lim (b + |
а(дс)) = b. |
||
x-*-oo |
|
x-*-oo |
|
|
Достаточность. Пусть существуют пределы (6.1), тогда, согласно следствию из теоремы 3.7. (см. § 3.4), из второго равенства (6.1) по лучаем f(x) = kx + b + а(дс) при дс—- ► + о о . Последнее представление функции ддс) означает, что прямая y = k x - \ - b является наклонной асимптотой графика функции y — f(x).
Итак, теорема доказана для случая дс-»- о о . Доказательство теоре мы для случая дс—*— оо производится аналогично. <
З а м е ч а н и е . При нахождении наклонных асимптот графика функции воз можны следующие случаи: 1 ) оба предела (6 . 1 ) существуют и не зависят от знака бесконечности, тогда прямая у = kx-\- b называется двусторонней асимптотой; 2 ) оба предела (6 . 1 ) существуют, но при х-*- ± оо они различны, тогда имеем две односто ронние наклонные асимптоты; 3) если хотя бы одни из пределов (6 .1 ) ие существует, то наклонных асимптот нет.
Пример |
6.10. Найти |
асимптоты |
линии |
у — -------- г- |
|
|
|||||
|
|
|
|
|
|
|
|
|
1 + х г |
|
|
Р е ш е н и е . |
Д ан н ая |
функция |
определена |
и непрерывна на R. Вертикальных |
|||||||
асимптот |
нет. Д л я нахождения |
невертикальиых |
асимптот вычисляем пределы (6.1): |
||||||||
lim |
X |
= |
lim |
—— 1 , |
■= |
0 , |
lim |
(f(x) — k x ) = lim |
■ |
....s = 0 . |
|
x-► dh oo |
|
x-*■± |
oo |
X |
)X |
x-*~zb oo |
x-*-± oo |
I |
-f- X |
Оба предела существуют и не зависят от знака оо. Следовательно, прямая у = 0 является двусторонней асимптотой (см. рис. 6 . 1 ).
150
|
Р и с . |
6.17 |
|
|
|
|
|
|
|
|
|
|
|
|
Пример 6.11. Найти |
асимптоты |
гиперболы х 2 — у 2 = |
|
1. |
|
|||||||||
Р е ш е н и е . Запишем |
функцию |
в явном виде: у = |
± |
~\Jx2— 1. Функция опреде |
||||||||||
лена и |
непрерывна иа |
множестве |
] — оо; |
— 1|U[1; |
°°[- Находим пределы (6 .1 ); |
|||||||||
|
|
lim |
№ = |
lim ( ± |
^ |
' ) = |
± 1 ^ / г = |
± 1 , |
||||||
|
|
дс—*оо |
X |
|
|
ДС—*■ ОО \ |
X |
/ |
|
|
|
|
|
|
|
lim (f(x) — kx) — |
lim (dz^Jx2 — 1 =F x) = |
± |
lim (л/х2 — 1 — дс) = |
||||||||||
|
X—►oo |
|
|
|
|
X—►oo |
_ I _ «2 |
|
X—►oo |
|
|
|||
|
|
|
|
|
|
|
|
0=>b = |
|
|
||||
|
|
|
|
= |
± |
lim |
|
- = |
0 . |
|
||||
|
|
|
|
|
|
|
л]х2 — 1 + * |
|
|
|
|
|
||
Асимптотами |
гиперболы |
являются прямые у = |
± х |
(рис. |
6.18). |
|||||||||
Пример 6.12. Найти |
асимптоты кривой у = л[х. |
|
|
|
|
|||||||||
Р е ш е н и е . |
Д анная |
функция определена и непрерывна на R+. Вертикальных |
||||||||||||
асимптот иет. Д л я |
определения наклонных |
|
|
|
|
|
|
|||||||
асимптот находим пределы (6 . 1 ); |
|
|
|
|
|
|
|
|||||||
* = Ш п Ж = И т |
|
|
|
|
|
|
|
|
|
|||||
|
х —►оо |
X |
х —►оо |
X |
|
|
|
|
|
|
|
|
||
6 = |
lim (f(x) — k x ) = |
lim |
- \Jx= oo. |
|
|
|
|
|
|
|||||
|
X-»*OO |
|
|
|
дс—►oo |
|
|
|
|
|
|
|
||
Второй предел бесконечен, следова- |
|
|
|
|
|
|
||||||||
тельно, кривая у = - 6 |
асимптот не |
имеет |
|
|
|
|
|
|
||||||
(рис. 6.19). |
|
|
|
|
|
|
|
|
|
|
|
Р и с . 6.19 |
6 .6 . ОБЩ АЯ СХЕМА И С СЛЕДОВА НИЯ ФУНКЦИИ
Исследование дважды дифференцируемой функции y = f(x) на D(f) (за исключением, быть может, конечного множества точек) и построение ее графика можно выполнять по приводимой ниже схеме.
1. Находим D(f), определяем точки разрыва, вертикальные асим тоты, нули, точки пересечения графика функции с осью Оу, пе риодичность, симметрию.
151
2. Находим невертикальные асимптоты графика функции (если они существуют).
3. С помощью первой производной функции определяем ста ционарные точки и интервалы монотонности.
4 . С помощью второй производной определяем интервалы вогну тости и выпуклости графика функции, точки перегиба.
5. Находим локальные экстремумы функции на D(f).
По результатам исследований строим график функции, для удоб ства сводя их в таблицу, построение которой покажем на примерах.
Пример 6.13. Исследовать |
функцию i/ = |
X '* |
построить |
ее график. |
|
|||||
--------- и |
|
|||||||||
|
|
|
|
3 — лг |
|
|
|
|
|
|
Р е ш е н и е . Д л я построения |
графика |
функции |
проведем |
ее |
исследование |
по |
||||
указанной выше схеме. |
|
|
|
|
|
|
|
|
|
|
1. Находим D(f). Д анная |
функция определена для х ф |
± -\/3: |
|
|
||||||
D(f) = |
] ~ оо; |
— л /з[и ] — V 3; |
У з [ и ] У з ; |
<*>[. |
|
|
||||
Функции непрерывна иа |
D(f)\ |
* = ± 3 — точки |
разрыва второго |
рода, * = ± 3 |
— |
вертикальные асимптоты графика функции. Функция у = 0 при х = 0, т. е. график функции пересекает координатные оси в начале координат. Функция иепериодичиа.
(— jc)3 —X3
Оиа нечетная, так как D(f) симметрична и f( — x) = — f(x), т. е. |
- ----- ---------------- г . |
|
|
3 — х |
3 — х |
Следовательно, график функции симметричен относительно начала координат и до
статочно исследовать функцию для |
|
х > |
0 . |
|
|
|
|
|
|
|
|
|
|
|
|||||
2. Д л я |
нахождения |
иевертикальиых асимптот |
вычислим |
пределы |
(6.1): |
|
|||||||||||||
|
|
fe= l i m |
Ж |
= |
„п, |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
х—►оо |
X х-*-оо (3 — X |
)Х |
|
|
|
|
|
|
|
|
|
|||||
Ь = |
lim (f(x ) — kx) = |
|
lim |
( — - — г- + |
лД = |
lim |
|
- Х |
. х— . = 0, |
|
|
||||||||
|
х- *■оо |
|
|
X * оо \ |
|
3 — |
X |
|
/ |
х—- оо |
|
3 — |
X |
|
|
|
|||
т. е. прямаи у = |
— х явлиетси наклонной |
асимптотой графика |
функции. |
|
|
||||||||||||||
3. Находим первую производную функции у = f(x): |
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
3х2(3 — я2) + 2х* |
|
х 2(9 — х 2) |
|
|
|
|
|
|
|
||||||
|
|
* |
= |
-----^ |
|
~ |
|
(З- |
, |
2)2 • |
|
|
|
|
|
|
|
||
Она определена иа D(f). В промежутке [0; |
-)- оо [ производная |
обращается |
в нуль |
в |
|||||||||||||||
точках xi = |
0 , Х2 = |
3. |
|
|
|
|
|
|
|
|
|
|
у' |
|
0 н у ' |
|
|
|
|
Определяем интервалы |
монотонности |
йз |
неравенств |
> |
< 0 |
V* > |
0. |
||||||||||||
Имеем |
х2( 9 - х 2) |
|
0=>9 — х 2 > |
0=>л:2 < |
9=>-0 < |
х < |
3, |
|
|
|
|||||||||
|
|
------щ - > |
|
|
|
||||||||||||||
|
(3 — х |
) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
т.е.функция возрастает иа ]0; |
|
|
|
3[. |
|
|
|
|
|
|
|
|
|
|
|
||||
Аналогично |
|
|
.y?\ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
х 2 < |
|
|
|
|
|
|
|
|
|
|||
|
|
" |
(з _ |
^ |
|
< |
0=>-9 |
0=>-л: > 3, |
|
|
|
|
|
|
|||||
т. е. функция убывает на ]3; |
оо [. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
4. Вычислием |
вторую производную функции y = f(x): |
|
|
|
|
|
|
|
|||||||||||
„ _ |
(18дс - |
4*3) (3 - |
х 2)2 - |
(9х2 - |
х4)2(3 - |
х2) ( - 2 х ) |
|
|
6х(9 + х 2) |
|
|
||||||||
У |
|
|
|
|
(3 - |
x 2f |
|
|
|
|
|
|
|
(3 — x 2f |
|
|
’ |
Оиа определена иа D(J). Находим интервалы вогнутости и выпуклости графика функции из неравенств у " > 0, у" < 0 Vx > 0. Тогда
152
6*(9 + х2)(х > |
0, |
(х > 0, |
„/— |
—7т------- м — > |
0 = Ч о |
I — -д/З < дс < |
?=>-]I—I—=> 0< |
(3— л:2)31 3 — дс2 > О |
^/3 |
т. е. кривая вогнута иа ]0; -\/з[. Аналогично
6дс(9 Ч- Д^2) . „ |
Г * > 0 |
Где |
> О, |
(3 — x 2f |
\3 — дс2 < 0= 1 3 < * 2 ~ |
|
|
|
|
|
|
=► / |
^ > |
°’ |
|
ф |
|
|
|
|
|
|
|
^ |
1 х < |
— ~^3, х > |
|||
т. е. кривая |
выпукла иа |
] л £ |
оо [. В |
точке |
|
||||||
дс = 0 |
у " = 0 |
и f " { x ) < 0 |
в |
0 « ( 0 — 0 ), |
|
||||||
a f i x ) > |
0 |
в 0 « ( 0 |
+ 0 ), |
т. е. точка кривой с |
|
||||||
абсциссой |
х — 0 |
отделяет |
интервал |
выпук |
|
||||||
лости кривой от ее интервала вогнутости. |
|
||||||||||
Тогда |
0 (0 ; |
0) |
является |
точкой |
перегиба |
|
|||||
кривой. |
|
|
|
|
|
|
|
|
|
|
|
5. |
|
Определяем |
с |
помощью второй |
про |
||||||
изводной f " ( х) локальные экстремумы. Так |
|
||||||||||
как f ” (3) = |
0 , точка At |
с абсциссой |
х = 3 |
|
|||||||
является |
точкой |
локального максимума. |
|
||||||||
В силу симметрии графика функции точка Л 2 |
|
||||||||||
с абсциссой х = |
— 3 является точкой локаль |
|
|||||||||
ного |
минимума. |
|
Итак, |
m ax fix) — |
—4,5, |
|
|||||
min |
f{x) = |
4,5. |
|
|
|
* 6 0 ,(3 ) |
|
|
|
||
|
|
|
|
|
|
|
|
||||
* 6 0 , ( - 3 ) |
|
|
|
|
|
|
функции у = |
|
|||
Результаты |
исследования |
|
|||||||||
= /(дс) иа |
[0 ; |
оо[ заносим в табл. 6 . 1 . |
|
|
|||||||
Исходя |
из |
результатов, |
содержащихся |
|
|||||||
в табл. 6 . 1 , строим график данной функции |
|
||||||||||
(рис. |
6 .2 0 ). |
|
|
|
|
|
|
|
|
|
|
X |
|
|
0 |
|
]0 ; |
|
|
|
|
|
|
у ' |
|
|
0 |
|
|
|
+ |
|
Не |
существует |
|
У" |
|
|
0 |
|
|
|
+ |
|
Не |
существует |
|
|
|
|
0 |
|
|
|
|
|
|
|
|
У |
|
|
точка |
|
|
|
|
Не |
существует |
||
|
|
перегиба |
|
|
|
|
|
|
|
>x>-\j3.
|
|
Таблица 6.1 |
h & 3 [ |
3 |
]3; °°[ |
+ |
0 |
— |
— |
— |
— |
|
- 4 , 5 |
ч |
|
ш ах |
|
6.7. И Н ТЕ РП О Л И РО В А Н И Е ФУНК ЦИЙ
Под интерполированием функции fix) на отрезке [а; Ь] понимают восстановление этой функции с заданной степенью точности по таблице п ее значений на [а; Ь]:
X |
Хо |
Xi |
Хг |
дся_ i |
Xn |
№ |
f i x) о |
f i x l) |
f i x 2) |
f i x * - 1) |
f i x . ) |
153
Точки Хо, х и |
..., |
хп (х0 = а, |
хп = Ь) называют узлами |
интерполи |
||
рования. |
|
|
|
|
|
|
Функцию f(x) приближенно |
заменяют многочленомРп(х) = а 0дсп + |
|||||
+ a t.x"~l + |
... + |
ап, |
значения |
которого |
совпадают со |
значениями |
f(x) в узлах |
интерполирования: |
|
|
|
||
|
|
|
P n (Xi) = |
f( X i) , i = |
О, п. |
(6.2) |
Многочлен Рп(х) называют интерполяционным многочленом степени
п |
функции f(x) на [а; Ь]. Д ля определения коэффициентов многочлена |
||||||||||||||
Р п(х) |
используем |
условия |
(6.2): |
|
|
|
|
|
|
|
|||||
|
|
|
|
ОоХо + |
а ,4 |
1 — ... — ап = fix о),' |
|
||||||||
|
|
|
|
а о х " aix" |
1 |
|
|
ап = |
f(xi), |
|
(6-3) |
||||
|
|
|
|
а0х" + |
а\хп |
1+ |
... + |
a„ = |
f( x n ) . , |
|
|||||
Получим |
систему |
п + |
1 |
уравнений |
с |
п + |
1 |
неизвестными. |
|||||||
|
Определитель системы |
(6.3) |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
ло |
„п — 1 |
.. |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
Хо |
|
|
|
|
|
|
|||
|
|
|
|
|
|
v n |
„ п — 1 |
.. |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
Х\ |
ДС1 |
|
Ф о, |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
X я |
vrt— 1 |
.. |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
если |
Хо, |
х\, |
х п |
различны. |
М ожно показать, что |
он равен |
|||||||||
|
П |
|
(xi — ^ . |
Действительно, |
определитель |
второго |
порядка |
||||||||
|
|
|
|
|
|
I * 0 |
1 I |
= х 0 — х й |
|
|
|
|
|||
|
|
|
|
|
|
I Xl |
1 I |
|
|
|
|
|
|
|
|
определитель третьего |
порядка |
|
|
|
|
|
|
|
|||||||
|
|
|
ДСо Хо |
1 |
|
|
ДСо |
|
ДСо |
|
1 |
= |
|
||
|
|
|
Х21 Xl |
1 |
= |
x i — Хо |
Х{ — Хо |
0 |
|
||||||
|
|
|
дс§ х 2 |
1 |
|
х \ — Хо |
ДС2 — |
ДС0 |
. 0 |
|
|
||||
= |
(x2t — X o ) { X 2 — Xo) — ( X 2 — |
X o ) ( X i — |
Xo) |
= |
( * 1 |
|
|
||||||||
и |
т. д. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Система (6.3) |
имеет |
единственное |
решение. Следовательно, су |
ществует единственный интерполяционный многочлен Р п{х), коэф фициенты которого определяются из системы (6.3). Однако опре деление коэффициентов интерполяционного многочлена из системы
(6.3) связано с громоздкими вычислениями. |
|
|
||
Будем искать интерполяционный |
многочлен Р п{х) в |
виде |
||
Рп(х) = а0(х — х ,) (* — х 2) • • • (дс — х п) + |
а , (* — х 0) {х — х 2) ■■■(х — хп) + |
|||
+ а2(х — хо) (дс — Jct) (х — х 3)---(х — х п) + ... + |
ап{х — х0) — (х — дс„_ 1). |
|||
П олагая в |
последнем равенстве |
x = xi, |
i — 0, п, |
и учитывая |
условия (6.2), |
находим коэффициенты Рп(х): |
|
: ________[ М ________
а о
(х0 — Xi) (х0 - х 2)---(х0 — хп) ’
154
а , = __________& й__________
( *1 — Хо) (Xl — Х2) ■■■(X! — Х„) '
а = _______________ м _______________
‘(х,- — Хо) • • ■(X/ — Xi_ 1) (X/ — Xi+ 1) •• • (Xi — Хп)
Интерполяционный многочлен примет вид
|
■ |
(х — х0)(х — х2)--- (х — Хп) f ( |
ч | |
, |
|
|
(Xl — Хо)(х, — Х2) ••• (х, — Хп) |
Т |
■" ' |
|
|
(х — х0)(х — X,) ••• (х — X„-i) |
Г/ |
ч |
|
|
(х„ — х0)(х„ — X!) ••• (х„ — x„_i) |
|
|
Тогда |
|
|
|
|
f(X) * |
Р„(*) = У |
f(Xi) |
|
. (6 .4) |
v |
L J |
(Xi-Xo)---(Xi— X/-l)(Xi — X,■+,)••.(Xi — X„) |
||
|
1= 0 |
|
|
|
Равенство (6.4) |
называется интерполяционной формулой Л а г |
ранжа.
Простейшим и широко используемым при работе с функциями, заданными таблично, является линейное интерполирование, т. е. интерполирование многочленом первой степени Pi(x) = а {х + а0. Линейное интерполирование удобно использовать, если известно, что на отрезке [а; Ь) данная функция мало отличается от линейной.
Пример 6.14. Функция y = f(x) задана таблично: |
|
|
||
X |
1 |
2,5 |
3 |
4 |
/ м |
2 |
1 |
2 , 2 |
3 |
Построить для этой функции интерполяционный многочлен Л агр ан ж а третьей степени
на отрезке |
[1; 4]. |
|
|
|
|
|
|
|
|
|
|
|
|
Р е ш е н и е . |
Воспользуемся |
формулой |
(6.4): |
|
|
|
|
|
|
||||
р,(х) = |
f (1) |
( * - 2 . 5 ) ( х - 3 ) ( х - 4 ) |
|
|
|
{х |
1) (х |
3) (х |
4) |
|
|||
|
Д ; |
( 1 - 2 , 5 ) ( 1 - 3 ) ( 1 - 4 ) |
|
|
' ’ ( 2 , 5 - 1) ( 2 , 5 - 3 ) ( 2 , 5 - 4 ) + |
||||||||
+ 1W |
(х |
1) (х |
2,5) (х |
4) |
|
(х |
1) |
(х |
2,5) |
(х |
3) |
р ( |
) _ |
(з _ |
1) ( 3 _ |
2,5) (3 — 4) |
) |
(4 - |
1) |
(4 - |
2,5) |
(4 - |
3) |
з( |
) |
= — -|-(х — 2,5)(х — 3) (х — 4) + -|-(х — I) (х — 3)(х — 4) —
— 2,2(х — 1) (х — 2,5) |
(х — 4) + |
у |
( х - 1) (х - |
2,5) (х - 3)=>-Рз(х) = |
|
13 |
|
. . . I . |
3 |
„ |
|
= ~ \ ь ^ |
+ 7 Т |
х |
- » 7 ТО* + 13- |
После построения интерполяционного многочлена Л агранж а для f(x) на [а; Ь] возникает вопрос об оценке точности приближения функции f(x) многочленом Р п(х).
Разность между функцией f(x) и ее интерполяционным много членом Р п(х), т. е.
155
Гп (х ) = f(x ) — Р П(Х),
называется остаточным членом интерполирования.
|
Если известно, что на [а; Ь) функция f(x) п + |
1 раз дифференци |
|||
руема, то |
|
|
|
|
|
|
Ы * ) \ |
< - $ $ $ ■ v * e [ a ; ь \ |
|
(6.5) |
|
где й = ш а х |/("+ |)(д:)|; |
Ц х ) = (х — х0)(.х — Х |)'"(* — •*«); |
^ — |
|||
|
______ *€ [о; 1>| |
|
|
|
|
= |
0 , п,— узлы интерполирования. |
|
|
|
|
|
Оценка остаточного члена интерполирования |
(6.5) |
требует опре |
||
деленной информации о поведении рассматриваемой |
функции |
на |
|||
[а; |
Ь], которая, как правило, отсутствует. |
|
|
|
Иногда, выбрав достаточно большое число узлов интерполиро вания, можно получить хорошее приближение интерполируемой функции. Известны примеры, когда последовательность интерполя ционных многочленов Л агр ан ж а функции f(x) с возрастанием числа узлов интерполирования не стремится к этой функции.
На практике для получения достаточно хорошего приближения функции вместо построения интерполяционного многочлена высокой степени используют кусочную интерполяцию многочленами более низких степеней, т. е. на каждом отрезке x<+i] строится свой многочлен. Часто используется кусочно-линейная интерполяция.
Получающиеся при этом кусочно-многочленные функции с однород ной структурой на каждом промежутке (многочлен одной и той же
степени) называются сплайн-функциями или просто |
сплайнами (от |
англ. spline — рей ка). Простейший пример сплайна |
— ломаная. |
В настоящее время в приложениях чаще всего используются приближения кубическими сплайнами, так как последовательность интерполяционных кубических сплайнов на равномерной сетке узлов всегда сходится к интерполируемой непрерывной функции, а если дифференциальные свойства функции улучшаются, то и сходимость улучшается. Кроме того, алгоритмы построения сплайнов достаточно
просто реализуются на |
ЭВМ. |
Пример 6.15. Построить |
линейный сплайн Si(x) н иайти S,(0,5) для функции |
у = f(x), заданной таблично: |
|
X |
0 |
0 , 2 |
0,4 |
0 , 6 |
0 , 8 |
1 |
1 ,2 |
m |
1 |
1,5 |
1,7 |
1,3 |
1,7 |
2 |
2 ,1 |
1,4 |
1 ,6 |
1 ,8 |
2 |
2,5 3 3,2 3,5
Р е ш е н и е . Линейный сплайн Si(x) задается уравнением прямой, проходящей через две точки:
У — Hi |
— — — |
_ |
о / ч |
X — X l , |
—Vi)- |
—------- |
—=>Si(x) = y = yt+ |
—------- —(«//+1 |
|||
i/ j 4 1 |
У i |
|
|
1 X i X t + i — Xi |
|
Вводя понятие ш ага разбиения A = |
xi+ i — Xi, имеем |
|
|
156
|
|
S, (x) = |
f(Xl) + ( x - X ,) |
K |
* + l)ft |
•^(* ) . |
|
|
Так |
как |
0,5 e [0,4; 0,61 |
( « = 0,4, |
xi + 1 |
= |
0,6), |
T O |
S,(0,5) = Д0,4) + (0,5 — |
— °-4) l , 3 0~ |
l : 7 = 1.5 (рис. 6 .2 1 ). |
|
|
|
|
|
||
|
|
6 .8 . П Р И Б Л И Ж Е Н Н О Е |
РЕ Ш Е Н И Е УРАВНЕНИЙ |
|||||
Постановка задачи. Рассмотрим задачу нахождения действи |
||||||||
тельных корней уравнения |
|
|
|
|
|
|||
|
|
|
fix) = 0, |
|
|
(6.6) |
||
где fix) |
— алгебраическая или трансцендентная |
функция. |
Известные точные методы решения уравнений пригодны только для узкого класса уравнений (квадратных, биквадратных, дву членных алгебраических, некоторых тригонометрических, показа тельных, логарифмических). Однако в общем случае решения урав нения (6 .6 ) находятся приближенно.
Нахождение приближенных значений действительных корней уравнения (6 .6 ) состоит из двух этапов: 1 ) отделение корней (опре деление малых отрезков [а; b], в которых находится только один действительный корень уравнения); 2 ) приближенное вычисление (уточнение) отделенного корня до заданной точности.
Отделение действительных корней. Под отделением действитель |
||
ного корня |
уравнения fix) = 0 |
понимают нахождение отрезка [а; b], |
в котором |
лежит только один |
корень данного уравнения. Этот от |
резок называют промежутком изоляции корня.
Известны различные графические и аналитические методы отде ления корней уравнения (6 .6 ).
Наиболее наглядный метод отделения корней уравнения fix) = = 0 состоит в определении координат точек пересечения графика функции у = fix) с осью абсцисс. Абсциссы x f точек пересечения
157