clear;
clc;
ag=0;
disp('请输入');
while ag==0
A=input('A=');
b=input('b=');
[m,n]=size(A);
if m~=n
disp('矩阵不是方阵，无法 Doolittle 分解,请重新输入 A，b');
else
C=zeros(n,1);
for i=1:n
C(i,1)=det(A(1:i,1:i));
end
if C(1:n,1)~=0
ag=1;
else
disp('无法 Doolittle 分解,请重新输入 A，b');
end
end
end
L=eye(n);
U=zeros(n);
Y=zeros(n,1);
X=zeros(n,1);
for j=1:n
U(1,j)=A(1,j);
end
for k=2:n
L(k,1)=A(k,1)/U(1,1);
end
for k=2:n
for j=k:n
U(k,j)=A(k,j)-sum(L(k,1:k-1)'.*U(1:k-1,j));
end
for i=k+1:n
L(i,k)=U(k,k)\(A(i,k)-sum(L(i,1:k-1)'.*U(1:k-1,k)));
end
end
Y(1)=b(1);
for i=2:n
Y(i)=b(i)-sum(L(i,1:i-1)'.*Y(1:i-1,1));
end
qq_43460878
- 粉丝: 0
- 资源: 1
我的内容管理
收起
- 我的资源
快来上传第一个资源
我的收益 登录查看自己的收益
我的积分
登录查看自己的积分
我的C币
登录后查看C币余额
我的收藏 
我的下载 
下载帮助
会员权益专享
最新资源
- zigbee-cluster-library-specification
- JSBSim Reference Manual
- c++校园超市商品信息管理系统课程设计说明书(含源代码) (2).pdf
- 建筑供配电系统相关课件.pptx
- 企业管理规章制度及管理模式.doc
- vb打开摄像头.doc
- 云计算-可信计算中认证协议改进方案.pdf
- [详细完整版]单片机编程4.ppt
- c语言常用算法.pdf
- c++经典程序代码大全.pdf
- 单片机数字时钟资料.doc
- 11项目管理前沿1.0.pptx
- 基于ssm的“魅力”繁峙宣传网站的设计与实现论文.doc
- 智慧交通综合解决方案.pptx
- 建筑防潮设计-PowerPointPresentati.pptx
- SPC统计过程控制程序.pptx
资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！
点击此处反馈
