计算广义积分$$\int_0^{+\infty}\cos x^p {\rm d}x,\int_0^{+\infty}\sin x^p {\rm d}x, p>1$$
${\bf 解:}$
在角状域$G=\{z\in\mathbb{C}|0<{\rm Arg}z<\frac{\pi}{2p}\}$上引入辅助函数$e^{iz^p}$, 其中$z^p=|z|^pe^{ip{\rm Arg}z}$,$0<{\rm Arg}z<\frac{\pi}{2p}$, 再设$0<\rho<R<+\infty$, 以及$\gamma_\rho=\partial B(0,\rho)\cap G$,$\gamma_R=\partial B(0,R)\cap G$, 逆时针为它们的正向. 由留数定理(或$Cauchy$积分公式), 得到
\begin{equation}\label{the1}
\int_\rho^Re^{ix^p}{\rm d}x+\int\limits_{\gamma_\rho}e^{iz^p}{\rm d}z+\int^\rho_R e^{x^p}e^{i\frac{\pi}{2p}}{\rm d}x-\int\limits_{\gamma_R}e^{iz^p}{\rm d}z=0
\end{equation}
下面证明$(\ref{the1})$中的第$2$,$4$项分别在$R\rightarrow+\infty,\rho\rightarrow 0^+$时趋向$0$.
当$R\rightarrow+\infty$时注意$e^{iz^p}=e^{iR^pe^{ip\theta}}=e^{R^p\cdot i(\cos p\theta+i \sin p\theta)}= e^{R^p\cdot (-\sin p\theta+i\cos p\theta)}$, 以及当$0<x<\frac{\pi}{2}$时成立$\sin x>\frac{2x}{\pi}$, 可得
\begin{align*} | \int\limits_{\gamma_\rho} e^{iz^p} {\rm d}z | &\leq \int \limits_{\gamma_\rho} |{e^{iz^p}|{\rm d}z } \\ &= \int^\frac{\pi}{2p}_0Re^{-R^p\sin p\theta}{\rm d}\theta\\ &\leq R\int^\frac{\pi}{2p}_0e^{-R^P\frac{2p\theta}{\pi}}{\rm d}\theta \\ &= -\frac{\pi R}{2pR^p}e^{-R^P\frac{2p\theta}{\pi}}|^{\frac{2p\theta}{\pi}}_0 \\ &= \frac{\pi R}{2pR^p}(1-e^{-R^p})\\ &\rightarrow 0(R\rightarrow +\infty)\end{align*}
当$\rho\rightarrow 0^+$时,
\begin{align*}|\int\limits_{\gamma_R}e^{iz^p}{\rm d}z|&=|\int^\frac{\pi}{2p}_0 e^{i\rho^pe^{ip\theta}} \rho e^{i\theta}i{\rm d}\theta|\\&\rightarrow 0(\rho\rightarrow 0^+)\end{align*}
于是可将$(\ref{the1})$化为
\begin{align*}\int_0^{+\infty} e^{ix^p}{\rm d}x &=e^{i\frac{\pi}{2p}}\int_0^{+\infty}e^{-x^p}\\ &=e^{i\frac{\pi}{2p}}\int_0^{+\infty}e^{-t}t^{\frac{1}{p}-1}{\rm d}t\\ &=\frac{1}{p}\Gamma(\frac{1}{p})e^{i\frac{\pi}{2p}}\end{align*}
故
\begin{align}\int_0^{+\infty}\cos x^p {\rm d}x &=\frac{1}{p}\Gamma(\frac{1}{p})\cos \frac{\pi}{2p}\\ \int_0^{+\infty}\sin x^p {\rm d}x &=\frac{1}{p}\Gamma(\frac{1}{p})\sin \frac{\pi}{2p}\end{align}
葛神给出了一个数学分析的做法:
\begin{align*} \int_0^\infty \sin \left( x^n\right)dx &= \frac{1}{n}\int_0^\infty x^{\frac{1}{n}-1} \sin(x) \ dx \quad (x^n \mapsto x) \\
&= \frac{1}{n \Gamma \left(
1-\frac{1}{n}\right)}\int_0^\infty \left(\int_0^\infty
u^{-\frac{1}{n}}e^{-xu}du\right) \sin(x) \ dx\\
&= \frac{1}{n
\Gamma \left( 1-\frac{1}{n}\right)} \int_0^\infty u^{-\frac{1}{n}}
\left( \int_0^\infty e^{-xu}\sin(x) \ dx\right)du\\
&= \frac{1}{n \Gamma \left( 1-\frac{1}{n}\right)} \int_0^\infty \frac{u^{-\frac{1}{n}}}{1+u^2}du \\
&=
\frac{1}{n \Gamma \left( 1-\frac{1}{n}\right)}
\int_0^{\frac{\pi}{2}}\tan^{-\frac{1}{n}}(\theta) d\theta \quad (u=\tan
\theta) \\
&= \frac{1}{n \Gamma \left(
1-\frac{1}{n}\right)}\int_0^{\frac{\pi}{2}}\sin^{-\frac{1}{n}}(\theta)
\cos^{\frac{1}{n}}(\theta) d\theta \\
&= \frac{1}{2n \Gamma \left( 1-\frac{1}{n}\right)} \mathrm{B} \left( \frac{1-n}{2},\frac{1+n}{2}\right) \\
&=
\frac{1}{2n \Gamma \left( 1-\frac{1}{n}\right)} \Gamma \left(
\frac{n-1}{2n}\right)\Gamma \left( \frac{n+1} {2n}\right) \\
&= \frac{\sin \left( \frac{\pi}{n}\right)}{2n\cos \left( \frac{\pi}{2n}\right)}\Gamma \left( \frac{1}{n}\right) \\
&= \frac{1}{n}\sin \left(\frac{\pi }{2n} \right)\Gamma \left( \frac{1}{n}\right)\end{align*}
随机推荐
- [置顶] IOS 基础入门教程
IOS 基础入门教程 教程列表: IOS 简介 IOS环境搭建 Objective C 基础知识 创建第一款iPhone应用程序 IOS操作(action)和输出口(Outlet) iOS - 委托( ...
- 提高你30%的设计效率的PPT快捷键
在编辑幻灯片的状态下: [Ctrl]+[A]选择全部对象或幻灯片 [Ctrl]+[B]应用(解除)文本加粗 [Ctrl]+[C]复制 [Ctrl]+[D]快速复制对象 [Ctrl]+[E]段落居中对齐 ...
- 构建ASP.NET MVC4+EF5+EasyUI+Unity2.x注入的后台管理系统(46)-工作流设计-设计分支
原文:构建ASP.NET MVC4+EF5+EasyUI+Unity2.x注入的后台管理系统(46)-工作流设计-设计分支 系列目录 步骤设置完毕之后,就要设置好流转了,比如财务申请大于50000元( ...
- XtraReport交叉表自适应行高及最佳列宽
1.自适应行头的行高,绑定CustomRowHeight事件,代码如下: private Graphics gr = Graphics.FromHwnd(IntPtr.Zero); private v ...
- Android 仿PhotoShop调色板应用(一)概述
版权声明:本文为博主原创文章,未经博主允许不得转载. 在前面的系列我已经将Android中颜色渲染的原理及使用做了一个整体上概述. 现在开始根据一个比较复杂的实现进行具体的分析,这就是PhotoSho ...
- hdu 4622 Reincarnation(后缀数组)
hdu 4622 Reincarnation 题意:还是比较容易理解,给出一个字符串,最长2000,q个询问,每次询问[l,r]区间内有多少个不同的字串. (为了与论文解释统一,这里解题思路里sa数组 ...
- Thread Runnable 多线程
1. 实现多线程的两种方法 a) 让这个类继承java.lang.Thread,然后重写run方法 b) 让这个类实现 java.lang.Runnable接口,实现r ...
- Image路径问题
BitmapImage image = new BitmapImage(new Uri("/Pictures/log.png", UriKind.Relative)); evaIm ...
- Codeforces 570D - Tree Requests【树形转线性,前缀和】
http://codeforces.com/contest/570/problem/D 给一棵有根树(50w个点)(指定根是1号节点),每个点上有一个小写字母,然后有最多50w个询问,每个询问给出x和 ...
- 关于mvc 分页的 这两个结合着用
http://www.cnblogs.com/JackFeng/archive/2010/01/25/JackFeng.html http://www.webdiyer.com/mvcpager/de ...