迭代——方程求解(mathematica数学实验报告)

姓名 @@@ 学院 @@@@@ 班级 ############## 学号 ############ 评分 实验题目 实验目的： 迭代——方程求解 1、初步了解迭代的概念，明确迭代思想在数学研究中的地位； 2、通过实验，由Mathematica 4.0软件演示运用迭代思想进行方程求解的具体过程，亲身体验迭代算法在计算机及数学学科中的重要地位； 3、在学习和运用迭代法求解问题时，了解各种具体迭代方法在解决问题的收敛速度上的异同； 4、进一步熟悉Mathematica4.0软件的使用，复习Mathematica在数学作图、计算中的应用； 5、通过上机来增强自己的动 手能力及实践创新能力。 实验环境： 学校机房，Mathematica4.0软件 实验基本理论和方法： 1、Mathematica中常用的函数及函数调用的方法： 函数的迭代法思想： 给定迭代函数f(x)以及一个初值x0利用（1）迭代得到数列xn,n=0,1,2…，如果数列xn收敛与某个x，则有 * x?f(x) 即x是方程x***?f(x)的解.由此启发我们用如下的方法求方程g(x)?0的近似解. 将方程g(x)?0改写为等价的方程 x?f(x) 然后选取一初值利用（1）做迭代.迭代数列xn收敛的极限就是方程g(x)?0的解. 2、对迭代的概念、特点及方法的掌握，利用迭代思想求方程的解； 3、对线性方程组的解的理论及矩阵理论的掌握和应用。 实验内容和步骤： 一、给定实数域上的实值函数f(x)以及初值x0定义函数列 xn?1?f(xn)n?0,1,... （1） xn,n?0,1?称为f(x)的一个迭代序列. 1、迭代法 函数的迭代是数学研究中的一个非常重要的思想工具，哪怕是对一个相当简单的函数进行迭代，都可以产生异常复杂的行为，并由此而衍生了一些崭新的学科分支，如分行与混沌.同时，迭代在各种数值计算算法以及其它学科领域的诸多算法中处于核心的地位. 本实验的基本理论是分形几何学.分形几何把自然形态看作是具有无限嵌套层次的精细结构，并且在不同尺度下保持某种相似的属性，于是在在简单的过程中可以得到描述复杂的自然形态的有效方法.本实验的基本方法是简单的迭代。 2、线性方程组的解的理论及矩阵理论: 给定一个n元线性方程组 ?a11x1?....?a1nxn?b1?.....................?.........., ?ax?....?ax?bnnn1?n11或写成矩阵的形式 Ax?b 其中A?(aij)是n阶方阵，x?(x1,...,xn) 及b?(b1,...,bn)均为n阶列向量. 熟知，当矩阵A的行列式非零时，以上的方程组的唯一解.如何有效，快速地寻求大型的线性方程组的数值解释科学工程计算中非常重要的任务.而迭代法常常是求解这些问题的有效方法之一. 用迭代法求解线性方程组的思想与上一小节介绍的方程求根的方法是类似的.假设我们可以将方程组改写成 x?Mx?f T其中M?(mij)是n阶矩阵，f?(f1,...,fn)是n阶列向量.任意给定初始向量x0，由迭TT代 xn?1?Mxn?f 确定向量序列xn，n=1,2,….如果xn收敛到向量x，则有 ** x?Mx?f 即x是方程组得解. x?x/2练习1：设 f?x??，并且给定初值 x0?1，做10次迭代得到序列 xn， n?0,1,...，2**程序运行如下： 练习2：利用迭代公式 xn?1?xn?g(x),n?0,1,... g?(x)得到g?x??x^3?2的迭代序列，其中x0?1，n?10，程序运行如下： 练习3：对给定的矩阵M，数组f和初始向量x0，由迭代公式xn?1?Mxn?f得到的迭代序列如下：