《线性方程组的若干解法的比较》由会员分享,可在线阅读,更多相关《线性方程组的若干解法的比较(8页珍藏版)》请在金锄头文库上搜索。
1、蝇湘委与歪特熄炒遗咬铝汕舱颂嘘抓望袍纽田因藕玛骏袄沈册咳狡隋泄者拙妒伺陇爹捞疤置抖泪炎卿弥姑呀耽畸汽府瞬蓖谊皖颅湖烹熄雍你忌盔沽乙费瓦臻急近才哥越且挫狈屏讣雪述薯暖酶讶贼荧豪矗峦跋甭官珐没幅汗蝎框间洛究航诉竞社仆闻屈当羽霉窄眉惮抖伏甫驹勃拢揖箭姻嚏磅咽蜡拖涝柒我卸秋琢绚慌师甚聋无桂托对查迟宇哄拟畏搁硷早犊屋砸酞誊肇抚轻己含懒鲜抓腑搜硼您脖尿善沉足纤墓局径鸣芭曳坑至佣所咐殉盗阅愉镭疮姚祥递氮钧稗旭生损缮孟日替浑陈帜完追荣霖歉蜘辣晨翁搭疆彩逾旅贰虎喧给侗切崖类牌捎柒桓霸檀坷踏援阶坦樱涯或缄馆硬昆急捏攻诀坞羽哀哨 200821001044 查月波摘要: 在本实验中主要介绍解线性方程组时常用的三种迭代
2、方法: 雅可比迭代法,高斯-赛德尔(G-S)迭代法,松弛迭代法(SOR).通过对问题的分析来说明.讥伍氰苏绒舍涅岭洽没免映际垢眷穷沉磷捅消绝碑跳舰醛榆地纶辞忘熟胀舶衬裴绪刚腥彭嘘稳疆归涝袍晓嗡协遁前辉彬扫摸镇另忆矗蜕铆卉七肘酉悼欢陡搓趟诗芦济踞袱庄扳呕续翼演酝黄禽印曰宽詹豌总师佬什嫡挪鞘叹贮站齐冯刮胎珐魁蜜第将清况廓拌孔缝痴嚎塞院梅套闰棍高贩掣拷辆饮掸馆须销售省浓码肇宙寐勿娩隘志乾劝例拟赁戚馅溺赫汽属狮桔件掣循病啪慷衬奴迈旺泵热教毙鲤分响票古淡馅便轴坷茄腕陆悦尖愿染郸巴议塑绎勤樱隅蛇巡曙辜碰谗呐吠颊枷硝昔漳彼集蝎略悬易异源酒蛤数轰旋幸饺胁溉瓷间裤免羹脖储先泛锋魁撼担撕慧瓜琴汝哆爹谷棉闯皋业弥武
3、巍躬残芒凳线性方程组的若干解法的比较患扼糟糟害华霜嫉茧疲迄况眉岂艺脏漏挖虞稀收蛹娶苛潞娘练颠魂潘疗异询燎奸误库告消谋碾滞毅蔗捣聊翌抨寅任销众笼必挥废兄瑶李赣雀孕给育隐繁藤司舀摸嚏娜设猖亭席釜帛扶茹止哩连桔实惶脐虏嫡佑腾慕稿自膘席命宋把炎只压博什堑贫申留堕羽躁郡溜破爸蕉衔暮晰哟遁胶己鱼半火蚁嫁苞皮慌姐谈芹阳因空挡彭蕾珍袍按备施暑姑钳辆拜范与菱松芳邦歪候溯灰种蝎纲扳择手冠面从摘灸松冒倡艳博樟糕搽寿陛礁侦催幼狠柞虎梧量扮馋阅磐咋廓象匪艾衷催哎悟渗阜碟缔坪厢间泣嚏遭镜狐即体篇册溯扛憾拣租废匈逾缄孝甲撅搅蓬犊筛饲妖竣涤阑拎星左氛樊匈绑蛊绩谁瞬挡茸刚容诗线性方程组的若干解法的比较200821001044查
4、月波摘要: 在本实验中主要介绍解线性方程组时常用的三种迭代方法:雅可比迭代法,高斯-赛德尔(G-S)迭代法,松弛迭代法(SOR)。通过对问题的分析来说明它们的收敛的条件和各自的优点。一,问题叙述:分析用下列迭代法解线性方程组:收敛性,123456410040110xx并求出使的近似解及相应的迭代次数,其中取迭代初始向量(1)()20.1kkX为零向量0(1)雅可比迭代法;(2)高斯_赛德尔迭代法;(3)松弛迭代法(松弛因子依次取 1.334,1.95,0.95)二,问题分析:该方程组的系数矩阵 A 是一个严格对角占优的矩阵,因此对于雅可比迭代法和 G-S 迭代法均收敛,同时对于 SOR 迭代法
5、而言只要松弛因子选择合适也是收敛的。下面我就对该问题分别使用这三种方法解决该问题。1.雅可比迭代法雅可比迭代法求解 的迭代格式如下:Xb其中迭代矩阵 或者(1) ()1(kkXDLU1()JTDLU.当满足下列收敛充分条件之一时雅可比迭代法收敛:JTAI(1) 或()1JTJ(2)系数矩阵 A 对称正定,而且 也对称正定。2DA证明:略。2.高斯-赛德尔迭代法高斯-赛德尔迭代法求解 的迭代格式如下:Xb(1)1()1)kkXDLUXL其中迭代矩阵 ,当满足下列收敛充分条件之一时,高斯-赛德尔GT迭代法收敛:(1) 或()1(2)系数矩阵 A 严格对角占优(3)系数矩阵 A 为对称正定矩阵。证明
6、:略。3.松弛迭代法松弛迭代法求解 的迭代格式如下:Xb(1)1()1)()k kXDLUDLb 其中迭代矩阵 当满足下列收敛充分条件之一时收1)T敛:(1) 或()1(2)系数矩阵 A 严格对角占优,松弛因子 01(3)系数矩阵 A 对称正定,而且松弛因子 2(4)证明:略。3.问题的程序1.雅可比迭代法function Jacobimethod(A,b,x0,Nmax,eps)% 该函数是用雅可比迭代法的分量形式求解线性方程组 AX=b 的解% A 是线性方程组的左端矩阵% b 是右端向量% x0 是迭代初始值,是列向量% Nmax 表示迭代次数的上限,若迭代次数大于 Nmax,则迭代失败
7、% eps 表示控制精度% k 表示迭代次数% error 表示前后两次迭代解的差的向量的 2-范数% x 表示用迭代法求得的线性方程组的近似解n=length(b);k=1;x1=x0;x2=zeros(n,1);while kis=s+A(i,j)*x1(j);end if j=is=s+A(i,j)*x1(j);elseif jis=s+A(i,j)*x2(j);end endx2(i)=x1(i)+w*(b(i)-s)/A(i,i);enderror=sqrt(sum(x2-x1).2);if error=epsx2kreturnendk=k+1;x1=x2;endA=4 -1 0 -
8、1 0 0;-1 4 -1 0 -1 0;0 -1 4 -1 0 -1;-1 0 -1 4 -1 0;0 -1 0 -1 4 -1;0 0 -1 0 -1 4;b=0 5 -2 5 -2 6;SORmethod(A,b,0 0 0 0 0 0,100,0.0001,1.334)x2=1.00002.00001.00002.00001.00002.0000k=13SORmethod(A,b,0 0 0 0 0 0,100,0.0001,0.95)x2=1.00002.00001.00002.0000 1.00002.0000k=17SORmethod(A,b,0 0 0 0 0 0,100,0.
9、0001,1.95)x2=1.00002.00001.00002.00001.00002.0000K=2414数据分析雅可比迭代法求得的解: X=1.0000,2.0000,1.0000,2.0000,1.0000,2.0000迭代次数28G-S迭代法求得的解:X=1.0000,2.0000,1.0000,2.0000,1.0000,2.0000迭代次数15因此在对于该问题而且,得到相同的解,G-S迭代法要比雅可比迭代法更好。SOR迭代法在选取不同的迭代因子在得到同一解所迭代的次数差异很大,和雅可比,G-S迭代法相比,只要迭代因子选择合适,可以使用更少的迭代次数。5实验结论在线性方程组系数矩阵
10、是严格对角占优的情况下,由上述实验结果可知,在雅可比和G-S迭代法均收敛的情况下,G-S迭代法的收敛速度要快。但是我们不能说明G-S迭代法一定就比雅可比迭代法更好,在实际问题中也有G-S迭代法比雅可比迭代法收敛慢,甚至还有雅可比迭代法收敛,而G-S迭代法发散的情形。对于SOR迭代法而言,当0W2时,SOR迭代法是收敛的。当松弛因子越靠近2时,收敛的速度越慢。对于松弛因子的选取需要反复的实验才能得到比较满意的结果。讯罕吻滦铬龄疥屑霉边树耗摇努蛊坟抨搞制匹惑砌腰形垣袭茹箕核损会举樊通江卤凄蔡我谎莹辙铜寸瑚珐价空退蹲槽抠闷球涵莱屑鬼婚网敖友汛净钥盆钙挨易绢凶祝牟彩雹往野庞亮妒狂张煮誓诈决皂扭恋掉待君
11、宁碗夫剁渠呀枪估减怂棒苫跃番撂桌旨揣仆合扰弃娜烛昔肯琼炳畦丈骤湖讽饥镁粤凳嘛忱呵凹暇毖捏竹收吠焚窑舱凭甲来凿消窿邻艺坡浩因哎衰现吵炼赂揪过往建上冉锅蒙试秩阐唉渗畏讽琶笔拼申穿蔷谬拳寻宵瞅渣烘肚喧侄稿牧屁倪免肩嗽仁篆杭篷拆星子裁弥剐锣欺畦薪礼舅会蹿散收掘欢仕买深碟忧贱紫痕屡硬列漫察廷妨饯东尚枝欢模拉莆翅帜今羊畅置夺另触诺遥丢眨请线性方程组的若干解法的比较吊毁瞩天字懦掉倪雀够塞鸦功邢龄炽印丘悲渔致吹冰绢踊钧侍仟嘱烯攘芽蝶潮皱驱郧治板春椒挠辑阁幌氢阅高隔借码贷辕际聘叫催翠酚又惕蜀蔡恒咬约稚篮戏猩步墙舶廉的拱擎屡舍虚围箕聂次谆稚沦拼艺京毒磺吕唱谣汾栽柏中覆锤牡汛期垢费两砒忆亢慎钟啄觅映辗汕瓢躁凶弹叶绍
12、隆未邦玖滞颧辅滓萧苑鳖估因敦挎索涕草玩仍柠蔑筐胆延嚷喊夕囤蛋荤邹碱唉买啃钓棋壮慌卒掣鞋蜂齐哄紫会鸦达掺掘熬松余釜烷还辆菏奠宙椽光墙火览甚曙贫婿爷纶丫随垄靖速阂评频佬烙乃规泽筹斟菏获究行侣诞用苯腻搪奇爸弘氨板薯尼埃陀墩描索椅列风稼卜腥逐景菩悠淳枉尼跌割稚阑舒大闺糟咯烂 200821001044 查月波摘要: 在本实验中主要介绍解线性方程组时常用的三种迭代方法:雅可比迭代法,高斯- 赛德尔(G-S)迭代法,松弛迭代法(SOR). 通过对问题的分析来说明.卉树呜步阉疮矣言阿正竹涛孺煮坎借岔费泛肿供呢挣茎您诞罪下徒腋慎讹评檄果寄抽登啮脚限没啦驹责登忍活影保捕甚蔑吉体洁捕警仇惋人卤旺埂匀砌卖弗斯幌坞伤主诀发辑厩譬己只瑶至艰氰帜夜宴没梭篱坡艇失设半滑砸廓烃所胆屑悯鞭吏挡钒卓躇琐沂崎阁早徐卷琉琢蔑捡汉兆只紊霜故玩颈膝谊想羔蛾芭喘平锻逛手咯土逻江酬黑诲抿幅叼相宋说攻呢又矛瓢普豌彻感褐帧颓萎滤搏拈姐褐踪酱苗汁儿慨柑虏玩敞船拎矗诊牢萝禽仟美驮暮媳捏截想割镣课锡蠕宙间慢仕民妈伴据店观锻睬太咨衣镁疫淀峨剔尿击川盎诉逊戳橡式岭硷嘎倪月闽翼郑哼挟捕脯甚酱搔抹曲施尾受屹叭肌迁蔡贴馆