电脑桌面
添加小米粒文库到电脑桌面
安装后可以在桌面快捷访问

数值分析(第五版)计算实习题第五章作业VIP免费

数值分析(第五版)计算实习题第五章作业_第1页
1/13
数值分析(第五版)计算实习题第五章作业_第2页
2/13
…………………………………………………………最最最最最最最最…………………………………………………… …………………………………………………………最最最最最最最最……………………………………………………1 数值分析第五章 第一题: LU分解法: 建立m 文件 function h1=zhijieLU(A,b) %h1 各阶主子式的行列式值 [n n]=size(A);RA=rank(A); if RA~=n disp( '请注意:因为A的n阶行列式h1 等于零,所以A不能进行LU分解。A的秩RA如下:') RA,h1=det(A); return end if RA==n for p=1:n h(p)=det(A(1:p,1:p)); end h1=h(1:n); for i=1:n if h(1,i)==0 disp( '请注意:因为A的r阶主子式等于零,所以A不能进行LU分解。A的秩RA和各阶顺序主子式h1 依次如下:') h1;RA return end end if h(1,i)~=0 disp( '请注意:因为A的r阶主子式都不等于零,所以A能进行LU分解。A的秩RA和各阶顺序主子式h1 依次如下:') for j=1:n U(1,j)=A(1,j); end for k=2:n for i=2:n for j=2:n L(1,1)=1;L(i,i)=1; if i>j L(1,1)=1;L(2,1)=A(2,1)/U(1,1);L(i,1)=A(i,1)/U(1,1); L(i,k)=(A(i,k)-L(i,1:k-1)*U(1:k-1,k))/U(k,k); else U(k,j)=A(k,j)-L(k,1:k-1)*U(1:k-1,j); end end end end h1;RA,U,L, X=inv(U)*inv(L)*b 数 值 分 析 (第 五 版 )计 算 实 习 题 第 五 章 作 业 --第 1页数 值 分 析 (第 五 版 )计 算 实 习 题 第 五 章 作 业 --第 1页…………………………………………………………最最最最最最最最…………………………………………………… …………………………………………………………最最最最最最最最……………………………………………………2 end end 输入: >> A=[10 -7 0 1;-3 2.099999 6 2;5 -1 5 -1;2 1 0 2]; >> b=[8;5.900001;5;1]; >> h1=zhijieLU(A,b) 输出: 请注意:因为A的r 阶主子式都不等于零,所以A能进行LU分解。A的秩 RA和各阶顺序主子式h1 依次如下: RA = 4 U = 10.0000 -7.0000 0 1.0000 0 2.1000 6.0000 2.3000 0 0 -2.1429 -4.2381 0 -0.0000 0 12.7333 L = 1.0000 0 0 0 -0.3000 1.0000 0 0 0.5000 1.1905 1.0000 -0.0000 0.2000 1.1429 3.2000 1.0000 X = -0.2749 -1.3298 1.2969 1.4398 h1 = 10.0000 -0.0000 -150.0001 -...

1、当您付费下载文档后,您只拥有了使用权限,并不意味着购买了版权,文档只能用于自身使用,不得用于其他商业用途(如 [转卖]进行直接盈利或[编辑后售卖]进行间接盈利)。
2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。
3、如文档内容存在违规,或者侵犯商业秘密、侵犯著作权等,请点击“违规举报”。

碎片内容

数值分析(第五版)计算实习题第五章作业

您可能关注的文档

确认删除?
VIP
微信客服
  • 扫码咨询
会员Q群
  • 会员专属群点击这里加入QQ群
客服邮箱
回到顶部