1、精選優(yōu)質(zhì)文檔-傾情為你奉上實(shí)驗(yàn)十一、編程并在計算機(jī)上調(diào)試修改運(yùn)行1. 利用公式計算 2. 利用公式計算 二、閱讀理解下列程序，并在計算機(jī)上運(yùn)行1. jzhexp.m %利用ex的展開式計算矩陣函數(shù)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為最大迭代次數(shù)function a,k=xjfdd(x0,x,w,n) a(1:2)=x0,x;for k=1:n if f

2、(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 endendif k=n disp('迭代失敗')endfunction y=f(x) %子函數(shù)y=x*exp(x)-1; 3. kxjfdd.m % x0,x為初值，w為精度，n為最大迭代次數(shù)function a,k=kxjfdd(x0,x1,w,n) a(1:2)=x0,x1;for k=1:n if f(x1)-f(x0)=0 disp('

3、分母為零') 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 endendif k=n disp('迭代失敗')endfunction y=f(x) %子函數(shù)y=x*exp(x)-1;三、選做題：修改單點(diǎn)弦截法程序dxjf.m，并在計算機(jī)上運(yùn)行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 a

4、bs(x2-x1)<w break end k=k+1; x1=x2;endif k=N x2='迭代失敗'endfunction y=f(x) %子函數(shù) y=4*cos(x)-exp(x) ;function y=g(x) %子函數(shù) y=(f(x)-f(x0)/(x-x0);一、1、function y=f(x,n)y=0;for k=0:ny=y+(xk)/(jiecheng(k);endfunction y1=jiecheng(n)y1=1;if n=0 y1=1;elsefor k=1:n y1=y1*k;endend>> f(1,10)ans = 2

5、.6385>> f(2,10)ans = 7.4708>> f(1,100)ans = 2.9046>> f(0,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,

6、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,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 thro

7、ugh 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.5570 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