Итак, по многочисленным просьбам в теме о распознавании речи
выделяю сабж в отдельный топик.
Для начала имеет смысл напомнить об общей процедуре классификации, приведенной http://forum.vingrad.ru/index.php?act=ST&f=13&t=5714.
Здесь только несколько слов о некоторых моментах практической реализации обсуждаемого. В основном будем придерживаться общей схемы процедуры классификации.
Начнем с главного. Теорию вейвлетов можно почитать в разных источниках. Наиболее понятные мне:
есть на сайте http://www.matlab.ru/wavelet/index.asp (ссылки внизу страницы). Очень толковые материалы здесь: http://www.autex.spb.ru/wavelet/, причем много в pdf.
Обращаю ваше внимание на то, что здесь используем именно вейвлет-пакет.
Ниже идет код функции, реализующей вейвлет-преобразование сигнала, записанного в массив sg. Для чего нужна переменная flag, скажу позже, там станет ясно. Физический смысл используемого преобразования состоит в том, что для разложения исходного сигнала по уровням используется его прореживание во временной области (для простоты скажем, разбиение на четные и нечетные отсчеты) плюс некоторые дополнительные преобразования. В частотной области это (опять же упрощенно) равносильно разбиению спектра на две половины. Применямый здесь вейвлет вряд ли вы где-либо найдете  . разве что только встретите моего коллегу, Сергея Агиевича (фамилия такая), с которым мы вместе работали. Он придумал этот вейвлет. Говорит, что построил его на основе слайна! Честно скажу, сам я глубоко эту тему не ковырял. Поэтому извиняйте, если не смогу на вопросы ответить. Но эта штука вполне работает, можете не сомневаться.
| Код |
(с) Агиевич С.Н., 1999-2000.
void wavelet(complex<double> *sg, int im, int flag)
{
// Вейвлет
int p1=log10l((long double)im)/log10l((long double)2); // p1 - степень двойки
// т.е. int im=pow(2,p1);
double a = 2*M_PI/(im/2);
double b1;
double y;
double *U7 = new double[(int)(im/2)];
double *beta = new double[(int)(im/2)];
complex<double> *sg1 = new complex<double>[im+1]; // массив входных данных
complex<double> *sg2 = new complex<double>[im+1]; // массив входных данных
complex<double> *sigS = new complex<double>[im+1]; // массив входных данных
complex<double> *vx = new complex<double>[im+1]; // ядро Фурье-преобразования
complex<double> *vyx = new complex<double>[im+1]; // ядро Фурье-преобразования (обратного)
complex<double> *wx = new complex<double>[im+1]; // ядро преобразования
complex<double> *wyx = new complex<double>[im+1]; // ядро преобразования
complex<double> *sigD = new complex<double>[im+1]; // ядро преобразования
for(int i=0; i < im/2+1; i++)
{
b1 = a*i;
vx[i] = complex<double>(cos(b1),-sin(b1));
vyx[i] = complex<double>(cos(b1),sin(b1));
}
//взятие половины
if(flag == 1) //1024
{
for(int i=0; i < im/2; i++)
{
sg1[i] = sg[i*2];
sg2[i] = sg[i*2+1];
}
}
else
{
for(int i=0; i < im/2; i++)
{
sg1[i] = sg[i];
sg2[i] = sg[i+im/2];
}
flag = 1;
}
fft(sg1, p1-1, vyx, im/2);
fft(sg2, p1-1, vyx, im/2);
for(int i=0; i < im/2+1; i++)
{
b1 = a*i/2; // для вейвлетов
wx[i] = complex<double>(cos(b1),-sin(b1));
wyx[i] = complex<double>(cos(b1),sin(b1));
}
// вспомогательные величины
for(int i=0; i < im/2; i++)
{
y=cos(M_PI*i/(im/2));
U7[i]=1-pow(1-y,4)*(pow(y,2)-28*y+61)/
(61+479*pow(y,2)+179*pow(y,4)+pow(y,6));
beta[i]= U7[i]/(1+pow(U7[i],2));
}
for(int i=0; i < im/2; i++)
{
sigD[i]=(sg2[i]-wx[i]*sg1[i]*U7[i])*(double)(im/2);
sigS[i]=(sg1[i]+beta[i]*wyx[i]*sigD[i]/(double)(im/2))*(double)(im/2);
}
// получено разрежение в частотной области
// переход во временную область
fft(sigD, p1-1, vx, im/2);
fft(sigS, p1-1, vx, im/2);
// возврат значений sg !!!
for(int i=0;i < im/2;i++)
{
sg[i] = sigS[i];
sg[i + (int)(im/2)] = sigD[i];
}
delete[] vyx; vyx = 0;
delete[] sg1; sg1 = 0;
delete[] vx; vx = 0;
delete[] sg2; sg2 = 0;
delete[] sigS; sigS = 0;
delete[] wyx; wyx = 0;
delete[] sigD; sigD = 0;
delete[] wx; wx = 0;
delete[] U7; U7 = 0;
delete[] beta; beta = 0;
} // конец wavelet
(с) Агиевич С.Н., 1999-2000.
|
Я здесь не привожу код алгоритма БПФ fft(...), т.к. он уже был показан, например, http://forum.vingrad.ru/index.php?act=ST&f=13&t=6185.
To be continued...
З.Ы. Просьба ждать терпеливо... Очень занят, понимаете ли. 
И напоминаю, что если вам нравится, как мы общаемся, то можете оценить свое настроение по 10-балльной шкале внизу страницы |
)
