Автор работы: Пользователь скрыл имя, 07 Января 2011 в 09:45, контрольная работа
6 задач.
Сравнивая полученные
суммы квадратов отклонений, делаем
вывод, что лучше в смысле МНК
является квадратическая зависимость.
Задача 3.
Методом Гаусса-Зейделя
для решения СЛАУ
По данным своего варианта записать систему трёх уравнений с тремя неизвестными, привести её к виду с диагональным преобладанием и решить методом Гаусса-Зейделя с точностью ε = 0.05. Сделать проверку.
Дано:
3.1x1+ 2.8x2+ 1.9x3= 0.2 (I)
1.9x1+ 3.1x2+ 2.1x3= 2.1 (II)
7.5x1+
3.8x2+ 4.8x3=
5.6 (III)
Выберем начальное приблежение:
Х1͑º͗ = Х2º͗ = Х3º͗ = 0.
Любое следующее
приблежение будем вычислять
по формулам:
Х1 = 1/5.9*(-16,8 + 2.2 x3 )
Х2 = 1/7.1*(9.7 - 0.2x1 -4.8 x3 )
Х3 = 1/6.0*(17 - 0.3x1 -5.6 x2 )
На каждом шаге необходимо проверять условие на окончание итерационного процесса:
|xi – xi |˂ 0.05, i=1,2,3
Результаты вычислений приведены в таблице:
Номер
итерации |
Х1 | Х2 | Х3 | max |xi – xi | |
0 | 0 | 0 | 0 | |
1 | -2,8475 | 1,3662 | 2,8333 | 2,8475 |
2 | -1,7910 | -0,4691 | 1,7005 | 1,8352 |
3 | -2,2133 | 0,2669 | 3,3607 | 1,6601 |
4 | -1,5943 | -0,8434 | 2,6948 | 1,1104 |
5 | -1,8426 | -0,4107 | 3,7003 | 1,0054 |
6 | -1,4677 | -1,0835 | 3,3088 | 0,6727 |
7 | -1,6137 | -0,8294 | 3,9179 | 0,6091 |
8 | -1,3865 | -1,2371 | 3,6881 | 0,4077 |
9 | -1,4722 | -1,0881 | 4,0573 | 0,3691 |
10 | -1,3346 | -1,3353 | 3,9225 | 0,2471 |
11 | -1,3848 | -1,2480 | 4,1463 | 0,2238 |
12 | -1,3014 | -1,3979 | 4,0674 | 0,1499 |
13 | -1,3308 | -1,3469 | 4,2031 | 0,1357 |
14 | -1,2802 | -1,4379 | 4,1570 | 0,0909 |
15 | -1,2974 | -1,4081 | 4,2393 | 0,0823 |
16 | -1,2667 | -1,4633 | 4,2124 | 0.0552 |
17 | -1,2767 | -1,4459 | 4,2624 | 0,0499 |
Таким образом, условие на окончание итерационного процесса выполнено после 17 шагов. Приближённое решение:
Х1 = -1,2767 | Х2 = -1,4459 | Х3 = 4,2624 |
Выполним проверку,
подставив полученное решение в
исходную систему:
3.1*(-1,2767)+ 2.8*(-1,4459)+ 1.9*4,2624= 0,0920 (I)
1.9*(-1,2767)+ 3.1*(-1,4459)+ 2.1*4,2624= 2,0428 (II)
7.5*(-1,2767)+ 3.8*(-1,4459)+ 4.8*4,2624= 5,3895 (III)
Задача 4. Решение нелинейных алгебраических уравнений методом простой итерации.
С точностью ε = 0.05 найти приближённое значение корня уравнения F(x) = 0 на интервале изоляции [a,b] методом простой итерации.
Уравнение f(x) = 0 имеет вид:
D1x³ = D2x² + D3x + D4 =0.
Дано нелинейное алгебраическое уравнение
F(x) = x³ - 10x – 5 = 0
И интервал изоляции корня [3,4] методом простой итерации.
Для применения метода простой итерации приведём исходное уравнение к виду:
X = φ(x),
Где φ(x) = x + c*(x³ - 10x – 5); c- const.
Подберём константу c так, чтобы выполнялось условие:
|φˈ(x)|˂1 для 3 ≤ x ≤ 4
φˈ(x) = 1+c*(2x²-10),
φˈ(3) = 1+ 8c, φˈ(4) = 1+ 14c.
Таким образом, для c можно взять значение -0.05. Убедимся в правильности построения φ(x).
max3 ≤ x ≤ 4|φˈ(x)| = max3 ≤ x ≤ 4|1-0.05(2x²-10)| = 0.1˂1
Итерациональная формула при выбранном c будет иметь вид
Xk+1=Xk – 0.05(Xk³ - 10Xk – 5) или
Xk+1= – 0.05Xk³ + 1.5Xk +
0.25
В качестве начального приближения выберем X0 = 3.5.
Результат вычислений оформим в виде таблицы
k | Xk | 0.05Xk³ | 1.5Xk | Xk+1 | |Xk+1 - Xk| |
0 | 3.5 | 2,14375 | 5.25 | 3,35625 | 0,14375 |
1 | 3,35625 | 1,89031 | 5,04375 | 3,39406 | 0,03781 |
2 | 3,39406 | 1,95491 | 5,09109 | 3,38618 | 0,00788 |
Ответ: Х ≈ 3,38618
Проверка: 3,38618³
- 10 * 3,38618 -5 = -0,03513.
Задача 5. Вычисление определённого интеграла по формулам трапеций и Симпсона.
а) по формуле трапеций;
б) по формуле Симпсона.
2. Оценить
погрешность вычисления
Дано
F(x) = 1/x² , количество узловых точек на [1,3] , n = 9
Решение
Вычислить интеграл ∫ 1/x² dx по формулам:
а) трапеций;
б) Симпсона.
Определим шаг интегрирования:
i | xi | 1/xi² |
0 | 1 | 1 |
1 | 1.25 | 0.64 |
2 | 1.5 | 0.444444 |
3 | 1.75 | 0.326531 |
4 | 2.0 | 0.25 |
5 | 2.25 | 0.197531 |
6 | 2.5 | 0.16 |
7 | 2.75 | 0.132231 |
8 | 3 | 0.111111 |
Вычисление интеграла по формуле трапеций:
I ≈ ∆x( (y0+y8)/2 + y1+ y2+ y3+ y4+ y5+ y6+ y7) =0,67657