inversion.py 文件源码

def _create_Toeplitz(data, npts_stf):
    # Desired dimensions: len(data) + npts_stf - 1 x npts_stf
    nrows = npts_stf + len(data) - 1
    data_col = data
    first_col = np.r_[data_col,
                      np.zeros(nrows - len(data_col))]

    ncols = npts_stf
    first_row = np.r_[data[0], np.zeros(ncols - 1)]
    return toeplitz(first_col, first_row)