计算方法实验报告
指导老师：易昆南 专业班级：统计0802班 学 号：1304080408 姓 名：梁远飞
学号 1304080408 实验题目 1、设计（实习）目的： 1.了解直接分解以及选主元分解法的原理 2.利用直接分解和选主元分解求解方程组 2、实验内容： A=[1,2,3;2,5,2;3,1,5]，b=[14,18,20]';A*x=b,用直接分解法和选主元分解法求解。 3、详细设计： 1. 直接分解法 function f=sanjiao(A,b) A=[1,2,3;2,5,2;3,1,5]; b=[14,18,20]'; [l,u,p]=lu(A) B=p*b; C=l\\B; x=u\\C 2．function F=LU_P(A) %选主元的三角分解法 %[L,A,P]=LU_P(A) % A 待分解阵 % L 返回单位下三角阵 % A 返回上三角阵,存储在原矩阵中 % P 返回排列阵 A=[ 1,2,3;2,5,2;3,1,5]; b=[14,18,20]'; n=size(A); if n(1)~=n(2) error('不是方阵!'); end n=n(1);flag=1:n; P=eye(n,n);L=P; for k=1:n if k~=1 A(k:n,k)=A(k:n,k)-A(k:n,1:k-1)*A(1:k-1,k); end t=find(abs(A(k:n,k))==max(abs(A(k:n,k)))); t=t(1)+k-1; %选主元 flag(k)=t; if t~=k p=A(k,:); A(k,:)=A(t,:); A(t,:)=p; end %交换行 班级 统计0802 姓名 梁远飞 指导教师 评 分 易昆南 三角分解法 if abs(A(k,k))
