第二章 机械工程学院 动力工程.docVIP

机械工程学院 动力工程 王加佳 4.用Doolittle分解法解方程组 解：由题意得，原方程组的系数矩阵与常数项分别为 求解，即 求解，即 已知，求. 解：由得： 由得： 由得： 如果P为正交矩阵，证明。 证：由题知P为正交矩阵； 又 21.设常数，求t的取值范围，使得解方程组的Jacobi迭代法收敛。 解：设则方程组可表示为；由定了2.24得当时收敛。 令，其中 则可改写为 存在，则 令 当时；当时； 综上当或者时Jacobi迭代法收敛。 编写LU分解法、改进平方根法、追赶法的Matlab程序，并进行相关数值实验。 Doolittle分解法解第3题： function [L,U]=Doolittle(A) b=size(A); n=b(1); if b(1)~=b(2) error(MATLAB:Doolittle:Input Matrix should be a Square matrix. See Doolittle.); end if n~=rank(A) error(MATLAB:Doolittle:Input Matrix should be FULL RANK. See Doolittle.); end L=zeros(n,n); U=zeros(n,n); for i=1:n L(i,i)=1; end for k=1:n for j=k:n temp_sum=0; for t=1:k-1 temp_sum=temp_sum+L(k,t)*U(t,j); end U(k,j)=A(k,j)-temp_sum; end for i=k+1:n temp_sum=0; for t=1:k-1 temp_sum=temp_sum+L(i,t)*U(t,k); end L(i,k)=(A(i,k)-temp_sum)/U(k,k); end end end 指令窗输入： A=[1 0 2 0;0 1 1 1;2 0 -1 1;0 0 1 1]; doolittle(A) 结果： ans = 1.0000 0 0 0 0 1.0000 0 0 2.0000 0 1.0000 0 0 0 -0.2000 1.0000 Crout分解法解第3题： function [L,U]=Crout(A) b=size(A); n=b(1); if b(1)~=b(2) error(MATLAB:Crout:Input Matrix should be a Square matrix. See Crout.); end if n~=rank(A) error(MATLAB:Crout:Input Matrix should be FULL RANK. See Crout.); end L=zeros(n,n); U=zeros(n,n); for i=1:n U(i,i)=1; end for k=1:n for i=k:n temp_sum=0; for t=1:k-1 temp_sum=temp_sum+L(i,t)*U(t,k); end L(i,k)=A(i,k)-temp_sum; end for j=k+1:n temp_sum=0; for t=1:k-1 temp_sum=temp_sum+L(k,t)*U(t,j)