山东大学数值计算试验报告9

1、山东大学软件学院数值计算课程实验报告学号：姓名：班级：5班实验题目：1、对教材例5.8的方程，用不动点迭代法，列出迭代过程，画图显示迭代过程（类似图5.5,不必画箭头）。迭代最大次数N=30,误差限E=10-2、用库函数求解（选做）实验学时：2实验日期：2017年12月5日实验目的：非线性方程求解硬件环境:PCW软件环境:MATLAB实验步骤与内容:等价不动点函数为g(x)=xA2-2代码如下：x=linspace(0,4);y=x;y1=x.A2-2;plot(x,y)%ffl出y=xholdonplot(x,y1)%®出g(x)holdonx1=2.1;P=；q=;fori=1:

2、3履代三次a=x1.A2-2;P=Px1;q=qa;P=Pa;q=qa;x1=a;endplot（p,q,'-bo'）holdon运行结果如下：4kFigure1X文件辑辑国直看国）mx（I）工具田桌面©国口他）帝助（H）Qdd、入踮0居0口日1412108602468101214等价不动点函数为g(x)=(x+2)A(1/2)代码如下：x=linspace(0,4);y=x;y1=(x+2).A(1/2);plot(x,y)%®出y=xholdonplot(x,y1)%®出g(x)holdonx1=1;p=；q=;fori=1:3履代三次a=(x

3、1+2).A(1/2);P=Pxl；q=qa；P=Pa；q=qa;x1=a;endplot(p,q,'-bo')holdon运行结果：，Figure1X文件由案第直看MKA(I)IRffi桌面©国口凶)帚助(H)3日UiR二，-J®0乂-息口明等价不动点函数为g(x)=1+2/x代码如下：x=linspace(0.1,4);y=x;y1=1+2./x;Plot(x,y)holdon%®出y=xplot(x,y1)%®出g(x)holdonx1=1;P=；q=;fori=1:10%S代十次a=1+2/x1;P=Px1;q=qa;P=Pa;q

4、=qa;x1=a;endplot(p,q,'-bo')holdon运行结果：KFigure1X文件(日编信查看M瓶入0)工具(D点面窗口也)帮助但)1.4等价不动点函数为g(x)=(xA2+2)/(2x-1)代码如下：x=linspace(0.9,4);y=x;y1=(x.A2+2)./(2*x-1);plot(x,y)%ffl出y=xholdonplot(x,y1)%®出g(x)holdonx1=1;P=；q=;fori=1:6%迭代六次a=(x1A2+2)/(2*x1-1);P=Px1;q=qa;P=Pa;q=qa;x1=a;endPlot(P,q,'-bo')holdon运行结果:AFigure1文件由强国堂盲N)插入5工具CD桌面电)帮助(H),0)运行结果:心口Rir.gcriteriadetails)ansfsolvecompletedbecausethevectoroffunctLonvaluesIsnearzeroasmeasuredbythedefaultvalueofthefunctisntdlEran七丘,andthe口reblrmappa:rsregularasmeasuredb