matlab程序设计实践解剖.doc

MATLAB程序设计实践 1、编程实现以下科学计算算法，并举一例应用之。（参考书籍《精 通ＭＡＬＡＢ科学计算》，王正林等著，电子工业出版社，２００９ 年） “里查森迭代法线性方程组求解” 解： 算法说明： 里查森迭代法是最简单的迭代法，它的迭代公式为：xk+1=（I-A）*xk+b；在MATLAB中编程实现的里查森迭代法函数为：richason。 功能：用里查森迭代法求线性方程组 调用格式：[x,n]=richason(A,b,x0,eps,M) 其中，A为线性方程组的系数矩阵； b为线性方程组的常数向量； x0为迭代初始向量； eps为解的精度控制（此参数可选）； M为迭代步数控制（此参数可选）； x为线性方程组的解； n为求出所需精度的解实际的迭代步数。 里查森迭代法的MATLAB程序代码如下： function [x,n] = richason(A,b,x0,eps,M) %采用里查森迭代法求线性方程组Ax=b的解 %线性方程组的系数矩阵：A %线性方程组的常数向量：b %迭代初始向量：x0 %解的精度控制：eps %迭代步数控制：M %线性方程组的解：x %求出所需精度的解实际的迭代步数：n if(nargin==3) eps=1.0e-6; %eps表示迭代精度 M=200;