计算广义积分$$\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*}
随机推荐
- 不要在精确计算中使用float和double类型
http://blog.csdn.net/androiddevelop/article/details/8478879 一 问题描述 float和double类型不能用于精确计算,其主要目的是为了科 ...
- android_小总结_方法过时的兼容处理
随着android系统的升级,有些过时的方法已经不再使用,但是又要兼容老的版本,所以这个时候可以使用反射来处理下 举个列子pull_to_refresh 中有个方法找不到==initializeScr ...
- [AngularJS] Sane, scalable Angular apps are tricky, but not impossible.
Read fromhttps://medium.com/@bluepnume/sane-scalable-angular-apps-are-tricky-but-not-impossible-less ...
- 虚拟化技术性能总结:Zones, KVM, Xen
[译]虚拟化技术性能总结:Zones, KVM, Xen 时间 2014-04-29 16:52:44 Babyfacer_陈晓炜 原文 http://blog.csdn.net/babyface ...
- html获取gps坐标
<script> function getLocation(){ var options={ enableHighAccuracy:true, maximumAge:1000 } if(n ...
- Canvas模糊化处理图片、毛玻璃处理图片之stackblur.js
Canvas实现毛玻璃效果解决方式1:使用stackblur.js 在Android系统中实现图片的毛玻璃效果比较好用的类库是:Android StackBlur 官方Git地址:https://gi ...
- HTML 5 Audio/Video DOM buffered 属性
1.实例1获取视频第一段缓冲范围部分,以秒计: myVid=document.getElementById("video1"); alert("Start: " ...
- jquery---点击弹出层
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/ ...
- ASP.NET如何发布更新
如果一个应用程序迭代开发后,我们如何来进行发布更新呢? 这里,我用Visual Studio 2008和SQL server来演示. 一.程序修改完毕后,我们要进行以下操作: 1.1 点击解决方案,重 ...
- 使用EMMET中的小坑
使用EMMET写HTML的时候,是一个非常爽的事情.但是今天我使用时,发现一个小坑.以前倒也没有注意,不过需要非常的小心. form[action="/process" metho ...