南通大学数学实验 非线性迭代.doc

南通大学实验报告 非线性迭代 学院： 理学院 班级： 数师153班 学号： 1502012072 姓名： 顾阳 一、实验解读 迭代是数学研究中的一个非常重要的工具，通过函数或向量函数由初始结点生成迭代结点列，也可通过函数或向量函数由初值（向量）生成迭代数列或向量列。 蛛网图也是一个有用的数学工具，可以帮助理解通过一元函数由初值生成的迭代数列的敛散性，也帮助理解平衡点（两平面曲线交点）的稳定性。 本实验在Mathematica平台上首先利用蛛网图和迭代数列研究不动点的类型；其次通过蛛网图和迭代数列研究Logistic映射，探索周期点的性质、认识混沌现象；第三通过迭代数列或向量列求解方程（组）而寻求有效的求解方法；最后，利用结点迭代探索分形的性质。 二、实验计划 1.迭代序列与不动点 1.1程序 对函数的迭代过程，我们可以用几何图象来直观地显示它——“蜘蛛网”。运行下列Mathematica程序： Clear[f] f[x_] := (25*x - 85)/(x + 3); Solve[f[x]==x , x] g1=Plot[f[x], {x, -10, 20}, PlotStyle - RGBColor[1, 0, 0], DisplayFunction - Identity]; g2=Plot[x, {x, -10, 10}, PlotStyle - RGBColor[0, 1, 0], DisplayFunction - Identity]; x0=5.5; r = {}; r0=Graphics[{RGBColor[0, 0, 1], Line[{{x0, 0}, {x0, x0}}]}]; For[i = 1, i = 100, i++, r=Append[r, Graphics[{RGBColor[0, 0, 1], Line[{{x0, x0}, {x0, f[x0]}, {f[x0], f[x0]}}]}]];x0=f[x0] ]; Show[g1, g2, r, r0, PlotRange - {-1, 20}, DisplayFunction - $DisplayFunction] x[0]=x0; x[i_]:=f[x[i-1]]; t=Table[x[i],{i,1,10}]//N ListPlot[t] 1.2实验思路 首先对函数研究不动点，需要 （1）对Plot中{x,-10,20}可改为{x,-50,50}；对PlotRange中{ -1,20}可改为{-50,50}； （2）x0=5.5中5.5分别改为-30，-20，-5，-1，0，1，5，5.1，5.001，16,17,18，20，30；（为了图像清楚，将t=Table[x[i],{i,1,10}]//N中10分别改为20） （3）对t=Table[x[i],{i,1,20}]//N中20分别改为100，200，500，1000； （4）对i=100中100分别改为200，500，1000。 运行程序后观察蛛网图与散点图！一看数列是否收敛，如收敛，极限是多少？收敛速度是快是慢？二看蛛网图中的轨道是否趋于平衡点，与平衡点处曲线的斜率有没有关系，三看初值对结果有没有影响。 其次，分别就，等函数迭代序列，观察蛛网图中的轨道是否趋于平衡点和序列的收敛性。 Logistic映射与混沌 2.1程序 从形如的二次函数开始做迭代 这里，是一个参数。对不同的系统地观察迭代的行为。 Mathematica程序： IterGeo[a_, x0_] := Module[ {p1, p2, i, pointlist = {}, v= x0, fv= a*x0*(1 - x0)}, p1=Plot[ {a*x*(1 - x), x}, {x, 0, 1}, DisplayFunction - Identity]; AppendTo[pointlist, {x0, 0}]; For[i = 1, i 20, i++, AppendTo[pointlist, {v, fv}];