matlab(实验5).doc

实验5
1：
左除命令：
A=[1 -9 -10;-9 1 5;8 7 1];
>> b=[-1 0 4]';
>> x=A\b
x =
-

-
Lu分解实验5
1：
左除命令：
A=[1 -9 -10;-9 1 5;8 7 1];
>> b=[-1 0 4]';
>> x=A\b
x =
-

-
Lu分解：
A=[1 -9 -10;-9 1 5;8 7 1];
b=[-1 0 4]';
[l,u]=lu(A);
y=l\b;
x=u\y
x =
-

-
迭代
雅可比迭代：
function x=yacbic(A,b)
eps=10^(-6);
n=1000;
D=diag(diag(A));
L=-triu(A,-1);
U=-tril(A,1);
B=inv(D)*(L+U);
f=inv(D)*b;
k=0;
x0=[1,1,1]';
while k<n
x1=B*x0+f;
k=k+1;
if norm(x1-x0,inf)<eps
disp(x1');
break
end
if k>n
end
x0=x1;
end
高斯—赛德尔迭代：
function x=agui(a,b)
n=length(b);
N=100;
e=1e-6;
x0=zeros(n,1);
x=x0;
x0=x+2*e;
k=0;
a1=tril(a);
a2=inv(a1);
while norm(x0-x,inf)>e&k<N
k=k+1;
x0=x;
x=-a2*(a-a1)*x0+a2*b;
format long
k
disp(x')
end
K
雅克比迭代
高斯--赛德尔迭代
0