matlab参考答案_第4568章.docVIP

第四章 1、 H=hilb(5); P=pascal(5); Dp=det(P); Dh=det(H); Kp=cond(P); Kh=cond(H) P更好，它的条件数更小 2、 A=[1,-1,2,3;0,9,3,3;7,-5,0,2;23,6,8,3] B=[3,(pi)/2,45;32,-76,sqrt(37);5,72,4.5e-4;exp(2),0,97] A1=diag(A); B1=diag(B); A2=triu(A); B2=triu(B); A3=tril(A); B3=tril(B); A4=rank(A); B4=rank(B); A5=norm(A); B5=norm(B); A6=cond(A); B6=cond(B); A7=trace(A); B7=trace(B); 3、 A=[31,1,0;-4,-1,0;4,-8,-2] [V,D]=eig(A) A*V V*D 4、 A=diag([-1,-1,-1,-1],-1)+diag([-1,-1,-1,-1],1)+diag([2,2,2,2,2]) b=[1,0,0,0,0]; x1=inv(A)*b; x2=A\b; [L,U]=lu(A); X3=U\(L\b); [Q,R]=qr(A); X4=R\(Q\b); R=chol(A); X5=R\(R\b) 5、 A=diag([-1,-1,-1,-1],-1)+diag([-1,-1,-1,-1],1)+diag([-2,-2,-2,-2,-2]) b=[1,0,0,0,0]; B=sparse(A); X1=inv(B)*b; X2=B\b; [L,U]=lu(B); X3=U\(L\b); X1=x1 X2=x2 X3=x3 第5章 1、 a=rand(1,30000); mean(a) std(a) max(a) min(a) length(find(a0.5))/30000 disp([num2str(ans*100) %]) 2 h=[466,715,950,1422,1635]; w=[7.04,4.28,3.40,2.52,2.13]; hh=[500,900,1500]; ww=interp1(h,w,hh,spline) 3 、 x=linspace(1,101,50); y=log10(x); f=polyfit(x,y,5); yy=polyval(f,x); polt(x,y,r-,x,yy,g.); 4、 N=64; T=5; t=linspace(0,5,64); x=exp(-t); dt=t(2)-t(1); f=1/dt; X=fft(x); F=X(1:N/2+1); f=f*(0:N/2)/N; plot(f,abs(F),m-*); xlabel(Frequency); ylabel(|F(k)|) 5. p1=[1 2 0 7];p2=[1 -2];p3=[1 0 5 1]; p=conv(p1,p2)+[0 p3]; roots(p); A=[-1,4,3;2,1,5;0,5,6]; Px=polyval(p,A) pxm=polyvalm(p,A) 6、 （1）z=fzero(3*x-sin(x)+1,0) (2)function F=byp(X) x=X(1);y=X(2); F(1)=x*x+y*y-9; F(2)=x+y-1; x=fsolve(@byp,[3,0],optimset(Display,off)); 第六章 1、 k=1; for x=1:0.01:3 fx(k)=det([x,x^2,x^3;1,2*x,3*x^2;0,2,6*x]); k=k+1; end F=diff (fx)/0.01 F(1) F(100) F(200) 2, x=0:0.01:1; y=x.^10+10.^x+1./log10(x); dy=diff(y)/0.01 x=0:0.01:1; y=log(1+x); dy=diff(y,2)/0.01; dy(end) 3. (1) fx=inline(x.^2.*sqrt(2*x.*x+3)); quad(fx,1,5) (2) fx=inline(x./sin(x).^2); quad(fx,pi/4,pi/3) (3) fx=inline(abs(cos(x+y))); dblquad(fx,0,pi,0,pi) syms x y; Int(int(x*y,x,y^2,y+2),x,-1,2) 4. x=0.3:0.2:1.5; y=[0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8325]; trapz(x,y) 5. （1）functi