计算方法实验10.docVIP

PAGE 2 实验十 一、编程并在计算机上调试修改运行 1. 利用公式计算 2. 利用公式计算 二、阅读理解下列程序，并在计算机上运行 1. jzhexp.m %利用ex的展开式计算矩阵函数eA的值 function [E,k]=jzhexp(A); %A是方阵 E=zeros(size(A));F=eye(size(A)); k=1; while norm(F,1)0 E=E+F; F=A*F/k; k=k+1; end 2. xjfdd.m % x0,x为初值，w为精度，n为最大迭代次数 function [a,k]=xjfdd(x0,x,w,n) a(1:2)=[x0,x]; for k=1:n if f(x)-f(x0)==0 disp(分母为零) break end y=x-(x-x0)*f(x)/(f(x)-f(x0));a(k+2)=y; if abs(y-x)w x=y; else break end end if k==n disp(迭代失败) end function y=f(x) %子函数 y=x*exp(x)-1; 3. kxjfdd.m % x0,x为初值，w为精度，n为最大迭代次数 function [a,k]=kxjfdd(x0,x1,w,n) a(1:2)=[x0,x1]; for k=1:n if f(x1)-f(x0)==0 disp(分母为零) break end y=x1-(x1-x0)*f(x1)/(f(x1)-f(x0)); a(k+2)=y; if abs(y-x1)w x0=x1;x1=y; else break end end if k==n disp(迭代失败) end function y=f(x) %子函数 y=x*exp(x)-1; 三、选做题：修改单点弦截法程序dxjf.m，并在计算机上运行 function x2=dxjf(x0,x1,w,N) k=1; while k~=N if g(x1)==0 x2=分母为零 end x2=x1-f(x1)/g(x1); if abs(x2-x1)w break end k=k+1; x1=x2; end if k==N x2=迭代失败 end function y=f(x) %子函数 y=4*cos(x)-exp(x) ; function y=g(x) %子函数 y=(f(x)-f(x0))/(x-x0); 一、 1、function y=f(x,n) y=0; for k=0:n y=y+(x^k)/(jiecheng(k)); end function y1=jiecheng(n) y1=1; if n==0 y1=1; else for k=1:n y1=y1*k; end end f(1,10) ans = 2.718281801146385 f(2,10) ans = 7.388994708994708 f(1,100) ans = 2.718281828459046 f(0,10) ans = 1 2、 function y=f(x,n) y=0; for k=1:n y=y+((-1)^(k-1))*x^(2*k-1)/(jiecheng(2*k-1)); end function y1=jiecheng(a) y1=1; if a==0 y1=1; else for k=1:a y1=y1*k; end end f(1,10) ans = 0.841470984807897 f(2,10) ans = 0.909297426825641 f(pi/2,100) ans = 1.000000000000000 二、 1、 A=[1 2;2 1]; [E,k]=jzhexp(A) E = 10.226708182179552 9.858828741008109 9.858828741008109 10.226708182179552 k = 225 A=[1 0;0 1]; [E,k]=jzhexp(A) E = 2.718281828