最新上海電力學(xué)院數(shù)值計(jì)算方法上機(jī)實(shí)習(xí)題

1、 2017 數(shù)值計(jì)算方法 上機(jī)實(shí)習(xí)報(bào)告 學(xué)院： 學(xué)院： 專業(yè)： 專業(yè)： 班級(jí)： 班級(jí)： 姓名： 姓名： 學(xué)號(hào)： 學(xué)號(hào)： 第一種算法為正向迭代算法， 每計(jì)算一步誤差增長 5 倍， 雖然所給的初始值很接近，隨著 n 的增大，

2、誤差也越來越大。 * n n n I I e ? ?n n n n n n n n n e I I n I n I e I I e 51 ) ( 51 ) 5151 ( 5151 * * *1 1 1 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?nn e e ) 51 ( 0 ?第二種算法為倒向迭代算法， 每計(jì)算一步誤差縮小 5 倍， 雖然所給的初始值之間差很多，隨著 n 的增大，誤差也越來越小。算法趨近穩(wěn)定，收斂，可以選

3、擇這種算法。 2． 求方程 0 2 10 ? ? ? x e x 的近似根，要求 41 10 5 ?? ? ? ? k k x x ，并比較計(jì)算量。 （1） 在[0，1]上用二分法； （2） 取初值 0 0 ? x ，并用迭代 1021xke x ? ? ? ； （3） 加速迭代的結(jié)果； （4） 取初值 0 0 ? x ，并用牛頓迭代法； （5） 分析絕對誤差。 解： （1）利用二分法，MATLAB 程序如下： %% 二分法程序 cle

4、ar all clc a=0;b=1; f=@(x)(exp(x)+10*x-2); %@是定義函數(shù)句柄的運(yùn)算符 c=(a+b)/2;%取區(qū)間中點(diǎn) i=0;%分割次數(shù) while abs(f(c))>5*10^(-4) %判斷f(x)的精度是否滿足要求 if f(a)*f(c)<0 b=c;c=(a+b)/2; elseif f(b)*f(c)<0