Информатика 2 семестр / obrazets_vypolnenia_laboratornykh_rabot_2_semestr
.pdfРеализация в программе MS Excel Построение интерполирующей функции с помощью тренда
Квадратичная интерполяция
1.Ввести исходные данные – значения Х разместить в 1- ой строке значения Y во второй строке.
2.Построить диаграмму по всем точкам (тип диаграммы – «Точечная» точки без соединения линиями).
3.Для построения полинома 2-й степени P2(x)= a0 + a1 x + a2 x2 выберем три точки (x0, y0), (x1, y1) и (x4, y4).
4.Добавим ряд из этих точек.
Построим для этого ряда полиномиальный тренд 2-й степени и показать уравнение на диаграмме.
51
Кусочно-линейная интерполяция (Метод неопределѐнных коэффициентов)
Постановка задачи: Дана таблично заданная функция {xi,yi}
i |
0 |
1 |
2 |
3 |
4 |
x |
0,2 |
0,4 |
0,7 |
0,85 |
1 |
y |
0,1 |
0,5 |
0,6 |
0,9 |
0,7 |
Выполнить кусочно – линейную интерполяцию.
Система для нахождения |
|
|
Ответ |
|
|||||
коэффициентов полинома на |
|
|
|
|
|
||||
каждом участке |
|
|
|
|
|
|
|||
1 участок |
2участок |
|
|
|
|
|
|||
a10 a11 x0 |
y0 |
a20 |
a21 x1 |
y1 |
a1 |
a1 |
x,приx |
x x |
|
|
y1 |
|
a21 x2 |
y2 |
|
0 |
1 |
0 |
1 |
a10 a11 x1 |
a20 |
a2 |
a2 |
x,приx |
x x |
||||
|
|
|
|
|
|
0 |
1 |
1 |
2 |
|
|
|
|
|
P1( x ) a3 |
a3 |
x,приx |
x x |
|
3 участок |
|
4участок |
|
|
0 |
1 |
2 |
3 |
|
|
|
a40 |
a41 x,приx3 x x4 |
||||||
a30 a31 x2 y2 |
a40 a41 x3 y3 |
||||||||
|
y3 |
|
a41 x4 |
y4 |
|
|
|
|
|
a30 a31 x3 |
a40 |
|
|
|
|
|
На интервале от 0,2 до 1 задано 5 точек, получаем 4 отрезка. Интерполируем (метод Неопределенных коэффициентов) полиномом 1 степени каждый отрезок. После нахождения каждого полинома запишем результат.
1 участок.
i |
0 |
1 |
Общий вид полинома P11( x ) a10 |
a11 x . По условию интерполяции полином |
|||||||
x |
0,2 |
0,4 |
|||||||||
должен проходить через точки, которые выбраны для построения, т.е. |
|||||||||||
y |
0,1 |
0,5 |
|||||||||
P11( x0 ) y0 |
a10 |
a11 x0 |
y0 |
|
|||||||
|
|
|
. Подставим значения |
||||||||
|
|
|
P11( x |
) y |
Следовательно a1 |
a1 x |
y |
||||
|
|
|
|
1 |
1 |
|
0 |
1 1 |
1 |
|
|
|
a10 a11 0,2 0,1 |
.Неизвестными в системе являются |
|||||
x0 ,x1, y0 , y1. В результате получаем |
a1 a1 0,4 0,5 |
|||||||
|
|
|
|
|
|
|||
|
|
0 |
1 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
0,2 |
0,1 |
|
|
|
|
a10 ,a11 .Решим систему методом Гаусса. |
|
. Приведем систему к треугольному виду, |
||||||
|
|
|
1 |
0,4 |
0,5 |
|
|
|
для этого 1 строку перепишем, из 2 строки вычтем 1 строку и результат запишем на место 2 |
|
|||||||
1 |
|
|
|
|
a10 a11 0,2 0,1 |
|
||
0,2 |
0,1 |
|
|
|
. Из |
|||
строки. 0 |
0,2 |
0,4 . Запишем полученную матрицу в виде системы. |
a1 0,2 0,4 |
|||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
2 уравнения найдем a1 0,4 |
a1 2 . Из 1 уравнения найдем a1 . |
a1 |
0,1 0,2 a1 |
|||||||||
|
|
|
|
1 |
0,2 |
1 |
|
|
0 |
0 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
||
a10 0,1 0,2 2 |
a10 0,3. Запишем найденное уравнение P11( x ) 0,3 2 x . Проверка. |
|||||||||||
Найденное уравнение должно проходить через точки x0 , y0 , x1, y1 . |
|
|||||||||||
P11( x0 ) 0,3 2 x0 |
|
|
|
|
|
|
|
|||||
P11( 0,2 ) 0,3 2 0,2 |
|
|
|
|
|
|
|
|||||
P11( 0,2 ) 0,1 |
|
|
|
|
|
|
|
|
||||
P11( x1 ) 0,3 2 x1 |
|
|
|
|
|
|
|
|||||
P11( 0,4 ) 0,3 2 0,4 |
|
|
|
|
|
|
|
|||||
P11( 0,4 ) 0,5 Следовательно прямая проходит через 0 и 1 точки. |
|
|
||||||||||
2 участок. |
|
|
|
|
|
|
|
|
|
|||
i |
1 |
2 |
Общий вид полинома P12( x ) a20 |
a21 x . По условию интерполяции полином |
||||||||
x |
0,4 |
0,7 |
||||||||||
должен проходить через точки, которые выбраны для построения, т.е. |
||||||||||||
y |
0,5 |
0,6 |
||||||||||
P12( x ) y |
a2 |
0 |
a2 x y |
|
|
|||||||
|
|
|
|
|
||||||||
|
|
|
|
1 |
1 |
Следовательно |
1 1 |
1 . Подставим значения |
||||
|
|
|
P12( x2 ) y2 |
a20 |
a21 x2 |
y2 |
|
|
x1,x2 , y1, y2 |
|
|
|
|
a20 a21 0,4 0,5 |
.Неизвестными в системе являются |
|||||||||
. В результате получаем |
|
|
|
||||||||||||
|
|
|
|
|
a20 a21 0,7 0,6 |
|
|
|
|
|
|
||||
|
|
|
|
|
1 |
0,4 |
|
0,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
a20 ,a21 .Решим систему методом Гаусса. |
|
|
. Приведем систему к треугольному виду, |
||||||||||||
|
|
|
|
|
1 |
0,7 |
|
0,6 |
|
|
|
|
|
|
|
для этого 1 строку перепишем, из 2 строки вычтем 1 строку, и результат запишем на место 2 |
|
||||||||||||||
1 |
0,4 |
|
0,5 |
|
|
|
|
|
a20 a21 0,4 0,5 |
. Из |
|||||
|
|
|
|
|
|
||||||||||
строки. 0 |
0,3 |
|
0,1 . Запишем полученную матрицу в виде системы. |
a2 |
0,3 0,1 |
||||||||||
|
|
|
|
0,1 |
|
|
|
|
|
|
1 |
|
|
||
|
|
|
|
|
|
|
|
|
|||||||
2 уравнения найдем a2 |
a2 0,333. Из 1 уравнения найдем a2 |
0 |
. a2 |
0 |
0,5 0,4 a2 |
||||||||||
|
|||||||||||||||
|
1 |
0,3 |
1 |
|
|
|
|
|
|
1 |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a20 0,5 0,4 0,333 a20 0,3668. Запишем найденное уравнение
P12( x ) 0,3668 0,333 x . Проверка. Найденное уравнение должно проходить через точки
x1, y1 , x2 , y2 .
P12( x1 ) 0,3668 0,333 x1 P12( 0,4 ) 0,3668 0,333 0,4
P12( 0,4 ) 0,5
P12( x2 ) 0,3668 0,333 x2 P12( 0,7 ) 0,3668 0,333 0,7
P12( 0,7 ) 0,6 Следовательно прямая проходит через 1-ю и 2-ю точки.
3 участок. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
i |
2 |
3 |
Общий вид полинома P13( x ) a30 |
a31 x . По условию интерполяции |
||||||||||||
x |
0,7 |
0,85 |
||||||||||||||
полином должен проходить через точки, которые выбраны для построения, т.е. |
||||||||||||||||
y |
0,6 |
0,9 |
||||||||||||||
P13( x2 ) y2 |
|
|
a30 |
a31 x2 |
y2 |
|
||||||||||
|
|
|
|
|
. Подставим значения |
|||||||||||
|
|
|
P13( x |
) y |
Следовательно a3 |
|
a3 x |
y |
||||||||
|
|
|
|
3 |
3 |
|
|
|
0 |
|
1 3 |
3 |
|
|||
|
|
|
|
|
|
|
|
a30 a31 0,7 0,6 |
.Неизвестными в системе |
|||||||
x2 ,x3 , y2 , y3 . В результате получаем |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
a30 a31 0,85 0,9 |
|
|
||||||
|
|
|
|
|
|
|
|
1 |
0,7 |
|
0,6 |
|
||||
|
|
|
|
|
|
|
|
|||||||||
являются a30 ,a31 .Решим систему методом Гаусса. |
|
|
|
|
. Приведем систему к |
|||||||||||
|
|
|
|
|
|
|
|
1 |
0,85 |
|
0,9 |
|
||||
треугольному виду, для этого 1 строку перепишем, из 2 строки вычтем 1 строку, и результат |
||||||||||||||||
|
|
|
|
|
1 |
0,7 |
|
0,6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
запишем на место 2 строки. |
|
|
. Запишем полученную матрицу в виде системы. |
|||||||||||||
|
|
|
|
|
0 |
0,15 |
|
0,3 |
|
|
|
|
|
|
|
a30 a31 0,7 0,6 |
. Из 2 уравнения найдем a31 |
0,3 |
|
2 . Из 1 уравнения найдем a30 . |
||
|
a31 0,15 0,3 |
|
a31 |
|||
0,15 |
||||||
|
|
|
|
a30 0,6 0,7 a31 a30 0,6 0,7 2 a30 0,8 . Запишем найденное уравнение P13( x ) 0,8 2 x . Проверка. Найденное уравнение должно проходить через точки
x2 , y2 , x3 , y3 .
P13( x2 ) 0,8 2 x2
P13( 0,7 ) 0,8 2 0,7 P13( 0,7 ) 0,6
P13( x3 ) 0,8 2 x3
P13( 0,85) 0,8 2 0,85
P13( 0,85) 0,9 .Следовательно, прямая проходит через 2-ю и 3-ю точки. 4 участок.
i 3 4
53
x |
0,85 |
1 |
Общий вид полинома |
P14( x ) a40 a41 x . По условию интерполяции |
|
||||||||||
y |
0,9 |
0,7 |
|
||||||||||||
полином должен проходить через точки, которые выбраны для построения, т.е. |
|||||||||||||||
|
|
|
|||||||||||||
P14( x ) y |
a4 |
0 |
a4 x |
3 |
y |
3 . Подставим значения x3 ,x4 , y3 , y4 |
. В |
||||||||
|
|
3 |
3 |
Следовательно |
|
1 |
|
|
|||||||
P14( x4 ) y4 |
a40 |
a41 x4 |
y4 |
|
|||||||||||
|
|
|
|
a40 a41 0,85 0,9 |
.Неизвестными в системе являются a40 ,a41 |
|
|||||||||
результате получаем |
1 0,7 |
|
|||||||||||||
|
|
|
|
a40 a41 |
|
|
|
|
|
||||||
|
|
|
|
1 |
0,85 |
|
0,9 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|||||||
.Решим систему методом Гаусса. |
|
|
|
|
|
. Приведем систему к треугольному виду, для |
|||||||||
|
|
|
|
1 |
|
|
1 |
|
0,7 |
|
|
|
|
|
этого 1 строку перепишем, из 2 строки вычтем 1 строку, и результат запишем на место 2 строки.
1 |
0,85 |
|
0,9 |
|
a40 a41 0,85 0,9 |
. Из 2 |
|||
|
|
||||||||
0 |
0,15 |
|
0,2 . Запишем полученную матрицу в виде системы. |
a4 |
0,15 0,2 |
||||
|
|
|
|
0,2 |
|
|
1 |
|
|
|
|
|
|
|
|||||
уравнения найдем a4 |
a4 1,333. Из 1 уравнения найдем a4 |
0 |
. |
|
|||||
|
1 |
0,15 |
1 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
a40 0,9 0,85 a41 |
a40 0,9 0,85 ( 1,333) a40 2,033. Запишем найденное |
|
уравнение P14( x ) 2,033 1,333 x . Проверка. Найденное уравнение должно проходить через
точки x3 , y3 , x4 , y4 .
P14( x3 ) 2,033 1,333 x3 P14( 0,85) 2,033 1,333 0,85
P14( 0,85) 0,9
P14( x4 ) 2,033 1,333 x4 P14(1) 2,033 1,333 1
P14(1) 0,7 .Следовательно, прямая проходит через 3-ю и 4-ю точки.
|
0,3 2 x,если0,2 x 0,4 |
|
|
|
|
Запишем ответ P1( x ) 0,3668 0,333 x,если0,4 x 0,7 |
||
|
0,8 2 x,если0,7 x 0,85 |
|
2,033 1,333 x,если0,85 x 1 |
||
|
Построим график |
. |
Реализация метода в Mcad
Метод неопределѐнных коэффициентов(кусочно-линейная интерполяция)
|
|
0. |
2 |
|
|
|
|
0.1 |
|||
|
|
0. |
4 |
|
|
|
|
|
0.5 |
|
|
|
|
|
|
|
|
|
|
||||
x |
0. |
7 |
|
y |
0.6 |
||||||
|
|
0.85 |
|
|
|
|
0.9 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
0.7 |
|||
1 участок |
|
|
|
|
|
||||||
|
|
1 |
x |
|
|
|
y |
0 |
|
|
|
C1 |
|
|
|
0 |
|
D1 |
|
|
|
||
|
|
|
|
|
|
||||||
|
1 |
x |
|
|
y |
|
|
|
|||
|
|
|
|
|
|||||||
|
|
|
|
1 |
|
|
1 |
|
|
a1 C1 1 D1
0.3
a1
2
P11(x) a10 a11 x
2 участок
|
1 |
x |
|
|
y |
1 |
|
|
C2 |
|
|
1 |
|
D2 |
|
|
|
|
|
|
||||||
|
1 |
x |
|
y |
|
|
||
|
|
|
||||||
|
|
|
2 |
|
|
|
2 |
a2 C2 1 D2
a2
0.3670.333
P12(x) a20 a21 x
Исходные данные
i 0 4
Матрицы C1 и D1 для системы линейных уравнений
Решение системы линейных уравнений нахождение коэффициентов а10 и а11
Интерполирующая функция 1 участка
Матрицы C2 и D2 для системы линейных уравнений
Решение системы линейных уравнений нахождение коэффициентов а20 и а21
Интерполирующая функция 2 участка
55
|
3 участок |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
1 |
x |
|
|
|
|
|
y |
|
|
|
|
Матрицы C3 и D3 для системы линейных |
|
|
|
|
|
|
|
2 |
|
|
уравнений |
|
|||||
|
C3 |
|
2 |
|
|
|
D3 |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
||||||
|
1 |
x |
|
|
|
|
y |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
3 |
|
|
|
|
3 |
|
|
|
|
|
||
|
a3 C3 1 D3 |
|
|
|
|
|
|
|
|
Решение системы линейных уравнений |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
нахождение коэффициентов а30 и а31 |
|
|
|
0.8 |
|
|
|
|
|
|
|
|
|
|
|||
|
a3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
P13(x) a30 a31 x |
|
|
|
|
|
Интерполирующая функция 3 участка |
||||||||
|
|
|
|
|
|
|
|
||||||||
|
4 участок |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
1 |
x |
|
|
|
|
y |
|
|
Матрицы C4 и D4 для системы линейных |
||
|
|
|
|
|
|
3 |
уравнений |
|
|||||||
|
C4 |
|
|
|
3 |
|
D4 |
|
|
|
|||||
|
|
|
|
|
|
|
|||||||||
|
|
1 |
x |
|
y |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
4 |
|
|
|
|
|
4 |
|
|
|
|
a4 C4 1 D4 |
|
|
|
|
|
|
Решение системы линейных уравнений |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
нахождение коэффициентов а40 и а41 |
|
|
|
|
2.033 |
|
|
|
|
|
|
|
|
||||
|
|
a4 |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
1.333 |
|
|
|
|
|
|
|
|||||
|
|
P14(x) |
a40 a41 x |
|
|
|
|
Интерполирующая функция 4 участка |
|||||||
|
|
|
|
|
|
|
|
||||||||
|
|
P1(t) if t x1 P11(t) if t x2 P12(t) if t x3 P13(t) P14(t) |
Ответ |
||||||||||||
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
P1 xi |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
y i |
|
0.5 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
0 0 |
|
|
|
|
0.5 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
xi |
|
|
Кусочно-параболическая интерполяция (Метод неопределѐнных |
|
||||||||||||||
коэффициентов) |
|
|
|
|
|
|
|
|
|
|
|||||
Постановка задачи: Дана таблично заданная функция {xi,yi} |
|
||||||||||||||
i |
0 |
|
|
|
1 |
|
|
2 |
|
|
|
3 |
4 |
|
|
x |
0,2 |
|
|
|
0,4 |
|
0,7 |
|
|
0,85 |
1 |
|
|||
y |
0,1 |
|
|
|
0,5 |
|
0,6 |
|
|
0,9 |
0,7 |
|
|||
|
Выполнить кусочно – параболическую интерполяцию. |
|
Система для нахождения |
|
|
Ответ |
|
|
|
|||||||||||
коэффициентов полинома |
|
|
|
|
|
|
|
||||||||||
|
на каждом участке |
|
|
|
|
|
|
|
|
||||||||
1 участок |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
a1 a1 x |
0 |
a1 |
2 |
x |
2 |
y |
0 |
|
|
|
|
|
|
|
|||
|
0 |
1 |
|
|
0 |
|
|
|
|
|
|
|
|
||||
|
|
a11 x1 |
|
|
|
2 |
y1 |
|
|
|
|
|
|
|
|||
a10 |
a12 x1 |
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
a10 a11 x |
2 |
a12 |
x |
2 |
y2 |
a10 |
a11 |
x a12 |
x |
, при x0 |
x x2 |
||||||
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
P2( x) |
a21 |
x a22 |
x 2 , при x2 |
x x4 |
|||
|
|
|
|
|
|
|
|
|
|
a20 |
|||||||
2 участок |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a2 |
0 |
a2 |
x |
2 |
a2 |
2 |
x2 |
y |
2 |
|
|
На |
|
1 |
|
|
2 |
|
|
|
|||||
a20 a21 x3 |
a22 |
x32 |
y3 |
|
|
инт |
||||||
a2 |
0 |
a2 |
x |
4 |
a2 |
2 |
x2 |
y |
4 |
|
|
ерв |
|
1 |
|
|
4 |
|
|
|
але |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
от |
0,2 до 1 задано 5 точек, разобьѐм его на 2 отрезка x0 x x2 |
и x2 x x4 . Интерполируем |
(метод Неопределенных коэффициентов) полиномом 1 степени каждый отрезок. После нахождения каждого полинома запишем результат.
1 участок.
i |
0 |
1 |
2 |
x |
0,2 |
0,4 |
0,7 |
y |
0,1 |
0,5 |
0,6 |
Общий вид полинома P21( x ) a10 a11 x a2 x2 . По условию интерполяции полином должен проходить через точки, которые выбраны для
P21( x ) y
0 0
построения, т.е. P21( x1 ) y1 Следовательно
P21( x ) y
2 2
a1 a1 x a1 x |
0 |
2 y |
0 |
|
|
|||||
|
0 |
1 |
|
0 |
2 |
|
|
|
||
|
|
|
|
a12 x12 y1 . Подставим значения x0 ,x1,x2 , y0 , y1, y2 . В результате |
||||||
a10 |
a11 |
x1 |
||||||||
a1 a1 x |
2 |
a1 x |
2 |
2 y |
2 |
|
|
|||
|
0 |
1 |
|
2 |
|
|
|
|||
|
|
a1 a1 0,2 a1 0,22 |
0,1 |
|
||||||
|
|
|
0 |
1 |
|
2 |
|
|
||
|
|
|
|
a11 0,4 |
a12 0,42 |
0,5 |
.Неизвестными в системе являются a10 ,a11,a12 |
|||
получаем a10 |
||||||||||
|
|
a1 a1 0,7 a1 0,72 |
0,6 |
|
||||||
|
|
|
0 |
1 |
|
2 |
|
|
|
|
|
0,2 |
2 |
0,1 |
|
|
|
|
1 0,2 |
2 |
|
|
||||
|
|
|
|
|
|
|
||
.Решим систему методом Гаусса. |
0,4 |
0,42 |
0,5 |
. Приведем систему к треугольному виду, |
||||
1 |
|
|||||||
|
|
|
|
|
0,6 |
|
|
|
|
1 |
0,7 |
0,7 |
|
|
|
|
|
|
|
|
|
для этого 1-ю строку перепишем, из 2-ой строки вычтем 1-ю строку, из 3-ей строки вычтем 1-ю
строку и результат запишем на место 3-ей строки.
2-ю строку поделим на 0,2, а 3-ю строку на 0,5.
1 |
0,2 |
0,04 |
|
0,1 |
|
||
|
|
||||||
|
0 |
0,2 |
0,12 |
|
0,4 |
|
. Перепишем 1-ю строку, а |
|
0 |
0,5 |
0,45 |
|
0,5 |
|
|
|
|
|
|
|
1 |
0,2 |
0,04 |
|
0,1 |
|
||
|
|
|
||||||
Система примет вид |
|
0 |
1 |
0,6 |
|
2 |
|
. Преобразуем систему перепишем 1-ю и 2-ю строки, а |
|
|
0 |
1 |
0,9 |
|
1 |
|
|
|
|
|
|
|
из 3-ей строки вычтем 2-ю строку, результат запишем на место 3-ей строки. Получим следующую
|
|
1 |
0,2 0,04 |
|
0,1 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|||||
систему |
|
0 |
1 |
0,6 |
|
2 |
. Запишем полученную матрицу в виде системы. |
|||||
|
|
|
0 |
0 |
0,3 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|||
a1 a1 0,2 a1 0,04 0,1 |
|
|
|
|
||||||||
|
0 |
|
1 |
|
2 |
|
|
|
|
1 |
|
a12 3,333. Из 2-го |
|
|
a11 |
a12 |
0,6 2 |
. Из 3 уравнения найдем a12 |
|
||||||
|
|
|||||||||||
|
|
a12 |
0,3 1 |
|
0,3 |
|
||||||
|
|
|
|
|
|
|
уравнения найдем a11. a11 2 0,6 a12 a21 2 0,6 3,333 a11 4 . Из 1-го уравнения
получим a10 a10 0,1 0,2 a11 0,04 a12 a10 0,1 0,2 4 0,04 ( 3,333)
a10 0,567.Запишем найденное уравнение P21( x ) 0,567 4 x 3,333 x2 . Проверка. Найденное уравнение должно проходить через точки x0 , y0 , x1, y1 , x2 , y2 .
57
P21( x0 ) 0,567 4 x0 3,333 x02 P21( 0,2 ) 0,567 4 0,2 3,333 0,22
P21( 0,2 ) 0,1
P21( x |
) 0,567 4 x |
|
3,333 x 2 |
P21( 0,4 ) 0,567 4 0,4 3,333 0,42 |
|||||||||||||||||||||||||||
|
|
|
1 |
|
|
|
|
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
P21( 0,4 ) 0,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
P21( x |
2 |
) 0,567 4 x |
2 |
3,333 x |
2 |
2 |
P21( 0,7 ) 0,567 4 0,7 3,333 0,72 |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
P21( 0,7 ) 0,6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Следовательно, полученный полином проходит через 0- ю, 1-ю и 2-ю точки. |
|||||||||||||||||||||||||||||||
2 участок. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
i |
2 |
|
|
|
3 |
|
4 |
|
|
Общий вид полинома P22( x ) a2 |
0 |
a2 |
x a2 |
2 |
x 2 . По условию |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
x |
0,7 |
|
|
0,85 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|||||
|
|
|
|
интерполяции полином должен проходить через точки, которые выбраны для |
|||||||||||||||||||||||||||
y |
0,6 |
|
|
0,9 |
0,7 |
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
P22(x ) y |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
построения, т.е. |
P22(x2 ) y2 |
Следовательно |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
3 |
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
P22(x4 ) y4 |
|
|
|
|
|
|
|
|
||||
a2 |
0 |
a2 x |
2 |
a2 |
2 |
x |
|
2 |
y |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
1 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
a22 x32 |
|
y3 . Подставим значения x2 ,x3 ,x4 , y2 , y3 , y4 . В результате |
|||||||||||||||||||||
a20 |
a21 x3 |
|
|||||||||||||||||||||||||||||
a2 |
0 |
a2 x |
4 |
a2 |
2 |
x |
|
2 |
y |
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
1 |
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
a2 |
0 |
a2 0,7 a2 |
2 |
0,72 0,6 |
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
a21 0,85 a22 |
0,852 |
0,9 .Неизвестными в системе являются |
||||||||||||||||||||||
получаем a20 |
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
0,7 |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
a20 a21 1 a22 1 |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,7 |
2 |
0,6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 0,7 |
|
|
|
|
||||
a20 ,a21,a22 .Решим систему методом Гаусса. |
1 |
0,85 |
0,852 |
0,9 |
. Приведем систему к |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
0,7 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
1 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
треугольному виду, для этого 1-ю строку перепишем, из 2-ой строки вычтем 1-ю строку, из 3-ей
|
1 |
0,7 |
0,49 |
|
0,6 |
||
|
|
||||||
строки вычтем 1-ю строку и результат запишем на место 3-ей строки. |
|
0 |
0,15 |
0,2325 |
|
0,3 |
. |
|
|
0 |
0,3 |
0,51 |
|
0,1 |
|
|
|
|
|
Перепишем 1-ю строку, а 2-ю строку поделим на 0,15, а 3-ю строку на 0,3.Система примет вид
1 |
0,7 0,49 |
|
|
|
0,6 |
|
|
|
|
|||||
|
|
|
|
|
|
|||||||||
|
0 |
|
1 |
1,55 |
|
|
|
2 |
|
. Преобразуем систему перепишем 1-ю и 2-ю строки, а из 3-ей строки |
||||
|
0 |
|
1 |
1,7 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,333 |
|
|
|
|
|||||
вычтем 2-ю строку, результат запишем на место 3-ей строки. Получим следующую систему |
||||||||||||||
1 |
0,7 0,49 |
|
0,6 |
|
|
|
||||||||
|
|
|
|
|||||||||||
|
0 |
|
1 |
1,55 |
|
2 |
|
|
. Запишем полученную матрицу в виде системы. |
|
||||
|
0 |
|
0 |
0,15 |
|
|
|
|
|
|
|
|
|
|
|
|
|
1,667 |
|
|
|
||||||||
a2 |
0 |
a2 0,7 a2 |
2 |
0,49 0,6 |
1,667 |
|
||||||||
|
|
1 |
|
|
|
|
|
|
|
a22 11,111. Из |
||||
|
|
|
a21 a22 |
1,55 |
2 |
. Из 3 уравнения найдем a22 0,15 |
||||||||
|
|
|
a22 0,15 1,667 |
|
|
|||||||||
|
|
|
|
|
2-го уравнения найдем a21 . a21 2 1,55 a22 a21 2 1,55 11,111 a21 19,222. Из 1-
го уравнения получим a20 a20 0,6 0,7 a21 0,49 a22
a21 0,6 0,7 19,222 0,49 ( 11,111) a20 7,411.Запишем найденное уравнение
P22(x) 7,411 19,222 x 11,111 x2 .
Проверка. Найденное уравнение должно проходить через точки x2 , y2 , x3 , y3 , x4 , y4 .
P22(x2 ) 7,411 19,222 x2 11,111 x22 P22(0,7) 7,411 19,222 0,7 11,111 0,72
P22(0,7) 0,6
P22(x3 ) 7,411 19,222 x3 11,111 x32
P22(0,85) 7,411 19,222 0,85 11,111 0,852
P22(0,85) 0,9
P22(x4 ) 7,411 19,222 x4 11,111 x42 P22(1) 7,411 19,222 1 11,111 12
P22(1) 0,7 Следовательно, полученный полином проходит через 2- ю, 3-ю и 4-ю точки.
0,567 4 x 3,333 x2 ,если 0,2 x 0,7
Запишем ответ P2( x)
7,411 19,222 x 11,111 x2 ,если 0,7 x 1
Построим график
Реализация метода в Mcad
59
Метод неопределѐнных коэффициентов(кусочно-параболическая интерполяция)
|
0.2 |
|
|
|
0.1 |
|
Исходные данные |
|
|
||||||
|
0.4 |
|
|
|
0.5 |
|
|
|
|
|
|
|
|
||
x |
0.7 |
|
y 0.6 |
i 0 4 |
|||
|
0.85 |
|
|
|
0.9 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
1 участок |
|
0.7 |
|
|
|
|
|
|
|
|
|
|
x0 |
x0 |
2 |
1 |
|
||
|
|
x1 |
|
C1 1 x1 |
2 |
||
|
|
x2 |
|
|
x2 |
2 |
|
1 |
|
a1 C1 1 D1
|
|
0.5667 |
||
a1 |
|
4 |
|
|
|
|
|||
|
|
|||
|
|
3.3333 |
y |
|
|
Матрицы C1 и D1 для системы линейных |
|
0 |
уравнений |
|||
|
|
|||
D1 y |
1 |
|
|
|
|
|
|
||
y |
|
|
|
|
|
2 |
|
|
Решение системы линейных уравнений нахождение коэффициентов а10 , а11 и а12
|
|
|
2 |
|
|
|
Интерполирующая функция 1 участок |
|
P21(x) a10 a11 x |
a12 x |
|
|
|
||||
|
|
|
|
|
||||
2 участок |
|
|
|
|
|
|
||
|
x2 |
2 |
|
|
y |
|
|
Матрицы C2 и D2 для системы линейных |
|
|
|
||||||
1 x2 |
|
|
|
2 |
уравнений |
|||
|
x3 |
|
|
|
|
|
|
|
C2 1 x3 |
2 |
D2 |
y3 |
|
|
|||
|
|
|
|
y |
|
|
|
|
|
x4 |
2 |
|
|
|
4 |
|
|
1 x4 |
|
|
|
|
|
|
|
|
a2 C2 1 D2 |
|
Решение системы линейных уравнений |
||||||
|
|
|
|
нахождение коэффициентов а20 , а21 и а22 |
7.4111 a2 19.222211.1111
|
2 |
Интерполирующая функция 2 участок |
||
P22(x) a20 a21 x a22 x |
|
|||
|
|
|
||
P2(t) |
if t x2 P21(t) P22(t) |
Ответ |
||
|
||||
xp 0 0.05 |
1 |
|
точки для построения полинома |
|
|
|
1 |
|
|
P2 ( xp ) |
|
|
0 |
|
|
y i |
|
|
1 0 |
0.5 |
1 |
|
xp xi |
|