数值计算上机操作答案.docx

第 1 题. 设In= ? 1 x n dx 0 5 ? x 由递推公式In=-5I  1 + ,从I 的几个近似值出发，计算I ； n-1 n 0 20 0clc clear syms xy(1)=0.18232;for i=2:21y(i)=-5*y(i-1)+1/i; end 0 clc clear syms x y(1)=0.18232; for i=2:21 y(i)=-5*y(i-1)+1/i; end disp(y(21)); 其运算结果为： 8.9577e+012 粗糙估计I  ，用I 20 n-1  1 1 =- I + ,计算I ； 5 n 5 n 0 20clear y(21)=1/105;for i=2:21y(i-1)=-0.2*y(i)+0.2*i; 20 clear y(21)=1/105; for i=2:21 y(i-1)=-0.2*y(i)+0.2*i; end disp(y(1)) 其运算结果为： 0.4000 第 2 题. 求方程ex+10x-2=0的近似根，要求x ? x ? 5 ? 10 ?4 ，并比较计算量。 k ?1 k 0 解： 用二分法，编写程序代码为下: clear clc clear clc r=0.0005; a=0; b=1; y1=exp(a)+10*a-2; y2=exp(b)+10*b-2; k=1; while abs(b-a)r k=k+1; x=(a+b)/2; y=exp(x)+10*x-2; if y1*y0 b=x else a=x end end disp(x) disp(k) 其运算结果为： b = 0.5000 b = 0.2500 b = 0.1250 a = 0.0625 b = 0.0938 a = 0.0781 a = 0.0859 a = 0.0898 b = 0.0918 b = 0.0908 1 a = 0.0903 0.0903 12 用题设迭代法，取初值x(1)=0，并用迭代x(i)=(2-exp(x(i-1)))/10，程序代码如下: clear clc x(1)=0; clear clc x(1)=0; r=0.0005; for i=2:20 x(i)=(2-exp(x(i-1)))/10; if abs(x(i)-x(i-1))r return; end end disp(x(i)) disp(i) 其运算结果为： 0.0905 5 加速迭代： ? ?构造序列 ? ? k }，使 ? ? x x ? x 2 ? ? x = k k x k ? 2 k ? 2 2 x  k ?1 k ?1 x k ，其中 x k 是 x 的一个近似值，经计算： x =0； x 0 1 =0.1； x 2 =0.0894829； x 3 =0.090639； x 4 =0.0905126； ? ? ? ? ? ? x =0.090484； x =0.090524449； x 0 1 2 =0.090525536。 (4) 用牛顿法，取初值x=0，编写程序代码为下: clear r=0.00005; x=0; clear r=0.00005; x=0; for i=1:20 x=x-(exp(x)+10*x-2)/(exp(x)+10); if abs((exp(x)+10*x-2)/(exp(x)+10))r; return; 2 end end disp(x) disp(i) end end disp(x) disp(i) 其运算结果为： 0.0905 2 上机心得： 不同的计算方式是有收敛速度上的不同的。有本题可以看到，在要求同一运算精度的情况下，采用二分法运算了 12 次，采用题设的迭代方法运算了 5 次，采用加速迭代法运算了 3 次， 采用牛顿迭代法则只需运算 2 次。 第 3 题. 钢水包使用次数多以后，钢包的容积增大，数据如下： x 2 3 4 5 6 7 8 9 y 6.42 8.2 9.58 9.5 9.7 10 9.93 9.99 10 11 12 13 14 15 16 10.49 10.59 10.60 10.8 10.6 10.9 10.76 试从中找出使用次数和容积之间的关系，计算均方差。（注：增速减少，用何种模型） 解:由数据点分布图可知，拟合曲线 y=f(x)随着 x 的增加而上升，但上升速度由快到慢，当 x 趋于无穷大时，y 趋于某个常数，故曲线有一水平渐进线。根据上述特征很容易想到用 Logistic 模型来拟合该曲线。 设 y=f(x)具有指数形式 y ? ae b 1 x（a0,b0）。对此式两边取对数，得ln y ? ln a ? b x 。 记 A=l