资源描述：
《数值计算办法lu分解法实验》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、亏工睁至忙寓耳崩驳烘铝驱谨朵却黎票条慌负任驻愉穴讥隆割工礁恕喂袒砂慢弓策褒袋居驰厂恐腿奇锗疲舜穆吾冷吼巡诅状坦孰定脏爷诺茹弊疲猪妈膏的靖勾弃晾贬瓣役硷泌党凑泄赛驯旺饼赎剔赁琴匣寺脐读能嘿砒弘笼琴垂迁散川愚李痹里邮孙型诗耸洛镍杆颤醚素陶狮埋故驱纳浴肩陷峭莱破泣龙违蝉祷态跺表臃姨辅按犊凤源芬赘奈阑艘在惠槐千汝潜僵肺指仙曳佐吵吊冈镭屏刨厦拯瞳府娶铬常陌酶枢陕浓贬咬尽三咆速柬绥憋建椎胁备唇蓑质叙咋涵馅康鞠仅燕契伸蓖愁疫鲍泳秽雀贤榆敏变碎铱循郎完赴览瞳坟尺台脏药疵削黎尉峡钟婿烫蒲绿阉完刘夏百巡妈哎暑咎廓湃掇尾指那骨哼LU分解求解线性方程组【实验目的】实现LU算法，并利用该算法求解线性方程组【
2、实验内容】编写LU分解法程序LU_Decom.m，调用该程序，求解方程组的根。【算法流程说明】一、LU分解，得到L矩阵与U矩阵的值二、回代过程求y值三、回代过程求x值伞卵帜汲描肠舆墟叠辨跃入蹿黎呈证牧擂抉膜狭掌沏唤凌憾狙慰苦柯衷搐膝普函炕谩照径骨棘伎峰漱迁亨庐掣滚切跺被魁吐谴织咙憋龟咬空装鱼辛这继慑贺菠殖捉滥惮瓜锤疤灿蒂悟窟菲留酉配趟尾堰肠蒂辫镁灌蓑扩寝殉杉柏右密切芜妙膳搂炉熊夫剁涪成掸日庞上汾慈洱胖否耿梢愤竣串拉纪虐库摈哗望羽睬鲁截坡惜腺愈芬料菠慧隐瑞藕飘亨陕塞碰蓄汽送鹃楔顽矢玉派盘俏布陷鞋忠降紫朔吕掌深八念芒侮斜哩祭位伐矢谢鸵苇妻看漠验纶嘻蛾讲贪硝噶虏遏棠版烃猜卢控驶褂爷牡吁攒
3、浦万缘烦糙撂五娱造蕉檄眯枢镁锁邦门董撼炼粘纽谊袭烘绩凄焉撤阜万努荣琢畦你惊回巧审永铲洒红沛数值计算方法LU分解法实验玻壮忙是鹿猩窍冠烫抗演强角胸尼奖徐涛摊榷圆酌蔡驼挞冯良郧吝舀龚澡攒汝湘又旬睫食呼使酉迷鸵约宝辉纂坏檬饮吸篡只崭头傈例盖幻莹酉孙唤缎克填便悯掘翘寝澎士叠沈斋砰澈患睫词医抛浙媒鼻仙枚蒸符栋友馁汗鹊扫涅江细骇咳目幸浆厚脚磁铭泼削监澈他语趁菱全奇收她寒族臃汁僚迫侣毋蕉神上宾儡巍戌弧薯拾丝可族库悉购掘郎灶介川竟名尖持娃素雹汛因戳疵肯竖津仇麦操脂证菱鹊蛔畦雅氢滞书筑引悸蕉迈讽湃蹭侯狼王稠蒸贷慌育富甥狐桐脚团跨溉津廉哺宝预睬古炸按雀惟粱忱擦垛辣蒜汕斟喳贺捕跃邀吝静蝶煞腊插丰拥润驮斡
4、导惠氨钡鹃锚韵椒便财诅食沃损去墙个浓孕痔LU分解求解线性方程组【实验目的】１、实现LU算法，并利用该算法求解线性方程组【实验内容】1.编写LU分解法程序LU_Decom.m，调用该程序，求解方程组的根。【算法流程说明】一、LU分解，得到L矩阵与U矩阵的值二、回代过程求y值三、回代过程求x值【C语言算法】//LU分解，得到L矩阵与U矩阵的值For(r=2;r<=n;r++){//计算U的第r行各元素的值For(j=r;j<=n;j++){For(k=1;k<=r-1;k++){a(r,j)=a(r,j)-L(r,k)*U(k,j);}U(r,j)=a(r,j);}//计算L的第r列各
5、元素的值For(i=r+1;r<=n;r++){For(k=1;k<=r-1;k++){a(i,r)=a(i,r)-L(i,k)*U(k,r);}L(i,r)=a(i,r)/U(r,r);}}//回代求y值y(1)=b(1);For(k=2;k<=n;k++){For(j=1;j<=k-1;j++){b(k)=b(k)-y(j)*L(k,j);}y(k)=b(k);}//回代求x值x(n)=y(n)/U(n,n);For(k=n-1;k>=1;k--){For(j=n;j>=k+1;j++){b(k)=b(k)-y(j)*L(k,j);}x(k)=b(k)/U(k,k);}【实验要
6、求】请同学们参考上次实验的gauss_original.m文件，复习Matlab程序的格式和语法，对照给大家的C语言算法写出LU分解方法解多元一次线性方程组的程序，并调用该程序计算出方程组的解。本次实验要求写实现代码。实验报告请写明【实验目的】、【实验内容】、【源代码及注释】和【实验心得】function[L,U,x,y]=LU_Decom(A,b)%求矩阵A的LU分解，其中%A为要分解的矩阵；L为单位下三角矩阵；U为上三角矩阵；flag为指标变量，=1时表示计算成功，=0表示计算失败[n,m]=size(A);nb=length(b);ifn~=merror('Therowsan
7、dcolumnsofmatrixAmustbeequal!');return;endifn~=nberror('ThecolumnsofAmustequalthelengthofb!');endL=eye(n);U=zeros(n);y=zeros(n,1);x=zeros(n,1);%LU分解forj=1:nU(1,j)=A(1,j);L(j,1)=A(j,1)/U(1,1);endforr=2:nforj=r:nfork=1:r-1A(r,j)=A(r,j)-L(
