3数值计算方法上机实习题.docVIP

上海电力学院 数值计算上机报告 课 程： 现代数值计算 题 目： 数值计算方法上机实习题报告 院 系： 自动化工程学院 专业年级： 电机与电器 学生姓名： 黄 丽 学号：ys1310401009 指导教师： 黄建雄 2013年12月20日 数值计算方法上机实习题 设， 由递推公式，从的几个近似值出发，计算； 粗糙估计，用，计算； 分析结果的可靠性及产生此现象的原因(重点分析原因)。 程序： （1）由递推公式，从的几个近似值出发，计算； I(0) =0.1820；I(1)= 0.0900；I(2)= 0.0500；I(3)= 0.0833；I(4)= -0.1667； I=0.182; for n=1:20 I=(-5)*I+1/n; end 故计算结果= -3.0666e+10 (2)粗糙估计，用，计算； I=0.008; for n=20:-1:1 I=(-1/5)*I+1/(5*n); end I = 0.1823 （3）分析结果的可靠性及产生此现象的原因(重点分析原因)。 假设的真值为，误差为，即。对于真值也有。综合2个递推等式，有，即意味着只要n足够大，按照这种每计算一步误差增长5倍的方式，所得的结果总是不可信的，因此整个算法是数值不稳定的。而第二种方式的误差会以每计算一步缩小到1/5的方式进行，这样的计算结果和实际是很相近的。 求方程的近似根，要求，并比较计算量。 在[0，1]上用二分法； 取初值，并用迭代； 加速迭代的结果； 取初值，并用牛顿迭代法； 分析绝对误差。 （1）在[0，1]上用二分法； 程序： a=0;b=1.0;i=0; while abs(b-a)5*1e-4 c=(b+a)/2; if exp(c)+10*c-20 b=c; else a=c; end i=i+1; end c 方程的近似根为：x*=（a+b）/2 = 0.0906 。步长为i=11 。 （2）取初值，并用迭代； 程序: x=0;y=0.1;i=0; while abs(y-x)5*1e-4 y=x; x=(2-exp(x))/10; i=i+1; end 方程的近似根为：x=0.0905 。步长为i=4 。 （3）加速迭代的结果； 程序：x=0;xx=1;i=0; while abs(xx-x)5*1e-4 y=exp(x)+10*x-2; z=exp(y)+10*y-2; xx=x; x=x-(y-x)^2/(z-2*y+x); i=i+1; end x 方程的近似根为：x= 0.0995。步长为i=3 。 取初值，并用牛顿迭代法； 程序：x=0;y=0.1;i=0; while abs(y-x)5*1e-4 y=x; x=x-(exp(x)+10*x-2)/(exp(x) + 10);%diff(exp(x)+10*x-2)=exp(x) + 10 i=i+1; end x i 方程的近似根为：x= 0.0905。步长为i=2 。 分析绝对误差。 solve(exp(x)+10*x-2=0) 方程的精确解为x= 0.0905。 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 试从中找出使用次数和容积之间的关系，计算均方差。（注：增速减少，用何种模型） 解：（1）设y=f(x)具有指数形式（a0,b0）。对此式两边取对数，得。 记A=lna，B=b x=[2 3 4 5 6 7 8 9 10 11 12 13 14 15 16]; for i=1:15 X(i)=1/x(i); end y=[6.42 8.2 9.58 9.5 9.7 10 9.93 9.99 10.49 10.59 10.60 10.8 10.6 10.9 10.76