第二单元用python学习微积分（十）曲线构图下和最值问题

本文链接：https://blog.csdn.net/bullseye/article/details/122385613

本文内容来自于学习麻省理工学院公开课：单变量微积分-最值问题-网易公开课

开发环境准备：CSDN

一、曲线构图例子

二、曲线构图总结：

1、描点

a 找到函数的不连续点（尤其是函数值趋向于无穷的点）

b 找到无限远端的点即x趋向于无穷（）的点

c 标出那些容易找到的点，比如在坐标轴上的....(optional)

4、根据2次导数检查（由于求二阶导数经常非常复杂，所以没有要求尽量避免这一步）

一、曲线构图例子

$f(x) = \frac{x+1}{x+2}$

$f'(x) = (\frac{x+2-1}{x+2})' = (1-\frac{1}{x+2})' =- \frac{1}{(x+2)^2}\ne0$ (看出原函数双曲线函数)

因为导数永不为0，所以这个函数没有驻点

画这个函数的简图，首先考虑非法点x=-2

当 $x=-2^+$ 则 $f(x) = \frac{-1}{-2^+ + 2} = \frac{-1}{0^+} = -\infty$

当 $x=-2^-$ 则 $f(x) = \frac{-1}{-2-+ + 2} = \frac{-1}{0^-} = \infty$

可以想象这个会是

这时考虑当x取 $\pm\infty$ 时 $f(x) = \frac{x+1}{x+2} = \frac{1+\frac{1}{x}}{1 + \frac{2}{x}}_{x\rightarrow \pm\infty} = 1$ ，同时由于f'(x) 永不为0，因此没有驻点，所以不会出现翻转（没有turning point，函数中没有一个点有一根平的切线）（非法点也已经考虑清楚了）。

检查：这双曲线在下-infinity&lt;x&lt;-2 -2&lt;x&lt;infinity 随x增长y增长，因为导数恒正

from sympy import *
import numpy as np 

import matplotlib.pyplot as plt 

fig = plt.figure()
ax = fig.add_subplot(1, 1, 1)
ax.spines['left'].set_position('zero')
ax.spines['bottom'].set_position('zero')
ax.spines['right'].set_color('none')
ax.spines['top'].set_color('none')
ax.xaxis.set_ticks_position('bottom')
ax.yaxis.set_ticks_position('left')
ax.set_aspect( 1 ) 

def DrawXY(xFrom,xTo,steps,expr,color,label,plt):
    yarr = []
    xarr = np.linspace(xFrom ,xTo, steps) 
    for xval in xarr:
        yval = expr.subs(x,xval)
        yarr.append(yval)
    y_nparr = np.array(yarr) 
    plt.plot(xarr, y_nparr, c=color, label=label)    

def TangentLine(exprY,x0Val,xVal):
    diffExpr = diff(exprY)
    x1,y1,xo,yo = symbols('x1 y1 xo yo')
    expr = (y1-yo)/(x1-xo) - diffExpr.subs(x,x0Val)
    eq = expr.subs(xo,x0Val).subs(x1,xVal).subs(yo,exprY.subs(x,x0Val))
    eq1 = Eq(eq,0)
    solveY = solve(eq1)
    return xVal,solveY

def DrawTangentLine(exprY, x0Val,xVal1, xVal2, clr, txt):
    x1,y1 = TangentLine(exprY, x0Val, xVal1)
    x2,y2 = TangentLine(exprY, x0Val, xVal2)
    if len(txt)>0:
        plt.plot([x1,x2],[y1,y2], color = clr, label=txt)
    else:
        plt.plot([x1,x2],[y1,y2], color = clr)
        
x= symbols('x')
y = (x+1)/(x+2)
DrawXY(-10,-2.1,100,y,'green','(x+1)/(x+2)',plt)
DrawXY(-1.9,8,100,y,'green','(x+1)/(x+2)',plt)

plt.plot([-12,10],[1,1], color = 'gray', label= 'y=1')

plt.legend(loc='upper right')
plt.show()

$f''(x) = -\frac{2}{(x+2)^3} (x\ne-2)$

$f''(x) >0 |_{ -\infty<x<-2}$ 上凹函数

$f''(x) <0 |_{ -2<x<\infty}$ 下凹函数

这两个二阶导数的结果告诉我们这个原函数的图像不是波浪形

二、曲线构图总结：

1、描点

a 找到函数的不连续点（尤其是函数值趋向于无穷的点）

b 找到无限远端的点即x趋向于无穷（ $\pm\infty$ ）的点

c 标出那些容易找到的点，比如在坐标轴上的....(optional)

2、计算

a 求导f'(x) = 0

b 求驻点和驻点的值

3、根据导数检查

在每个以驻点或不连续点为端点的区间根据f'(x)检查是否递增或递减

4、根据2次导数检查（由于求二阶导数经常非常复杂，所以没有要求尽量避免这一步）

在每个以驻点或不连续点为端点的区间根据f''(x) >或< 0, 检查是否符合上凹(concave up)或下凹(concave down)

f''(x) = 0 求出拐点( inflection point )

5、综合以上

三、例子

$f(x) = \frac{x}{ln(x)}$ (x>0 , 对数的定义域没有负）

$x\rightarrow1, f(x)\rightarrow\infty$

1、描点

a $f(1^+) = \frac{1}{ln1^+} = \frac{1}{0^+} = \infty$

$f(1^-) =- \infty$

b 远端点

$f(0^+) = \frac{0^+}{ln0^+} = \frac{0^+}{\infty} = 0$

$f(10^{10}) = \frac{10^{10}}{ln10^{10}} = \frac{10^{10}}{10ln10}$ (2<ln10 <3) 因此 $f(\infty) = \infty$

2、计算

$f'(x) = \frac{x'ln(x) - ln(x)'x}{(ln(x))^2} = \frac{ln(x)-1}{ln(x)^2}|_{x=e} = 0$

$f(e) = \frac{e}{lne} = e$

驻点为：(e,e)

x= symbols('x')
y = x/ln(x)
DrawXY(0,0.9,100,y,'r','',plt)
DrawXY(1.1,10,100,y,'r','',plt)
plt.plot([np.e,np.e],[-10,12.5], color = 'gray', label= 'x=e')
plt.plot([1,1],[-10,12.5], color = 'g', label= 'x=1')

3、根据导数检查

$f'(x) = \frac{x'ln(x) - ln(x)'x}{(ln(x))^2} = \frac{ln(x)-1}{ln(x)^2}$

0<x<1时 ln(x) < 0 , ln(x)*ln(x) > 0 因此此时f'(x) < 0，f(x) 递减

1<x<e时 0<ln(x) < 1 , ln(x)*ln(x) > 0 因此此时f'(x) < 0，f(x) 递减

e=x 时 ln(x) -1 = 0 , 驻点

x>e 时ln(x) - 1 > 0 , ln(x)*ln(x) > 0 因此此时f'(x) > 0，f(x) 递增

符合图像

4、根据2次导数检查

$f'(x) = \frac{ln(x)-1}{ln(x)^2}$

$f''(x) = (\frac{1}{ln(x)} - \frac{1}{ln^2(x)} )' = -\frac{1}{ln^2(x)}\frac{1}{x} + 2\frac{1}{ln^3(x)}\frac{1}{x} = \frac{2-ln(x) }{ln^3(x)x}$ （从这能看出其实x=0点处的f'(x) = 0）