public static double[] inverse(int overlap, double[][] specto) {
DoubleFFT_1D fft = new DoubleFFT_1D(specto[0].length / 2);
Window w = new BlackmanWindow();
double[][] spectoNew = specto.clone();
double[][] spectoN = new double[spectoNew.length][spectoNew[0].length * 2];
for (int i = 0; i < spectoN.length; i++) {
for (int j = 0; j < spectoN[0].length; j++) {
if (j < spectoNew[0].length) spectoN[i][j] = spectoNew[i][j];
else spectoN[i][j] = spectoNew[i][spectoNew[i].length - 1 - (j - spectoNew[i].length)];
}
}
for (int i = 0; i < spectoN.length; i++) {
spectoN[i] = w.inverseW(specto[i]);
fft.realInverseFull(spectoN[i], true);
}
int window = spectoN[0].length * 2;
double[] signal = new double[spectoN.length * (window - overlap)];
for (int i = 0; i < spectoN.length; i++) {
int length = spectoN[0].length;
if ((window - overlap) * i + length > signal.length)
length = signal.length - (window - overlap) * i;
System.arraycopy(spectoN[i], 0, signal, (window - overlap) * i, length);
}
return signal;
}
public static double[][] create(int window, double[] signal, int overlap) {
Window w = new BlackmanWindow();
int partitions = signal.length / (window - overlap); // window;
double[][] specto = new double[partitions][window / 2];
System.out.println("Partitions: " + partitions);
for (int i = 0; i < partitions; i++) {
int start = i * (window - overlap);
int length = window * 2;
if (start + window > signal.length) {
break;
}
// System.out.println("[" + i + "] Start at: " + start + " end at: " + (start + length));
double[] fSignals = new double[length];
System.arraycopy(signal, start, fSignals, 0, length / 2);
DoubleFFT_1D fft = new DoubleFFT_1D(length / 2);
// fSignals = DFT.DiscreteFourier(fSignals);
fSignals = w.timesW(fSignals); // Apply the window-function
fft.complexForward(fSignals);
double[] outSignals = new double[fSignals.length / 4];
for (int j = 0; j < outSignals.length / 2; j++) {
outSignals[j] =
Math.sqrt(
fSignals[2 * j] * fSignals[2 * j] + fSignals[2 * j + 1] * fSignals[2 * j + 1]);
}
// fSignals = DFT.forward(fSignals);
System.arraycopy(outSignals, 0, specto[i], 0, outSignals.length);
// System.arraycopy(fSignals, 0, specto[i], 0, length/2);
}
return specto;
}
