第三章 线性代数方程组的直接解法2.pptVIP

* ?平方根法(Cholesky乔列斯基分解方法) 当 为实对称正定矩阵时，三角分解法的变形。 ?实对称正定矩阵的几个重要性质 ? A?1 亦对称正定，且 aii 0 ? A 的顺序主子阵 Ak 亦对称正定 ? A 的特征值 ?i 0 ? A 的全部顺序主子式 det ( Ak ) 0（充要条件） §3.3 平方根法 疆绊蛆终钮冲疵纵禾歼肇挪坷巩郡魔延黍帛蛮尹昧速收狄匡蝗煎阮练添缎第三章 线性代数方程组的直接解法2第三章 线性代数方程组的直接解法2 （Cholesky分解） 如果 是正定矩阵，则存在一个对角元素为正数的 下三角矩阵 ，使得 。 证明： 设 则 的所有顺序主子式为正 矩阵 存在Doolittle分解 易证 其中 为 的主对角元素，且有 泞又谁痊桌宵床拈梳膀枢先根悍潭功悬羞太隐愿痒攘热歼脾霞羊柑借遗黎第三章 线性代数方程组的直接解法2第三章 线性代数方程组的直接解法2 单位上三角 记 其中 存董匡舱孝智洋畜硷豫识盒爹粮潍笛穴畜段蓖洽搁绞绷释沥侣呆眠笆宽香第三章 线性代数方程组的直接解法2第三章 线性代数方程组的直接解法2 思想 ?Cholesky分解的计算公式 设 由 对应元素相等得 铀奇缠田鲜螺拔垛元羌拴蔽主肯秘畴赶芋蛹讳钨梗酝隐媒沤感膘括委糯敬第三章 线性代数方程组的直接解法2第三章 线性代数方程组的直接解法2 Cholesky分解公式 具酪燃汹畜浓尧切财氨碴倾篓码默氛工史怖砍螟美柞峙暮钞副底炽惨癣翅第三章 线性代数方程组的直接解法2第三章 线性代数方程组的直接解法2 因对称性无需存储 Step1 Step2 Step3 Stepn 的计算过程: 逐 列 计 算 元素 仍然存放在矩阵 的相应位置上 它杨纠龄绊谁挤脱掸桓饱睛莹读或就央伞喘卑信蹬芒火籽胯惋瘦乳蘸帆镰第三章 线性代数方程组的直接解法2第三章 线性代数方程组的直接解法2 矩阵 分解的实际计算公式(算法3.3.1)： for for for for 咽值效蚜酶厉聘近半革刃俞剁骏磷样浙拎测值驼蒲市揉熄咯猾铱坪拨岁吁第三章 线性代数方程组的直接解法2第三章 线性代数方程组的直接解法2 例5：用Cholesky分 解法求解下列方程组 解： 系数矩阵为 Step1 Step2 Step3 恋谎卿俐爆凛喳乞薪炬沦威叠坷烽奏帧栈惕汗灯藤溃娃抚麦吴署陆抠涧泻第三章 线性代数方程组的直接解法2第三章 线性代数方程组的直接解法2 求解方程组 求解方程组 典缩凹钉博迁赢仟哄淡吉指姜腥哭乍撬骡金肾陨艳醋掷袒截锗证褪器面逝第三章 线性代数方程组的直接解法2第三章 线性代数方程组的直接解法2 ? Cholesky分解法求解方程组中需说明的几个问题 ?工作量：约为 分解的一半； ?不必选主元： 的正定性和算法的稳定性 ?稳定性：是数值稳定的； ?缺陷：存在开平方运算。 ?改进方法： 分解 改进的平方根法 兑蒙属吻扯阁唤佃低霜剂钓蒙蒸汇腊问淖骋曙镊勉漏问浑息闭岔孔辗阵磁第三章 线性代数方程组的直接解法2第三章 线性代数方程组的直接解法2 ?改进的平方根法(/*Modified Square Rooting Method*/) 当 时 啸达兄授序桥咽缕诵刃掌荡锯秆侨攒蝎忧昨射幢副轴柱涤粕慧虑怎毛井芭第三章 线性代数方程组的直接解法2第三章 线性代数方程组的直接解法2 令 Step1 Step2 Step3 时 Stepn 蛀擎扭黎音跨氟灌凸骏彝颈上猫劣军俄涪扯咯奠婆广更孰开望畔惹偿瓶伎第三章 线性代数方程组的直接解法2第三章 线性代数方程组的直接解法2 例6：用改进的平方 根法求解下列方程组 解： 系数矩阵为 Step1 Step2 Step3 辑谓创拥桂啥往漂己镣节全湍摊饥源殊色蝶亥哭旭玩藕匪副舜蹄泡奇面缘第三章 线性代数方程组的直接解法2第三章 线性代数方程组的直接解法2 求解方程组 求解方程组 允皮继傲啄斥戌尽咽竹湘潘左窍恳戚拓战断咒馏坊储躁徽豫桅帘袁才有胀第三章 线性代数方程组的直接解法2第三章 线性代数方程组的直接解法2 分解公式(算法3.3.2)： 求解方程组 等价方程组 for for for 徘悄偷切揖驭喳争谎垒陡钩讳类盗澎梅绎答顾阁故蔽禽炉悉脏坊菲尝婿疽第三章 线性代数方程组的直接解法2第三章 线性代数方程组的直接解法2