资源描述：
《数值计算方法I实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、重庆交通大学学生实验报告实验课名称计算方法A开课实验室数学实验室学院年级专业班级学生姓名学号开课时间至学年第学期评分细则评分报告表述的清晰程度和完整性（20分）程序设计的正确性（40分）实验结果的分析（30分）实验方法的创新性（10分）总成绩教师签名实验一误差分析实验目的：掌握Newton迭代法，了解拟Newton迭代法，培养编程能力和用迭代法解非线性方程(组)的实践技能。实验环境：MATLAB实验1.1：用二分法求非线性方程f(x)=x3+x2−3x−3=0在1.5附近的根。实验代码：bisecMethod：functionbisecMethod(f_name,a
2、,b)fprintf('二分法求非线性方程的根');tolerance=1e-6;it_limit=30;fprintf('It.acbfa=f(a)');fprintf('fb=f(b)abs(fb-fa)');it=0;Y_a=feval(f_name,a);Y_b=feval(f_name,b);if(Y_a*Y_b>0)fprintf('停止，因为f(a)f(b)>0');fprintf('重新给出端点a或者b.');elsewhile1it=it+1;c=(a+b)/2;Y_c=feval(f_name,c);fprintf
3、('%3.0f%10.6f%10.6f',it,a,c);fprintf('%10.6f%10.6f%10.6f%12.3e',b,Y_a,Y_b,abs(Y_b-Y_a));%fprintf('%12.3e',abs(Y_b-Y_a))if(abs(b-a)<=tolerance)fprintf('结果满足允许误差tolerance=1e-6.');breakendif(it>it_limit)fprintf('迭代超过限定次数.');breakendif(Y_a*Y_c<=0)b=c;Y_b=Y_c;elsea=c;Y_a=Y_c;endendf
4、printf('最终结果：x*=%12.6f',c);endchap4_fun1.mfunctionf=chap1_fun1(x)f=x^3+x^2-3*x-3;在命令窗口调用M文件bisecMethod(f_name’,a,b)>>bisecMethod('chap4_fun1',1,2)二分法求非线性方程的根It.acbfa=f(a)fb=f(b)abs(fb-fa)11.1.2.-4.3.7.000e+0021.1.2.-1.3.4.875e+0031.1.1.-1.0.2.047e+0041.1.1.-0.0.1.115e+0051.1.1.-0.0.5
5、.813e-0161.1.1.-0.0.2.967e-0171.1.1.-0.0.1.468e-0181.1.1.-0.0.7.379e-0291.1.1.-0.0.3.699e-02101.1.1.-0.0.1.847e-02111.1.1.-0.0.9.241e-03121.1.1.-0.0.4.622e-03131.1.1.-0.0.2.311e-03141.1.1.-0.0.1.155e-03151.1.1.-0.0.5.776e-04161.1.1.-0.0.2.888e-04171.1.1.-0.0.1.444e-04181.1.1.-0.0.7.221
6、e-05191.1.1.-0.0.3.610e-05201.1.1.-0.0.1.805e-05211.1.1.-0.0.9.026e-06结果满足允许误差tolerance=1e-6.最终结果：x*=1.实验1.2：公元1225年，Lenardo宣布他求得方x3+2x2+10x−20=0实验代码：functionx=Newton(fname,dname,x0,e,N)ifnargin<5,N=500;endifnargin<4,e=1e-4;endx=x0;x0=x+2*e;k=0;fprintf('It.no=%2dx%[2d]=%12.9f',k,k,x)
7、whileabs(x0-x)>e&k>fun=inline('x^3+2*x^2+10*x-20');>>dfun=inline('3*x^2+4*x+10');>>Newton(fun,dfun,1.5,0.5e-6);It.no=0xIt.no=1x
