期末微积分考试试题求解 ：利用python求解

在 今年期末微积分考试试题：看看你能够在两个小时内做对几道题？ 搜集到了一份期末微积分考试试题。为了对其内容进行进一步分析，对其内容进行整理如下。

 

§01 填空题

---

每个空3分，共10题

1. 求解常微分方程通解

$$y^{\prime }=1+2x+y^2+2x{y}^2$$

◎ 求解：

$$2x{y}^2\left(x\right)+2x+y^2\left(x\right)+1=\frac{d}{dx}y\left(x\right)$$$$y\left(x\right)=\tan \left(C_1+x^2+x\right)$$

from headm import *                 # =

from sympy                  import symbols,simplify,expand,print_latex
from sympy                  import *

x,z = symbols('x z')
y,f,g = symbols('y f g', cls=Function)

diffeq = Eq(1 + 2*x + y(x)**2 + 2*x*y(x)**2, y(x).diff(x))
print_latex(diffeq)
print_latex(dsolve(diffeq, y(x)))

_ = tspexecutepythoncmd("msg2latex")
clipboard.copy(str(result))

2. 求解常微分方程通解

$$y^{\prime \prime }-2y^{\prime }+y=2$$

◎ 求解：

$$y\left(x\right)-2\frac{d}{dx}y\left(x\right)+\frac{d^2}{d{x}^2}y\left(x\right)=2$$$$y\left(x\right)=\left(C_1+C_{2}x\right)e^x+2$$

diffeq = Eq(y(x).diff(x,x)-2*y(x).diff(x)+y(x),2)

print_latex(diffeq)
print_latex(dsolve(diffeq, y(x)))

3. 求解级数极限

$$\underset{n\to +\infty }{\lim }\sum _{k=1}^n\frac{1}{n+3k}$$

◎ 数值求解：

▲ 数值求解n→无穷大对应的累加和

当 $n=10000$时对应的取值为：

sumn(10000): 0.46206062271704684

数值求解代码：

from headm import *                 # =

def sumn(n):
    return sum([1/(n+3*k) for k in range(1,n+1)])

ndim = [sumn(i+1) for i in range(1000)]

plt.clf()
plt.figure(figsize=(10,6))
plt.plot(ndim)
plt.xlabel("n")
plt.ylabel("Sum(n)")
plt.grid(True)
plt.tight_layout()
plt.savefig(r"d:\temp\figure1.jpg")
tspshowimage(image=r"d:\temp\figure1.jpg")

printt(sumn(10000)|)

根据 Euler–Maclaurin formula ，可以将有限长累加和转换成对应的Euler-Maclauring公式。

$$\sum _{k=1}^n\frac{1}{n+3k}=\left(\frac{\log \left(4n\right)}{3}-\frac{\log \left(n+3\right)}{3}+\frac{1}{2\left(n+3\right)}+\frac{1}{8n},\left|\frac{1}{4{\left(n+3\right)}^2}-\frac{1}{64n^2}\right|\right)$$

sumf = Sum(1/(n+3*k),(k,1,n))
result=sumf.euler_maclaurin()

所以：$$\underset{n\to +\infty }{\lim }\sum _{k=1}^n\frac{1}{n+3k}=\underset{n\to +\infty }{\lim }\left(\frac{\log \left(4n\right)}{3}-\frac{\log \left(n+3\right)}{3}+\frac{1}{2\left(n+3\right)}+\frac{1}{8n}\right)=\frac{1}{3}\log \left(4\right)$$

这个取值与前面数值计算符合的。

• sympy文献参考：https://docs.sympy.org/latest/modules/concrete.html

4. 定积分求解

$${\int }_0^2\left|1-x\right|dx$$

◎ 求解：

题目所求的定积分实际上是下面函数与x轴包含的面积，根据给定的参数可以知道这个面积为1，所以：$${\int }_0^2\left|1-x\right|dx=1$$

▲ 图1.4.1 定积分图像

5. 求函数高阶导数

$$f\left(x\right)=\sin \left(x^3\right),f^{\left(15\right)}\left(0\right)=?$$

◎ 求解：

$\sin \left(x^3\right)$进行台劳级数展开：
$$x^3-\frac{x^9}{6}+\frac{x^{15}}{120}+O\left(x^{20}\right)$$

所以：$$f^{\left(15\right)}\left(0\right)=\frac{1}{120}$$

result = sin(x**3).series(x, 0, 20)
print_latex(result)

n,k,x = symbols('n,k,x')
result = integrate(abs(1-x), (x,0,2))

6. 求解积分函数导数

$$\frac{d}{dx}{\int }_{x^2}^{x^3}\frac{\sin t}{t}dt$$

◎ 求解：

$$-\frac{2\sin \left(x^2\right)}{x}+\frac{3\sin \left(x^3\right)}{x}$$

n,k,x,t = symbols('n,k,x,t')
result = integrate(sin(t)/t, (t,x**2,x**3))
result = result.diff(x)

7. 定积分求解

$${\int }_0^{\pi }x\cdot {\left(\sin x\right)}^{2}dx$$

◎ 求解：

n,k,x,t = symbols('n,k,x,t')
result = integrate(x*(sin(x)**2),(x,0,pi))

$${\int }_0^{\pi }x\cdot {\left(\sin x\right)}^{2}dx=\frac{{\pi }^2}{4}$$

8. 求解微分方程解

  常微分方程 $y^{\prime }+y=e^{-x}$满足 $y\left(0\right)=0$的解$y=y\left(x\right)$的拐点的横坐标为        。

◎ 求解：

常微分方程的通解为：$$y\left(x\right)=\left(C_1+x\right)e^{-x}$$

n,k,x,t = symbols('n,k,x,t')
y = Function('y')
diffeq = Eq(y(x).diff(x)+y(x), exp(-x))
result = dsolve(diffeq,y(x))

根据 $y\left(0\right)=0$，可以得到：$0=\left(C_1+0\right)\cdot e^{-0}=C_1$，所以常微分方程的完全解为：$$y\left(x\right)=x\cdot e^{-x}$$

改函数对应的拐点$x_0$满足：$y^{\prime \prime }\left(x_0\right)=0$。$$y^{\prime \prime }\left(x\right)=\left(x-2\right)e^{-x}$$

可以得到：$x_0=2$。

9. 求解曲线弧长

  求解下面曲线段的弧长：

$$y=2x^{\frac{3}{2}},\left(0\le x\le 1\right)$$

◎ 求解：

▲ 图1.9.1 曲线函数波形

曲线的$y=y\left(x\right)$对应的弧长公示为：$$L_{0,1}={\int }_0^1\sqrt{1+y^{\prime }{\left(x\right)}^2}dx$$

上面公示的不定积分为：

$$L\left(x\right)=\frac{x\sqrt{9x^2+1}}{2}+\frac{\log \left(x+\sqrt{x^2+1}\right)}{6},x\ge 0$$

定积分：$$L\left(1\right)=\frac{\mathrm{asinh}\left(3\right)}{6}+\frac{\sqrt{10}}{2}=\frac{\log \left(3+\sqrt{10}\right)}{6}+\frac{\sqrt{10}}{2}$$

n,k,x,t = symbols('n,k,x,t')
result = integrate(sqrt(1+(3*x)**2),(x,0,1))

10. 分析级数无穷小阶次

  设当 $x\to 0$时，下面函数为 $p$阶无穷小，则 $p=$        。

$${\left(\frac{\sin x}{x}\right)}^{\frac{1}{3}}-e^{-\frac{x^2}{3}}$$

◎ 求解：

分别将表达式的两项在0点进行台劳级数展开：$${\left(\frac{\sin x}{x}\right)}^{\frac{1}{3}}=1-0.000308641975308642x^{4.0}-0.0555555555555556x^{2.0}+O\left(x^5\right)$$

$$e^{-\frac{x^2}{3}}=1-\frac{x^2}{3}+\frac{x^4}{18}+O\left(x^5\right)$$

它们之间的差对应的无穷小量的阶次 $p=2$。

 

§02 解答题

---

共8题，写出详细的计算过程和必要的根据！

11. （10分）

  讨论 $p$ 取何值时，下面广义积分是收敛的。$${\int }_0^{+\infty }\frac{x^p\ln x}{{\left(1+x^2\right)}^2}dx$$

◎ 求解：

▲ p值不同对应的积分

12.（10分）

  求数列$\left\{n^{1/n}\right\},\left(n=1,2,3\cdots \right)$的最大项的值。

◎ 求解

数量$\left\{n^{1/n}\right\},\left(n=1,2,3\cdots \right)$ 的取值变化如下图所示：

▲ N值对应函数曲线

$$\frac{d}{dn}\left(n^{1/n}\right)=n^{\frac{1}{n}}\left(-\frac{\log \left(n\right)}{n^2}+\frac{1}{n^2}\right)$$

可以看到 $n=e$附近，$n^{1/n}$的取值最大值，综上判断，可以知道序列的最大值出现在$n=3$，对应的最大值为：$3^{1/3}\approx 1.4422$。

13. （13分）

  函数 $$f\left(x\right)=\{\begin{array}{c}{e}^{\frac{1}{x}},x\ne 0\\ 0,x=0\end{array}$$  讨论函数 $f\left(x\right)$的连续性，并求 $f\left(x\right)$的单调区间、极值点与极值、凸性区间、拐点和渐进线。

● 求解：

$f\left(x\right)$在0附近函数图像，可以知道改函数在0点附近不连续。

▲ f(x)在0附近的函数图像

下图是$1/x$在0附近的函数图像。

▲ 1/x在0附近的函数图像

Ⅰ.区间 0 ∞

$$f^{\prime }\left(x\right)=-\frac{e^{\frac{1}{x}}}{x^2}$$

$$f^{\prime \prime }\left(x\right)=\frac{\left(2+\frac{1}{x}\right)e^{\frac{1}{x}}}{x^3}$$

在 $x>0$，$x<0$区间内，函数是单调下降，下凹，渐进性是x轴。

14. （12分）

  设曲线段 $\Gamma$ 为圆心在点（0,1）的单位圆周位于正方形 $0\le x\le 1,0\le y\le 1$的部分，平面区域 $D$为由 $\Gamma$，$x$轴以及直线 $x=1$ 围城的有界区域。

▲ 图2.4.1 曲线段Gama

  （I） 求区域 $D$ 饶 $x$轴旋转一周所产生的旋转体体积；
  （II）求曲线段 $\Gamma$ 绕 $x$ 轴旋转一周所产生的旋转面面积；

● 求解：

Ⅰ.求解旋转体体积

$$V={\int }_0^{\frac{\pi }{2}}\pi {\left(1-\cos \theta \right)}^2\cdot \cos \theta \cdot d\theta =\pi \left(-\frac{\pi }{2}+\frac{5}{3}\right)$$

▲ 图2.4.2 另外一种分析方法

$$V={\int }_0^1\pi {\left[1-{\left(1-x^2\right)}^{0.5}\right]}^{2}dx=\pi \left(-\frac{\pi }{2}+\frac{5}{3}\right)$$

Ⅱ.求解旋转体面积

$$S={\int }_0^{\frac{\pi }{2}}2\pi \cdot \left(1-\cos \left(\theta \right)\right)\cdot dx=2\pi \left(-1+\frac{\pi }{2}\right)$$

$$S={\int }_0^1\frac{2\pi \left[1-\sqrt{1-x^2}\right]}{\sqrt{1-x^2}}dx=2\pi \left(-1+\frac{\pi }{2}\right)$$

15. （10分）

  求常微分方程的初值问题的解 （$x<1$）。

$$\{\begin{array}{c}\sqrt{1+{\left(y^{\prime }\right)}^2}=\left(1-x\right)\cdot y^{\prime \prime }\\ y\left(0\right)=0\\ y^{\prime }\left(0\right)=0\end{array}$$

● 求解：

16. （5分）

  设 $f\in C\left(0,+\infty \right)$， 并且 $\forall a>0,b>1$，都有积分${\int }_a^{ab}f\left(x\right)dx$ 与 $a$ 无关。
  求证：存在常数 $C$，使得 $$f\left(x\right)=\frac{C}{x},x\in \left(0,+\infty \right)$$

17. （5分）

  设 $f\left(x\right)$在 $\left[0,1\right]$ 上非负连续， 且满足 ：$$f^2\left(x\right)\le 1+2{\int }_0^{x}f\left(t\right)dt,x\in \left[0,1\right]$$  证明：$$f\left(x\right)\le 1+x,x\in \left[0,1\right]$$

18. （5分）

  设：$p\left(x\right)=x^n+a_{n-1}{x}^{n-1}+\cdots +a_{1}x+a_0$ 为实系数 $n$ 次多项式。若 $p\left(x\right)\ge 0,x\in \left(-\infty ,+\infty \right)$

  证明：$p\left(x\right)+p^{\prime }\left(x\right)+\cdots +p^{\left(n\right)}\left(x\right)\ge 0,x\in \left(-\infty ,+\infty \right)$

  这里 $p^{\prime }\left(x\right),p^{\prime \prime }\left(x\right),\cdots ,p^{\left(n\right)}\left(x\right)$表示 $p\left(x\right)$的一阶、二阶以及 $n$阶导数。

● 证明：

证明过程请参见： 正值实系数多项式函数所对应的导数累加和是否非负？

 

§03 附加题

---

本题全对才给分，其分数不计入总评，仅用于评判（A+）

1. 附加题内容

  设 $h>0$， $f\left(x\right)$为闭区间$\left[-h,h\right]$上的无穷可导函数，且 $\forall x\in \left[0,h\right]$，以及任意的非负整数 $n$，都有 $f^{\left(n\right)}\left(x\right)\ge 0$。 记 $$r_n\left(x\right)=\frac{1}{n!}{\int }_0^x{\left(x-t\right)}^n{f}^{\left(n+1\right)}\left(t\right)dt$$  求证：$\forall x\in \left(0,h\right)$，均有$\underset{n\to +\infty }{\lim }{r}_n\left(x\right)=0$

● 证明：

本题的证明过程参见： 期末考试中微积分的证明题的分析：Taylor级数展开 。