matlab第四章课后答案.docVIP

matlab第四章课后答案 %Exercise 1(1) roots([1 1 1]) %Exercise 1(2) roots([3 0 -4 0 2 -1]) %Exercise 1(3) p=zeros(1,24); p([1 17 18 22])=[5 -6 8 -5]; roots(p) %Exercise 1(4) p1=[2 3]; p2=conv(p1, p1); p3=conv(p1, p2); p3(end)=p3(end)-4; %原p3最后一个分量-4 roots(p3) %Exercise 2 fun=inline(x*log(sqrt(x^2-1)+x)-sqrt(x^2-1)-0.5*x); fzero(fun,2) %Exercise 3 fun=inline(x^4-2^x); fplot(fun,[-2 2]);grid on; fzero(fun,-1),fzero(fun,1),fminbnd(fun,0.5,1.5) %Exercise 4 fun=inline(x*sin(1/x),x); fplot(fun, [-0.1 0.1]); x=zeros(1,10);for i=1:10, x(i)=fzero(fun,(i-0.5)*0.01);end; x=[x,-x] %Exercise 5 fun=inline([9*x(1)^2+36*x(2)^2+4*x(3)^2-36;x(1)^2-2*x(2)^2-20*x(3);16*x(1)-x(1)^3-2*x(2)^2-16*x(3)^2],x); [a,b,c]=fsolve(fun,[0 0 0]) %Exercise 6 fun=@(x)[x(1)-0.7*sin(x(1))-0.2*cos(x(2)),x(2)-0.7*cos(x(1))+0.2*sin(x(2))]; [a,b,c]=fsolve(fun,[0.5 0.5]) %Exercise 7 clear; close; t=0:pi/100:2*pi; x1=2+sqrt(5)*cos(t); y1=3-2*x1+sqrt(5)*sin(t); x2=3+sqrt(2)*cos(t); y2=6*sin(t); plot(x1,y1,x2,y2); grid on; %作图发现4个解的大致位置，然后分别求解 y1=fsolve([(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4],[1.5,2]) y2=fsolve([(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4],[1.8,-2]) y3=fsolve([(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4],[3.5,-5]) y4=fsolve([(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4],[4,-4]) %Exercise 8(1) clear; fun=inline(x.^2.*(x.^2-x-2)); fplot(fun,[-2 2]);grid on; %作图观察 x(1)=-2; x(2)=fminbnd(fun,-1,-0.5); x(4)=fminbnd(fun,1,2); fun2=inline(-x.^2.*(x.^2-x-2)); x(3)=fminbnd(fun2,-0.5,0.5); x(5)=2 feval(fun,x) %答案: 以上x(2)(4)是局部极小，x(1)(3)(5)是局部极大，从最后一句知道x(1)全局最大， x(4)最小。 %Exercise 8(2) clear; fun=inline(3*x.^5-20*x.^3+10); fplot(fun,[-3 3]);grid on;%作图观察 x(1)=-3; x(3)=fminsearch(fun,2.5); fun2=inline(-(3*x.^5-20*x.^3+10)); x(2)=fminsearch(fun2,-2.5); x(4)=3; feval(fun,x) %Exercise 8(3) fun=inline(abs(x^3-x^2-x-2)); fplot(fun,[0 3]);grid on;%作图观察 fminbnd(fun,1.5,2.5) fun2=inline(-abs(x^3-x^2-x-2)); fminbnd(fun2,0.5,1.5) %Exercise 9 close; x=-2:0.1:1;y=-7:0.1:1