Java 用梳状滤波器检测基频

我需要使用梳状滤波器/变换来检测java中wav文件的基频。 我还需要实现ZCR，但这很容易

现在我有这个：

            int best = 0, best_step = 0;
            for (int step = 3; step < 400; ++step) {
                int sum = 0;
                for (i = 1; i < 10 && i * step < spectrum.length; ++i) {
                    for (int di = 0; di < i; ++di) {
                        sum += spectrum[i * step + di] / i;
                    }
                }


                sum *= 100;
                comb.add(sum);
            }
            int sum = 0;

            for (i = 0; i < comb.size(); ++i) {
                sum = comb.get(i); // 3 * comb[i] - comb[i-1] - comb[i+1];
                 System.out.println(i + " - " + sum);
                 if (sum > best) {
                    best_step = i;
                    best = sum;
                }
            }

int最佳=0，最佳步长=0；
对于（int步数=3；步数<400；++步数）{
整数和=0；
对于（i=1；i<10&&i*步长最佳值）{
最佳步长=i；
最佳=总和；
}
}

我的问题是这个代码检测到了错误的频率；（我搜索了算法/实现（任何语言），但没有找到任何东西

注意，我不能使用自动关联等。它必须是梳状过滤

编辑：关于我的代码的更多解释：

我加载一个wav文件，将帧放到数组帧中，然后对其进行fft，得到复数数组（名为widmo）（处理复数的简单结构）

现在我将复数的绝对值放入数组谱中：

double[] spectrum = new double[widmo.length];

            for (i = 0; i + 1 < widmo.length; ++i) {
                spectrum[i] = widmo[i].abs();
            }


            ArrayList<Integer> comb = new ArrayList<Integer>();
            int best = 0, best_step = 0;
            for (int step = 3; step < 400; ++step) {
                int sum = 0;
                for (i = 1; i < 10 && i * step < spectrum.length; ++i) {
                    for (int di = 0; di < i; ++di) {
                        sum += spectrum[i * step + di] / i;
                    }
                }

                // sum /= step + 100; // ta linijka pozwala usunąć sporo
                // niespodziewanych skoków częstotliwości

                sum *= 100;
                comb.add(sum);
            }
            int sum = 0;

            for (i = 0; i < comb.size(); ++i) {
                sum = comb.get(i); // 3 * comb[i] - comb[i-1] - comb[i+1];
                // ctx.fillRect(i, canvas.height, 1, -sum);
                System.out.println(i + " - " + sum);
                // tmp.add(new freqTime(sum,));
                if (sum > best) {
                    best_step = i;
                    best = sum;
                }
            }
            System.out.println();
            System.out.println(best_step);
            System.out.println(4 * 44100);
            System.out.println((frames.length / numChanels));
            System.out.println(best_step * 44100
                    / (frames.length / numChanels));

double[]光谱=新的double[widmo.length]；
对于（i=0；i+1最佳值）{
最佳步长=i；
最佳=总和；
}
}
System.out.println（）；
系统输出打印LN（最佳步骤）；
系统输出打印项次（4*44100）；
System.out.println（（frames.length/numChanels））；
System.out.println（最佳步骤*44100
/（帧长度/单位长度）；

最后一个

println

应该显示我的基本频率，但它没有

这是我朋友用javascript写的。

我想出来了（我想：p）。最后一次println给出了基本频率：）。也许有人会需要它，甚至会想要它：）

ArrayList okna=new ArrayList（）；
ArrayList lista=新的ArrayList（）；
intpo2=（int）Math.pow（2，
parseInt（（字符串）potega2Input.getText（））；
po2/=2；
双[]三角形=新的双[po2]；
双maxWykres=0；
int-licznik=0；
int-licznik2=0；
int T=帧长度；
布尔wykresFlaga=false；
对于（int k=0；kmaxF）{
maxF=sumT；
maxFI=a；
}
}
//
//maxFI/=2；
//
如果（wykresFlaga==false）{
maxWykres=maxH；
}
对于（int i=0；iArrayList<double[]> okna = new ArrayList<>();
            ArrayList<freqTime> lista = new ArrayList<freqTime>();

            int po2 = (int) Math.pow(2,
                    Integer.parseInt((String) potega2Input.getText()));

            po2 /= 2;
            double[] triangles = new double[po2];
            double maxWykres = 0;
            int licznik = 0;
            int licznik2 = 0;
            int T = frames.length;

            boolean wykresFlaga = false;

            for (int k = 0; k < T; k += po2) {

                licznik = 0;
                licznik2 = 0;
                double[] tmp = new double[po2];

                Complex[] zespolone = new Complex[po2];

                for (int i = k; i < k + po2; i++) {

                    if (i < T) {
                        tmp[licznik] = frames[i];
                        zespolone[licznik] = new Complex(frames[i], 0);
                        licznik2 = licznik;
                    } else {
                        tmp[licznik] = frames[licznik2];
                        zespolone[licznik] = zespolone[licznik2];

                    }
                    licznik++;
                }

                okna.add(tmp);

                FFT fft = new FFT();

                zespolone = fft.fft(zespolone);

                double maxF = 0;
                double maxFI = 0;

                double maxH = findMaxComp(zespolone);
                double[] doWykresu = new double[zespolone.length];
                for (int a = 2; a < 100; a++) {

                    for (int i = 0; i < po2; i++) {
                        doWykresu[i] = zespolone[i].abs();
                        triangles[i] =  Math.abs(i % (2 * a) - a)
                                * (maxH) / a;

                        // triangles[i] = Math.abs(i % (2 * a) - a) * (maxH)
                        // / a;

                    }

                    double sumT = 0;
                    for (int i = 0; i < po2 / 2; i++) {

                        sumT += triangles[i] * doWykresu[i];
                    }

                    if (sumT > maxF) {
                        maxF = sumT;
                        maxFI = a;
                    }

                }
                //
                // maxFI /= 2;
                //
                if (wykresFlaga == false) {
                    maxWykres = maxH;
                }
                for (int i = 0; i < po2; i++) {
                    doWykresu[i] = zespolone[i].abs();
                    triangles[i] =  Math.abs(i % (2 * maxFI) - maxFI)
                            * (maxWykres) / maxFI;

                }
                if (wykresFlaga == false) {
                    System.out.println("Max w widmie: " + maxWykres);
                    new Wykres(doWykresu, 1, triangles);
                    wykresFlaga = true;
                }
                // System.out.println((2 * 44100 / po2) * maxFI);
                System.out.println((float) (44100 / (float) po2)
                        * 2*(float) maxFI + " Znalzione a: " + maxFI);

                lista.add(new freqTime(
                        (int) ((float) (44100 / (float) po2) *2* (float) maxFI),
                        (double) ((double) po2 / 44100)));

                /*
                 * System.out.println((44100 / po2) * maxFI + " " + maxFI +
                 * " " + maxFI / 44100 + " " + 44100 / (maxFI / po2 * 44100)
                 * + " " + 44100 * maxFI / T);
                 */
                // System.out.println(zespolone[(int) maxFI/2].abs());
            }