中科院矩阵分析与应用大作业.docxVIP

中科院矩阵分析与应用大作业 中科院矩阵分析与应用大作业 PAGE 中科院矩阵分析与应用大作业 中科院矩阵分析与应用大作业 实现LU分解 QR分解 Householder reduction、Givens reduction Matlab 代码： function [] =juzhendazuoye A=input(请输入一个矩阵A=); x=input(请输入序号 1 LU分解 2 Gram-Schmidt分解 3 Householder reduction 4 Givens reduction: ); if(x==1) %%*************LU分解*****************%% disp(PA=LU) m=size(A,1); % m等于矩阵A的行数 n=size(A,2); % n等于矩阵A的列数 if(m==n) % 判断矩阵A是不是方阵 % 如果矩阵A不是方阵那么就输出“error” U=A; % 把矩阵A赋值给矩阵U L=zeros(n); % 先将L设为单位阵 P=eye(n); % 首先将交换矩阵P设为单位矩阵 for j=1:n-1 for i=j+1:n if (U(j,j)~=0) %判断主元元素是否不为0 L(i,j)=U(i,j)/U(j,j); U(i,:)=U(i,:)-U(j,:)*U(i,j)/U(j,j); % U(j,j)为主元元素 else a=j+1; % 令a等于j+1 while((U(a,j)==0)(an)) % 判断主元元素所对的下一行元素是不是0，a是否小于n a=a+1; % 寻找下一个元素 end temp=U(j,:); % 判断主元元素所在列（除主元元素外）第一个不为零的元素的所在行与主元元素所在行进行行交换 U(j,:)=U(a,:); % U两行交换位置 U(a,:)=temp ; m=L(j,:); L(j,:)=L(a,:); % L矩阵两行交换位置 L(a,:)=m; q=P(j,:); P(j,:)=P(a,:); % 交换矩阵的两行交换 P(a,:)=q; L(i,j)=U(i,j)/U(j,j); U(i,:)=U(i,:)-U(j,:)*U(i,j)/U(j,j); end end end for k=1:n L(k,k)=1; % 把L矩阵的对角线赋值为1 end L % 输出下三角矩阵L U % 输出上三角矩阵U P % 输出交换矩阵P A=inv(P)*L*U else disp(error) end end if(x==2) %% 判断如果x=2，那么将执行schmid分解 %%**************Gram-Schmidt正交分解*****************%% disp(A=Q*R) Q=zeros(size(A,1),size(A,2)); %% 先把Q设为全零矩阵 R=zeros(size(A,2)); %% R设置为全零矩阵 a=A(:,1);