2
:
˳ .
3-1: ..
:
.
: . m - MatLab.
ճ
1. :
. ( ) . , . ().
() . ( ), . : ( : , , .), , LU-. , , 103.
. , , . ε 0 (01, 02, 03, ). : , , , 䳺 . 106 .
(2.3) .
1) ( = (ij) ). .
2) , . , n - , n-1, n-2, n-33, 2, 1 .
2. :
atLAB, ( ) LU- ( ). , . , , ().
|
|
5
3. ˳ :
function X=Gauss(A,B,n)
A=[5 2 1;1 3 7;4 6 -1];
B=[3 14 21 ];
n=3;
X=zeros(n,1);
for i=1:1:n-1
for t=i+1:1:n
k=A(i,i)/A(t,i);
for j=i:1:n
A(t,j)=k*A(t,j);
A(t,j)=A(t,j)-A(i,j);
end
end
end
A
for i=n:-1:1
X(i)=B(i)/A(i,i);
for j=n:-1:1
A(j,i)=A(j,i)*X(i);
B(j)=B(j)-A(j,i);
end
end
4.
: MatLab.