Автор работы: Пользователь скрыл имя, 07 Января 2011 в 09:45, контрольная работа
6 задач.
Контрольная работа
№3
Задача 1. Аппроксимация функций. Интерполирование.
Дана таблица значений функции y=f(x) и два значения аргумента x=a и x=b.
x | 0,6 | 0,7 | 0,8 | 0,9 | 1,0 | a = 0,88 |
y | 0,5646 | 0,6442 | 0,7174 | 0,7833 | 0,8415 | b = 0,62 |
Требуется с
помощью полиномов Ньютона третьей
степени вычислить приближенное значение
функции y=f(x) в точках x=a и x=b, т.е. f(a) и f(b).
Решение
Построим таблицу конечных разностей:
i | xi | yi | ∆yi | ∆²yi | ∆³yi |
0 | 0,6 | 0,5646 | 0,0796 | -0,0064 | 0,0009 |
1 | 0,7 | 0,6442 | 0,0732 | -0,0073 | 0,0004 |
2 | 0,8 | 0,7174 | 0,0659 | -0,0077 | |
3 | 0,9 | 0,7833 | 0,0582 | ||
4 | 1,0 | 0,8415 |
Поскольку точка b= 0,62 находится ближе к началу таблицы, а точка a=0,88 – ближе к концу таблицы, то для нахождения значения f(b) будем использовать первый полином Ньютона, а для нахождения f(a) – второй. Ближайшее меньшее к точке b значение узла таблицы – 0.6, а ближайшее большее к точке a – 1,0.
Выпишем первый и второй полиномы третьей степени:
P (x) =y0 + (∆y0/1!)*q + (∆²y0/2!)*q*(q-1) + (∆³y0/3!)*q*(q-1)*(q-2),
q= (0.62-0.6)/0.1=0.2,
P (x) =y4 + (∆y3/1!)*q + (∆²y2/2!)*q*(q+1) + (∆³y1/3!)*q*(q+1)*(q+2),
q= (0.88-1.0)/0.1= -1.2,
f(0.62) ≈ P (0.62) =0,5646+
(0.0796/1!)*0.2 + (-0.0064/2!)*0.2*(0.2-1) + (0.0009/3!)*0.2*(0.2-1)*(0.2-
= 0.5811
f (0.88) ≈ P (0.88) = 0,8415 + (0.0582/1!)*(-1.2) + (-0.0077/2!)*(-1.2)*(-1.2+1) +
+(0.0004/3!)*(-1.2)*(-1.2+1)*(
Задача 2. Аппроксимация функций. Подбор эмпирической зависимости.
1. Используя систему нормальных уравнений, найти параметры a0, a1 для линейной ỹ=a0+a1x и a0,a1,a2 для параболической ỹ=a0+a1x+a2x² эмпирических зависимостей. Записать эти зависимости.
2. Построить
на одном чертеже графики
3. Определить, какая из двух эмпирических зависимостей (линейная или параболическая) лучше в смысле МНК.
Дана таблица значений функций y=f(x)
x | 0.45 | 0.58 | 0.89 | 0.97 | 1.15 | 1.35 | 1.43 | 1.55 | 1.65 | 1.71 |
y | 1.87 | 2.05 | 2.66 | 2.85 | 3.36 | 4.03 | 4.52 | 4.80 | 5.23 | 5.50 |
Найти эмпирические зависимости для двух случаев:
ỹ=a0+a1x и
ỹ=a0+a1x+a2x²
и определить, какая из них лучше описывает табличные данные.
Составим таблицу промежуточных вычислений:
i | xi | yi | xiyi | xi² | xi³ | xi | xi²yi |
1 | 0.45 | 1.87 | 0.841 | 0.202 | 0.091 | 0.041 | 0.377 |
2 | 0.58 | 2.05 | 1.189 | 0.336 | 0.195 | 0.113 | 0.689 |
3 | 0.89 | 2.66 | 2.367 | 0.792 | 0.704 | 0.627 | 2.106 |
4 | 0.97 | 2.85 | 2.764 | 0.941 | 0.912 | 0.885 | 2.682 |
5 | 1.15 | 3.36 | 3.864 | 1.322 | 1.521 | 1.749 | 4.442 |
6 | 1.35 | 4.03 | 5.440 | 1.822 | 2.460 | 3.321 | 7.342 |
7 | 1.43 | 4.52 | 6.463 | 2.045 | 2.924 | 4.181 | 9.243 |
8 | 1.55 | 4.80 | 7.440 | 2.402 | 3.723 | 5.772 | 11.529 |
9 | 1.65 | 5.23 | 8.629 | 2.722 | 4.492 | 7.412 | 14.236 |
10 | 1.71 | 5.50 | 9.405 | 2.924 | 5.000 | 8.550 | 16.082 |
Ʃ | 13.44 | 36.87 | 48.402 | 15.508 | 22.022 | 32.651 | 68.728 |
Для определения параметров a0 и a1 линейной эмпирической зависимости подставим необходимые значения найденных сумм в систему нормальных уравнений.
Получим
10a0 + 13.440a1 = 36.870
13.44a0 + 15.508a1 = 48.402
Решив любым методом, например методом последовательного исключения неизвестных, систему, получим a0 = 3.081 и a1 = 0.451. Следовательно, линейная эмпирическая зависимость имеет вид:
ỹ=3.081+0.451x.
Поступая аналогичным образом, запишем систему нормальных уравнений для определения a0,a1,a2 квадратической эмпирической зависимости:
10a0 + 13.44a1 + 15.508a2 = 36.84
13.44a0 + 15.508a1 + 22.022a2 = 48.402
15.508a0 + 22.022a1 + 32.651a2 = 68.728
Решив эту систему, получим:
a0 =0.481, a1 = 0.982, a2 = 1.214
ỹ = 1.214x² + 0.982x+ 0.481
Построим график
f(x)=1.214x^2+0.982x+0.481
x f(x) |
-4 15.977 |
-3 8.461 |
-1 0.713 |
0 0.481 |
1 2.677 |
2 7.301 |
3 14.353 |
f(x)= 0.451x + 3.081
x f(x) |
-4 3.081 |
3 3.081 |
Чтобы определить,
какая из двух эмпирических зависимостей,
линейная или квадратическая, лучше
описывает табличные данные, вычислим
суммы квадратов отклонений значений
каждой из указанных эмпирических зависимостей
от табличных:
i | xi | yi | [yi +3.081+0.451xi)]^2 | [yi-(1.214xi^2+0.982xi+0.481)] |
1 | 0,45 | 1,87 | 1,999255 | 0,491773 |
2 | 0,58 | 2,05 | 1,670763 | 0,349341 |
3 | 0,89 | 2,66 | 0,676325 | 0,117931 |
4 | 0,97 | 2,85 | 0,446852 | 0,07519 |
5 | 1,15 | 3,36 | 0,057432 | 0,020789 |
6 | 1,35 | 4,03 | 0,115702 | 0,000116 |
7 | 1,43 | 4,52 | 0,630547 | 0,023174 |
8 | 1,55 | 4,8 | 1,040298 | 0,014336 |
9 | 1,65 | 5,23 | 1,973604 | 0,031122 |
10 | 1,71 | 5,5 | 2,715212 | 0,044133 |
Ʃ | 11,32599 | 1,167905 |