유한차분 옵션 가격결정: 부호, 안정성, 장애조건과 수렴

시간 부호, 이산 연산자, 경계·사건 조건, 미국형 상보성, 감쇠, 가격과 그릭 수렴을 일관되게 적용해 감사 가능한 유한차분 솔버를 구축합니다.

교육 목적으로만 제공되며 투자 조언이 아닙니다. 투자로 손실이 발생할 수 있습니다.

직접 답변

유한차분 옵션 가격결정은 연속 가치평가 방정식을 가격과 시간 노드로 바꿉니다. 유효한 구현에는 촘촘한 격자만 필요한 것이 아닙니다. 달력시간과 잔존만기의 부호, 공간 연산자, 종단·경계조건, 사건 점프, 시간 이산화 방식, 미국형 장애문제 솔버, 잔차와 수렴 통제가 모두 일관되어야 합니다.

격자 세분화가 줄이는 것은 수치오차뿐입니다. 잘못된 지급함수, 확률모형, 배당 처리, 결제 권리, 시장 입력 또는 체결 가정을 고칠 수는 없습니다. 솔버가 틀린 답에 정밀하게 수렴할 수도 있습니다.

코딩 전에 격자부터 도출하기

계약상 권리, 지급함수 Phi(S), 행사 유형, 결제, 금리, 주식 차입비용, 배당, 변동성 모형과 좌표를 확정합니다. 달력시간 t를 쓰면 Black-Scholes PDE는 V_t + L V = 0이고, 잔존만기 tau = T - t를 쓰면 V_tau = L V이며 V(S, 0) = Phi(S)에서 시작합니다.
가격 또는 로그가격 격자, 절단 영역, 시간 격자를 선택하고 현물, 행사가, 배리어, 배당과 관찰시점을 정렬합니다. N이 구간 수인지 노드 수인지 명시하고 S_max는 Delta S와 별도로 검사합니다.
균일 격자 S_i = i Delta S에서 a_i = 0.5[sigma^2 i^2 - (r - q)i], b_i = -(sigma^2 i^2 + r), c_i = 0.5[sigma^2 i^2 + (r - q)i]를 도출해 (L_h V)_i = a_i V_(i-1) + b_i V_i + c_i V_(i+1)로 둡니다.
시간 단계를 정확히 정의합니다. 명시적 방식은 V^(n+1) = (I + Delta tau L_h)V^n, 암시적 방식은 (I - Delta tau L_h)V^(n+1) = V^n, Crank-Nicolson 방식은 (I - 0.5 Delta tau L_h)V^(n+1) = (I + 0.5 Delta tau L_h)V^n입니다. 경계항을 올바른 우변과 시간층에 넣습니다.
전체 영역에서 계수의 부호, CFL형 제한, 양수성, 단조성, 행렬 특성과 진동을 검사합니다. 명시적 Euler와 후향 Euler는 보통 모두 시간 1차입니다. 암시적 방식은 소산성이 더 강하지만 항상 정확도가 낮은 것은 아닙니다. Crank-Nicolson의 선형 안정성은 단조성이나 꺾임 부근의 무진동성을 보장하지 않습니다. 각 불연속점 주변에 Rannacher 감쇠를 적용했는지와 적용 방법을 기록합니다.
미국형 행사에는 PSOR, 벌점법 또는 문서화된 연산자분할처럼 이름이 명시된 이산 선형상보 기법을 사용합니다. V - Phi >= 0, A V - b >= 0, (V - Phi) elementwise_times (A V - b) = 0을 요구하고 완화계수, 반복 한도와 상보성 잔차를 보고합니다.
현물에서 가격과 그릭을 보간한 뒤 영역, 공간, 시간, 사건 정렬과 솔버 허용오차를 세분화합니다. 유럽형 사례는 공식과, 미국형 사례는 깊은 이항트리 또는 독립 솔버와 비교하고 가격, Delta, Gamma, 자유경계, 무차익 검사와 잔차를 기록합니다.

경계조건은 계약에 속합니다. 표준 가정에서 유럽형 풋은 V(0, tau) = K exp(-r tau)를 사용할 수 있고 S가 커지면 0에 가까워집니다. 미국형 풋은 일반적으로 V(0, tau) = K를 사용합니다. 배당이 있는 유럽형 콜의 고가격 점근식은 다릅니다. 이산 현금배당에는 V(S, t-) = V(max(S - D, 0), t+) 같은 점프와 보간이 필요하며, 이후 올바른 순서로 행사 여부를 판단합니다.

시간 단계가 작아도 중앙차분은 음의 이류계수를 만들 수 있습니다. 풍상, 적합 또는 로그가격 이산화가 필요할 수 있습니다. 확산 PDE에 가격 점프가 자동으로 생기는 것도 아닙니다. 모형, 매개변수, 보정, 수치와 체결 가능 가격 오차를 분리해야 합니다.

네 가지 계산 예시

시계를 바꾸면 시간 부호도 바뀝니다. V_t + L V = 0과 tau = T - t에서 연쇄법칙으로 V_tau = L V를 얻습니다. 따라서 tau = 0의 만기 지급액에서 더 큰 tau로 진행할 때는 양의 L 연산자를 사용합니다. 달력시간을 뒤로 진행할 때의 음의 부호를 그대로 복사하면 다른 문제를 풀게 됩니다.
한 노드만 보면 국소적으로 안전해 보일 수 있습니다. sigma = 25%, r = 4%, q = 0, i = 95, Delta tau = 0.0005로 두면 a_i = 280.13125, b_i = -564.1025, c_i = 283.93125입니다. 명시적 가중치는 0.140065625, 0.71794875, 0.141965625이고 합은 0.99998입니다. 이웃 노드값 6, 5, 4는 4.998로 갱신됩니다. 장애값이 5인 미국형 풋에서는 명시적 투영이 5를 유지하지만, 암시적 또는 Crank-Nicolson 시스템은 여전히 결합 LCP 잔차를 확인해야 합니다.
제한은 전체 영역에서 충족되어야 합니다. 같은 입력에서 i = 200이면 명시적 가중치는 0.623, -0.25002, 0.627입니다. 현물 부근 노드가 괜찮아 보여도 음의 중앙 가중치는 단조성을 무너뜨립니다. 단계를 Delta tau = 0.0002로 줄이면 가중치는 0.2492, 0.499992, 0.2508이 되어 모두 음이 아니고 합은 0.999992입니다.
가격과 그릭은 독립적으로 통제해야 합니다. 유럽형 풋의 S = $95, K = $100, T = 0.5, sigma = 25%, r = 4%, q = 0일 때 Black-Scholes는 d_1 = -0.088633257689, d_2 = -0.265409952985, 가격 $8.413258404601, Delta -0.535313311861, Gamma 0.023662224614를 제시합니다. 격자가 세 출력에 각각 수렴한 뒤에만 미국형 프리미엄이나 자유경계를 신뢰할 수 있습니다.

7단계 통제의 실패 위험

달력시간 t와 잔존만기 tau를 일관되지 않은 부호로 진행할 수 있습니다.
일·년, 금리, 변동성, 배당 단위를 혼용할 수 있습니다.
지급함수, 행사 유형, 결제 권리 또는 종단층을 잘못 구현할 수 있습니다.
연속 배당수익률과 이산 현금배당을 중복 반영하거나 순서를 잘못 적용할 수 있습니다.
S_max가 현물에 너무 가까울 수 있으며, 이는 Delta S가 작아 보여도 마찬가지입니다.
저가격·고가격 경계식이 다른 권리나 보유비용 가정에 해당할 수 있습니다.
배당, 배리어, 고정값 또는 관찰시점이 시간 노드 사이에 놓여 잘못 보간될 수 있습니다.
노드 수와 구간 수 또는 비균일 격자 미분식에서 오프바이원 오류가 날 수 있습니다.
드리프트가 중앙 공간계수를 음수로 만들어 단조성을 훼손할 수 있습니다.
명시적 가중치가 현물에서 먼 곳에서 전체 안정성이나 양수성을 위반할 수 있습니다.
Crank-Nicolson이 지급액 꺾임, 배리어 또는 사건 점프 주변에서 진동할 수 있습니다.
Rannacher 또는 다른 감쇠를 빠뜨리거나 잘못된 사건에 적용하거나 문서화하지 않을 수 있습니다.
가치표면이 매끄러워 보여도 선형계 잔차가 클 수 있습니다.
PSOR 완화계수, 정지 허용오차, 초깃값 또는 반복 한도가 불안정할 수 있습니다.
한 번의 무제약 풀이 뒤 노드별 최댓값을 취하는 것을 정확한 암시적 또는 Crank-Nicolson LCP 풀이로 오인할 수 있습니다.
현물에서의 보간이 가치, Delta, Gamma 또는 자유경계를 편향시킬 수 있습니다.
가격이 수렴해 보여도 Gamma와 행사경계에는 잡음이 남을 수 있습니다.
매끄러운 점근 영역에 들어가기 전에 관측된 수렴차수나 Richardson 외삽을 사용할 수 있습니다.
확산, 로컬 변동성 보간, 점프, 매개변수 또는 보정이 구조적으로 틀릴 수 있습니다.
이론가치를 체결 가능한 매수호가, 매도호가, 헤지 또는 차익거래로 제시할 수 있습니다.

흔한 오해

“격자를 더 세분화하면 가격이 맞다는 뜻이다.” 계약이나 모형의 정확성을 입증하지 못합니다.
“암시적 방식은 본질적으로 명시적 방식보다 부정확하다.” 두 Euler 방식 모두 보통 시간 1차이며 안정성과 오차 특성이 다릅니다.
“Crank-Nicolson은 무조건 안정적이므로 단조적이고 진동이 없다.” 이들은 서로 다른 성질입니다.
“미국형 가치는 한 번의 무제약 풀이 후 점별 최댓값을 취하면 된다.” 결합 암시적 시스템에는 LCP 기법 또는 정량화된 근사가 필요합니다.
“하나의 호가와 일치하면 솔버와 그릭이 검증된다.” 보정과 시장 잡음이 수치 또는 구조 오류를 감출 수 있습니다.