Курсовая работа: Метод квадратных корней для симметричной матрицы при решении систем линейных алгебраических уравнений
for k=(i+1):n
c=c+t(i,k)*x(k);
end
x(i)=(y(i)-c)/t(i,i);
end
else
error('Внимание! Ошибка! Размерность вектора F не соответствует размерности матрицы А');
end
else
error('Внимание! Ошибка! Определитель матрицы А равен 0');
end
end
3. Описание тестовых задач
После того, как функция была разработана, для ее отладки была составлена программа, где задавались матрица А, вектор fи откуда вызывалась написанная функция.
Программа имеет вид:
a=[1 0 0; 0 1 0; 0 0 1];
f=[7;8;9];
[e,x]=mkk(a,f)
Решение для данной программы выдано такое:
e =
0
0
0
x =
7
8
9
Как видим, решение правильное.
Начнем исследование метода квадратных корней. Для начала исследуем влияние мерности матрицы А на точность решения.
Для этого будем последовательно решать СЛАУ, каждый раз увеличивая мерность А. Для этого составим такую программу, которая