기본 콘텐츠로 건너뛰기

[matplotlib] 등고선(Contour)

삼각함수의 미분

내용

삼각함수(trigonometric functions)의 미분

sin, cos, tan의 미분

일반적으로 각을 표시하기 위해 그리스 문자 θ으로 사용합니다. 다음 함수를 고려해 봅니다.

y = sin(θ)

이 함수에서 조사할 것은 그 함수의 변화량 (미분 계수) d(sin(θ))dθ 입니다. 즉, 각도 θ와 sin(θ)의 변화 사이의 관계를 찾는 것입니다. 특히 증가가 무한히 작을 경우가 주된 관심 사항입니다. 이러한 관계를 그림 1에 나타내었습니다.

import numpy as np
import pandas as pd
from sympy import *
import matplotlib.pyplot as plt
plt.figure(figsize=(5, 5))
a=plt.axes(xlim=(-1.5, 2.5), ylim=(-1.5, 2.5))
r=plt.Circle((0, 0), 1, facecolor="none", edgecolor="navy", linewidth=2, label="radius=r")
a.add_patch(r)
plt.arrow(0, 0, 0.8, np.sqrt(1-(-0.8)**2), color="green", label=r"$\mathbf{degree=\theta}}$")
plt.arrow(0, 0, 0.7, np.sqrt(1-(-0.7)**2), color="red", label=r"$\mathbf{degree=\theta+\Delta\theta}}$")
plt.hlines(0, -1.5, 1.5, color="black")
plt.vlines(0, -1.5, 1.5, color="black")
plt.vlines(0.7, 0,np.sqrt(1-(-0.7)**2), linestyle="--", color="red", label=r"$\mathbf{y+\Delta y}}$")
plt.vlines(0.8,0, np.sqrt(1-(-0.8)**2), linestyle="--", color="green", label=r"$\mathbf{y}}$")
plt.legend(loc="best", prop={"weight":"bold"})
plt.text(0.5, -0.2, "X", size=12, weight="bold")
plt.text(-0.2, 0.5,"Y", size=12, weight="bold")
plt.xticks([]) 
plt.yticks([])
plt.show()
그림 1. 원에서 각도의 변화.

그림 1은 반지름 r의 원에서 y는 높이, θ 각을 나타낼 수 있으며 r = 1인 경우 높이와 각 사이 y= sin(θ) 관계가 성립합니다. θ 변화에 따른 sin(θ) 즉, 높이의 변화를 dy(Δ y)라 할 수 있습니다. 새로운 높이는 y+dy가 되고 각은 θ+dθ 될 것입니다.

y+dy=sin(θ+dθ)dy=sin(θ+dθ)sin(θ)

위 식의 우항은 삼각함수의 합법칙을 적용하여 식 1과 같이 전개할 수 있습니다.

(1)sin(m)sin(n)=2cos(m+n2)sin(mn2)m=θ+dθ,n=θdy=2cos(θ+dθ+θ2)sin(θ+dθθ2)=2cos(θ+12dθ)sin(12dθ)

θ 무한히 작다면 다음 코드의 결과와 같이 그 값의 크기는 무시 가능하므로 위 식에서 θ+12dθθ로 간주할 수 있으며 또한 각(θ)과 sin(θ) 값이 거의 일치하기 때문에 sin(12θ)12θ이 성립합니다.

th=[0.001,0.0001, 0.00001, 0.000001]
rad=np.deg2rad(th)
np.around(rad, 8)
array([1.745e-05, 1.750e-06, 1.700e-07, 2.000e-08])
val=[N(sin(i), 8) for i in rad]
pd.DataFrame([rad, val], index=["radian", "sin(x)"])
0 1 2 3
radian 0.000017 0.000002 0.0 0.0
sin(x) 1.74e-5 1.74e-6 1.74e-7 1.74e-8

그러므로 위 식은 다음과 같이 정리 됩니다.

dy=2cos(θ)2dθdydθ=cos(θ)

위 함수 y=sin(θ)와 미분 함수인 cos(θ)는 각각 그림 2와 같습니다.

그림 2. sin(θ)와 미분 곡선.

식 1과 동일한 과정으로 y=cos(θ)와 y=tan(θ)의 미분을 유도하면 각각 식2와 3과 같습니다.

(2)cos(θ)=sin(π2θ)y=cos(θ)dy=d(sin(π2θ))=cos(π2θ)d(θ)=cos(π2θ)d(θ)dydθ=cos(π2θ)=sin(θ) (3)y=tan(θ)=sin(θ)cos(θ)dydθ=d(sin(θ))cos(θ)sin(θ)d(cos(θ))cos2(θ)=cos(θ)cos(θ)+sin(θ)sin(θ)cos2(θ)=cos2(θ)+sin2(θ)cos2(θ)=1cos2(θ)=sec2(θ)

simpy 함수 diff()를 삼각함수의 미분에 적용할 수 있습니다.

theta=symbols("theta")
diff(sin(theta))
cos(θ)
diff(cos(theta))
sin(θ)
diff(tan(theta))
tan2(θ)+1

위 결과 tan2(θ)+1를 정리하면 다음과 같습니다.

tan2(θ)+1=sin2(θ)cos2(θ)+1=1cos2(θ)cos2(θ)+1=1cos2(θ)=sec2(θ)
simplify(diff(tan(theta)))
1cos2(θ)

원의 총 각도는 2π 또는 360도 입니다. 어떤 총 기간을 T로 할 때 일정한 시간 t에 따라 각도의 변화는 다음과 같이 나타낼 수 있습니다.

θ=2πtT

초당 주기(T)의 수를 주파수(frequency, n)라고 하면 n=1T의 관계를 갖습니다. 이 관계를 위 식에 대입하면 다음과 같이 표현됩니다.

θ=2πnt

예)
  함수 y=sin(2πnt)를 미분 합니다.

이 식의 변수는 θ가 아닌 시간 t에 대해 미분하기 위해 다음 관계를 적용합니다.

dydt=dydθdθdty=sin(θ),θ=2πntdydt=cos(θ)2nπ=2nπcos(2πnt)
n, t=symbols("n, t")
y=cos(2*pi*n*t)
y.diff(t)
2πnsin(2πnt)

2차미분과 역함수 미분

삼각함수의 2차미분

sin과 cos의 미분은 각각 상호간의 함수가 전환됩니다. 이 관계를 사용하여 각 함수의 2차 미분을 계산할 수 있습니다.

(4)d(sin(θ))dθ=cos(θ)d2(sin(θ))dθ2=d(cos(θ))dθ=sin(θ) (5)d(cos(θ))dθ=sin(θ)d2(cos(θ))dθ2=d(sin(θ))dθ=cos(θ) (6)d2(tan(θ))dθ=d(1cos2(θ))dθ=2cos(θ)sin(θ)cos4(θ)=2tan(θ)cos2(θ)=2tan(θ)sec2(θ)
theta=symbols("theta")
diff(sin(theta), theta, 2)
sin(θ)
diff(cos(theta), theta, 2)
cos(θ)
simplify(diff(tan(theta), theta, 2))
2tan(θ)cos2(θ)

sin(x), cos(x)함수의 2차 미분계수는 원함수와 부호만 반대인 결과를 갖습니다.

표 1. sin, cos, tan 함수의 미분규칙
함수(y)1차 미분(dydθ)2차 미분(d2ydθ2)
sin(θ) cos(θ) -sin(θ)
cos(θ) -sin(θ) -cos(θ)
tan(θ) sec2(θ) 2tan(θ)sec2(θ)

역함수의 미분

식 7은 sin, cos, tan의 역함수를 나타낸 것입니다.

(7)y=arcsin(x)y=sin1(x)x=sin(y)y=arccos(x)y=cos1(x)x=cos(y)y=arctan(x)y=tan1(x)x=tan(y)

예를 들어 y=sin1(x)는 x = sin(y)를 의미합니다. 이 삼각함수의 역함수의 미분은 다음과 같습니다.

y=sin1(x)x=sin(y)dxdy=cos(y)dydx=1cos(y)

삼각함수의 합공식 sin2(y) + cos2(y) = 1을 적용하여 위 결과를 정리할 수 있습니다.

cos(y)=1sin2(y)=1x2dydx=1cos(y)=11x2
x, y=symbols("x, y")
eq=y-asin(x) #원 식 
eq_x=solve(eq, x) #원 식을 x에 대해 정리(역함수) 
eq_x
[sin(y)]
eq_y=solve(eq, y) #원 식을 y에 대해 정리
eq_y
[asin(x)]
dxdy=eq_x[0].diff(y) # 역함수 미분 
dxdy
cos(y)
(1/dxdy).subs(y, eq_y[0])
11x2
eq_y[0].diff(x)#원 함수 미분
11x2

예)
  cos-1(x)의 미분합니다.

이 함수 역시 위의 sin-1(x)와 동일한 방식으로 계산할 수 있습니다.

y=cos1(x)x=cos(y)dxdy=sin(y)=1cos2(y)dydx=1sin(y)=11cos2(y)=11x2
x, y=symbols("x, y") 
eq=y-acos(x) 
eq_x=solve(eq, x) 
eq_x
[cos(y)]
eq_y=solve(eq, y) #원 식을 y에 대해 정리 
eq_y
[acos(x)]
dxdy=eq_x[0].diff(y) # 역함수 미분  
dxdy
sin(y)
(1/dxdy).subs(y, eq_y[0])
11x2
eq_y[0].diff(x)#원 함수 미분
11x2

예)
  tan-1(x)의 미분합니다.

y=tan1(x)x=tan(y)=sin(y)cos(y)dxdy=1cos2(y)=sin2(y)+cos2(y)cos2(y)=1+tan2(y)=1+x2dydx=11+x2
x= symbols("x")
diff(atan(x))
1x2+1
표 2. 역삼각함수의 미분규칙
함수(y)1차 미분(dydx)
sin-1(x) 11x2
cos-1(x)11x2
tan-1(x)11+x2

예)
  y = cos3(θ)을 미분해 봅니다.

이 식에서 cos3(θ) = (cos(θ))3이 성립하므로 cos(θ) = u로 치환하여 미분할 수 있습니다.

theta, u=symbols("theta u")
u1=cos(theta)
y=u**3
dydu=y.diff(u)
dydu
3u2
dudth=u1.diff(theta)
dudth
sin(θ)
dydth=dydu*dudth
dydth
3u2sin(θ)
dydth.subs(u, u1)
3sin(θ)cos2(θ)
diff(cos(theta)**3, theta) #원 함수를 직접미분
3sin(θ)cos2(θ)

예)
 y = sin(x+a)를 미분합니다.

이 함수 역시 치환 방법을 적용합니다. 이 문제에 대한 미분 과정은 다음과 같습니다.

u=x+ay=sin(u)dudx=1dydu=cos(u)dydududx=cos(u)dudxdydx=cos(x+a)
a, x=symbols("a x")
y=sin(x+a)
y.diff(x)
cos(a+x)

결과적으로 치환한 부분의 미분과 전체 미분을 곱한 것과 같습니다.

예)
  y = log(sin(θ))를 미분합니다. sin(θ) = u로 치환하면 y = log(u)이 됩니다.

theta, u=symbols("theta u")
y=log(sin(theta))#원 함수 
y.diff(theta)
cos(θ)sin(θ)
u1=sin(theta) #원함수의 부분을 u(=u1)로 치환 
y1=log(u)
dydu=y1.diff(u)
dydu
1u
dudth=u1.diff(theta)
dudth
cos(θ)
dydth=dydu*dudth
dydth
cos(θ)u
dydth.subs(u, u1)
cos(θ)sin(θ)

예)
 y = cot(θ)를 미분합니다.

y=cot(θ)=cos(θ)sin(θ)dydx=sin2(θ)cos2(θ)sin2(θ)=1sin2(θ)=csc2(θ)
theta=symbols("theta")
y=cot(theta)
dy=diff(y, theta)
simplify(dy)
1sin2(θ)

예)
 y=tan(3θ)를 미분합니다.

theta=symbols("theta")
y=tan(3*theta)
dy=diff(y, theta)
simplify(dy)
3cos2(3θ)

예)
y=1+3tan2(θ)를 미분합니다.

이 함수에서 3tan2(θ)=u로 치환할 수 있습니다.

theta, u=symbols("theta u")
y=sqrt(1+3*tan(theta)**2)
y.diff(theta)
3(2tan2(θ)+2)tan(θ)23tan2(θ)+1
u1=3*tan(theta)**2
y1=sqrt(1+u)
dydu=y1.diff(u)
dydu
12u+1
dudth=u1.diff(theta)
dudth
3(2tan2(θ)+2)tan(θ)
dydth=dydu*dudth
dydth
3(2tan2(θ)+2)tan(θ)2u+1
dydth.subs(u, u1)
3(2tan2(θ)+2)tan(θ)23tan2(θ)+1

예)
 y = sin(x)cos(x)를 미분합니다.

이 함수의 미분은 곱법칙을 적용합니다.

y=A(x)B(x)dydx=A(x)B(x)+A(x)B(x)dydx=d(sin(x))cos(x)+sin(x)d(cos(x))=cos2(x)sin2(x)
x=symbols("x")
y=sin(x)*cos(x)
dy=diff(y, x)
dy
sin2(x)+cos2(x)

댓글

이 블로그의 인기 게시물

[Linear Algebra] 유사변환(Similarity transformation)

유사변환(Similarity transformation) n×n 차원의 정방 행렬 A, B 그리고 가역 행렬 P 사이에 식 1의 관계가 성립하면 행렬 A와 B는 유사행렬(similarity matrix)이 되며 행렬 A를 가역행렬 P와 B로 분해하는 것을 유사 변환(similarity transformation) 이라고 합니다. (1)A=PBP1P1AP=B 식 2는 식 1의 양변에 B의 고유값을 고려한 것입니다. (식 2)BλI=P1APλP1P=P1(APλP)=P1(AλI)P 식 2의 행렬식은 식 3과 같이 정리됩니다. det(BλI)=det(P1(APλP))=det(P1)det((AλI))det(P)=det(P1)det(P)det((AλI))=det(AλI)det(P1)det(P)=det(P1P)=det(I) 유사행렬의 특성 유사행렬인 두 정방행렬 A와 B는 'A ~ B' 와 같...

[sympy] Sympy객체의 표현을 위한 함수들

Sympy객체의 표현을 위한 함수들 General simplify(x): 식 x(sympy 객체)를 간단히 정리 합니다. import numpy as np from sympy import * x=symbols("x") a=sin(x)**2+cos(x)**2 a sin2(x)+cos2(x) simplify(a) 1 simplify(b) x3+x2x1x2+2x+1 simplify(b) x - 1 c=gamma(x)/gamma(x-2) c Γ(x)Γ(x2) simplify(c) (x2)(x1) 위의 예들 중 객체 c의 감마함수(gamma(x))는 확률분포 등 여러 부분에서 사용되는 표현식으로 다음과 같이 정의 됩니다. 감마함수는 음이 아닌 정수를 제외한 모든 수에서 정의됩니다. 식 1과 같이 자연수에서 감마함수는 factorial(!), 부동소수(양의 실수)인 경우 적분을 적용하여 계산합니다. (식 1)Γ(n)={(n1)!n:자연수0xn1exdxn:부동소수 x=symbols('x') gamma(x).subs(x,4) 6 factorial 계산은 math.factorial() 함수를 사용할 수 있습니다. import math math.factorial(3) 6 a=gamma(x).subs(x,4.5) a.evalf(3) 11.6 simpilfy() 함수의 알고리즘은 식에서 공통사항을 찾아 정리하...

sympy.solvers로 방정식해 구하기

sympy.solvers로 방정식해 구하기 대수 방정식을 해를 계산하기 위해 다음 함수를 사용합니다. sympy.solvers.solve(f, *symbols, **flags) f=0, 즉 동차방정식에 대해 지정한 변수의 해를 계산 f : 식 또는 함수 symbols: 식의 해를 계산하기 위한 변수, 변수가 하나인 경우는 생략가능(자동으로 인식) flags: 계산 또는 결과의 방식을 지정하기 위한 인수들 dict=True: {x:3, y:1}같이 사전형식, 기본값 = False set=True :{(x,3),(y,1)}같이 집합형식, 기본값 = False ratioal=True : 실수를 유리수로 반환, 기본값 = False positive=True: 해들 중에 양수만을 반환, 기본값 = False 예 x2=1의 해를 결정합니다. solve() 함수에 적용하기 위해서는 다음과 같이 식의 한쪽이 0이 되는 형태인 동차식으로 구성되어야 합니다. x21=0 import numpy as np from sympy import * x = symbols('x') solve(x**2-1, x) [-1, 1] 위 식은 계산 과정은 다음과 같습니다. x21=0(x+1)(x1)=0x=1or1x4=1의 해를 결정합니다. solve() 함수의 인수 set=True를 지정하였으므로 결과는 집합(set)형으로 반환됩니다. eq=x**4-1 solve(eq, set=True) ([x], {(-1,), (-I,), (1,), (I,)}) 위의 경우 I는 복소수입니다.즉 위 결과의 과정은 다음과 같습니다. x41=(x2+1)(x+1)(x1)=0x=±1,±1=±i,±1 실수...