Как вычислить спектральную производную с помощью scipy.fftpack (DST, DCT) ⇐ Python

[Anonymous]

Пакет scipy.fftpack предоставляет большое количество процедур, связанных с дискретными преобразованиями Фурье. Мне нужно вычислить 1-ю и 2-ю производную функции, используя только летнее время (DCT). Однако пакет содержит процедуру diff, которая возвращает kth производную с использованием FFT

Кто-нибудь знает, как получить 1-е и Например, 2-е производные с использованием летнего времени? Вот мой черновик:

Код: Выделить всё

import numpy as np
import matplotlib.pyplot as plt
from scipy.fftpack import dst, idst, dct, idct

L   = 10
N   = 100
a   = 0.4
x   = np.linspace(0,L,N)
# function
u   = np.sin(2*np.pi*x/L)*np.exp(-a*x)

# exact 1st derivative
du  = np.exp(-a*x)*(-a*np.sin(2*np.pi*x/L) + np.cos(2*np.pi*x/L)*2*np.pi/L)

# get 1st derivative
dufft = idst(-dct(u))

plt.figure()
plt.plot(x,u)
plt.plot(x,du)
plt.figure()
plt.plot(x,dufft)
plt.show()

Подробнее здесь: https://stackoverflow.com/questions/576 ... ck-dst-dct