数值计算方上机实习题.docxVIP

数值计算方法上机实习题设，由递推公式，从, 出发，计算；，, 用，计算；(1)由I0计算I20递推子程序：function f=fib(n,i)if n=1 f=fib(n-1,i)*(-5)+(1/(n));elseif n==0 f=i;end计算和显示程序：I=0.1824;I1=0.1823;fib1=fib(20,I);fib2=fib(20,I1);fprintf(I_0=0.1824时，I_20=%d\n,fib1);fprintf(I_0=0.1823时，I_20=%d\n,fib2);计算结果显示：I_0=0.1824时，I_20=7.480927e+09I_0=0.1823时，I_20=-2.055816e+09(2)由I20计算I0程序：n=21;i1=0;i2=10000;f1=i1;f2=i2;while n~=0; f1=f1*(-1/5)+(1/(5*n)); f2=f2*(-1/5)+(1/(5*n)); n=n-1;endfprintf(I_20=0 时，I_0=%4.8f\n,f1);fprintf(I_20=10000时，I_0=%4.8f\n,f2);计算结果显示：I_20=0 时，I_0=0_20=10000时，I_0=0析结果的可靠性及产生此现象的原因(重点分析原因)。答：第一个算法可得出易知第一个算法每一步计算都把误差放大了5倍，n次计算后更是放大了5n倍，可靠性低。第二个算法可得出可以看出第二个算法每一步计算就把误差缩小5倍，n次后缩小了5n倍，可靠性高。求方程的近似根，要求，并比较计算量。在[0，1]上用二分法；[0，1]上的二分法二分法子程序：function [root,n]=EFF3(f,x1,x2)%第二题(1)二分法f1=subs(f,symvar(f),x1);%函数在x=x1的值f2=subs(f,symvar(f),x2);%x=x2n=0;%步数er=5*10^-4;%误差if(f1==0) root=x1; return;elseif(f2==0) root=x2; return;elseif(f1*f20) disp(两端点函数值乘积大于0！); return;else while(abs(x1-x2)er)%循环 x3=(x1+x2)/2; f3=subs(f,symvar(f),x3); n=n+1; if(f3==0) root=x3; break; elseif(f1*f30) x1=x3; else x2=x3; end end root=(x1+x2)/2;%while循环少一步需加上end计算根与步数程序：fplot(@(x) exp(x)+10*x-2,[0,1]);grid on;syms x;f=exp(x)+10*x-2;[root,n]=EFF3(f,0,1);fprintf(root=%6.8f ,n=%d \n,root,n);计算结果显示：root=0,n=11 取初值，并用迭代；初值x0=0迭代迭代法子程序：function [root,n]=DDF(g,x0,err,max)（接下页）%root根，n+1步数，g函数，x0初值，err误差，max最大迭代次数X(1)=x0;for n=2:max X(n)=subs(g,symvar(g),X(n-1)); c=abs(X(n)-X(n-1)); root=X(n); if(cerr) break; end if n==max disp(超出预设迭代次数); endendn=n-1;%步数减1计算根与步数程序：syms x;f=(2-exp(x))/10;（接下页）x0=0; err=5*10^(-4);max=100;[root,n]=DDF(f,x0,err,max);fprintf(root=%6.8f ,n=%d \n,root,n);计算结果显示：root=0,n=4加速迭代的结果;加速迭代加速迭代计算程序：x0=0;err=5*10^(-4);max=100;syms x;g=x-(f(x)-x)^2/(f(f(x))-2*f(x)-x);[root,n]=DDF(g,x0,err,max);fprintf(root=%6.8f,