单位圆的双重积分极坐标计算

爱代码的小黄人

于 2024-11-01 10:42:51 发布

阅读量1.3k

点赞数 11

分类专栏： MATLAB 文章标签： matlab

本文链接：https://blog.csdn.net/weixin_44114030/article/details/143425967

版权

MATLAB 专栏收录该内容

59 篇文章

订阅专栏

在这里插入图片描述

要计算如下的双重积分：

$\iint_D e^{-(x^2+y^2)}\,dx\,dy$

其中，区域 $D$ 是单位圆，定义为 $x^2+y^2\leq1$ 。

求解步骤

将积分转换为极坐标

在极坐标中，我们有：
$x=r\cos\theta,\quad y=r\sin\theta$
因此， $x^2+y^2=r^2$ ，并且面积元素 $dx\,dy$ 转换为 $r\,dr\,d\theta$ 。

于是积分可以写成：
$\iint_D e^{-(x^2+y^2)}\,dx\,dy=\int_0^{2\pi}\int_0^1 e^{-r^2}\cdot r\,dr\,d\theta$
对 $r$ 积分

现在我们先对 $r$ 积分：
$\int_0^1 e^{-r^2}\cdot r\,dr$
令 $u=r^2$ ，则 $du=2r\,dr$ ，因此 $r\,dr=\frac{1}{2}\,du$ 。

当 $r = 0$ 时， $u = 0$ ；当 $r = 1$ 时， $u = 1$ 。

于是积分变为：
$\int_0^1 e^{-r^2}\cdot r\,dr=\int_0^1 e^{-u}\cdot\frac{1}{2}\,du=\frac{1}{2}\int_0^1 e^{-u}\,du$
计算 $u$ 的积分

计算 $\int_0^1 e^{-u}\,du$ ：
$\int_0^1 e^{-u}\,du=\left[-e^{-u}\right]_0^1=-(e^{-1}-1)=1-\frac{1}{e}$

因此：
$\int_0^1 e^{-r^2}\cdot r\,dr=\frac{1}{2}\left(1-\frac{1}{e}\right)$
对 $\theta$ 积分

现在对 $\theta$ 从 $0$ 到 $2\pi$ 积分：
$\int_0^{2\pi}\int_0^1 e^{-r^2}\cdot r\,dr\,d\theta=\int_0^{2\pi}\frac{1}{2}\left(1-\frac{1}{e}\right)\,d\theta$

因为 $\frac{1}{2}\left(1-\frac{1}{e}\right)$ 是常数，可以提到积分号外，得到：
$=\frac{1}{2}\left(1-\frac{1}{e}\right)\int_0^{2\pi}d\theta=\frac{1}{2}\left(1-\frac{1}{e}\right)\cdot2\pi$
计算结果

$=\pi\left(1-\frac{1}{e}\right)$

最终答案

$\iint_D e^{-(x^2+y^2)}\,dx\,dy=\pi\left(1-\frac{1}{e}\right)$

MATLAB代码

在 MATLAB 中，我们可以使用数值积分的方式或者直接使用符号积分计算这个积分。以下是使用符号计算的代码示例：

syms r theta  % 定义符号变量

% 定义积分表达式
f = exp(-r^2) * r;

% 先对 r 从 0 到 1 积分
inner_integral = int(f, r, 0, 1);

% 再对 theta 从 0 到 2*pi 积分
result = int(inner_integral, theta, 0, 2*pi);

% 显示结果
disp(result);