迭代数学实验..docVIP

北京建筑工程学院 理学院 实验报告 课程名称 数学实验 实验名称 迭代（一）—方程求解 实验地点 基c425 日期 2015-7-1 姓名 周盈海 班级 信131 学号 201307010137 指导教师 高雁飞 成 绩 【实验目的】 （1）.利用迭代法求解方程 （2）.线性方程组的迭代求解 【实验要求】 1.（1）.在同一直角坐标系中画出及的图像。从图上观察，方程的解的个数。 （2）.利用迭代法，对作迭代，取极限，得到方程的解。 2. 利用迭代法，即，M是n阶矩阵，f是n维列向量，对作迭代，取极限，得到线性方程组的解。 【实验内容】（主要包含问题分析、计算过程、实验结果等，按课程要求完成） 问题的分析 1.（1）.在同一直角坐标系中画出及的图像。由理论可知，方程有3个解。 （2）. 取初值为，对作迭代，迭代序列应当收敛于0。但对初值取不同的点，对作迭代，迭代序列会出现不收敛的现象。 2. 给定实数域上光滑的实值函数f(x)以及初值x0定义数列 称为f(x)的一个迭代序列. 给定迭代函数f(x)以及一个初值x0利用上式迭代得到数列 如果数列收敛于某个,则有 .即是方程x=f(x)的解.由此启发我们用如下的方法求方程g(x)=0的近似解. 将方程g(x)=0改写为等价的方程 x=f(x),然后选取一初值利用上式做迭代.迭代数列收敛的极限就是方程g(x)=0的解. 给定一个n元线形方程组或写成矩阵的形式 Ax=b, 其中A=()是n阶方阵,x=均为n维列向量.任意给定初值向量，由迭代确定向量序列收敛到向量，则有即是方程组的解。 计算过程 1、 f[x_]:=2Sin[x] y[x_]:=x Plot[{f[t],y[t]},{t,-Pi,Pi}] （2） Iterate[f_, x0_, n_Integer] := Module[{t = {}, i, temp = x0}, AppendTo[t, temp]; For[i = 1, i = n, i++, temp = f[temp]; AppendTo[t, temp]];t] f[x_] := 2Sin[x] - x; Iterate[f, 0.5, 10] Iterate[f_, x0_, n_Integer] := Module[{t = {}, i, temp = x0}, AppendTo[t, temp]; For[i = 1, i = n, i++, temp = f[temp]; AppendTo[t, temp]];t] f[x_] := 2Sin[x] - x; Iterate[f, 10., 10] 2、 LSIterate[m_,f_List,f0_List,n_Integer]:= Module[ { i,var=f0,t=Table[{},{i,n}]}, For[i=1,i=n,i++,t[[i]]=var;var=m.var+f];t] m={{0.1,0},{0,0.2}};f={1,1};f0={0,0}; LSIterate[m,f,f0,20] LSIterate[m_,f_List,f0_List,n_Integer]:= Module[ { i,var=f0,t=Table[{},{i,n}]}, For[i=1,i=n,i++,t[[i]]=var;var=m.var+f];t] m={{0.1,0},{0,0.2}};f={2,2};f0={0,0}; LSIterate[m,f,f0,20] LSIterate[m_,f_List,f0_List,n_Integer]:= Module[ { i,var=f0,t=Table[{},{i,n}]}, For[i=1,i=n,i++,t[[i]]=var;var=m.var+f];t] m={{0.1,0},{0,0.2}};f={1,1};f0={2,2}; LSIterate[m,f,f0,20] LSIterate[m_,f_List,f0_List,n_Integer]:= Module[ { i,var=f0,t=Table[{},{i,n}]}, For[i=1,i=n,i++,t[[i]]=var;var=m.var+f];t] m={{0.2,0},{0,0.4}};f={1,1};f0={0,0}; LSIterate[m,f,f0,20] LSIterate[m_,f_List,f0_List,n_I