东南大学数值分析上机题matlab(前三章).docxVIP

数值分析上机题第一章（17题）(1)从2依次累加到N的程序functionsn = sum1( n )sn=0;sn=single(sn);fori=2:nai=1/(i^2-1);sn=sn+ai;endend（2）从N依次累加到2的程序functionsn = sum2( n )sn=0;sn=single(sn);fori=n:-1:2ai=1/(i^2-1);sn=sn+ai;endend（3）编制求精确值的求和函数sum0functionsn = sum0( n )sn=0;sn=single(sn);sn=1/2*(3/2-1/n-1/(n+1));end按第一种顺序得到的值及有效位数如下：N=100时sn0=sum0(100);sn=sum1(100)n=fix(-log10(2*abs(sn-sn0)))得到：sn =0.7400495n =7N=10e4时sn0=sum0(10e4);sn=sum1(10e4)n=fix(-log10(2*abs(sn-sn0)))得到：sn =0.7498521n =3N=10e6时sn0=sum0(10e6);sn=sum1(10e6)n=fix(-log10(2*abs(sn-sn0)))得到：sn =0.7498521n =3按第二种顺序得到的值及有效位数如下：N=100时sn0=sum0(100);sn=sum2(100)n=fix(-log10(2*abs(sn-sn0)))得到：sn =0.7400495n =7N=10e4时sn0=sum0(10e4);sn=sum2(10e4)n=fix(-log10(2*abs(sn-sn0)))得到：sn =0.7499900n =7N=10e6时sn0=sum0(10e6);sn=sum2(10e6)n=fix(-log10(2*abs(sn-sn0)))得到：sn =0.7499999n =7（4）通过这道上机题，我明白了，应用计算机进行求和运算时，求和的顺序不同对结果的精度是有影响的。具体而言，小数先求和，大数后求和算得的结果更精确；反之，则会出现“大数吃小数”的现象，从而精度会降低。第二章（20题）（1）用牛顿法解f（x）=0的根的通用程序如下：function x = newton( f,df,x0,ep )%用牛顿法解f（x）=0的根 %f为待求解函数，df为f的导函数，x0为迭代初值，ep为允许误差epi=100;%为了保证能进入下列循环，epi取一个大于ep的初值，取epi=100forepiep x=x0-feval(f,x0)/feval(df,x0);epi=abs(x-x0);x0=x;endend（2）1）编程算法如下：显然x0=0时，迭代序列收敛于根x2。设置步长b=0.00001，从x0=0开始每次加一个步长b,然后判断迭代序列是否收敛于根x2,若收敛于根x2，则继续加步长b；若不收敛与根x2，则退出循环，此时的x0即为题中所求的尽可能大的δ。程序如下：正向找：f=@(x)x^3/3-x;df=@(x)x^2-1;x0=0;x=0;while x==0x=newton(f,df,x0,1e-10);x0=x0+0.00001;endx0得到：x0 =0.77461负向找f=@(x)x^3/3-x;df=@(x)x^2-1;x0=0;x=0;while x==0x=newton(f,df,x0,1e-10);x0=x0-0.00001;endx0得到x0 =-0.77461故δ=0.774612）取x0∈(-∞,-1),x0=-100f=@(x)x^3/3-x;df=@(x)x^2-1;x=newton(f,df,-100,1e-10)收敛于：x =-1.7321取x0∈(-1,-0.77461),x0=-0.79f=@(x)x^3/3-x;df=@(x)x^2-1;x=newton(f,df,-0.79,1e-10)收敛于：x = -1.7321取x0∈(-0.77461, 0.77461),x0=-0.76f=@(x)x^3/3-x;df=@(x)x^2-1;x=newton(f,df,-0.76,1e-10)收敛于：x =0取x0∈(-0.77461, 0.77461),x0=0.77f=@(x)x^3/3-x;df=@(x)x^2-1;x=newton(f,df,0.77,1e-10)收敛于：x =0取x0∈(0.77461, 1),x0=0.79f=@(x)x^3/3-x;df=@(x)x^2-1;x=newton(f,df,0.79,1e-10)收敛于：x =1.7321取x0∈(1, +∞),x0=100f=@(x)x^3/3-x;df=@(x)x^2-1;x=newton(f