今天,收到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. 安装 Cacti 监控

    简介:                Cacti是一套基于PHP,MySQL,SNMP及RRDTool开发的网络流量监测图形分析工具.         Cacti是通过 snmpget来获取数据,使用 ...

  2. 用bootstrap来放置天气和图标的位置 自适应

    今天写了个关于天气的页面,他的摆放位置有点难,花了一两个小时用bootstrap来摆放,但是感觉bug很多 所以今天写下自己的心得,放上代码,以后这种就知道怎么写了 <div class=&qu ...

  3. Magicodes.IE基础教程之导出Pdf

    原文作者:hueifeng 说明 本教程主要说明如何使用Magicodes.IE.Pdf完成Pdf收据导出 要点 导出PDF数据 自定义PDF模板 导出单据 如何批量导出单据 导出特性说明 PdfEx ...

  4. search(0)- 企业搜索,写在前面

    计划研究一下搜索search,然后写个学习过程系列博客.开动之前先说说学习搜索的目的:不是想开发个什么搜索引擎,而是想用现成的搜索引擎在传统信息系统中引进搜索的概念和方法.对我来说,传统的管理系统le ...

  5. 单源最短路 Bellman-Ford算法(有向图)

    // 单源最短路问题 // Bellman-Ford算法 // 复杂度O(V*E) //! 可以判断负圈 #include <cstdio> #include <iostream&g ...

  6. 10.python内置模块之os模块

    os模块的作用:os 模块提供了非常丰富的方法用来处理文件和目录(管理和维护目录以及文件). os.path模块的作用:主要用于获取文件的属性(管理路径的(物理地址)). 小生总结了一些平时常用到的属 ...

  7. java.math包下计算浮点数和整数的类

    (1)BigIntege:实现任意精度的整数运算.(2)BigDecimal:实现任意精度的浮点运算. 例如: 使用BigDecimal进行浮点数比较 import java.math.BigDeci ...

  8. watch实现监听Vuex状态监听(利用computed)

    Vuex 通过 store 选项,提供了一种机制将状态从根组件"注入"到每一个子组件中(需调用 Vue.use(Vuex)):通过在根实例中注册 store 选项,该 store ...

  9. 【redis】基于redis实现分布式并发锁

    基于redis实现分布式并发锁(注解实现) 说明 前提, 应用服务是分布式或多服务, 而这些"多"有共同的"redis"; (2017-12-04) 笑哭, 写 ...

  10. 不用循环游标,一句update代码实现滚动计算

    发现一段经典SQL,不用循环游标,一句update代码实现滚动计算结存.为方便理解,结合实例测试之 --1,源数据#t1,jcshl初值为每个sid的当前库存数量,要实现的效果:每个sid的后一结存数 ...