MATLAB 应用 第四章 导数应用[精品].ppt

输出： c1=0 c2=0 h1=-1 f1=1 即函数有唯一驻点x=0，且g’(0)=f’’(0)=-10，f(0)=1。由于在驻点x=0处函数的二阶导数小于0，该点是极大值点。又因为这是唯一驻点，所以它是函数的最大值点。于是， 在x=0取最大值1。因此当x1且x≠0时 ，有 * 【例7】证明 ，其中x0。 作图，输入： ezplot(atan(x)+1/x,[4,20]) 输出如图4-7所示。 当x趋向于正无穷时，直线 是 的渐近线。这也是用单调性证明的不等式。 * 图4-7 * 输入： syms x y=atan(x)+1/x; y1=diff(y,x) limit(y,x,inf) 输出： y1=1/(1+x^2)-1/x^2 ans=1/2*pi * 再研究单调性，输入 c=fzero(1/(1+x^2)-1/x^2,[1,inf]) 输出： c=Inf 即当x0时，函数y=arctan(x)+1/x无驻点。 再输入： x=2; y2=eval(1/(1+x^2)-1/x^2) 输出： y2=-0.0500 即当x0时，y1(x)=y’(x)0。于是y=arctanx+1/x单调减少，趋向于 。因此当x0时，有 。 * 【例8】 证明 ，当 时。 作图，输入： x=0:0.1:pi/2; y1=sin(x); y2=2*x/pi; plot(x,y1,b-,x,y2,r) 输出如图4-8所示。 * 图4-8 * 当x=0或x= 时，两曲线相交。在两点之间差距较大。这是需要用凹凸性证明的不等式的特征。 输入： x=0; f1=eval(sin(x)-2*x/pi) x=pi/2; f2=eval(sin(x)-2*x/pi ) 输出： f1=0 f2=0 即在区间端点，函数值等于零。 * 再输入： syms x f=sin(x)-2*x/pi; g=diff(f,x,2) 输出： g=-sin(x) 即g(x)=f’’(x)=-sin(x)。注意在区间(0, ) 内f’’(x)=-sin(x)0，函数f(x)的曲线上凸。又在区间的两个端点都等于零，因此在区间内部f(x)恒大于零，即当 时，有 。 * MATLAB 高等数学实验 实验四 导数应用 实验目的 理解并掌握用函数的导数确定函数的单调区间、凹凸区间和函数的极值的方法。进一步熟悉和掌握用MATLAB做平面图形的方法和技巧，掌握用MATLAB求方程的根（包括近似根）和求函数极值（包括近似极值）的方法。 * 4.1 学习MATLAB命令 * 4.1.1 求多项式方程近似根的命令 用MATLAB求多项式方程 的解的命令是roots，使用方法为： roots（c） 其中c是上述方程左端多项式的系数向量 * 4.1.2 求方程f(x)=0近似根的命令 命令的一般形式如下： (1)建立函数： f=inline(表达式) (2)求函数零点： c=fzero(f,[a,b]) %求函数f(x)在区间[a,b]内的零点c c=fzero(f,x0) %求函数f(x)在x0附近的零点c * 4.1.3 求非线性函数f(x)的极小值 用MATLAB求一元函数极小值命令是fminbnd，格式如下： (1) x=fminbnd(fun,x1,x2) (2) [x,fav1]=fminbnd(fun,x1,x2) (3) [x,fav1,exitflag,output]=fminbnd(fun,x1,x2) 其中： x=fminbnd(fun,x1,x2)是求(x1,x2)上fun函数的最小值x。 [x,fav1]=fminbnd(fun,x1,x2)返回解x处目标函数的值。 [x,fav1,exitflag,output]=fminbnd(fun,x1,x2)返回包含优化信息的结构输出。 * 注(1) 函数fminbnd的算法基于黄金分割法和二次插值法，要求目标函数必须是连续函数，此命令可能给出局部最优值。 注(2) 命令fminbnd是求函数f(x)的极小值，若要求函数f(x)的极大值，只需求-f(x)的极小值即可。 * 4.2 实验内容 * 4.2.1 求函数的单调区间 【例1】求 的单调区间。 输入： s