편미분(Partial differential)

내용

• 편미분
• 결합함수의 편미분
• 편미분에 의한 극대와 극소

편미분

편미분

다음과 같이 하나 이상의 독립변수가 존재하는 경우의 종속변수는 모든 독립변수들 각각에 영향을 받습니다. 이러한 함수에서 종속변수의 변화는 각 독립변수의 변화에 대한 합으로 나타냅니다. 두 변수에 의한 함수의 경우 식 1과 같이 표현합니다.

$$\begin{array}{}\text{(1)}& y=f(u,v)\end{array}$$

식 1과 같이 두 독립변수로 구성된 함수의 미분은 각 독립변수에 대한 종속변수의 미분(변화율)으로 구성되며 이러한 미분을 편미분(partial differentiation)이라 합니다. 편미분은 특정한 독립변수에 대해 실행할 경우 다른 독립변수를 상수로 간주하는 방식으로 이루어집니다. 그러므로 식1의 미분은 식 2와 같습니다.

$$\begin{array}{r}\text{(2)}& d{y}_v=v\cdot dud{y}_u=u\cdot dv\end{array}$$

식 2의 각 결과에서 아래첨자는 그 시행에서 상수로 간주되는 변수를 의미합니다. 이러한 편미분은 식 3과 같이 그리스 문자 $\partial$를 사용하여 다음과 같이 나타냅니다.

$$\begin{array}{rl}\text{(3)}& \frac{\partial y}{\partial u}=v& \to d{y}_v=\frac{\partial y}{\partial u}\cdot du\frac{\partial y}{\partial v}=u& \to d{y}_u=\frac{\partial y}{\partial v}\cdot dv\end{array}$$

함수 y의 전체적인 미분은 각 편미분의 합으로 나타냅니다.(식 4)

$$\begin{array}{rl}\text{(4)}& dy& =\frac{\partial y}{\partial u}du+\frac{\partial y}{\partial v}dv& =vdu+udv\end{array}$$

import numpy as np
import pandas as pd
from sympy import *
import matplotlib.pyplot as plt

u, v= symbols('u, v')
y=u*v
dydu=diff(y, u) #=∂y/∂u
dydu 

v

dydv=diff(y, v) #=∂y/∂v
dydv

u

예)
  w = 2ax² + 3bxy + 4cy³의 각 독립변수에 대한 편미분계수를 계산해 봅니다.

$$\begin{array}{rl}\frac{\partial w}{\partial x}& =4ax+3by\\ \frac{\partial w}{\partial y}& =3bx+12c{y}^2\end{array}$$

a, b, c, x, y=symbols("a, b, c, x, y")
w=2*a*x**2+3*b*x*y+4*c*y**3
diff(w, x)

 4 a x + 3 b y

diff(w, y)

$3bx+12c{y}^2$

예)
  z = x^y의 x와 y에 대한 편미분계수와 총 미분 dz을 계산해 봅니다.

x, y=symbols("x, y")
z=x**y
dzdx=z.diff(x)
dzdx

$\frac{x^{y}y}{x}$

dzdy=z.diff(y)
dzdy

$x^y\log \left(x\right)$

dz은 결과의 합입니다.

dz=dzdx+dzdy
dz

$x^y\log \left(x\right)+\frac{x^{y}y}{x}$

예)
  높이가 h, 밑면의 반지름이 r인 원뿔의 부피 ${\displaystyle V=\frac{1}{3}\pi r^{2}h}$입니다. 높이가 고정된 상태에서 반지름의 변화 동안 V의 변화율과 반지름이 고정되고 h가 변하는 동안의 V의 변화율?

이 문제의 각각의 부피 변화율은 다음과 같이 나타낼 수 있습니다.

r, h=symbols("r, h")
V=Rational('1/3')*pi*r**2*h
dVdr=diff(V, r)
dVdr

$\frac{2\pi hr}{3}$

dVdh=diff(V, h)
dVdh

$\frac{\pi r^2}{3}$

그러므로 원뿔의 부피변화율은 다음과 같이 나타낼수 있습니다.

$$dv=\frac{2}{3}\pi hr\cdot dr+\frac{1}{3}\pi r^2\cdot dr$$

결합함수의 편미분

삼각 함수 또는 지수 함수 등과 같이 다른 종류의 함수들 결합으로 이루어진 결합함수의 편미분은 치환 방법을 적용할 수 있습니다. 예를 들어 변수 x와 t의 결합된 형태라고 가정하면 그 편미분은 다음과 같이 치환방법을 적용하여 전개할 수 있습니다.

1단계

$$\begin{array}{rl}y& =F(x+at)+f(x-at)\\ w& =x+at\\ v& =x-at\end{array}$$

2단계

$$\begin{array}{rl}\frac{\partial y}{\partial x}& =\frac{\partial F(w)}{\partial w}\frac{\partial w}{\partial x}+\frac{\partial f(v)}{\partial v}\frac{\partial v}{\partial x}\\ & =F^{\prime }(w)\cdot 1+f^{\prime }(v)\cdot 1\end{array}$$

3단계

$$\frac{{\partial }^{2}y}{\partial x^2}=F^{\prime \prime }(w)+f^{\prime \prime }(v)$$

4단계

$$\begin{array}{rl}\frac{\partial y}{\partial t}& =\frac{\partial F(w)}{\partial w}\frac{\partial w}{\partial t}+\frac{\partial f(v)}{\partial v}\frac{\partial v}{\partial t}\\ & =F^{\prime }(w)\cdot a+f^{\prime }(v)\cdot a\end{array}$$

5단계

$$\begin{array}{rl}\frac{{\partial }^{2}y}{\partial t^2}& =F^{\prime \prime }(w)a^2+f^{\prime \prime }(v)a^2\\ & =a^2\frac{{\partial }^{2}y}{\partial x^2}\end{array}$$

편미분에 의한 극대와 극소

여러 개의 독립변수들을 고려하는 점 이외에 편미분 역시 미분과 같은 개념으로 극값을 계산하는 과정 역시 동일합니다. 몇몇의 예제들을 통해 극값을 계산하는 과정을 알아봅니다.

예)
  끈으로 삼각형을 만드는 경우 x, y를 두 부분의 길이라고 합시다. 세 번째 부분의 길이는 30 - (x + y)이고 면적은 다음과 같이 계산됩니다.

$$A=\sqrt{(s-x)(s-y)(s-30+x+z)]}$$

s는 둘레의 절반인 15이므로 $A=\sqrt{15P}$라고 하면 다음 식이 성립합니다. 이러한 관계에서 최대 면적을 계산해 봅니다.

$$\begin{array}{rl}P& =(15-x)(15-y)(x+y-15)\\ & =x{y}^2+x^{2}y-15x^2-15y^2-45xy+450x+450y-3375\end{array}$$

면적 A는 P에 의존되므로 최대 면적은 P가 극대인 경우에 해당하므로 P의 미분 계수가 0인 지점을 계산하는 것입니다. 이 경우는 두 개의 독립변수가 존재이므로 각 변수의 편미분 계수가 0인 지점을 계산해야 합니다. 즉, 다음을 계산합니다.

$$\frac{\partial P}{\partial x}=0,\frac{\partial P}{\partial y}=0$$

x,y=symbols("x, y")
P=(15-x)*(15-y)*(x+y-15)
dPdx=diff(P, x)
dPdx

$(15-x)(15-y)+(y-15)(x+y-15)$

dPdy=diff(P, y)
dPdy

$(15-x)(15-y)+(x-15)(x+y-15)$

sol=solve(Eq(dPdx, dPdy))
sol

[{x: y}, {x: 15 - y}]

위 결과에 의하면 각 편미분 계수가 0이 되는 지점은 x = y 또는 15 - y입니다. 이 값을 원 함수에 대입하여 극값을 계산할 수 있습니다. 그러나 x = 15-y인 경우는 원 식은 0이 되므로 x = y인 경우만을 고려하여 봅니다.

# x = y
P1=P.subs(y, x)
P1

${(15-x)}^2(2x-15)$

dP1=diff(P1, x)
dP1

$2{(15-x)}^2+(2x-30)(2x-15)$

dP1_sol=solve(dP1, x)
dP1_sol

[10, 15]

for i in dP1_sol:
    print(f"{i} 에서 이차미분계수는 {diff(P1, x, 2).subs(x, i)} 입니다.")

10 에서 이차미분계수는 -30 입니다.
    15 에서 이차미분계수는 30 입니다.

위 결과에 의하면 변환된 식의 극점은 x = 10, 15 두 개이고 각 지점에서의 이차 미분계수는 각각 -30, 30 입니다. 그러므로 극대값은 x=10일 때 발생합니다.

예)
  주어진 부피에서 측면과 밑면의 면적이 최소가 되도록 직사각형태의 열차 석탄 트럭의 치수를 계산합니다. 길이가 x, 너비를 y라고 하면 부피 V와 h 그리고 표면적 S는 다음과 같이 계산할 수 있습니다.

$$\begin{array}{rl}& V=x·y·h\\ & h=\frac{V}{xy}\\ & S=xy+\frac{2V}{x}+\frac{2V}{y}\end{array}$$

측면과 바닥의 면적인 최소인 지점으로 표면적 S의 극소점을 결정합니다. V가 주어진 조건이므로 상수입니다.

x,y, V=symbols("x, y, V", real=True, positive=True)
S=x*y+2*V/x+2*V/y
dSdx=diff(S, x)
dSdx

$-\frac{2V}{x^2}+y$

dSdy=diff(S, y)
dSdy

$-\frac{2V}{y^2}+x$

solx=solve(dSdx, x)
solx

[sqrt(2)*sqrt(V)/sqrt(y)]

soly=solve(dSdy, y)
soly

[sqrt(2)*sqrt(V)/sqrt(x)]

solve(Eq(solx[0], soly[0]))

[{x: y}]

위 결과에 의하면 x = y인 지점에서 극점이 형성됩니다.

S1=S.subs(y, x)
S1

$\frac{4V}{x}+x^2$

dS1=diff(S1, x)
dS1

$-\frac{4V}{x^2}+2x$

solS1=solve(dS1, x)
solS1

[2**(1/3)*V**(1/3)]

diff(S1, x, 2).subs(x, solS1[0])

6

극점에 대응하는 극값은 $x=\sqrt[3]{2V}$이며 이 값에서 이차미분은 6으로 양수 입니다. 즉, 이 값은 극소지점임을 의미합니다.