资源描述:
《matlab练习第四章 (5)》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、第四章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)=;[
2、0,3]解:(1)clear;fun=inline('x.^2.*sin(x.^2-x-2)');fplot(fun,[-22]);gridon;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)=2feval(fun,x)(2)clear;fun=i
3、nline('3*x.^5-20*x.^3+10');fplot(fun,[-33]);gridon;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,[03]);gridon;fminbnd(fun,1.5,2.5)fun2=in
4、line('-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;,(1)作出f(x,y)在-2<x<1,-
5、7<y<1的图,观察极值点位置,(2)用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);gridon;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,[00])fun2=inline
6、('-(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(
7、t+)clear;t=0:24;c=[15141414141516182022232528...313231292725242220181716];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,[00],t,c)f=feval(fun,a,t)norm(f-c)plot(t,c,t,f)fun2=inline('b(1)*sin(pi/12*t+
8、b(2))+20','b','t');b=lsqcurvefit(fun2,[00],t,c)figuref2=feval(fun2,b,t)norm(f2-c)plot(t,c,t,f2)答:(1)-0.0936x+2.5943x+8.4157;(2)-0.0080x+0.1931x-0.1022x+13.2513;(3)a=27.9628,b=-0.0058;(4)r=-6.9721,.