个人收集整理 仅供参考学习
1 / 8
%高斯赛德尔迭代解线
性方程组
a=input('请输入系数矩阵
a:') %输入系数矩阵 a
disp(a)
[m,n]=size(a)
if m~=n %若 a 不是方阵
则显示错误
disp('false')
end
de=det(a);
if de==0 %判断矩阵 a 是
否为奇异矩阵
disp('矩阵奇异')
end
b=input('请输入 b 矩阵:')
%输入矩阵 b
disp(b)
x=input('请输入迭代初始
值:')
disp(x)
k=x;
jindu=10e-7 %迭代精度
p=jindu+1;
while p>=jindu %迭代次
数控制条件
s=zeros(m,1);
te=zeros(m,1);
for i=1:m %迭代计算
for j=1:i-1
s(i)=s(i)+a(i,j)*k(j);
end
for j=i+1:n
te(i)=te(i)+a(i,j)*x(j);
x=k;
end
k(i)=(b(i)-s(i)-te(i))/a(i,i)
end
if
norm(k(i)-x(i),inf)>jindu
资源评论
SlumberingPerson
- 粉丝: 3
最新资源
- 清水池结构计算书.doc
- 墩台身施工工艺流程图.doc
- 住宅小区保温工程标准化施工工艺(PPT-图文并茂).ppt
- vcam-ios-main.zip
- (我设计)校园网络安全问题分析与对策.doc
- 湘潭万达广场项目C组团15#塔吊基础施工方案.doc
- 通风系统风道的设计.ppt
- 土建造价员考试模拟试题(含答案).doc
- S7-1200-控制器变更固件包含哪些附加信息？.pdf
- 幼儿园小班健康活动：自己穿袜子.doc
- 消防水泵扬程的计算.ppt
- 分部分项工程量清单.doc
- H-reflex-FCR.ppt
- 屋面防水工程-合成高分子卷材屋面防水技术交底模板.doc
- 基坑降水及工程案例分析.ppt
- 物流信息化解决方案怎么做ppt.pptx
资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！
点击此处反馈
