北京科技大学计算方法大作业概要.docVIP

PAGE PAGE 1 计算方法大作业 机械电子工程系 老师：廖福成 注：本文本只有程序题，证明题全部在手写已交到理化楼204了。 2. 证明方程 在上有一实根，并用二分法求这个根。要求。请给出程序和运行结果。 证明： 设f(x)=x3-x-1 则f(1)= -1， f(2)= 5，f(1)*f(2)= -50 因此，方程在[1,2]上必有一实根。 二分法求解程序： %预先定义homework2.m文件如下： function lc=homework2(x) lc=x^3-x-1; 在MALAB窗口运行： clear a=1;b=2;tol=10^(-3);N=10000; k=0; fa=homework2(a); % f 需事先定义 for k=1:N p=(a+b)/2; fp=homework2(p); if( fp==0 || (b-a)/2tol) break end if fa*fp0 b=p; else a=p; end end k,p 程序运行结果： k = 10 p = 1.325195312500000 3. 用Newton迭代法求方程 的一个正根，计算结果精确到7位有效数字. 要求给出程序和运行结果. 解： 取迭代初值 ,并设,则． 牛顿迭代函数为 牛顿迭格式为： Matlab程序如下： %定义zuoye3