跳转至

微积分

alt text

积分

几何意义

如果在 [a,b](ab)[a,b](a\neq b) 上函数 f(x)f(x) 连续且恒有 f(x)0f(x)\geq 0,那么定积分 abf(x)dx\displaystyle\int_{a}^{b}f(x)dx 表示由曲线 y=f(x)y=f(x) 以及直线 x=a,x=b,y=0x=a,x=b,y=0 围成的曲边梯形的面积。

f(x)0f(x)\leq 0,那么定积分 abf(x)dx\displaystyle\int_{a}^{b}f(x)dx 表示由曲线 y=f(x)y=f(x) 以及直线 x=a,x=b,y=0x=a,x=b,y=0 围成的曲边梯形的面积的负值

如果我们把 xx 轴上方的面积赋予正号,下方的面积赋予负号,那么在一般情形下,定积分 abf(x)dx\displaystyle\int_{a}^{b}f(x)dx 表示由曲线 y=f(x)y=f(x) 以及直线 x=a,x=b,y=0x=a,x=b,y=0 围成的各部分图形面积的代数和。

微积分基本定理 / NewtonLeibniz\mathrm{Newton-Leibniz} 牛顿 - 莱布尼兹公式

以下设 CC 是一个常数。

如果 f(x)f(x)[a,b][a,b] 上的连续函数且 F(x)=f(x)F'(x)=f(x),那么

abf(x) dx=F(x)ab=F(b)F(a)\boxed{\int_{a}^{b}f(x)\ dx=F(x)|_{a}^{b}=F(b)-F(a) }

我们称 F(x)F(x)f(x)f(x) 的原函数。因为 [F(x)+C]=f(x)[F(x)+C]'=f(x),所以 F(x)+CF(x)+C 也是 f(x)f(x) 的原函数。

常用定积分公式:

abC dx=Cxab=CbCa=C(ba)\int_{a}^{b}C\ dx=Cx|_{a}^{b}=Cb-Ca=C(b-a)
abxn dx=1n+1xn+1ab=1n+1bn+11n+1an+1=bn+1an+1n+1\int_{a}^{b}x^n\ dx=\frac{1}{n+1}x^{n+1}|_{a}^{b}=\frac{1}{n+1}b^{n+1}-\frac{1}{n+1}a^{n+1}=\frac{b^{n+1}-a^{n+1} }{n+1}
absinx dx=(cosx)ab=cosb+cosa\int_a^b\sin x\ dx=(-\cos x)|_a^b=-\cos b+\cos a
abcosx dx=sinxab=sinbsina\int_a^b\cos x\ dx=\sin x|_a^b=\sin b-\sin a
ab1x dx=lnxab=lnblna=lnba\int_a^b\frac{1}{x}\ dx=\ln x|_a^b=\ln b-\ln a=\ln\frac{b}{a}
abex dx=exab=ebea\int_a^b e^x\ dx=e^x|_a^b=e^b-e^a
abnx dx=nxlnnab=nblnbnalna\int_a^b n^x\ dx=\frac{n^x}{\ln n}|_a^b=\frac{n^b}{\ln b}-\frac{n^a}{\ln a}
LRcx+dax+b dx=cxa+(dabca2)×lnax+b LR\int_L^R \frac{cx+d}{ax+b}\ dx=\frac{cx}{a}+(\frac{d}{a}-\frac{bc}{a^2})\times\ln|ax+b|\ |_L^R

定积分的基本性质

  1. abCf(x) dx=Cabf(x) dx\int_a^b Cf(x)\ dx=C\int_a^b f(x)\ dx
  2. ab[f(x)±g(x)] dx=abf(x) dx±abg(x) dx\int_a^b[f(x)\pm g(x)]\ dx=\int_a^b f(x)\ dx\pm\int_a^b g(x)\ dx
  3. abf(x) dx=acf(x) dx+cbf(x) dx\int_a^b f(x)\ dx=\int_a^c f(x)\ dx+\int_c^b f(x)\ dx
  4. 在区间 [a,b][a,b] 上满足 f(x)0f(x)\geq 0,则 $abf(x) dx0\int_a^b f(x)\ dx\geq 0$

  5. 在区间 [a,b][a,b] 上满足 f(x)g(x)f(x)\leq g(x),则 $abf(x) dxabg(x) dx\int_a^b f(x)\ dx\leq\int_a^b g(x)\ dx$

  6. abf(x) dxabf(x) dx\left|\int_a^bf(x)\ dx\right |\leq\int_a^b|f(x)|\ dx
  7. f(x)f(x) 是偶函数,且在 [a,a][-a,a] 上连续,则 $aaf(x) dx=20af(x) dx\int_{-a}^af(x)\ dx=2\int_0^af(x)\ dx$

  8. f(x)f(x) 是奇函数,且在 [a,a][-a,a] 上连续,则 $aaf(x) dx=0\int_{-a}^af(x)\ dx=0$

不定积分

我们现在需要一种简单的表示反导数的方式。根据微积分基本定理,我们可以用 f(x) dx\displaystyle\int f(x)\ dx 表示“函数 ff 的反导数的集合”,注意任何可积函数都有无数个反导数,唯一不同的是常数部分。例如,

x2 dx=x33+C\int x^2\ dx=\frac{x^3}{3}+C

对于任意常数 CC 都成立。也就是说,若 F(x)=f(x)F'(x)=f(x),则

f(x) dx=F(x)+C\int f(x)\ dx=F(x)+C

不定积分的性质同定积分的性质。

换元法

  • 第一换元积分法:令 y=g(x)y=g(x) 有:
abf(g(x))g(x) dx=g(a)g(b)f(y) dy\red{\boxed{\int_a^bf(g(x))g'(x)\ dx=\int_{g(a)}^{g(b)}f(y)\ dy}}
  1. 2sin2x dx=u=2xsinu du=cosu+C=cos2x+C    (du=2dx)\int 2\sin 2x\ dx\xlongequal{u=2x}\int\sin u\ du=-\cos u+C=-\cos 2x+C\ \ \ \ (du=2dx)
  2. (2x+1)5 dx=u=2x+1u5 du=112u6+C=112(2x+1)6+C    (du=2dx)\int (2x+1)^5\ dx\xlongequal{u=2x+1}\int u^5\ du=\frac{1}{12}u^6+C=\frac{1}{12}(2x+1)^6+C\ \ \ \ (du=2dx)
  3. 2x(x2+1)3 dx=(x2+1)32x dx=u=x2+1u3 du=12u2+C=12(x2+1)+C    (du=2xdx)\begin{aligned}\int\frac{2x}{(x^2+1)^3}\ dx&=\int(x^2+1)^{-3}\cdot 2x\ dx\xlongequal{u=x^2+1}\int u^{-3}\ du\\&=-\frac{1}{2}u^{-2}+C=-\frac{1}{2(x^2+1)}+C\ \ \ \ (du=2xdx)\end{aligned}
  4. x1x2 dx=121x2(2x dx)=u=1x212u12 du=13u32+C=13(1x2)32+C    (du=2dx)\begin{aligned}\int x\sqrt{1-x^2}\ dx&=-\frac{1}{2}\int\sqrt{1-x^2}(-2x\ dx)\xlongequal{u=1-x^2}-\frac{1}{2}\int u^{\frac{1}{2}}\ du\\&=-\frac{1}{3}u^{\frac{3}{2}}+C=-\frac{1}{3}(1-x^2)^{\frac{3}{2}}+C\ \ \ \ (du=-2dx)\end{aligned}
  5. 1x2+a2 dx=1a21x2a2+1 dx=1a11+(xa)2 d(xa)=1aarctanxa+C\int\frac{1}{x^2+a^2}\ dx=\frac{1}{a^2}\int\frac{1}{\frac{x^2}{a^2}+1}\ dx=\frac{1}{a}\int\frac{1}{1+(\frac{x}{a})^2}\ d(\frac{x}{a})=\frac{1}{a}\arctan\frac{x}{a}+C
  6. 01e2x dx=0212ey dy=12(e2e0)=12(e21)\int_0^1e^{2x}\ dx=\int_0^2\frac{1}{2}e^y\ dy=\frac{1}{2}(e^2-e^0)=\frac{1}{2}(e^2-1)
  7. 第二换元积分法:令 x=g(t),dx=g(t)dtx=g(t),dx=g'(t)dt 有:

f(x) dx=f[g(t)]g(t) dt\red{\boxed{\int f(x)\ dx=\int f[g(t)]g'(t)\ dt}}
  1. 11+x dx=x=t211+t2t dt=2(111+t) dt=2t2ln1+t+C=2x2ln1+x+C\begin{aligned}\int\frac{1}{1+\sqrt x}\ dx&\xlongequal{x=t^2}\int\frac{1}{1+t}\cdot 2t\ dt=2\int(1-\frac{1}{1+t})\ dt=2t-2\ln|1+t|+C\\&=2\sqrt{x}-2\ln|1+\sqrt{x}|+C\end{aligned}
  2. a2x2 dx=x=asintacostacost dt=a2cos2t dt=a21+cos2t2dt=a2(12t+14sin2t)+C=a2(12arcsinxa+142a2xa2x2)+C=a22arcsinxa+x2a2x2+C\begin{aligned}\int\sqrt{a^2-x^2}\ dx&\xlongequal{x=a\sin t}\int a\cos t\cdot a\cos t\ dt=a^2\int \cos^2 t\ dt\\&=a^2\int\frac{1+\cos 2t}{2}dt=a^2(\frac{1}{2}t+\frac{1}{4}\sin 2t)+C\\&=a^2(\frac{1}{2}\arcsin\frac{x}{a}+\frac{1}{4}\frac{2}{a^2}x\sqrt{a^2-x^2})+C\\&=\frac{a^2}{2}\arcsin\frac{x}{a}+\frac{x}{2}\sqrt{a^2-x^2}+C\end{aligned}

分部积分法

(uv)=uv+uv    uv=(uv)vu    uv dx=uvvu dx(uv)'=u'v+uv'\implies uv'=(uv)'-vu'\implies \int uv'\ dx=uv-\int vu'\ dx
u dv=uvv du\red{\boxed{\int u\ dv=uv-\int v\ du}}
  1. 计算 xcosx dx\displaystyle\int x\cos x\ dx,设 u=x,dv=cosx dxu=x,dv=\cos x\ dx,则 du=dx,v=sinxdu=dx,v=\sin x $xcosx dx=xsinxsinx dx=xsinx+cosx+C\int x\cos x\ dx=x\sin x-\int \sin x\ dx=x\sin x+\cos x+C$

  2. 计算 arctanx dx\displaystyle\int\arctan x\ dx,设 u=arctanx,dv=dxu=\arctan x,dv=dx,则 du=11+x2 dx,v=x\displaystyle du=\frac{1}{1+x^2}\ dx,v=x $arctanx=xarctanx11+x2 dx=xarctanx1211+x2d(1+x2)=xarctanx12ln(1+x2)+C\begin{aligned}\int \arctan x&=x\arctan x-\int\frac{1}{1+x^2}\ dx=x\arctan x-\frac{1}{2}\int\frac{1}{1+x^2}d(1+x^2)\\&=x\arctan x-\frac{1}{2}\ln(1+x^2)+C\end{aligned}$

  3. 计算 xlnx dx\displaystyle\int x\ln x\ dx,设 u=lnx,dv=x dxu=\ln x,dv=x\ dx,则 du=1x dx,v=12x2\displaystyle du=\frac{1}{x}\ dx,v=\frac{1}{2}x^2 $xlnx dx=12x2lnx1x12x2 dx=12x2lnx14x2+C\int x\ln x\ dx=\frac{1}{2}x^2\ln x-\int\frac{1}{x}\cdot\frac{1}{2}x^2\ dx=\frac{1}{2}x^2\ln x-\frac{1}{4}x^2+C$

Gamma 函数(第二类欧拉积分)

欧拉发现当 nn 为整数时,0+tnet dt=n!\displaystyle\int_0^{+\infty}t^ne^{-t}\ dt=n!,于是我们把 nn 替换成实数,就可变为

Γ(s)=0+ts1et dt\Gamma(s)=\int_0^{+\infty}t^{s-1}e^{-t}\ dt
  1. ss 为非正整数的复数 ss 时,Γ(s+1)=sΓ(s)\Gamma(s+1)=s\Gamma(s)
  2. nn 为整数时,Γ(n+1)=n!\Gamma(n+1)=n!
  3. 余元公式(欧拉反射原理):s(0,1)s\in(0,1) 时,Γ(s)Γ(1s)=πsinπs\displaystyle\Gamma(s)\Gamma(1-s)=\frac{\pi}{\sin\pi s},得 Γ(12)=π\displaystyle\Gamma\left(\frac{1}{2}\right)=\sqrt\pi

应用

  1. 计算 y=exy=e^xx=0x=0x=1x=1 之间与 xx 轴围成的曲边梯形的面积。

只需计算 01ex dx=e1e0=e1\displaystyle\int_0^1e^x\ dx=e^1-e^0=e-1

  1. 计算 x=y2x=y^2x=1x=1 之间围成的图形的面积。

只需计算 01(x(x)) dx=01(2x) dx=43x301=430=43\displaystyle\int_0^1(\sqrt{x}-(-\sqrt{x}))\ dx=\int_0^1(2\sqrt{x})\ dx=\frac{4}{3}\sqrt{x^3}|_0^1=\frac{4}{3}-0=\frac{4}{3}

  1. 【2025 GD 一模 T19】如果函数 F(x)F(x) 的导数为 F(x)=f(x)F'(x)=f(x),可记为 f(x) dx=F(x)\displaystyle\int f(x)\ dx=F(x),若 f(x)0f(x)\geq 0,则 abf(x) dx=F(b)F(a)\displaystyle\int_a^bf(x)\ dx=F(b)-F(a) 表示曲线 y=f(x)y=f(x),直线 x=a,x=bx=a,x=b 以及 xx 轴围成的“曲边梯形”的面积。如:2x dx=x2+C\displaystyle\int 2x\ dx=x^2+C,其中 CC 为常数;022x dx=(22+C)(0+C)=4\displaystyle\int_0^2 2x\ dx=(2^2+C)-(0+C)=4,则表示 x=0,x=1,y=2x+Cx=0,x=1,y=2x+C 以及 xx 轴围成的面积为 44. (1)(1)f(x)=(ex+1) dx,f(0)=2f(x)=\displaystyle\int(e^x+1)\ dx,f(0)=2,求 f(x)f(x) 的表达式。 (2)(2) 求曲线 y=x2y=x^2 与直线 y=x+6y=-x+6 所围成图形的面积。 (3)(3)f(x)=ex12mx,x[0,+)f(x)=e^x-1-2mx,x\in[0,+\infty),其中 mR,a,b[0,+)m\in\R,\forall a,b\in[0,+\infty),若 a>ba>b,都满足 0af(x) dx>0bf(x) dx\displaystyle\int_0^a f(x)\ dx>\int_0^b f(x)\ dx,求 mm 取值范围。

解:(1) f(x)=ex+x+1(1)\ f(x)=e^x+x+1.

(2)(2) 先解方程 x2=x+6x^2=-x+6 得到 x=3x=-322.

答案即为 32(x+6x2) dx=(12x2+6x13x3)32=223(272)=1256\displaystyle\int_{-3}^2 (-x+6-x^2)\ dx=(-\frac{1}{2}x^2+6x-\frac{1}{3}x^3)|_{-3}^2=\frac{22}{3}-(-\frac{27}{2})=\frac{125}{6}.

(3)(3) 由题意可知,a,b[0,+),a>b\forall a,b\in[0,+\infty),a>b 满足 F(a)F(0)>F(b)F(0)F(a)-F(0)>F(b)-F(0),即 F(a)F(a)\uparrow,进而 f(x)0f(x)\geq 0[0,+)[0,+\infty) 恒成立。

接下来就是常规导数了。分离参数 mex12x=g(x),g(x)=(x1)ex+12x2m\leq\frac{e^x-1}{2x}=g(x),g'(x)=\frac{(x-1)e^x+1}{2x^2},令 h(x)=(x1)ex,h(x)=xex0h(x)=(x-1)e^x,h'(x)=xe^x\geq 0,即 h(x),g(x),g(x)h(x)\uparrow,g'(x)\uparrow,g(x)\uparrow,由洛必达法则,mlimx0ex12x=limx0ex2=12\displaystyle m\leq\lim_{x\to 0}\frac{e^x-1}{2x}=\lim_{x\to 0}\frac{e^x}{2}=\frac{1}{2}.

微积分在数列中的应用

秒杀神秘放缩的终极方法。但是需要一些高等知识,比如说解简单的常微分方程。但是一般的方程一眼看出来就可以了,这里给出通法:

对于一阶齐次线性微分方程:

y+P(x)y=0y'+P(x)y=0

解为

y=CeP(x)dxy=Ce^{-\int P(x)dx}

对于一阶非齐次线性微分方程:

y+P(x)y=Q(x)y'+P(x)y=Q(x)

先解 y+P(x)y=0y'+P(x)y=0 得到 y=CeP(x)dxy=Ce^{-\int P(x)dx},再令 C=u(x)C=u(x)

y=u(x)eP(x)dxy=u(x)e^{-\int P(x) dx}

代入原方程得

u(x)=Q(x)eP(x)dxu'(x)=\frac{Q(x)}{e^{-\int P(x)dx}}

u(x)u'(x) 积分带入得到通解:

y=CeP(x)dx+eP(x)dxQ(x)eP(x)dxdxy=Ce^{-\int P(x)dx}+e^{-\int P(x)dx}\int Q(x)e^{\int P(x)dx}dx

以上的 CC 都为常数,由初始条件决定。


a1=1,an+1=an+1an\displaystyle a_1=1,a_{n+1}=a_n+\frac{1}{a_n},求 a2024a_{2024} 的整数部分?

ana_n 近似为一个连续函数 f(x)f(x),考虑 an+1ana_{n+1}-a_n1an\frac{1}{a_n} 的关系。对于 n+n\to +\infty,有

f(x)=1f(x)    df(x)dx=1f(x)    f(x)df(x)=dx    f(x)df(x)=dx    f2(x)2=x+Cf'(x)=\frac{1}{f(x)}\implies\frac{\mathrm{d}f(x)}{\mathrm{d}x}=\frac{1}{f(x)} \implies f(x)\mathrm{d}f(x)=\mathrm{d}x\\\implies\int f(x)\mathrm{d}f(x)=\int\mathrm{d}x\implies \frac{f^2(x)}{2}=x+C

f(1)=1f(1)=1 得到 C=12\displaystyle C=-\frac{1}{2},微分方程解得 f(x)=2x1f(x)=\sqrt{2x-1},得到 f(2024)=404763f(2024)=\sqrt{4047}\approx 63


已知 a1=1,an+1=an13an2\displaystyle a_1=1,a_{n+1}=a_n-\frac{1}{3}a_n^2,则 100a100100a_{100} 的范围是哪个:(2,52),(52,3),(3,72)\displaystyle\left(2,\frac{5}{2}\right),\left(\frac{5}{2},3\right), \left(3,\frac{7}{2}\right)

1an+1=3an(3an)=1an+13an>1an+13\frac{1}{a_{n+1}}=\frac{3}{a_n(3-a_n)}=\frac{1}{a_n}+\frac{1}{3-a_n}>\frac{1}{a_n}+\frac{1}{3}
1an>a1+n13=n+23    an<3n+2    100a100<3\frac{1}{a_n}>a_1+\frac{n-1}{3}=\frac{n+2}{3}\implies a_n<\frac{3}{n+2}\implies 100a_{100}<3
1an+1<1an+133n+2=1an+n+23(n+1)=1an+13+13(n+1)\frac{1}{a_{n+1}}<\frac{1}{a_n}+\frac{1}{3-\frac{3}{n+2}}=\frac{1}{a_n}+\frac{n+2}{3(n+1)}=\frac{1}{a_n}+\frac{1}{3}+\frac{1}{3(n+1)}
1a100<1a1+993+13(12+13++1100)=34+13(12+13++1100)\frac{1}{a_{100}}<\frac{1}{a_1}+\frac{99}{3}+\frac{1}{3}\left(\frac{1}{2}+\frac{1}{3}+\dots+\frac{1}{100}\right)=34+\frac{1}{3}\left(\frac{1}{2}+\frac{1}{3}+\dots+\frac{1}{100}\right)

注意到

12+13++1100<12×4+16×96=18\frac{1}{2}+\frac{1}{3}+\dots+\frac{1}{100}<\frac{1}{2}\times 4+\frac{1}{6}\times 96=18
1a100<140    100a100>52\frac{1}{a_{100}}<\frac{1}{40}\implies 100a_{100}>\frac{5}{2}

这道题是我们班的补充题,当时无人证出 >52\displaystyle >\frac{5}{2},我们老师上课思考了不到 1515 分钟讲了这种方法(他从不备课),注意力简直惊为天人。

积分法:f(x)=an,f(x)=df(x)dx=13f2(x)    f2(x)df(x)=13dxf1(x)=13x+C,f(x)=3x+C    3x+2\displaystyle f(x)=a_n,f'(x)=\frac{df(x)}{dx}=-\frac{1}{3}f^2(x)\implies \int f^{-2}(x)df(x)=\int-\frac{1}{3}dx\\-f^{-1}(x)=-\frac{1}{3}x+C,f(x)=\frac{3}{x+C}\implies\frac{3}{x+2},显然 100a100(52,3)\displaystyle 100a_{100}\in\left(\frac{5}{2},3\right). 可以代入 a2,a3a_2,a_3 获得更高的精度。


已知 a1=1,an+1=an1+an3\displaystyle a_1=1,a_{n+1}=\frac{a_n}{1+\sqrt[3]{a_n}},求 S100S_{100} 的整数部分。

答案:准确值是 2.382.38.

1an+1=1an+1an23=bn=1anbn+1=bn+bn23=f(x)=bnf(x)=f23(x)\frac{1}{a_{n+1}}=\frac{1}{a_n}+\frac{1}{a_n^{\frac{2}{3}}}\xlongequal{b_n=\frac{1}{a_n}}b_{n+1}=b_n+b_n^{\frac{2}{3}}\xlongequal{f(x)=b_n}f'(x)=f^{\frac{2}{3}}(x)
f23(x)df(x)=dx    3f13(x)=x+C    f(x)=(x+23)3=(x2+3233)3\int f^{-\frac{2}{3}}(x) df(x)=\int dx\implies 3f^{\frac{1}{3}}(x)=x+C\implies f(x)=\left(\frac{x+2}{3}\right)^3=\left(\frac{x-2+3\sqrt[3]{2}}{3}\right)^3
an=27(x+2)3,ζ(3)=n=11n31.202,S=270.077=2.079a_n=\frac{27}{(x+2)^3},\zeta(3)=\sum_{n=1}^{\infty}\frac{1}{n^3}\approx 1.202,S_{\infty}=27\cdot 0.077=2.079

另外,欧拉常数 γ=limn[(k=1n1k)lnn]=1(1x1x)0.57721566490153286060\displaystyle\gamma=\lim_{n\to\infty}\left[\left(\sum_{k=1}^n\frac{1}{k}\right)-\ln n\right]=\int_1^{\infty}\left(\frac{1}{\lfloor x\rfloor}-\frac{1}{x}\right)\approx 0.57721 56649 01532 86060,于是可以估计 k=1n1klnn+γ\displaystyle\sum_{k=1}^n\frac{1}{k}\approx\ln n+\gamma.


也可以用一个结论(Basel 问题,知乎专栏 ):n=1+1n2=π261.645\red{\boxed{\sum_{n=1}^{+\infty}\frac{1}{n^2}=\frac{\pi^2}{6}\approx 1.645}}

相关证明:设黎曼 Zeta 函数 ζ(s)=n=1+1ns (s>1)\zeta(s)=\displaystyle\sum_{n=1}^{+\infty}\frac{1}{n^s}\ (s>1),显然 lims1+ζ(s)=+,lims+ζ(s)=1\displaystyle\lim_{s\to 1^+}\zeta(s)=+\infty,\lim_{s\to+\infty}\zeta(s)=1.

由泰勒展开得 sinx=xx36+x5120+=k=0+(1)k(2k+1)!x2k+1\displaystyle\sin x=x-\frac{x^3}{6}+\frac{x^5}{120}+\dots=\sum_{k=0}^{+\infty}\frac{(-1)^k}{(2k+1)!}x^{2k+1}

因此 sinxx=1x23!+x45!+=n=0+(1)nx2n(2n+1)!\displaystyle\frac{\sin x}{x}=1-\frac{x^2}{3!}+\frac{x^4}{5!}+\dots=\sum_{n=0}^{+\infty}\frac{(-1)^nx^{2n}}{(2n+1)!}①sinxx\displaystyle\frac{\sin x}{x} 的无穷级数展开)

熟知 sinxx\displaystyle\frac{\sin x}{x} 的无穷乘积展开 sinxx=n=1+(1x2n2π2)=(1x2π2)(1x24π2)\displaystyle\frac{\sin x}{x}=\prod_{n=1}^{+\infty}\left(1-\frac{x^2}{n^2\pi^2}\right)=\left(1-\frac{x^2}{\pi^2}\right)\left(1-\frac{x^2}{4\pi^2}\right)\cdots②

Basel 问题即 ζ(2)=n=1+1n2\displaystyle\zeta(2)=\sum_{n=1}^{+\infty}\frac{1}{n^2},比较 ①②①② 式中二次项系数得 13!=1π2122π2\displaystyle -\frac{1}{3!}=-\frac{1}{\pi^2}-\frac{1}{2^2\pi^2}-\dots

由此可得 ζ(2)=π26\displaystyle\zeta(2)=\frac{\pi^2}{6},进一步地我们有 ζ(4)=π490,ζ(6)=π6945\displaystyle\zeta(4)=\frac{\pi^4}{90},\zeta(6)=\frac{\pi^6}{945}.

ζ(2n)=[k=1n1(1)k+1π2kζ(2n2k)(2k+1)!]+(1)n+1nπ2n(2n+1)!      nN\zeta(2n)=\left[\sum_{k=1}^{n-1}\frac{(-1)^{k+1}\pi^{2k}\zeta(2n-2k)}{(2k+1)!}\right]+\frac{(-1)^{n+1}n\pi^{2n}}{(2n+1)!}\ \ \ \ \ \ n\in\N^*

上述思路可总结出一个结论:若 k=0+akxk=k=1+(1+bkx)\displaystyle\sum_{k=0}^{+\infty}a_kx^k=\prod_{k=1}^{+\infty}(1+b_kx),设 Sn=k=1+bknS_n=\displaystyle\sum_{k=1}^{+\infty}b_k^n,则 $S1=a1,   Sn=k=1n1(1)k+1akSnk+(1)n+1nanS_1=a_1,\ \ \ S_n=\displaystyle\sum_{k=1}^{n-1}(-1)^{k+1}a_kS_{n-k}+(-1)^{n+1}na_n$

积分放缩在数列中的应用

我们熟知的积分放缩一般是这样的(证明:几何意义)

1n+1f(x) dx<k=1nf(k)<0nf(x) dx        f(x)\int_1^{n+1}f(x)\ dx <\sum_{k=1}^nf(k)<\int_0^nf(x)\ dx\ \ \ \ \ \ \ \ f(x)\downarrow
0nf(x) dx<k=1nf(k)<1n+1f(x) dx        f(x)\int_0 ^{n}f(x)\ dx<\sum_{k=1}^nf(k)<\int_1^{n+1}f(x)\ dx\ \ \ \ \ \ \ \ f(x)\uparrow

例如:

k=1n1k>1n+11x dx=2x1n+1=2(n+11)<0n1x dx=2x0n=2n\begin{aligned}\sum_{k=1}^n\frac{1}{\sqrt{k} }&>\int_1^{n+1}\frac{1}{\sqrt{x} }\ dx=2\sqrt{x}|_1^{n+1}=2\left(\sqrt{n+1}-1\right)\\&<\int_0^n\frac{1}{\sqrt{x} }\ dx=2\sqrt{x}|_0^n=2\sqrt{n}\end{aligned}

练习 1:证明 k=2n1kk<2    \displaystyle\sum_{k=2}^n\frac{1}{k\sqrt{k}}<\sqrt{2}\ \ \ \ \text{}    23n32k=1nk23[(n+1)32n]\ \ \ \ \displaystyle\frac{2}{3}n^\frac{3}{2}\leq\sum_{k=1}^n\sqrt{k}\leq\frac{2}{3}[(n+1)^\frac{3}{2}-n]

练习 2:证明 k=2n+11k<ln(n+1)<k=1n1k    \displaystyle\sum_{k=2}^{n+1}\frac{1}{k}<\ln(n+1)<\sum_{k=1}^n\frac{1}{k}\ \ \ \ \text{}    47<k=n+12n1k<ln2<2536<22\ \ \ \ \displaystyle\frac{4}{7}<\sum_{k=n+1}^{2n}\frac{1}{k}<\ln 2<\frac{25}{36}<\frac{\sqrt{2}}{2}

练习 3:证明 {k=1n1kp<pp1p>1k=1n1kpk=1n1k>lnnp1\displaystyle\begin{cases}\displaystyle\sum_{k=1}^n\frac{1}{k^p}<\frac{p}{p-1}&p>1\\\displaystyle\sum_{k=1}^n\frac{1}{k^p}\geq\displaystyle\sum_{k=1}^n\frac{1}{k}>\ln n&p\leq 1\end{cases}

我们还可以根据 f(x)f(x) 的凹凸性进一步放缩。

f(x)f(x)\downarrow 且下凸时,

1n+1f(x) dx+f(1)f(n+1)2<k=1nf(k)<0nf(x) dx\int_1^{n+1}f(x)\ dx+\frac{f(1)-f(n+1)}{2}<\sum_{k=1}^nf(k)<\int_0^nf(x)\ dx

f(x)f(x)\downarrow 且上凸时,

1n+1f(x) dx<k=1nf(k)<0nf(x) dxf(0)f(n)2\int_1^{n+1}f(x)\ dx<\sum_{k=1}^nf(k)<\int_0^nf(x)\ dx-\frac{f(0)-f(n)}{2}

f(x)f(x)\uparrow 且下凸时,

0nf(x) dx+f(n)f(0)2<k=1nf(k)<1n+1f(x) dx\int_0^nf(x)\ dx+\frac{f(n)-f(0)}{2}<\sum_{k=1}^nf(k)<\int_1^{n+1}f(x)\ dx

f(x)f(x)\uparrow 且上凸时,

0nf(x) dx<k=1nf(k)<1n+1f(x) dxf(n+1)f(1)2\int_0^nf(x)\ dx<\sum_{k=1}^nf(k)<\int_1^{n+1}f(x)\ dx-\frac{f(n+1)-f(1)}{2}