楕円の面積と楕円体の体積の求め方をまとめてみた

#目的

「楕円の面積」や「楕円体の体積」の求め方を紹介します。

##理解のためのステップ
【ステップ】
ステップとして下記のステップを踏んで「4.楕円体の体積」を求めたいと思います。

1.円の面積
2.楕円の面積
3.球の体積
4.楕円体の体積

【解法】
A.直接積分する
B.微小面積（体積）を幾何学的に計算して積分する方法
C.ヤコビ行列を使用する方法

![1.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/3b470747-a53c-b293-4c67-ed70ba0e3da6.png)



チェックを入れた方法(AとBとCの方法）で計算して、公式と一致しているかどうかを確認しようと思います。

ここでは、「(1-B)について説明する」と書けば、「1.円の面積」を「B.微小面積（体積）を幾何学的に計算して積分する方法」で計算する方法を説明すると理解してください。




#円の面積

![2.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/bcd6eb7c-dc2e-bbce-1144-1985ae1c69fc.png)

$x^2+y^2=R^2$の内部の面積を求めることになります。

つまり、
$x^2+y^2 \leq R^2$である領域の面積を求める。
変数$x,y$の領域は、
$-R \leq x \leq R$
$-R \leq y \leq R$
です。

>**円の面積の公式**
$\pi R^2$

##(1-A)直接的な解法

![5.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/fe7dbe0b-dae0-6ae9-f9ba-5146e149964a.png)


上記の(1-A)をやっていくことになります。


まず、求めるべく円の面積を$S$とします。

直接的な方法というのは、円の方程式$x^2+y^2=R^2$を関数として扱って、下記の図の微小面積を$x$について足し合わせることを意味しています。
$$S=\int_{-R}^{R} y dx ・・・(1)$$


![3.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/a49283f7-aa9e-a631-a292-5aedf25fce13.png)

ここで注意しなければならないのが、**「円の方程式」とは言っても「円の関数」**とは言わないということです。

関数、「xの値をひとつ決めるとyの値がただいひとつに決まる」ものです（そのように中学生の時に習いました）。

円の方程式は、「xの値をひとつ決めるとyの値が2つ存在することになる」ので関数とは言えません。
確かに、
$$y=\pm\sqrt{R^2-x^2}・・・(2)$$
とすると$y$の値が2つ出てきます。
これでは積分が実行できないので、+記号だけ使います。

$$y=\sqrt{R^2-x^2}・・・(3)$$

そうすると関数として扱うことができるので積分が実行できます。



求める面積は図のように上半分の面積を求めて2倍すれば良いことになります。


$$S=2\int_{-R}^{R} \sqrt{R^2-x^2} dx ・・・(4)$$

さてこれを実行すれば良いのですが、ルートが存在すると積分できないので、置換積分しましょう。（変数を置き換えることになります）

例えば、この場合であると、
$$x=R\sin\theta・・・(5)$$
と置くと、良いです。

>まず$dx$を求めます。
(5)式を使って、
$$\frac{dx}{d\theta}=R\cos\theta$$
となるので、
$$dx=R\cos\theta d\theta・・・(6)$$
と決まります。

次に・・・・
>$$\sqrt{R^2-x^2} =R\sqrt{1-\sin^2\theta}=R\cos\theta・・(6)$$
となります。

さらに・・・
>積分区間は$x$から$\theta$に変わります。
![4.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/523db77e-f54e-91db-29a8-bc096c898a8e.png)


よって、(4)式は、
$$S=2\int_{-\frac{\pi}{2}}^{\frac{\pi}{2}}R^2\cos^2\theta d\theta・・・・(7)$$
となりますので、これを求めれば良いことになります。

半角の公式$\cos^2\theta=\frac{1+\cos2\theta}{2}$を使えば、(7)式は、

$$S=2\int_{-\frac{\pi}{2}}^{\frac{\pi}{2}}R^2\cos^2\theta d\theta=2R^2\int_{-\frac{\pi}{2}}^{\frac{\pi}{2}}\frac{1+\cos2\theta}{2}\theta d\theta=2R^2\bigg[\frac{\theta}{2}+\frac{1}{4}\sin2\theta\bigg]_{-\frac{\pi}{2}}^{\frac{\pi}{2}}=\pi R^2・・・(8)$$
となります。

よって、公式通りになります。

複雑な関数や3次元の体積を求める方法の基本的な積分はこのような手法であると思います。

しかし、何もかもデカルト座標$(x,y,z)$として積分を実行することが賢いやり方だとは思わないので、場合によっては別の座標系に置き換えたりして上手く積分が簡単に実行できる形にすれば良いと思っています。


このような円の面積などのような場合には次のような方法も有効的です。


##(1-B)微小面積（体積）を幾何学的に計算して積分する方法

![6.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/65e8aeec-8319-8f51-ac63-bc255e41dabc.png)

求めたい面積$S$は、変数を$x,y$とすると、
$$S=\int\int dxdy$$
として微小面積について$x$と$y$について足し合わせればよい。
![8.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/7b82ba8d-33c3-e234-faec-1ba63855d39b.png)
この時の微小面積(ABCD)は、$dxdy$と書けます（デカルト座標での表記）が、これが唯一の表記方法ではないはずです。


下記のように、変数を$r$と$\theta$にした極座標表示で書いても良いです。

![7.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/1a8473bb-cbba-c2f2-9486-38862d82f9dc.png)

この時の微小面積（ABCD）は、$rd\theta dr$になります。
>※念のためですが、微小面積は少し四角形とは形が崩れて見えています。長さでいうと$AB \neq CD$です。
しかし、CDは厳密には$(r+dr)d\theta$と書けますが、微小量の2次以上は消去するとすると、結局$rd\theta+drd\theta\simeq rd\theta$となるので、AB=CDとなります。


このように簡単な場合は微小面積は幾何学的に求めることができるので、求めたい面積$S$というのは、
$$S=\int\int rdrd\theta・・・(9)$$となります。

ここで、変数$r$と$\theta$の範囲は、
$0\leq r\leq R$
$0\leq \theta \leq 2\pi$
なので、(9)式の積分すると、
$$S=\int_{0}^{R}\int_{0}^{2\pi} rdrd\theta=2\pi\bigg[\frac{r^2}{2}\bigg]_{0}^{R}=\pi R^2・・・(10)$$

となります。

面積公式と同じ値になりました。



##(1-C)ヤコビ行列を使用する方法

![9.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/8039c808-feb7-ee25-f666-f12c2c1f83c3.png)


(1-B)の方法と似ているのですが、座標変換をして面積を求める方法を紹介します。

円の面積を求めるので、何もデカルト座標で頑張って面積を求める必要はないです。
下記のように極座標表示に変数を変換して、面積を求めても良いです。
むしろ面積が求めやすいように自由に変数を変換して計算した方が計算が楽です。

(1-B)の方法は、微小面積が視覚的にわかりやすいので、微小面積を計算してそこから面積$S$を計算しましたが、今から紹介する「変数変換による面積を求める方法」は、より一般的な変数変換の面積計算にも使えます。


![10.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/4ef2c1c8-6265-08d7-4cd3-b688a34515ec.png)

このように座標をデカルト座標$(x,y)$から極座標$(r,\theta)$に変換します。

$$x(r,\theta)=r\cos\theta・・・(11)$$
$$y(r,\theta)=r\sin\theta・・・(12)$$

$x,y$は$r,\theta$を独立な変数として持ちます。
逆に、
$r,\theta$は、$x,y$を独立な変数として持ちます。
実際、(11)(12)から、
$$r(x,y)=\sqrt{x^2+y^2}・・・(13)$$
$$\theta(x,y)=\tan^{-1}\frac{y}{x}・・・(14)$$
と書けます。

求めたい面積は、
$$S=\int\int dxdy$$
なのですから、$dx$,$dy$が必要なります。

全微分公式より、
$$dx=\frac{\partial x}{\partial r}dr+\frac{\partial x}{\partial \theta}d\theta・・・(15)$$
$$dy=\frac{\partial y}{\partial r}dr+\frac{\partial y}{\partial \theta}d\theta・・・(16)$$

これを求めるのに、
$\frac{\partial x}{\partial r}=\cos\theta$
$\frac{\partial x}{\partial \theta}=\sin\theta$
$\frac{\partial x}{\partial r}=-r\sin\theta$
$\frac{\partial x}{\partial \theta}=r\cos\theta$

だから、(15)(16)式は、
$$dx=\cos\theta dr-r\sin\theta d\theta・・・(17)$$
$$dy=\sin\theta dr+r\cos\theta d\theta・・・(18)$$

これを行列式にすると、

```math
\begin{pmatrix}
\ dx \\ 
\ dy 
\end{pmatrix}=
\begin{pmatrix}
\cos\theta & -r\sin\theta\\ 
\sin\theta & r\cos\theta
\end{pmatrix}
\begin{pmatrix}
\ dr \\ 
\ d\theta 
\end{pmatrix}・・・(19)
```

ここで、

```math
J=
\begin{pmatrix}
\cos\theta & -r\sin\theta\\ 
\sin\theta & r\cos\theta
\end{pmatrix}
```

この行列を「ヤコビ行列」といいます。


なので、(19)式は


```math
\begin{pmatrix}
\ dx \\ 
\ dy 
\end{pmatrix}=
J
\begin{pmatrix}
\ dr \\ 
\ d\theta 
\end{pmatrix}・・・(20)
```

面積は外積の絶対値で求められるので、

$$dxdy=|J|drd\theta ・・・(21)$$


よって、

```math
|J|=\begin{vmatrix}
\cos\theta & -r\sin\theta\\ 
\sin\theta & r\cos\theta
\end{vmatrix}
```
を求めれば(21)式より微小面積の変数変換が完了します。

$$J=r\cos^{2}\theta+r\sin^{2}\theta=r・・・(22)$$

よって、(21)式は、
$$dxdy=rdr d\theta・・・(23)$$

だから求めたい面積$S$は、
$$S=\int\int dxdy=\int\int r dr d\theta・・・(23)$$
となります。

(1-B)で示した(9)式と同じ結果になりました。


このようにして変数変換を行って、微小面積は「ヤコビ行列」で変数変換前と後を結びつける（（21)式）ことで、求めたい面積が簡単に求まります。

ではこの(1-C)の方法を使って、円の面積ではなく楕円の面積を求めてみましょう。

#楕円の面積

![11.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/fb8365fc-5253-94dd-0fb9-a85bd34b8b8b.png)

$\frac{x^2}{a^2}+\frac{y^2}{b^2}=1$の内部の面積を求めることになります。

つまり、
$\frac{x^2}{a^2}+\frac{y^2}{b^2} \leq 1$である領域の面積を求める。
変数$x,y$の領域は、
$-a \leq x \leq a$
$-b \leq y \leq b$
です。

>**円の面積の公式**
$\pi ab$

もちろん(1-A)の方法で計算することもできますが、面倒なので(1-C)の方法を駆使して楕円の面積を求めたいと思います。

##(2-C)ヤコビ行列を使用する方法
![12.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/9e9af2bd-d574-9b2a-05ab-b49155920216.png)


まず、楕円のまま面積を計算するのはめんどうなので、楕円を円に変換してやりましょう。

それは、
$x=au$ $(-1 \leq u \leq1)$
$y=bv$ $(-1 \leq v \leq1)$
とすれば良いです。

このとき、楕円の面積$\frac{x^2}{a^2}+\frac{y^2}{b^2}=1$が、
$$u^2+v^2=1・・・(24)$$
になります。
![13.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/be08698b-2825-3dbe-962d-b5e442cc0a7c.png)
では、求めたい面積、
$$S=\int\int dxdy$$
を計算するにあたって、
ヤコビ行列を用いて微小面積を変換します。

$$dxdy=|J|dudv・・・(25)$$
であり、

```math
J=
\begin{vmatrix}
\ \frac{\partial x}{\partial u} & \frac{\partial x}{\partial v} \\
\ \frac{\partial y}{\partial u} & \frac{\partial y}{\partial v}
\end{vmatrix}
=
\begin{vmatrix}
\ a & 0 \\
\ 0 & b
\end{vmatrix}
=ab
```

よって、(25)式に代入して、
$$dxdy=ab dudv・・・(25)$$

だから求めたい面積は、
$$S=\int\int dxdy=ab\int\int dudv$$

もうここまで来たら、$dudv$を積分した値は「半径1の円の面積と同じ($\pi 1^2$)」だから、
$$S=\pi ab・・・(26)$$

楕円の面積の公式通りの値になりましたね。

結果は、公式通りになったので公式を覚えておけば良いわけですが考え方は非常に重要です。



#球の体積

球の体積
$x^2+y^2+z^2\leq R^2$
$-R\leq x\leq R$
$-R\leq y\leq R$
$-R\leq z\leq R$
の領域の体積を計算することになります。
![15.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/b2f58acf-d648-0859-ae1b-23b7861a4c9f.png)

>球の体積
$\frac{4}{3}\pi R^3$

ちなみにこれを$R$で微分すると、表面積の$4\pi R^2$になります。
それは後ほど明らかになるでしょう。


##(3-B)微小面積（体積）を幾何学的に計算して積分する方法

![14.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/7351e22b-2f0c-3ab6-a186-851a17e049c7.png)

繰り返しになりますが、級の体積であったとしても(A)の方法でももちろん計算は可能です。
しかし、(A)の方法はあくまでデカルト座標での積分の方法であって球を表現するのに必ずしもデカルト座標で考えることが良いとは限らないので積分の方法もそれに適した方法を選択する方が良いでしょう。

まずは、微小面積を幾何学的に計算してみましょう。
![16.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/2aa08e97-783e-c923-a7d5-25335eb72aee.png)



微小体積は絵の通り、
$dV=r^2\sin\theta  dr d\phi d\theta$

だから求める体積は、
$$V=\int\int\int r^2\sin\theta  dr d\phi d\theta・・・(27) $$
積分区間は、
$0\leq r\leq R$
$0\leq \phi\leq 2\pi$
$0\leq \theta\leq \pi$

だから、

$$V=\int_{0}^{R}\int_{0}^{2\pi}\int_{0}^{\pi} r^2\sin\theta  dr d\phi d\theta=\bigg[\frac{1}{3}r^3\bigg]_{0}^{R}\bigg[\phi \bigg]\_{0}\^{2\pi}\bigg[-\cos\theta\bigg]\_{0}\^{\pi}=\frac{4}{3}\pi R^3・・・(28) $$

となり、公式通りになります。

>ちなみに・・・表面積と体積の関係性について・・・
![17.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/9b4f67b6-e994-f731-a7ca-c23b22edf1b0.png)
このように表面積を仮に$4\pi r^2$と知っていたとします。
球の体積は表面積を$r$方向に積算していった値であると考えれば、
$$V=\int_{0}^{R}4\pi r^2 dr=\frac{4}{3}\pi R^3$$
となり、球の表面積と体積は微分積分のの関係にあることが理解できます。
もちろん(27)式のr以外の積分を先にやってしまっても同じ結果が出てきます。


##(3-C)ヤコビ行列を使用する方法
![18.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/f00aff55-779a-0379-9df1-102116ae065f.png)


3変数になるので少々複雑になりますが、基本的な手続きは同じです。
求めたい体積は、
$$V=\int\int\int dxdydz$$
です。

ここで、変数変換をしましょう。
$x=r\sin\phi\cos\theta$
$y=r\sin\phi\sin\theta$
$x=rcos\theta$

それぞれの変数の範囲は、
$0\leq r\leq R$
$0\leq \phi\leq 2\pi$
$0\leq \theta\leq \pi$

では、微小体積はどのようになるかと言いますと、

$$dV=dxdydz=|J|dr d\phi d\theta$$


ここでヤコビ行列式

```math
J=
\begin{vmatrix}
\ \frac{\partial x}{\partial r} & \frac{\partial x}{\partial \phi} & \frac{\partial x}{\partial \theta} \\
\ \frac{\partial y}{\partial r} & \frac{\partial y}{\partial \phi} & \frac{\partial y}{\partial \theta}\\
\ \frac{\partial z}{\partial r} & \frac{\partial z}{\partial \phi} & \frac{\partial z}{\partial \theta}
\end{vmatrix}
=
\begin{vmatrix}
\ \sin\phi \cos\theta & r\cos\phi \cos\theta & -r\sin\phi \sin\theta\\
\ \sin\phi \sin\theta & r\cos\phi \sin\theta & r\sin\phi \cos\theta\\
\ \cos\theta & 0 & -r\cos\theta
\end{vmatrix}
=r^2 \sin\theta 
・・・(29)
```
ゆえに、求めたい体積は、
$$V=\int\int\int dxdydz=\int\int\int |J|dr d\phi d \theta=\frac{4}{3}\pi R^3・・・(30)$$

※計算過程は(28)式と同じなので省略

よって公式通りの結果になりました。


では最後に楕円の体積を求めてみましょう。
これは楕円の面積を求めた方法と手順は同じで、もう既に手法は紹介し尽くしているので結構簡単に結果を導くことができます。


#楕円体の体積


楕円体
$\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}\leq 1$
$-a\leq x\leq a$
$-b\leq y\leq b$
$-c\leq z\leq c$
の領域の体積を計算することになります。

![21.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/a96e9a9f-856a-5ac6-3af3-f20e488499eb.png)


>楕円体の体積の公式
$\frac{4}{3}\pi abc$


##(3-C)ヤコビ行列を使用する方法
![20.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/21723ac7-1c98-8899-b2ca-748781891eec.png)

まず、楕円体のまま体積を計算するのはめんどうなので、楕円体を球体に変換してやりましょう。

それは、
$x=au$ $(-1 \leq u \leq1)$
$y=bv$ $(-1 \leq v \leq1)$
$z=cw$ $(-1 \leq w \leq1)$

とすれば良いです。

このとき、楕円体の体積$\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}=1$が、
$$u^2+v^2+w^2=1・・・(31)$$
になります。

![22.PNG](https://qiita-image-store.s3.amazonaws.com/0/220427/2c82c2ee-f774-4bbf-6c1a-af903c559432.png)


では、求めたい体積、
$$V=\int\int dxdydz$$
を計算するにあたって、
ヤコビ行列を用いて微小体積を変換します。

$$dxdydz=|J|dudvdw・・・(32)$$
であり、

```math
J=
\begin{vmatrix}
\ \frac{\partial x}{\partial u} & \frac{\partial x}{\partial v}  & \frac{\partial x}{\partial w} \\
\ \frac{\partial y}{\partial u} & \frac{\partial y}{\partial v}  & \frac{\partial y}{\partial w} \\
\ \frac{\partial z}{\partial u} & \frac{\partial z}{\partial v}  & \frac{\partial z}{\partial w} 
\end{vmatrix}
=
\begin{vmatrix}
\ a & 0 & 0\\
\ 0 & b & 0\\
\ 0 & 0 & c
\end{vmatrix}
=abc
```

よって、(32)式に代入して、
$$dxdy=abc dudv・・・(33)$$

だから求めたい体積は、
$$V=\int\int dxdydz=abc\int\int dudvdw$$

もうここまで来たら、$dudvdw$を積分した値は「半径1の球の体積と同じ($\frac{4}{3}\pi 1^3$)」だから、
$$V=\frac{4}{3}\pi abc・・・(34)$$

と、公式通りになりました。