索引器错误：索引666超出了大小为501的轴1的界限

Question

import math
import numpy as np

def ExplicitMethod(S0, K, r, q, T, Sigma, M, N, Option):

    M = int(M)
    N = int(N)
    dt = T / N
    K = float(K)
    Smax = 2 * K
    dS = Smax / N

    FGrid = np.zeros(shape=(N+1, M+1))

    if Option == 'Call':
        FGrid[-1, :] = np.maximum(np.arange(0, M+1) * dS - K, 0)  
    elif Option == 'Put':
        FGrid[-1, :] = np.maximum(K - np.arange(0, M+1) * dS, 0)

    A = np.zeros(shape=(M+1, M+1))
    A[0,0] = 1
    A[-1,-1] = 1
    for j in range(1, M-1):
        A[j, (j-1,j, j+1)] = \
        [0.5 * dt * (Sigma ** 2 * j ** 2 - (r - q) * j  ), \
         1 - dt * (Sigma ** 2 * j ** 2 + r), \
         0.5 * dt * (Sigma **2 * j ** 2 + (r - q) * j)]


    for i in range(N):
        Fhat = FGrid [i, ]
        FGrid[i, :] = np.dot(A , Fhat)
        Fhat[0] = 0
        Fhat[-1] = Smax - K * math.exp(-r * (N - i) * dt)

    k = math.floor(S0 / dS)

    V = FGrid[0, k] + (FGrid[0, k] - FGrid[0, k]) / dS * [(S0 - k * dS)]

    print (V)

ExplicitMethod (50, 30, 0.1, 0.05, 2, 0.3, 400, 800, Option = 'Call')

我得到了一个错误：

IndexError Traceback (most recent call last) in () 37 print (V) 38 ---> 39 ExplicitMethod (50, 30, 0.1, 0.05, 2, 0.3, 400, 800, Option = 'Call')

in ExplicitMethod(S0, K, r, q, T, Sigma, M, N, Option) 33 k = math.floor(S0 / dS) 34 ---> 35 V = FGrid[0, k] + (FGrid[0, k] - FGrid[0, k]) / dS * [(S0 - k * dS)] 36 37 print (V)

IndexError: index 666 is out of bounds for axis 1 with size 401

作为参考，我正在编写期权定价显式方法：

请帮忙，谢谢！在

Answer 1

您已经将FGrid初始化为（801401）。它告诉我们，出于某种原因，k=666

k = math.floor(S0 / dS)
V = FGrid[0, k] + (FGrid[0, k] - FGrid[0, k]) / dS * [(S0 - k * dS)]

您需要优化k的设置方式。也许数学是错的。至少你需要确保它不会超过400。在

^{pr2}$

这是你的666。在