计算方法实验

1、实验十一、编程并在计算机上调试修改运行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;end2. 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(

2、分母为零) 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 endendif k=n disp(迭代失败)endfunction 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

3、);a(k+2)=y; if abs(y-x1)w x0=x1;x1=y; else break endendif k=n disp(迭代失败)endfunction 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) f(1,10)ans = 2.6385 f(2,10)ans = 7.4708 f(1,100)ans = 2.9046 f(0,

4、10)ans = 12、function y=f(x,n)y=0;for k=1:ny=y+(-1)(k-1)*x(2*k-1)/(jiecheng(2*k-1);endfunction y1=jiecheng(a)y1=1;if a=0 y1=1;elsefor k=1:a y1=y1*k;endend f(1,10)ans = 0.7897 f(2,10)ans = 0.5641 f(pi/2,100)ans = 1.0000二、1、 A=1 2;2 1; E,k=jzhexp(A)E = 10.9552 9.8109 9.8109 10.9552k = 225 A=1 0;0 1; E,

5、k=jzhexp(A)E = 2.9046 0 0 2.9046k = 1792、 a,k=xjfdd(0,1,1e-4,100) a = Columns 1 through 3 0 1.0000 0.1442 Columns 4 through 6 0.5347 0.3637 0.5597 Columns 7 through 9 0.5027 0.3379 0.1089 Columns 10 through 12 0.0177 0.1049 0.9061 Columns 13 through 15 0.6883 0.8284 0.1495 Columns 16 through 18 0.55

6、70 0.0788 0.9357 Column 19 0.7215k =17 a,k=xjfdd(0.5,1,1e-4,100)a = Columns 1 through 3 0.0000 1.0000 0.7742 Columns 4 through 6 0.5088 0.2704 0.8608k = 43、 a,k=kxjfdd(0.5,1,1e-4,100) a = Columns 1 through 3 0.0000 1.0000 0.7742 Columns 4 through 6 0.1843 0.4186 0.8679 Column 7 0.9005k = 5 a,k=kxjfdd(0.5,6,1e-4,100)a = Columns 1 through 3 0.000