计算广义积分$$\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*}
随机推荐
- dom0的cpu hotplug【续】
上一篇说到,手动xm vcpu-pin住,在hotplug就好了. 本质上,还是因为代码有bug,导致vcpu offline的时候,信息没有清理干净,有残留,当vcpu online的时候,如果调度 ...
- Python:urllib和urllib2的区别(转)
原文链接:http://www.cnblogs.com/yuxc/ 作为一个Python菜鸟,之前一直懵懂于urllib和urllib2,以为2是1的升级版.今天看到老外写的一篇<Python: ...
- dreamvc框架(一)ioc容器的集成
我的dreamvc框架最终写得差点儿相同了,借鉴了非常多开源框架,SpringMVC.Struts2等,眼下放在github上面.地址请猛戳我 写得差点儿相同了,是要写一个总结,把自己当时的思路记录下 ...
- react native mac install
Mac上使用react native tips: 1. 安装Homebrew ruby -e "$(curl -fsSL https://raw.githubusercontent.com/ ...
- Java源码中的发现:快速判断一个int值是几位数
判断一个int值是几位数,要是我自己实现,估计又会想到除法和模运算了,偶然在java标准API源码中发现的写法,很强大. public class Test { final static int[] ...
- ip、数字的互转
# ip ==> 数字 >>> ip2num = lambda x:sum([256**j*int(i) for j,i in enumerate(x.split('.')[: ...
- Android中dip,dp,sp,pt和px的区别
dip: device independent pixels(设备独立像素). 不同设备有不同的显示效果,这个和设备硬件有关,一般我们为了支持WVGA.HVGA和QVGA 推荐使用这个,不依赖像素. ...
- jQuery实现页面滚动时顶部动态显示隐藏
http://www.jqcool.net/jquery-scroll.html 另外headroom.js也行:http://www.bootcss.com/p/headroom.js/
- sid超过8个字符处理步骤
服务端配置如下: [oracle@p3 admin]$ cat listener.ora # listener.ora Network Configuration File: /home/oracle ...
- oracle使用LEFT JOIN关联产生的问题在查询结果中使用CASE WHEN 无法判断
oracle使用LEFT JOIN关联产生的问题在查询结果中使用CASE WHEN 无法判断 查询方式一: SELECT CASE WHEN (SELECT CAST(SUM(CASE ) THEN ...