다변수 함수의 적분을 정의하고, 적분 가능성, 반복적분, 수치적분, determinant와 volume, change of variables, Lebesgue integration의 핵심 개념을 정리합니다.
4.1 Defining the Integral
1) Integration
Integration: 작은 영역에서의 값을 합하여 전체량을 계산하는 과정.
예를 들어 density가
ρ(x)
이면
∫ρ(x)dx
는 전체 mass를 나타낸다.
다변수에서는
f:Rn→R
을 작은 n-dimensional cubes로 나누어 같은 개념을 확장한다.
2) Indicator Function
집합 A⊂Rn의 indicator function은
1A(x)={1,0,x∈A,x∈/A.
이다.
Meaning: 영역 자체를 함수로 표현한다.
영역 A 위에서만 적분하고 싶으면
∫Af=∫Rnf(x)1A(x)∣dnx∣
로 생각할 수 있다.
3) Dyadic Cube
level N의 dyadic cube는 각 변의 길이가
2N1
인 cube이다.
n-dimensional volume은
volnC=2nN1.
4) Dyadic Paving
Dyadic Paving: 같은 level N에 있는 모든 dyadic cubes의 collection.
표기:
DN(Rn).
Meaning: 공간을 같은 크기의 작은 cubes로 분할한다.
5) Upper and Lower Sums
cube C에서
MC(f)=x∈Csupf(x)
mC(f)=x∈Cinff(x)
라 하자.
Nth upper sum:
UN(f)=C∈DN∑MC(f)volnC.
Nth lower sum:
LN(f)=C∈DN∑mC(f)volnC.
resolution이 증가하면
UN(f)
는 nonincreasing하고
LN(f)
는 nondecreasing한다.
6) Upper and Lower Integrals
U(f)=N→∞limUN(f)
L(f)=N→∞limLN(f).
7) Riemann Integrable Function
bounded이고 bounded support를 가지는 함수 f가
U(f)=L(f)
를 만족하면 integrable이다.
그 공통값을 integral로 정의한다.
∫Rnf(x)∣dnx∣=U(f)=L(f)
8) Riemann Sum
f가 integrable이면 각 cube C에서 임의의 sample point
xC∈C
를 선택하여
C∑f(xC)volnC
을 계산할 수 있다.
paving이 충분히 fine해지면
C∑f(xC)volnC→∫Rnf(x)∣dnx∣.
Meaning: 실제 계산에서는 supremum과 infimum 대신 sample values를 사용할 수 있다.
9) Volume
집합 A의 volume은 indicator function의 integral로 정의한다.
volnA=∫Rn1A(x)∣dnx∣
10) Pavable Set
Pavable Set: indicator function이 integrable하여 volume이 잘 정의되는 집합.
A pavable⟺1A integrable.
4.2 Probability and Centers of Gravity
1) Mass
영역 A의 density가
μ:A→R
이면 mass는
M(A)=∫Aμ(x)∣dnx∣
이다.
2) Center of Gravity
center of gravity의 i번째 coordinate는
xˉi=M(A)∫Axiμ(x)∣dnx∣
이다.
Meaning: density로 weighted된 coordinate의 average이다.
density가 constant이면
xˉi=volnA∫Axi∣dnx∣.
3) Probability Density
Probability Density: nonnegative integrable function
μ:Rn→R
가 모든 event A에 대해
P(A)=∫Aμ(x)∣dnx∣
를 만족하면 μ를 probability density라고 한다.
따라서
μ(x)≥0
이고
∫Rnμ(x)∣dnx∣=1
이다.
4) Expected Value
random variable f의 expected value는
E(f)=∫Sf(x)μ(x)∣dnx∣
이다.
5) Variance
Var(f)=E(f2)−E(f)2
6) Standard Deviation
σ(f)=Var(f)
7) Covariance
cov(f,g)=E(fg)−E(f)E(g).
8) Normal Distribution
standard normal probability density는
μ(x)=2π1e−x2/2
이다.
책에서 기억할 대표적인 범위는 대략
- 68%: mean에서 1σ 이내
- 95%: mean에서 2σ 이내
- 99%: mean에서 약 2.5σ 이내
이다.
4.3 What Functions Can Be Integrated?
1) Integrability
모든 함수가 Riemann integrable한 것은 아니다.
Riemann integral에서는 기본적으로
- bounded function
- bounded support
를 먼저 가정한다.
2) Continuous Functions
f:Rn→R이 continuous이고 bounded support를 가지면 integrable이다.
f continuous with bounded support⇒f integrable
3) Discontinuities
함수가 모든 곳에서 continuous일 필요는 없다.
discontinuity가 충분히 작은 집합에만 존재하면 integrable할 수 있다.
특히 bounded support를 가지며 volume 0인 집합을 제외하고 continuous이면 integrable이다.
4) Oscillation
집합 C 위에서 함수의 oscillation은
oscC(f)=Csupf−Cinff.
Meaning: 해당 영역에서 함수가 얼마나 많이 변하는지를 나타낸다.
integrability는 충분히 작은 cubes에서 전체 oscillation의 영향이 작아지는 것과 연결된다.
4.4 Measure Zero
1) Measure Zero
집합
X⊂Rn
이 Measure Zero라는 것은 모든
ϵ>0
에 대해 open boxes
B1,B2,…
가 존재하여
X⊂i=1⋃∞Bi
이고
i=1∑∞voln(Bi)≤ϵ
이 되도록 할 수 있다는 뜻이다.
Meaning: 집합을 임의로 작은 total volume의 boxes로 덮을 수 있다.
2) Almost Everywhere
Almost Everywhere, abbreviated a.e., 는
except on a set of measure zero
라는 뜻이다.
3) Countable Union
measure-zero sets
A1,A2,…
의 countable union도 measure zero이다.
Ai measure zero⇒i=1⋃∞Ai measure zero.
4) Riemann Integrability Criterion
bounded이고 bounded support를 가진 함수 f에 대해
f is Riemann integrable⟺f is continuous almost everywhere
이다.
즉, discontinuity set이 measure zero이면 된다.
Meaning: Riemann integral은 measure-zero 집합 위의 discontinuities에는 영향을 받지 않는다.
4.5 Fubini's Theorem and Iterated Integrals
1) Fubini's Theorem
multiple integral을 one-dimensional integrals의 반복으로 계산할 수 있다.
예를 들어
f:Rn→R
에 대해 적절한 조건 아래
∫Rnf(x)∣dnx∣=∫−∞∞⋯∫−∞∞f(x1,…,xn)dx1⋯dxn
이다.
2) Two Dimensions
∬Af(x,y)dxdy
를
∫(∫f(x,y)dx)dy
또는
∫(∫f(x,y)dy)dx
형태로 계산할 수 있다.
3) Order of Integration
적절한 조건에서는 integration order를 바꿀 수 있다.
∫(∫f(x,y)dx)dy=∫(∫f(x,y)dy)dx.
Meaning: 계산하기 쉬운 순서를 선택할 수 있다.
4) Setting Up Iterated Integrals
계산은 inside-out으로 하지만 integration region을 설정할 때는 outside-in으로 생각하는 것이 유용하다.
예를 들어 영역
A={(x,y)∣a≤y≤b,g1(y)≤x≤g2(y)}
에서는
∫ab∫g1(y)g2(y)f(x,y)dxdy
로 표현한다.
4.6 Numerical Methods of Integration
1) Numerical Integration
antiderivative를 정확히 구하기 어려운 경우 function values를 sampling하여 integral을 approximation한다.
일반적인 형태는
∫abf(x)dx≈i∑wif(pi)
이다.
여기서
- pi: sample points
- wi: weights
이다.
2) Simpson's Rule
interval [a,b]를 2n개의 equal intervals로 나누고
x0,…,x2n
에서 function을 평가한다.
Simpson approximation은
Sn(f)=6nb−a[f(x0)+4f(x1)+2f(x2)+⋯+4f(x2n−1)+f(x2n)]
이다.
weights는
1,4,2,4,…,2,4,1
의 pattern을 가진다.
Simpson's rule은 cubic polynomial을 exactly integrate하도록 만들어져 있다.
3) Product Rule
one-dimensional numerical integration rule
∫abf(x)dx≈i∑wif(pi)
은 higher dimensions에서 product rule로 확장된다.
예를 들어
[a,b]n
에서는
∫[a,b]nf(x)∣dnx∣≈i1,…,in∑wi1⋯winfpi1⋮pin
이다.
4) Curse of Dimensionality
각 coordinate에서 k개의 points를 사용하면 n dimensions에서는
kn
개의 function evaluations가 필요하다.
Meaning: deterministic product rules는 dimension이 증가하면 계산량이 매우 빠르게 증가한다.
5) Monte Carlo Method
high-dimensional integral에서는 random sampling을 사용할 수 있다.
영역 A에서 uniformly random points
x1,…,xN
을 선택하고
ai=f(xi)
를 계산한다.
sample mean:
aˉ=N1i=1∑Nai
sample variance:
s2=N1i=1∑N(ai−aˉ)2
이다.
uniform probability에서는
E(f)=volnA∫Af(x)∣dnx∣.
따라서
∫Af(x)∣dnx∣≈voln(A)aˉ
이다.
Meaning: Monte Carlo는 높은 dimension에서 product rule보다 유리할 수 있다.
4.7 Other Pavings
1) Limitation of Dyadic Pavings
dyadic paving은 이론을 만들기에는 편리하지만 모든 cubes가 같은 구조로 계속 나뉜다.
따라서 함수가 빠르게 변하는 곳에 더 작은 pieces를 배치하고, 천천히 변하는 곳에는 큰 pieces를 사용하는 것이 어렵다.
2) General Pavings
integration은 dyadic cubes에만 의존하지 않는다.
적절한 general paving을 사용하여 domain을 분할할 수 있다.
Meaning: integral 자체는 특정 grid의 선택에 의존하지 않는다.
3) Adaptive Idea
함수가 빠르게 변화하는 영역에서는 small paving pieces를 사용하고, 변화가 적은 영역에서는 larger pieces를 사용할 수 있다.
이 flexibility는 이후 change of variables를 이해하는 데 중요하다.
4.8 Determinants
1) Determinant
n개의 vectors
v1,…,vn∈Rn
을 입력받는 determinant는 다음 세 properties로 특징지어진다.
Multilinearity
각 argument에 대해 linear하다.
Antisymmetry
두 vectors를 교환하면 sign이 바뀐다.
det(…,vi,…,vj,…)=−det(…,vj,…,vi,…).
Normalization
det(e1,…,en)=1
이다.
Matrix
A=[v1,…,vn]
에 대해
detA=det(v1,…,vn).
2) Equal Columns
두 columns가 같으면
detA=0.
3) Product Rule
det(AB)=det(A)det(B)
이다.
4) Invertibility
A invertible⟺detA=0
이다.
5) Transpose
det(AT)=detA
이다.
6) Triangular Matrix
triangular matrix에서는 determinant가 diagonal entries의 product이다.
detA=a11a22⋯ann
7) Orthogonal Matrix
orthogonal matrix Q에서는
QTQ=I.
따라서
(detQ)2=1
이고
detQ=±1
이다.
8) Determinant and Eigenvalues
matrix A의 eigenvalues를 multiplicity를 포함하여
λ1,…,λn
이라고 하면
detA=λ1⋯λn
이다.
9) Trace
trace는 diagonal entries의 sum이다.
trA=i=1∑naii
책에서는 determinant의 derivative와 trace의 관계도 다룬다.
특히
[Ddet(I)]B=trB.
4.9 Volumes and Determinants
1) Determinant as Volume Scaling
Linear transformation
T:Rn→Rn
에 대해
voln(T(A))=∣det[T]∣voln(A)
이다.
Meaning: ∣detT∣는 linear transformation에 의한 volume scaling factor이다.
2) Geometric Interpretation
∣detT∣>1
이면 volume이 확대되고,
0<∣detT∣<1
이면 volume이 축소된다.
detT=0
이면 dimension이 collapse되어 n-dimensional volume은 0이 된다.
orthogonal transformation에서는
∣detQ∣=1.
따라서
voln(Q(A))=voln(A)
이다.
Meaning: rotation과 reflection은 volume을 보존한다.
4) Determinant and Orientation
absolute value
∣detT∣
는 volume scaling을 나타내며 determinant의 sign은 orientation과 연결된다.
이 개념은 Chapter 6의 differential forms에서 다시 사용된다.
1) Change of Variables
복잡한 영역의 integral을 더 단순한 coordinate system에서 계산한다.
mapping
ϕ:U→V
를 사용한다.
2) Local Linear Approximation
작은 영역에서
ϕ(x+h)≈ϕ(x)+Dϕ(x)h.
따라서 작은 volume은 derivative에 의해 approximately
∣detDϕ(x)∣
배 scale된다.
3) Jacobian Determinant
Jacobian Determinant:
Jϕ(x)=detDϕ(x)
Meaning: coordinate transformation이 local volume을 얼마나 확대 또는 축소시키는지를 나타낸다.
적절한 regularity와 invertibility 조건 아래
∫Vf(y)∣dny∣=∫Uf(ϕ(x))∣detDϕ(x)∣∣dnx∣
이다.
5) Why the Absolute Value?
determinant의 sign은 orientation을 나타내지만 ordinary volume은 nonnegative이므로
∣detDϕ∣
를 사용한다.
6) Polar Coordinates
ϕ(r,θ)=(rcosθrsinθ).
Derivative는
Dϕ=[cosθsinθ−rsinθrcosθ].
따라서
∣detDϕ∣=r
이다.
그래서
dxdy=rdrdθ
가 된다.
Meaning: polar-coordinate integral에 나타나는 r은 임의로 붙는 것이 아니라 local area scaling factor이다.
4.11 Lebesgue Integrals
1) Limitation of the Riemann Integral
이 장 앞부분의 Riemann integral은 기본적으로
- bounded functions
- bounded support
를 중심으로 정의된다.
그러나 실제 문제에서는
- unbounded domain
- unbounded function
- limits of functions
을 다루어야 한다.
이를 위해 Lebesgue integration을 도입한다.
2) Lebesgue Equality
두 functions가 measure-zero set을 제외하고 같으면 Lebesgue의 의미에서 같은 것으로 취급한다.
f=ga.e.
Meaning: measure-zero set에서의 값은 Lebesgue integral에 영향을 주지 않는다.
3) L-Integrable Function
책에서는 Lebesgue-integrable function을 Riemann-integrable functions의 absolutely integrable series로 구성한다.
즉,
f=k=1∑∞fka.e.
이고 각 fk가 Riemann integrable이며
k=1∑∞∫Rn∣fk(x)∣∣dnx∣<∞
이면 f를 L-integrable하게 다룬다.
Lebesgue integral은
∫f=k=1∑∞∫fk
로 정의된다.
4) Relation to Riemann Integration
Riemann integrable function은 모두 Lebesgue integrable이며 두 integrals는 같다.
f R-integrable⇒f L-integrable
그리고
∫Rf=∫Lf.
5) Measure-Zero Sets
measure-zero set에서 function 값을 변경해도 Lebesgue integral은 변하지 않는다.
예를 들어 rational numbers의 indicator function은
1Q(x)={1,0,x∈Q,x∈/Q
이지만
1Q=0a.e.
이므로 Lebesgue integral은
0
이다.
6) Infinite Series and Integration
적절한 absolute integrability 조건 아래
∫k=1∑∞fk=k=1∑∞∫fk
이다.
7) Monotone Convergence
nonnegative L-integrable functions가 monotone하게 limit function으로 증가하는 경우 적절한 조건 아래
∫k→∞limfk=k→∞lim∫fk
로 limit과 integral을 교환할 수 있다.
8) Dominated Convergence Theorem
L-integrable functions fk가
fk(x)→f(x)
almost everywhere이고 어떤 L-integrable function F가 존재하여
∣fk(x)∣≤F(x)
almost everywhere이면
f
도 L-integrable이고
∫f=k→∞lim∫fk
이다.
Meaning: 적절한 dominating function이 있으면 limit과 integral의 순서를 바꿀 수 있다.
9) Fubini for Lebesgue Integrals
L-integrable function에 대해서도 iterated integration이 가능하다.
∫Rn+mf(x,y)∣dnx∣∣dmy∣=∫Rm(∫Rnf(x,y)∣dnx∣)∣dmy∣
이다.
10) Change of Variables for Lebesgue Integrals
Lebesgue integral에서도
∫Vf(v)∣dnv∣=∫Uf(ϕ(u))∣detDϕ(u)∣∣dnu∣
형태의 change of variables formula가 성립한다.
Essential Study Checklist
반드시 이해해야 할 내용:
1A
- Dyadic Paving
- Upper Sum / Lower Sum
UN(f),LN(f)
U(f)=L(f)
- Riemann Sum
- Volume as an Integral
volnA=∫1A∣dnx∣
M(A)=∫Aμ(x)∣dnx∣
- Center of Gravity
- Probability Density
P(A)=∫Aμ(x)∣dnx∣
E(f)=∫fμ
- Measure Zero
- Almost Everywhere
- Riemann Integrability Criterion
f integrable⟺f continuous a.e.
for bounded functions with bounded support.
- Fubini's Theorem
- Iterated Integrals
- Simpson's Rule
- Product Rules
- Monte Carlo Method
- Determinant의 세 핵심 properties
- Multilinear
- Antisymmetric
- Normalized
- Determinant Product Rule
det(AB)=detAdetB
A−1 exists⟺detA=0
- Determinant as Volume Scaling
vol(T(A))=∣detT∣vol(A)
Jϕ=detDϕ
∫Vf=∫U(f∘ϕ)∣detDϕ∣
- Polar Coordinate Jacobian
dxdy=rdrdθ
- Riemann Integral과 Lebesgue Integral의 차이
- Almost Everywhere Equality
- Lebesgue Integrability
- Dominated Convergence Theorem
∫limfk=lim∫fk
under the domination condition.