数学实验综合实验报告.docxVIP

、实验目的: 1、 初步认识迭代，体会迭代思想的重要性。 2、 通过在mathematica环境下编写程序，利用迭代的方法求解方程的根、线性 方程组的解、非线性方程组的解。 3、 了解分形的的基本特性及利用 mathematica编程生成分形图形的基本方法， 在欣赏由mathematica生成的美丽的分形图案的同时对分形几何这门学科有一 个直观的了解。从哲理的高度理解这门学科诞生的必然性， 激发读者探寻科学真 理的兴趣。 4、 从一个简单的二次函数的迭代出发，利用mathematica认识混沌现象及其所 蕴涵的规律。 5、 进一步熟悉Mathematic软件的使用，复习总结Mathematic在数学作图中的 应用，为便于研究数学图像问题提供方便，使我们从一个新的视角去理解数学问 题以及问题的实际意义。 6在学习和运用迭代法求解过程中，体会各种迭代方法在解决问题的收敛速度 上的异同点。 二、 实验的环境： 学校机房，mathematica4环境 三、 实验的基本理论和方法： 1、迭代(一)一方程求解 函数的迭代法思想： 给定实数域上光滑的实值函数f(x)以及初值x0定义数列 Xn 1 f (Xn)，n 0,1,23 ， (1) Xn , n 0,123,，称为f(x)的一个迭代序列。 方程求根 给定迭代函数f (x)以及初值X0利用(1)迭代得到数列Xn , n 0,123,. 如果数列收敛到某个x*，则有 X* f (X*). (2) 即X*是方程X f(X)的解。由此启发我们用如下的方法求方程 g(X) 0的近似 解。 将方程g(X) 0改写为等价的方程 X f(X)， ⑶ 然后选取一初值利用(1)做迭代。迭代数列Xn收敛的极限就是方程g(X) 0的 解。 为了使得迭代序列收敛并尽快收敛到方程 g(X) 0的某一解的条件是迭代 函数f(X)在解的附近的导数将的绝对值尽量小，因此迭代方程修订成 x h(X) f (x) (1 )x (4) 选取 使得山(x)|在解的附近尽量小.为此，我们可以令 h (x) f (x) 1 0, 得 1 1 f (x) 于是 h(x) xf(x) x h(x) x f(x) x f (x) 1 特别地，如果取f(x) g(x) X,则可得到迭代公式 Xn 1Xn Xn 1 Xn 響,n 0,1, g (Xn) Xn)T及b (bib，,bn)T均为n维列向量. 给定一个n元线性方程组 a〔1 x〔 a1n xn b1, a1 nx1 annxn bn 或写成矩阵的形式 Ax b, (2)线性方程组的数值解的迭代求解理论与矩阵理论 其中A (aj是n阶方阵,x (Xi，X2, 熟知，当矩阵A的行列式非零时，以上的方程组有唯一解.如何有效，快速 地寻求大型的线性方程组的数值解释科学工程计算中非常重要的任务 .而迭代法 常常是求解这些问题的有效方法之一。 用迭代法求解线性方程组的思想与上一小节介绍的方程求根的方法是类似 的。将方程组(7)改写成 TOC \o 1-5 \h \z x Mx f, (8) 其中M (mj)是n阶矩阵,f (fi， ,f n)T是n维列向量.任意给定初试向量x0, 由迭代 xn 1 Mxn f (9) 确定向量序列xn,n 0,1,.如果xn收敛到向量x*,则有 x* Mx* f, 则x*为方程组(7)的解. 假设矩阵A的对角元素aj 0,i 1,2,L n。令D diag(a11,L ,ann)，贝賊们 可以将方程(7)改写成 Dx (D A)x b x (I D 1A)x D 1b (10) 由上式即可确定一种迭代格式。 如果即将矩阵(I D 1A)分解为U L ,其中L,U分别为下三角阵与上三 角阵，则(10)可以进一步改成 (I L)x Ux D 1b 或 x (I L) 1Ux (I L) 1D 1b ( 11) 上式又可确定另一种迭代格式。 (3)非线性方程组的迭代求解理论 类似于单变量的方程组及线性方程组的求解，用迭代方法可以求更加复杂的 非线性方程组的解， 给定非线性方程组 TOC \o 1-5 \h \z f,X1丄 X) 0, LLLLLLLL, (12) fn(X1,L X) 0. 将它改写为等价的方程组 X1 g1(X1,L ,Xn), LLLLLLLL, Xn gn (x1, L , Xn )? 或 x g(x) (13) 其中，X为n维列向量X (X1,...,Xn)T,g ?(x), gn(x))T为n维列向量函 数，由上式即确定了一种迭代格式xn 1 g（xn）,n 0,1L . 由于非线性方程组可能有许多解（甚至有无穷多个解），因此对它的求借比 线性方程组的求解要面临更多的挑战。 2、迭代（二）一分形 分形几何一描述自然界的几何形态