用梳状滤波器检测基频

时间:2012-06-10 11:19:00

标签: java audio

我需要使用梳状滤波器/变换来检测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;
                }
            }

我的问题是这段代码检测到错误的频率。 ;(我搜索了一种算法/实现(用任何语言)但没有找到任何东西。

注意,我不能使用autocorelation等。它必须是梳状滤波。

修改:对我的代码进行更多解释:

我加载一个wav文件并将帧放到数组帧中。然后我在它上面做了fft并且有一个Complex(名为widmo)的数组(处理复数的简单结构)。

现在我将复数的abs放入阵列谱中:

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));

最后println应该向我展示我的基本频率,但事实并非如此。

我在朋友的javascript中得到了这个。

1 个答案:

答案 0 :(得分:1)

我弄明白了(我想:P)。最后println给出了基本频率:)。也许有人会需要它甚至改进它:)

ArrayList<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());
            }