Matlab线性方程组的迭代解法Gauss-Seidel迭代法.doc
- 文档编号:1399158
- 上传时间:2022-10-22
- 格式:DOC
- 页数:1
- 大小:22KB
Matlab线性方程组的迭代解法Gauss-Seidel迭代法.doc
《Matlab线性方程组的迭代解法Gauss-Seidel迭代法.doc》由会员分享,可在线阅读,更多相关《Matlab线性方程组的迭代解法Gauss-Seidel迭代法.doc(1页珍藏版)》请在冰豆网上搜索。
Matlab线性方程组的迭代解法Gauss-Seidel迭代法实验报告
1.熟悉Gauss-Seidel迭代法,并编写Matlab程序
function[y,n]=gauss_seidel(A,b,x0,eps)
%gaussseideliterativemethodtosolveAx=b
ifnargin==3
eps=1e-6;
elseifnargin<3
error;
return
end
D=diag(diag(A));
L=-tril(A,-1);
U=-triu(A,1);
G=(D-L)\L;
f=(D-L)\b;
y=G*x0+f;
n=1;
whilenorm(y-x0)>=eps
x0=y;
y=G*x0+f;
n=n+1;
end
2.
用Gauss-Seidel迭代法求解线性方程组。
输入:
>>A=[430;33-1;0-14];
>>b=[24;30;-24];
>>x0=[0;0;0];
>>[y,n]=gauss_seidel(A,b,x0)
输出:
y=
6
-2
-7
n=
4
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Matlab 线性方程组 解法 Gauss Seidel 迭代法
![提示](https://static.bdocx.com/images/bang_tan.gif)