MATLAB练习第四章.docVIP

第四章 2.求方程xln(+x)--0.5x=0的正根； fun=inline(x*log(sqrt(x^2-1)+x)-sqrt(x^2-1)-0.5*x); fzero(fun,2) 答：所求正根为2.1155。 8.作出下列函数图形，观察所有局部极大，局部极小和全局最大，全局最小值点的粗略位置，并用MATLB函数fminbnd和fminsearch求各极值点的确切位置。 （1）f(x)=x2sin(x2-x-2); [-2,2] (2) f(x)=3x5-20x3+10; [-3,3] (3) f(x)= ; [0,3] 解：（1）clear; fun=inline(x.^2.*sin(x.^2-x-2)); fplot(fun,[-2 2]);grid on; x(1)=-2; x(3)=fminbnd(fun,-1,-0.5); x(5)=fminbnd(fun,1,2); fun2=inline(-x.^2.*sin(x.^2-x-2)); x(2)=fminbnd(fun2,-2,-1); x(4)=fminbnd(fun2,-0.5,0.5); x(6)=2 feval(fun,x) （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) （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) 答：局部极大值点,-1.5326，0，2；局部极小值点-0.7315,1.5951，-2；全局极小值为-2,；全局极大值-1.5326.（2）局部极大值点-2,3；局部极小值点-3,2；全局极小值为-3；全局极大值3；（3）局部极大值点1,3；局部极小值点0,2；全局极小值为2；全局极大值3. 9.考虑函数 f(x,y)=y3/9+3x2y+9x2+y2+xy+9;, 作出f(x,y)在-2＜x＜1，-7＜y＜1的图，观察极值点位置， 用MATLB函数fminsearch求极值点和极值。 解：close; x=-2:0.1:1;y=-7:0.1:1; [x,y]=meshgrid(x,y); z=y.^3/9+3*x.^2.*y+9*x.^2+y.^2+x.*y+9; mesh(x,y,z);grid on; fun=inline(x(2)^3/9+3*x(1)^2*x(2)+9*x(1)^2+x(2)^2+x(1)*x(2)+9); x=fminsearch(fun,[0 0]) fun2=inline(-(x(2)^3/9+3*x(1)^2*x(2)+9*x(1)^2+x(2)^2+x(1)*x(2)+9)); x=fminsearch(fun2,[0 -5]) 答：极大值点（-1/3，-6），极小值点（0,0），鞍点（-7/6,-7/2）,(5/6,-5/2); 10.假定某天的气温变化记录如第二章习题5，试用最小二乘方法找出这一天的气温规律，考虑下列类型函数，作图比较效果，并计算均方误差。 （1）二次函数， （2）三次函数； （3）钟形函数f(x)=aeb(t-14) （4）函数f(x)=rsin(t+) clear;t=0:24; c=[15 14 14 14 14 15 16 18 20 22 23 25 28 ... 31 32 31 29 27 25 24 22 20 18 17 16]; p2=polyfit(t,c,2) p3=polyfit(t,c,3) fun=inline(a(1)*exp(a(2)*(t-14).^2),a,t); a=lsqcurvefit(fun,[0 0],t,c) f=feval(fun, a,t) norm(f-c) plot(t,c,t,f) fun2=inline(b(1)*sin(pi/12*t+b(2))+20,b,t); b=lsqcurvefit(fun2,[0 0],t,c) figure f2=feval(fun2, b,t) norm(f2-c) plot(t,c,t,f2) 答：（1）-0.0936x+2.5943x+8.4157;(2)-0