图

误差函数

预备知识 高斯积分

   误差函数的定义为

\begin{equation} \erf(x) = \frac{1}{\sqrt{\pi}} \int_{-x}^x \E^{-t^2} \dd{t} = \frac{2}{\sqrt{\pi}} \int_{0}^x \E^{-t^2} \dd{t} \end{equation}

   由高斯积分中的结论可得

\begin{equation} \int_{-\infty}^\infty \E^{-t^2} \dd{t} = \sqrt{\pi} \end{equation}
所以有
\begin{equation} \erf(\pm\infty) = \pm 1 \end{equation}

   由式 13 , 误差函数的导数为

\begin{equation} \dv{x} \erf(x) = \frac{2}{\sqrt{\pi}} \E^{-x^2} \end{equation}
即 $\E^{-x^2}$ 的不定积分(原函数)为
\begin{equation} \int \E^{-x^2} \dd{x} = \frac{\sqrt{\pi}}{2} \erf(x) + C \end{equation}

级数展开

预备知识 泰勒级数

   由指数函数的级数展开得

\begin{equation} \E^{-x^2} = \sum_{n=0}^\infty \frac{1}{n!} (-x^2)^n = \sum_n \frac{(-1)^n}{n!} x^{2n} \end{equation}
对各项做不定积分代入式 5 可得误差函数的级数展开为
\begin{equation} \erf(x) = \frac{2}{\sqrt{\pi}} \sum_{n=0}^\infty \frac{(-1)^n}{(2n+1)n!} x^{2n+1} = \frac{2}{\sqrt{\pi}}\qtyRound{x - \frac{x^3}{3} + \frac{x^5}{10} - \frac{x^7}{42} + \frac{x^9}{216} \dots} \end{equation}

例1 

   令 $a$ 为实数且 $a > 0$, 计算不定积分

\begin{equation} \int \expRound{-ax^2 + bx} \dd{x} \end{equation}

   我们可以先将指数部分凑平方得

\begin{equation} -ax^2 + bx = -t^2 + \frac{b^2}{4a} \end{equation}
其中
\begin{equation} t = \sqrt{a} x - \frac{b}{2\sqrt{a}} \qquad \end{equation}
使用上式进行换元积分得
\begin{equation} \int \expRound{-ax^2 + bx} \dd{x} = \frac{1}{\sqrt{a}} \E^{\frac{b^2}{4a}} \int \E^{-t^2} \dd{t} = \frac12 \sqrt{\frac{\pi}{a}} \E^{\frac{b^2}{4a}} \erf\qtyRound{\sqrt{a} x - \frac{b}{2\sqrt{a}}} \end{equation}

例2 

   计算无穷积分

\begin{equation} g(k) = \int_{-\infty}^{+\infty} \E^{-a x^2} \E^{-\I kx} \dd{x} \end{equation}
令 $b = -\I k$, 使用例 1 的结论有
\begin{equation} g(k) = \frac12 \sqrt{\frac{\pi}{a}} \E^{-\frac{k^2}{4a}} \eval{\erf\qtyRound{\sqrt{a} x + \frac{\I k}{2\sqrt{a}}}}_{-\infty}^{+\infty} = \sqrt{\frac{\pi}{a}} \E^{-\frac{k^2}{4a}} \end{equation}

   该积分在高斯波包的傅里叶变换中需要使用.

致读者: 小时物理百科一直以来坚持所有内容免费且不做广告,这导致我们处于日渐严重的亏损状态。长此以往很可能会最终导致我们不得不选择商业化,例如大量广告,内容付费,会员制,甚至被收购。因此,我们鼓起勇气在此请求广大读者热心捐款,使网站得以健康发展。如果看到这条信息的每位读者能慷慨捐助 10 元,我们几天内就能脱离亏损状态,并保证网站能在接下来的一整年里向所有读者继续免费提供优质内容。感谢您的支持。
—— 小时(项目创始人)

编辑词条 返回目录 返回主页 捐助项目 © 小时物理百科 保留一切权利