Exp05姜启源数学建模课件第五章：线性代数方程组的数值解法.pptVIP

大学数学实验; 许多实际问题归结为线性（代数）方程组;3. 线性方程组数值解法的MATLAB实现 ;线性方程组的一般形式、两类解法;直接法 - 高斯消元法;直接法 - 高斯消元法的矩阵表示;直接法 - 高斯消元法的矩阵表示;直 接 法 - 矩 阵 LU 分 解;若A可逆，则存在交换阵 P 使 PA=LU L为单位下三角阵，U为上三角阵。;直接法 - 对称正定矩阵的分解;直接法 - 三对角矩阵的LU分解 ;线性方程组Ax=f可通过等价的两个三角形方程组Ly=f和Ux=y求解如下 :;方程组的病态性;向量和矩阵的范数;1）设b有扰动 ?b，分析x的变化 ?x;x的(相对)变化不超过b的(相对)变化的Cond(A)倍, 也大致上是A的(相对)变化的Cond(A)倍。;迭代法 - 一个例子;迭代法 - 雅可比迭代;迭代法 – 高斯-塞德尔迭代;迭代法的收敛性;迭代法的收敛性;;;超定线性代数方程组的最小二乘解 ;+; 先选定一组函数 r1(x), r2(x), …rm(x), mn, 令 f(x)=a1r1(x)+a2r2(x)+ …+amrm(x) （1） 其中 a1,a2, …am 为待定系数。;线性最小二乘法的求解;当 RTR 可逆时(4)有唯一解;为什么要规定mn?若m=n或mn,会如何?;基函数 {r1(x), …rm(x)}的选取 ; 线性方程组数值解法的MATLAB实现 ; 线性方程组数值解法的MATLAB实现 ;例. 解;当n很大时Hilbert矩阵呈病态;H=hilb(5), h=rats(H), b=ones(5,1); x=H\b; b(5)=1.1; x1=H\b; [x,x1], n1=cond(H), n2=rcond(H), ;1. 提取（产生）对角阵;例. 用迭代法解;MATLAB对稀疏矩阵的处理: 进行大规模计算的优点;n=500;b=[1:n]; a1=sparse(1:n,1:n,4,n,n); a2=sparse(2:n,1:n-1,1,n,n); a=a1+a2+a2; tic;x=a\b;t1=toc aa=full(a); tic;xx=aa\b;t2=toc y=sum(x) yy=sum(xx);用MATLAB作线性最小二乘拟合;实例1 投入产出模型;4）如果对于任意给定的、非负的外部需求，都能得到非负的总产出，模型就称为可行的。问为使模型可行，投入系数应满足什么条件？;1）基本模型;基本模型;若?d=(1,0,0)T, 即农业外部需求增加1单位时，三部门总产出应分别增加1.3459，0.5634，0.4382单位。 即C的第1列。 C的第2，3列给出了什么？;4）如果对于任意给定的、非负的外部需求，都能得到非负的总产出，模型就称为可行的。问为使模型可行，投入系数应满足什么条件？;模型可行;实验2给出的模型为 ; 开始有100棵植物，要求50年后有1000棵植物 ;实例3 汽车刹车距离(续);布置实验