1 #ifndef _NUMcomplex_h_ 2 #define _NUMcomplex_h_ 3 /* NUMcomplex.h 4 * 5 * Copyright (C) 2017-2020 David Weenink 6 * 7 * This code is free software; you can redistribute it and/or modify 8 * it under the terms of the GNU General Public License as published by 9 * the Free Software Foundation; either version 2 of the License, or (at 10 * your option) any later version. 11 * 12 * This code is distributed in the hope that it will be useful, but 13 * WITHOUT ANY WARRANTY; without even the implied warranty of 14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 15 * General Public License for more details. 16 * 17 * You should have received a copy of the GNU General Public License 18 * along with this work. If not, see <http://www.gnu.org/licenses/>. 19 */ 20 21 #include "NUM2.h" 22 23 /* Gamma [alpha,x] = integral{x, infty, t^(alpha-1)exp(-t)dt}, Gamma [alpha]= Gamma [alpha,0] 24 * alpha and x are complex numbers with Re(alpha) > 0 25 */ 26 dcomplex NUMincompleteGammaFunction (const dcomplex alpha, const dcomplex x); 27 28 /* 29 Get response of a truncated gammatone filter defined on the interval [0, truncationTime] at the resonance frequency 30 */ 31 dcomplex gammaToneFilterResponseAtCentreFrequency (double centre_frequency, double bandwidth, double gamma, double initialPhase, double truncationTime); 32 33 #endif 34 35