p>
Функція Зейделя: [z1, z2] = zeide1 (a, b, eps)% векторне привласнення
a = [-1 4 серпня 4 червень; -5 2 -8 5 квітня; 4 -5 1 5 10; -4 1 Березня -7 5, 5 7 жовтня -4 -3]
b = [6; -5; 3, 5, 3]
eps = 10 ^ (-6); = size (a, 1)% повертає кількість рядків масиву А = a '* a = a' * b% приведення системи до нормального вигляду
for i = 1: N% розмірність задачі (i) = D (i)/C (i, i);
end
% D11 = D1 ';% D11-початкове наближення (стовпець правій частині)
d1 = D1; i = 1: Nj = 1: N
if i == j% перевірка на рівність
C1 (i, j) = 0; (i, j) =-C (i, j)/C (i, i);
end
% рішення системи методом Зейделя
R1 = d1; = 0; k == 0i = 1: N = C1 (i, 1: N);% вирізка = dot (v ', d1);% скалярний твір
d1 (i) = A + D1 (i); = d1; = max (abs (R2-R1)); S
Одержаний результат:
a =
8 Квітень 6 Квітня
2 -8 4 травня
-5 1 5 10
3 січня -7 5
10 липня -4 -3
b =
N =
5 =
4 67 -3 -19
154 84 -52 -41
84179 -22 -6
-52 -22 151 71
-41 -6 71166 =
= 1 through 4
-0.22784165101666 -0.11317732900043 0.77418152802107 -0.24929480082327
Column 5
.38178538703810
3. Інтерполяції ФУНКЦІЙ
.1 Завдання
З першої частини курсової роботи взяти функцію, яка має більше вигинів. На проміжку (-4?; 4?) Взяти рівномірну сітку з десяти значень (x). Обчислити значення функції в цих точках. Сформувати сітку c кроком 1/100. p> Виконати інтерполяцію 3 методами.
.2 Інтерполяція заданої функції
Часто для відновлення залежностей; порівняння значень двох або більше функцій, заданих на різній сітці; підготовки якісних графічних матеріалів потрібно побудувати інтерполюються функцію f (xi) = yi, яку називають інтерполянтой. Інтерполянту шукають у вигляді многочлена:
В
де-невизначені коефіцієнти.
Більш зручним є представлення полінома у вигляді полінома Лагранжа ступеня n:
Формула Лагранжа має вигляд:
.
Формат функції в Matlab:
= interp (x, y, xi, метод ). ...
