Методы решения задач математической физики. Рындин Е.А. - 48 стр.

UptoLike

Составители: 

48
% B - вектор-столбец свободных членов;
% X0 - вектор-столбец начального приближения;
% delta - допустимая относительная погрешность;
% Imax - максимальное число итераций;
% X - приближенное решение линейной системы AX = B.
n=length(B);
err=5*delta;
X=X0;
ct=0;
while err>delta
Xp=X;
for i=1:n
X(i)=(B(i)-A(i,[1:i-1,i+1:n])*…
Xp([1:i-1,i+1:n]))/A(i,i);
end
if max(abs(X0))==0
error('Следует изменить начальное приближение')
break
end
err=max(abs(X-Xp))/max(abs(X0));
ct=ct+1;
Xe(ct,:)=[X' err];
if ct>Imax
error('Требуемое число итераций слишком велико')
break
end
end
Xe
ct
Получить информацию о функциях и синтаксисе написания программ в
системе инженерных и научных расчетов MATLAB можно в [10, 11].
Для запуска описанного выше итерационного процесса с допустимой от-
носительной погрешностью результата 10
-5
и максимальным числом итераций
не более 100, необходимо войти в систему MATLAB, создать файл функции и
сохранить его под именем yakobi.m в рабочем каталоге /WORK , находящемся в
корневом каталоге системы MATLAB, после чего в командной строке последо-
вательно набрать команды, приведенные ниже.
A=[ 2 -1 1; 4 -6 1; -2 1 8];
B=[2; -4; 16];
X0=[1; 2; 2];
X=yakobi(A,B,X0,1e-5,1e2)