matlab之UL分解实验报告.docVIP

计算机上机实验报告 课程名称：MATLAB数值计算 教师姓名：谭小球 日期：2016年9月12日 学号：150607116 学生姓名：李妙锋 班级：Z15计算机 院系：数理与信息学院 实验题目：UL分解实验报告 1、实验名称 实现LU算法，并利用该算法求解线性方程组 2、实验目的 了解如何用LU三角分解法解线性方程组，利用LU三角分解法解线性方程组 3、实验原理 设无行交换变换的高斯消去法可求解一般线性方程组AX=B，则矩阵A可分解为一个下三角矩阵L和一个上三角矩阵U乘积： A=LU 而且L的对角线元素为1，U的对角线元素非零。得到L和U后，可通过如下步骤得到X： 利用向前替换法对方程组LY=B求解Y。 利用回带法对方程组UX=Y求解X。 4、实验步骤 输入矩阵A LU分解A，得到L矩阵与U矩阵的值 [L U]=LU_1(A) 输入矩阵B，利用向前回带法求出Y值 Y=upsub(L,B) 利用回带发求出X值 [X]=backsub(U,Y) 5、实验程序 LU分解 代码： function [L U]=LU_1(A) n=length(A(1,:)); L=eye(n); U=zeros(n); for j=1:n U(1,j)=A(1,j); end for i=2:n L(i,1)=A(i,1)/U(1,1); end for k=2:n for j=k:n U(k,j)=A(k,j)-L(k,1:k-1)*U(1:k-1,j); end for i=k+1:n L(i,k)=(A(i,k)-L(i,1:k-1)*U(1:k-1,k))/U(k,k); end end 结果： 向前回带法 代码： %向前代入法 function Y=upsub(A,B) n=length(B); Y=zeros(n,1); Y(1)=B(1)/A(1,1); for k=2:n Y(k)=(B(k)-A(k,1:k-1)*Y(1:k-1))/A(k,k); end 结果： 回带法 代码： %回代法 function [X]=backsub(A,B) n=length(B); X=zeros(n,1); X(n)=B(n)/A(n,n); for k=n-1:(-1):1 X(k)=(B(k)-A(k,k+1:n)*X(k+1:n))/A(k,k); end 结果： 6、实验分析 LU分解法比较简便迅速，当解多个系数矩阵为A的线性方程做时，LU分解法就显得特别优越，只要对系数矩阵做一次LU分解，以后只要解三角形方程即可。也可以根据系数矩阵的形状来设计算法。