void StarSeemTest::star_external_ffem() { QFETCH(QString, data); AudioBuffer c; QVERIFY2(c.read(data), ("Failed reading test wave; " + data).toUtf8().data()); /*=== WORLD による分析 ===*/ int timeLength = GetSamplesForDIO(c.format().sampleRate(), c.length(), msFramePeriod); int fftLength = GetFFTSizeForStar(c.format().sampleRate()); double *f0 = new double[timeLength]; double *t = new double[timeLength]; double **specgram = new double*[timeLength]; specgram[0] = new double[timeLength * (fftLength / 2 + 1)]; for(int i = 1; i < timeLength; i++) { specgram[i] = specgram[0] + i * (fftLength / 2 + 1); } Dio(c.data()[0], c.length(), c.format().sampleRate(), msFramePeriod, t, f0); Star(c.data()[0], c.length(), c.format().sampleRate(), msFramePeriod, f0, timeLength, specgram); /*=== WORLD による分析ここまで ===*/ /*=== StarSeem による分析 === */ Envelope *e = new Envelope; QVERIFY2(DioFfem().estimate(e, c.data()[0], c.length(), c.format().sampleRate(), msFramePeriod), "Failed dio FFEM"); // Envelope が正しいか確認する for(int i = 0; i < e->size(); i++) { QVERIFY2(e->value(i) == f0[i], "Error; invalid f0 envelope"); } Specgram *testset = new Specgram; QVERIFY2(StarSeem(e).estimate(testset, c.data()[0], c.length(), fftLength, c.format().sampleRate(), msFramePeriod), "Failed STAR SEEM"); /*=== StarSeem による分析ここまで === */ // 比較開始 for(int t = 0; t < timeLength; t++) { for(int f = 0; f <= fftLength / 2; f++) { if(testset->value(t, f) != specgram[t][f]) { QString s; s = "t = " + QString::number(t) + ", f = " + QString::number(f); s += " actual :" + QString::number(testset->value(t, f)) + " , expceted :" + QString::number(specgram[t][f]); // WORLD 自体は同じだから値が全く同じでないとおかしい。 QFAIL(("Error ;" + s).toUtf8().data()); } } } delete[] specgram[0]; delete[] specgram; delete[] t; delete[] f0; }