计算广义积分$$\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*}
随机推荐
- 用ChooseALicense帮自己选一个开源license,然后用AddALicense给自己的github自动加上license文件
在我之前的一篇博客里面介绍过tl;drLegal ——开源软件license的搜索引擎,可以很方便的查询各种license,并且给出了很简洁的解释.今天又发现了另外一个帮助你选择你的开源软件licen ...
- 关于无光盘无u盘状态下该如何安装系统
看到好东西,跟大家分享一下,需要装系统的可以作为参考资料 无光盘无u盘状态下该如何安装系统 重点 : 安装虚拟光驱(用来打开镜像文件) 一个你要安装的系统文件的iso镜像文件 http://www. ...
- springframework hibernate Transaction not successfully started
先贴出错误:org.springframework.transaction.TransactionSystemException: Could not commit Hibernate transac ...
- const char*, char const* and char *const 分类: C/C++ OpenCV 2014-11-08 18:10 114人阅读 评论(0) 收藏
const char*, char const*, char*const的区别问题几乎是C++面试中每次都会有的题目. 事实上这个概念谁都有只是三种声明方式非常相似很容易记混. Bjarne在他的 ...
- [Redux] Navigating with React Router <Link>
We will learn how to change the address bar using a component from React Router. In Root.js: We need ...
- NSNotificationCenter消息机制的介绍
转载自http://www.cnblogs.com/pengyingh/articles/2367374.html NSNotificationCenter的作用是专门提供程序中不同类之间的消息通讯而 ...
- Robolectric Test-Drive Your Android Code
Running tests on an Android emulator or device is slow! Building, deploying, and launching the app o ...
- 关于ADO.NET的一些知识整理
ADO.NET是什么 虽然我们都知道ADO.NET是对数据库的操作,但是要真的说出ADO.NET的具体含义还不是很容易. ADO.NET是ActiveX Data Objects的缩写,它是一个COM ...
- Linux命令行编辑快捷键
Linux命令行编辑快捷键: history 显示命令历史列表 ↑(Ctrl+p) 显示上一条命令 ↓(Ctrl+n) 显示下一条命令 !num 执行命令历史列表的第num条命令 !! 执行上一条命令 ...
- oracle 里面定时执行任务,比如存储过程内容等。
DECLARE job_no_ NUMBER; BEGIN DBMS_JOB.SUBMIT(job_no_, 'proc_qszx_dw_sc(' ...