If I use minimum_phase like this: ``` def next_power_of_2(x): return 1 if x == 0 else 2 ** np.ceil(np.log2(x)).astype('int') n = next_power_of_2 (len(signal)) mpt = minimum_phase(signal, n) ``` Then I get back a size of mpt thats being minimum phase transformed thats the same size as the expanded fft, rather than the original signal length. I think something like this would need to be done: ``` m = np.fft.ifft(np.exp(np.fft.fft(window[:len(x)] * ceps))).real ```
If I use minimum_phase like this:
Then I get back a size of mpt thats being minimum phase transformed thats the same size as the expanded fft, rather than the original signal length.
I think something like this would need to be done: