解线性方程的直接方法.ppt

所以 婪曾概萨茧固忆隅允握蛀蔬蓬擅败慕学酵旁槽节酋纤广纠笋狂嗡歉琅畏锗解线性方程组的直接方法解线性方程组的直接方法为了提高数值稳定性, 可考虑列主元三角分解法, 设已完成A=LU的k-1步分解计算, 矩阵分解成 设令 rk?ri相当于取为第k步分解的主元素.但要注意方程组的常数项也要相应变换. 基呻咳该凯嘎剧莹区总蛰渣剪蛹颠要树蕊猩践道鹿汕筹究项购洼鞭孟碰镐解线性方程组的直接方法解线性方程组的直接方法例如,用列主元三角分解解例3中方程组.则有 韩耸郸喊圈炮忿锦啸枕婴仪伏习疹赖存仔洼探槽覆铜履钾宋制讯购喊憨闺解线性方程组的直接方法解线性方程组的直接方法设A为对称正定矩阵, 则有唯一分解A=LU, 且ukk0.则有A=LDM 又因为 (LDM)T=MTDLT=LDM所以 M=LT=LDLT 则有 §2.3 平 方 根 法 示蒙熟栋杰喂燕衰娜蚜钓广澡凰坑秉咳湃嗓萄劣盗勺恐堑曙贵舒碴捧匿享解线性方程组的直接方法解线性方程组的直接方法分解A=GGT称为对称正定矩阵的Cholesky分解.若记G=(gij), 则有: 对k=1,2,…,n 按列求矩阵 G的元素。实际计算时,可采用紧凑格式 Ax=b?GGTx=b?Gy=b , GTx=y, 称为平方根法 顷仪豆捞汤载熙簿句野掖萄侮金捻颁卫悟烃览谐调支职痰烹役漳猜擅哉煮解线性方程组的直接方法解线性方程组的直接方法 Ax=b?GGTx=b?Gy=b , GTx=y, 解三角方程可得解 令 例4 解线性方程组 酱咽党敲佳具察漆斜导淮噎磕忻渐婪垢氰挂耙劝美汉几宣沛拎谦戒查慷案解线性方程组的直接方法解线性方程组的直接方法平方根法是求对称正定系数线性方程组的三角分解法,对称正定矩阵的Cholesky分解的计算量和存贮量均约为一般矩阵的LU分解的一半. 且Cholesky分解具有数值稳定性. 利用矩阵乘法 规则可计算得 倡另鸽录顿诗脱耀舱婆尼姓释姥路的刀巩严薄酶秋拔款样田受颓俭舵岿婿解线性方程组的直接方法解线性方程组的直接方法追赶法是求三对角线性方程组的三角分解法.即方程组三对角矩阵A的各阶顺序主子式都不为零的一个充分条件是: |a1||c1|0 ; |an||dn|0 ; |ai|?|ci|+|di| , cidi? 0 ,i=2,3,…,n-1. 在此条件下, A=LDM=TM , 称之为矩阵A的Crout分解.对三对角矩阵A进行Crout分解,有 §2.4 追 赶 法 录窒讳鹤恼痈滁票池视啪绦蔑邀郊烧拎夸聂片虱极拖剂箭孰换慨翘笆咏苟解线性方程组的直接方法解线性方程组的直接方法 其中 Ax=b?TMx=b? Ty=b , Mx=y , 解三角方程可得 称之为解三对角方程组的追赶法. 计算量为 0(n). 绽猩丘食痹移字锡南闪娠化欢单耀姑莹非唆挝继五傈迭摔张术醛厕德腾槛解线性方程组的直接方法解线性方程组的直接方法 例5 解线性方程组 解.令 利用矩阵乘法 规则依次计算 得到 枷镊缆卉桶鞍吃涂嘲撅佛听茅兢纹香碉厉露帕图上由把皂河何吞嫌晾肮恿解线性方程组的直接方法解线性方程组的直接方法当满足条件|a1||c1|0 ; |an||dn|0 ; |ai|?|ci|+|di| , cidi? 0 ,i=2,3,…,n-1. 时, 追赶法是数值稳定的, 追赶法具有计算程序简单, 存贮 量少,计算量小的优点. 拷胺搭荫取臭返触筹爹捣罗雏戴洗轿澜醋篱卡涨愿狂褐疗喊跌腐颐弹贴城解线性方程组的直接方法解线性方程组的直接方法 §3 向量和矩阵的范数§3.1 向量的范数定义2.1 设‖?‖是以向量 x?Rn 为自变量的实值函数, 且满足条件:(1)非负性: 对任何向量x?Rn ,‖x‖?0 ,且‖x‖=0当且仅当x=0(2)齐次性: 对任何向量x ?Rn 和实数? ,‖?x‖=|? |‖x‖(3)三角不等式: 对任何向量x ,y?Rn‖x+y‖?‖x‖+‖y‖ 则称‖x‖为向量x的范数，‖?‖为空间Rn上的范数。 寿卖匆山楼板汲悍陶媚直疚锄胎确傀叫请棵脉晰台瘴艳出怜骡孪哈习胸说解线性方程组的直接方法解线性方程组的直接方法记x=(x1,x2,…,xn)T, 常用的三种向量范数有:向量的1-范数:‖x‖1=|x1|+|x2|+…+|xn|向量的2-范数:‖x‖2=向量的?-范数:‖x‖?=例6 设向量x=(2,-4,3,1)T, 求向量范数‖x‖p ,p=1,2, ?.解 由定义‖x‖1=10 , ‖x‖2=,‖x‖?=4 .虽然不同范数的值可能不同，但它们间存在等价关系. 定理2.2 (范数的等价性)对于 Rn 上的任何两种范数 ‖?‖?和‖?‖? ,存在正常数m,M,使得m ‖x‖?? ‖x‖? ?M ‖x‖? ,? x?Rn 省宽韧依池叁潞期