2011数值分析作业.docVIP

第一章作业 求积分 function qiujifen = qiuzhi() y = 1.^9.*exp(1-1); y1 = 0.5*y; y2 = y1./2+1./2.*0.5.^9.*exp(0.5-1); S1 = 4.*y2./3-y1./3; for k = 2:inf z = 0; for i = 2^k for j = 1:2:i-1 z = z+j./i.^9.*exp(j./i-1); end y1 = y2; y2 = y1./2+1./i.*z; S2 = 4.*y2./3-y1./3; end if max(abs(S2-S1))15*0.5e-6 %zhong zhi tiaojian break else S1 = S2; end end disp(结果为) y = y2 Mesh命令画出二元函数 function mesh = mesh(x,y,z) x = linspace(-10,10,100); y = x; [a,b] = meshgrid(x, y); z = exp(-abs(a))+cos(a+b)+1./(a.^2+b.^2+1); surf(x,y,z); 第二章作业 列主元三角分解，求逆并验证 function pascal = pascal() A = [1 1 1 1 1;1 2 3 4 5;1 3 6 10 15;1 4 10 20 35;1 5 15 35 70]; n = length(A); P = eye(n); L = eye(n); U = zeros(n,n); for j = 1:n max = j; for i = j+1:n B = zeros(1,n-j+1); %%%%%%%求最大列主元所在行数%%%%%% if j == 1 abs(A(i,j))abs(A(max,j)) max = i; elseif j ~= 1 w = abs(A(i,j)-A(i,1:j-1)*A(1:j-1,j)); v = abs(A(max,j)-A(max,1:j-1)*A(1:j-1,j)) if wv max = i; end end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% end B = P(j,:); P(j,:) = P(max,:); P (max,:) = B; temp = A(j,j:n); A(j,j:n) = A(max,j:n); A(max,j:n) = temp; if j == 1 A(j,j) = A(j,j) else if j ~= 1 A(j,j) = A(j,j)-A(j,1:j-1)*A(1:j-1,j); end end U(j,j) = A(j,j); for i = j+1:n if j == 1 A(i,j) = A(i,j)./A(j,j); else A(j,i) = A(j,i)-A(j,1:j-1)*A(1:j-1,i); A(i,j) = (A(i,j)-A(i,1:j-1)*A(1:j-1,j))./A(j,j); end U(j,i) = A(j,i); L(i,j) = A(i,j); end end disp(列主元三角分解为：) P = P L = L U = U V = U\eye(n); V = V*(L\eye(n)); V = V*P; disp(矩阵A的逆矩阵:) V = V A = [1 1 1 1 1;1 2 3 4 5;1 3 6 10 15;1 4 10 20 35;1 5 15 35 70]; disp(左除命令验证:) V = A\eye(n) 运行结果： 列主元三角分解为： 可以验证列主元三角分解法正确 求逆矩阵函数 function ni = nijuzhen()