内容发布更新时间 : 2024/12/23 0:04:04星期一 下面是文章的全部内容请认真阅读。

数值分析作业

第二章

1、用Gauss消元法求解下列方程组： 2x1-x2+3x3=1, (1) 4x1+2x2+5x3=4, x1+2x2 =7; (2) 解：

A=[2 -1 3 1;4 2 5 4;1 2 0 7] n=size(A,1);x=zeros(n,1);flag=1; % 消元过程 for k=1:n-1 for i=k+1:n if abs(A(k,k))>eps

A(i,k+1:n+1)= A(i,k+1:n+1)-A(k,k+1:n+1)*A(i,k)/A(k,k); else flag=0; return end end end % 回代过程 if abs(A(n,n))>eps x(n)=A(n,n+1)/A(n,n); else

flag=0; return end for i=n-1:-1:1 x(i)=(A(i,n+1)-A(i,i+1:n)*x(i+1:n))/A(i,i); end return x

A = 2 -1 3 1 4 2 5 4 1 2 0 7

x = 9

-1 -6

11x1-3x2-2x3=3, (2) -23x1+11x2+1x3=0,

x1+2x2+2x3=-1;

(2) 解：

A=[11 -3 -2 3;-23 11 1 0;1 2 2 -1] n=size(A,1);x=zeros(n,1);flag=1; % 消元过程 for k=1:n-1 for i=k+1:n

if abs(A(k,k))>eps A(i,k+1:n+1)=

A(i,k+1:n+1)-A(k,k+1:n+1)*A(i,k)/A(k,k); else flag=0; return end end end

% 回代过程

if abs(A(n,n))>eps x(n)=A(n,n+1)/A(n,n); else flag=0; return end

for i=n-1:-1:1

x(i)=(A(i,n+1)-A(i,i+1:n)*x(i+1:n))/A(i,i); end return x

A = 11 -3 -2 3

-23 11 1 0 1 2 2 -1 x = 0.2124 0.5492

-1.1554

4、用Cholesky分解法解方程组

3 2 3 x1 5 2 2 0 x2 3 3 0 12 x3 7

解： . A=[3 2 3;2 2 0;3 0 12]; b=[5 3 7]; lambda=eig(A); if lambda>eps&isequal(A,A') [n,n]=size(A); R=chol(A); %解R'y=b y(1)=b(1)/R(1,1); if n>1 for i=2:n y(i)=(b(i)-R(1:i-1,i)'*y(1:i-1)')/R(i,i); end end

%解Rx=y

x(n)=y(n)/R(n,n); if n>1

for i=n-1:-1:1

x(i)=(y(i)-R(i,i+1:n)*x(i+1:n)')/R(i,i); end end x=x'; else x=[];

disp('该方法只适用于对称正定的系数矩阵!'); end

R= 1.7321 1.1547 1.7321 0 0.8165 -2.4495 0 0 1.7321

y= 2.8868 -0.4082 0.5774

x= 1.0000 0.5000 0.3333

5. 用列主元Doolittle分解法解方程组

解：A=[3 4 5; -1 3 4; -2 3 -5;]; 3 4 5 X1 2 b=[2,-2 6]'; -1 3 4 X2 -2 [L,U,pv]=luex(A); -2 3 -5 X3 6