今天,收到G.Han的提问,第一个是计算积分

\[\int_0^{\infty}{\frac{\ln x}{(x^2+1)^n}dx}\]顿时不明觉厉,然后在宝典《Table of Integrals, Series, and Products》上找到一个更一般的结果:

\[\int_0^\infty  {\ln x\frac{{dx}}{{{{\left( {{a^2} + {b^2}{x^2}} \right)}^n}}}}  = \frac{{\Gamma \left( {n - \frac{1}{2}} \right)\sqrt \pi  }}{{4\left( {n - 1} \right)!{a^{2n - 1}}b}}\left[ {2\ln \frac{a}{{2b}} - \gamma - \psi \left( {n - \frac{1}{2}} \right)} \right]\qquad \text{$a>0,b>0$}.\]

其中$\gamma$为Euler-Mascheroni常数,$\psi(x)$为Digamma 函数,有:

\begin{align}\psi\left(n+\frac{1}{2} \right)&=-\gamma-2\ln2+\sum\limits_{k=1}^n{\frac{2}{2k-1}}.\\\psi\left(\frac{1}{2}\right)&=-\gamma-2\ln2\end{align}

解:

\begin{align*}{I_n} &= \int_0^{ + \infty } {\frac{{\ln x}}{{{{\left( {{x^2} + 1} \right)}^n}}}dx}  = \int_0^{ + \infty } {\frac{1}{{{{\left( {{x^2} + 1} \right)}^n}}}d\left( {x\ln x - x} \right)}\\&= \left[ {\frac{{x\ln x - x}}{{{{\left( {{x^2} + 1} \right)}^n}}}} \right]_0^{ + \infty } + 2n\int_0^{ + \infty }{\frac{{{x^2}\ln x - {x^2}}}{{{{\left( {{x^2} + 1} \right)}^{n + 1}}}}dx}  \\&= 2n\int_0^{ + \infty } {\frac{{\left( {{x^2} + 1} \right)\ln x - \ln x - {x^2}}}{{{{\left( {{x^2} + 1} \right)}^{n + 1}}}}dx}  \\&= 2n\int_0^{ + \infty } {\frac{{\ln x}}{{{{\left( {{x^2} + 1} \right)}^n}}}dx}  - 2n\int_0^{ + \infty } {\frac{{\ln x}}{{{{\left( {{x^2} + 1} \right)}^{n + 1}}}}dx} - 2n\int_0^{ + \infty } {\frac{{{x^2}}}{{{{\left( {{x^2} + 1} \right)}^{n + 1}}}}dx}\end{align*}
\begin{align*}\Rightarrow 2n{I_{n + 1}} &= \left( {2n - 1} \right){I_n} - 2n\int_0^{ + \infty } {\frac{{{x^2}}}{{{{\left( {{x^2} + 1} \right)}^{n + 1}}}}dx}  \\&= \left( {2n - 1} \right){I_n} - 2n\int_0^{ + \infty } {\frac{{{x^2} + 1 - 1}}{{{{\left( {{x^2} + 1} \right)}^{n + 1}}}}dx}  \\&= \left( {2n - 1} \right){I_n} - 2n\int_0^{ + \infty } {\frac{1}{{{{\left( {{x^2} + 1} \right)}^n}}}dx}  + 2n\int_0^{ + \infty } {\frac{1}{{{{\left( {{x^2} + 1} \right)}^{n + 1}}}}dx} .\end{align*}
\[2n{I_{n + 1}} = \left( {2n - 1} \right){I_n} - 2n\int_0^{ + \infty } {\frac{1}{{{{\left( {{x^2} + 1}\right)}^n}}}dx}  + 2n\int_0^{ + \infty } {\frac{1}{{{{\left( {{x^2} + 1} \right)}^{n + 1}}}}dx} .\]
先证明一个引理
\begin{align*}\int_0^{ + \infty } {\frac{1}{{{{\left( {{x^2} + 1} \right)}^n}}}dx} &= \frac{1}{2}\mathbf{B}\left( {n - \frac{1}{2},\frac{1}{2}} \right) = \frac{{\sqrt \pi  }}{2}\frac{{\mathbf{\Gamma} \left( {n - \frac{1}{2}} \right)}}{{\left( {n - 1} \right)!}} \\&= \left\{ \begin{array}{l}\frac{{\sqrt \pi  }}{{2\left( {n - 1} \right)!}}\left[ {\sqrt \pi  .\frac{{\left( {2n - 3} \right)!!}}{{{2^{n - 1}}}}} \right] = \frac{{\left( {2n - 3} \right)!!}}{{{2^n}\left( {n - 1} \right)!}}\pi \qquad n \ge 2\\\frac{\pi }{2} \qquad n = 1\end{array} \right..\end{align*}
引理的证明
\begin{align*}\int_0^{ + \infty } {\frac{1}{{{{\left( {{x^2} + 1} \right)}^n}}}dx} &\underline{\underline {\text{令}t = \frac{1}{{{x^2} + 1}}}} \int_1^0 {{t^n}d\sqrt {\frac{1}{t} - 1} }  = \frac{1}{2}\int_0^1 {{t^{n - \frac{3}{2}}}{{\left( {1 - t} \right)}^{ - \frac{1}{2}}}dt}  \\&= \frac{1}{2}\mathbf{B}\left( {n - \frac{1}{2},\frac{1}{2}} \right) = \frac{{\sqrt \pi  }}{2}\frac{{\mathbf{\Gamma} \left( {n - \frac{1}{2}} \right)}}{{\left( {n - 1} \right)!}} \\&= \left\{ \begin{array}{l}\frac{{\sqrt \pi  }}{{2\left( {n - 1} \right)!}}\left[ {\sqrt \pi  .\frac{{\left( {2n - 3} \right)!!}}{{{2^{n - 1}}}}} \right] = \frac{{\left( {2n - 3} \right)!!}}{{{2^n}\left( {n - 1} \right)!}}\pi \qquad n \ge 2\\\frac{\pi }{2} \qquad n = 1\end{array} \right..\end{align*}
回到原题
\begin{align*}2n{I_{n + 1}} &= \left( {2n - 1} \right){I_n} - 2n\int_0^{ + \infty } {\frac{1}{{{{\left( {{x^2} + 1} \right)}^n}}}dx}  + 2n\int_0^{ + \infty } {\frac{1}{{{{\left( {{x^2} + 1} \right)}^{n + 1}}}}dx}  \\&= \left( {2n - 1} \right){I_n} - 2n\frac{{\left( {2n - 3} \right)!!}}{{{2^n}\left( {n - 1} \right)!}}\pi  + 2n\frac{{\left( {2n - 1} \right)!!}}{{{2^{n + 1}}n!}}\pi \\&= \left( {2n - 1} \right){I_n} - \frac{\pi }{{\left( {n - 1} \right)!}}\frac{{\left( {2n - 3} \right)!!}}{{{2^n}}}\\&\Rightarrow \frac{{\left( {2n} \right)!!}}{{\left( {2n - 1} \right)!!}}{I_{n + 1}} = \frac{{\left( {2n - 2} \right)!!}}{{\left( {2n - 3} \right)!!}}{I_n} - \frac{\pi }{2} \cdot \frac{1}{{2n - 1}} \\&\Rightarrow \frac{{\left( {2n - 2} \right)!!}}{{\left( {2n - 3} \right)!!}}{I_n} = \frac{{2!!}}{{1!!}}{I_2} - \sum\limits_{k = 2}^{n - 1} {\frac{\pi }{{2k - 1}}}  =  - \frac{\pi}{2}\sum\limits_{k = 1}^{n - 1} {\frac{1 }{{2k - 1}}} \\&\Rightarrow {I_n} = \left\{ \begin{array}{l}0 \qquad n = 1\\- \frac{\pi}{2}\frac{{\left( {2n - 3} \right)!!}}{{\left( {2n - 2} \right)!!}}\sum\limits_{k = 1}^{n - 1} {\frac{1}{{2k - 1}}} \qquad n \ge 2\end{array} .\right.\end{align*}
 
 

第二题是个重要的Steffensen积分不等式

设$f,g\in \mathcal{R}[a,b]$,且$f$在$[a,b]$单减,$0<g(x)\leq1$,求证:

\[\int_{b-\lambda}^b{f(x)dx}\leq\int_a^b{f(x)g(x)dx}\leq\int_a^{a+\lambda}{f(x)dx}.\]

其中\[\lambda=\int_a^b{g(x)dx}.\]

证:先证明右边不等式

\begin{align*}&\int_a^{a + \lambda } {f\left( x \right)dx} - \int_a^b {f\left( x \right)g\left( x \right)dx} = \int_a^{a + \lambda } {f\left( x \right)\left[ {1 - g\left( x \right)} \right]dx} - \int_{a + \lambda }^b {f\left( x \right)g\left( x \right)dx}\\ &\ge f\left( {a + \lambda } \right)\int_a^{a + \lambda } {\left[ {1 - g\left( x \right)} \right]dx} = f\left( {a + \lambda } \right)\left( {\lambda - \int_a^{a + \lambda } {g\left( x \right)dx} } \right) - \int_{a + \lambda }^b {f\left( x \right)g\left( x \right)dx} \\&= f\left( {a + \lambda } \right)\int_{a + \lambda }^b {g\left( x \right)dx} - \int_{a + \lambda }^b {f\left( x \right)g\left( x \right)dx} = \int_{a + \lambda }^b {\left[ {f\left( {a + \lambda } \right) - f\left( x \right)} \right]g\left( x \right)dx} \ge 0\end{align*}

左边不等式同理可证.

积分题1之来自G.Han的一道积分题的更多相关文章

  1. 玉伯的一道课后题题解(关于 IEEE 754 双精度浮点型精度损失)

    前文 的最后给出了玉伯的一道课后题,今天我们来讲讲这题的思路. 题目是这样的: Number.MAX_VALUE + 1 == Number.MAX_VALUE; Number.MAX_VALUE + ...

  2. 又一道简单题&&Ladygod(两道思维水题)

    Ladygod Time Limit: 3000/1000MS (Java/Others)     Memory Limit: 65535/65535KB (Java/Others) Submit S ...

  3. 一道js题

    <script> var a = 5; function test(){ this.a = 10; a = 15 this.func = function(){ var a = 20 ; ...

  4. 每天一道Java题[3]

    问题 为什么在重写equals()方法的同时,必须重写hashCode()方法? 解答 在<每天一道Java题[2]>中,已经对hashCode()能否判断两个对象是否相等做出了解释.eq ...

  5. 每天一道Java题[11]

    题目 synchronized怎么实现线程同步?请修改<每天一道Java题[10]>中的MyRunnableThread类以解决三个线程都获取到10的问题. 解答 方法一: 采用synch ...

  6. 第三届山西省赛1004 一道大水题(scanf)

    一道大水题 时间限制: C/C++ 2000ms; Java 4000ms 内存限制: 65535KB 通过次数: 44 总提交次数: 1020 问题描述 Dr. Pan作为上兰帝国ACM的总负责人, ...

  7. ny525 一道水题

    一道水题时间限制:1000 ms  |  内存限制:65535 KB 难度:2描述 今天LZQ在玩一种小游戏,但是这游戏数有一点点的大,他一个人玩的累,想多拉一些人进来帮帮他,你能写一个程序帮帮他吗? ...

  8. 【数据结构】 最小生成树(四)——利用kruskal算法搞定例题×3+变形+一道大水题

    在这一专辑(最小生成树)中的上一期讲到了prim算法,但是prim算法比较难懂,为了避免看不懂,就先用kruskal算法写题吧,下面将会将三道例题,加一道变形,以及一道大水题,水到不用高级数据结构,建 ...

  9. 简单tarjan》一道裸题(BZOJ1051)(easy)

    这是一道水题,实际考察的是你会不会写强连通分量...(在BZOJ上又水了一道题) Description 每一头牛的愿望就是变成一头最受欢迎的牛.现在有N头牛,给你M对整数(A,B),表示牛A认为牛B ...

随机推荐

  1. DOM - Document Object Model

    Document Object Model

  2. css 纯css自适应页面 示例

    <!doctype html> <html> <head> <meta charset="utf-8"> <title> ...

  3. OptaPlanner 7.32.0.Final版本彩蛋 - SolverManager之批量求解

    上一篇介绍了OptaPlanner 7.32.0.Final版本中的SolverManager接口可以实现异步求解功能.本篇将继续介绍SolverManager的另一大特性 - 批量求解. 适用场景 ...

  4. python新手如何编写一个猜数字小游戏

    此文章只针对新手,希望大家勿喷,感谢!话不多说先上代码: import random if __name__ == '__main__': yourname = input("你好! 你的名 ...

  5. js this是什么?[多次书写]

    前言 以前的时候,我写了一个关于js this的博客,写的非常复杂,分析了各种情况. 现在我想简化. 如果你有后台基础,专门去理解过this,那么请忘记. 这东西是有口诀的: 在方法中,this 表示 ...

  6. StarUML之四、StarUML的Diagrams(图)与Elements(元素)及相关属性

    Diagrams(图)可以理解为画布  1:创建图 在右侧的Model Explorer管理界面的第一个节点右键,或者选择菜单中Model | Add Diagram | [DiagramType]都 ...

  7. Linux 任务计划,周期性任务执行

    未来的某个时间点执行一次某任务的命令:at,batch 命令执行的结果,会发送到任务发起者的邮箱 周期性运行某任务:crontab 命令执行的结果,会发送到任务发起者的邮箱 本地邮件服务 发邮件的协议 ...

  8. 展讯sprd_battery.c 充电驱动

    sprd_battery.c 是充电驱动,这个是充电功能的核心内容,电量显示策略.温度检测策略.充电保护机制等功能在这里实现,功能实现与硬件细节剥离,调用通用接口实现逻辑控制: 1 sprdbat_p ...

  9. vue中报错Do not use built-in or reserved HTML elements as component id details

    原因是定义了一个叫做details的comonent 跟现有的html网页中的标签重合 export default { name: 'details', data () { return { equ ...

  10. Android中获取目标布局文件中的组件

    方法如下: LayoutInflater flater= LayoutInflater.from(getContext()); //R.layout.title处填写目标布局 final View v ...