Обработка WAV файла

[Karateka]

Задача состоит в следующем: есть WAV файл (преподаватель сказал записать свой голос)
1) нужно удостовериться что это WAV файл - тут как я понял это можно посмотреть в заголовках самого файла http://audiocoding.ru/%D1%81%D1%82%D...%BB%D0%B0.html
2) Нужно провести фильтрацию сигнала сохранив звук только в диапазоне от 300Гц до 4кГц.

Как это сделать?!

Почитав некоторые статьи, я понял что нужно разложить в ряд фурье (http://ru.wikipedia.org/wiki/%D0%91%...80%D1%8C%D0%B5) потом обрезать частоты и все это собрать. Есть у кого какие мысли по этому поводу?

На вики есть алгоритм быстрого преобразования в Ряд Фурье, а дальше как отсечь не нужные высоты частот и потом все это собрать в wav файл?

[Karateka]

Еще хотелось бы узнать есть какие-нибудь библиотеки которые все это делают?

[raxp]

1- да, именно в сигнатуре файла
2- провести разложение последовательно каждого сэмпла прямым преобразованием Фурье в спектр, откидыванием частот ниже 300 и выше 4000 Гц и последующим обратным преобразованием Фурье спектра в отсчеты.
3- библиотеки FFT/iFFT - FFTW или моя FFTR.DLL
4- также есть пакеты обработки ACS.

Замечу, что программное создание ФНЧ для обрезания полосы, скажем сверху, совсем необязательно осуществлять через Фурье, если вспомнить про частотную характеристику фильтра "скользящего среднего". Если опустить некоторые выкладки из теории ЦОС (см. мою публикацию "ПЛИС. Часть 3 или... Цифровая фильтрация. Сглаживание и борьба с импульсными помехами. - Радиолюбитель, 2014, №3, с.32"), допустим, имеем сигнал-синус (набор данных, сэмплов WAV PCM) с частотой дискретизации 44100 Гц и разрешением 32 бита на сэмпл. Среднее из скольки сэмплов нужно взять, чтобы остались частоты ниже 4000 Гц? То есть стоит задача отфильтровать высокие частоты.

Обратимся к готовой формуле (ее вывод см. у меня) и подставим исходные данные:

Nопт = round(Fd/2/Fs) = round(44100/2/4000) = round(5.51) ~ 6

Как видим, порядок фильтра или длина окна должны быть равны 6 сэмплам. Таким образом, если нам необходимо убрать все частоты выше 4000 Гц или что эквивалентно введению окна ФНЧ, то мы должны взять один сэмпл, сформированный как среднее из 6-ти сэмплов аудиосигнала.

[Karateka]

Сообщение от raxp

1- да, именно в сигнатуре файла
2- провести разложение последовательно каждого сэмпла прямым преобразованием Фурье в спектр, откидыванием частот ниже 300 и выше 4000 Гц и последующим обратным преобразованием Фурье спектра в отсчеты.
3- библиотеки FFT/iFFT - FFTW или моя FFTR.DLL
4- также есть пакеты обработки ACS.

Есть файл WAV с записью моего голоса. Извлек его шапку. Получил следующие:

0-3 байт RIFF
4-7 байт 926234
8-11 байт WAVE
12-15 байт fmt 
16-19 байт 18
20-21 байт 1
22-23 байт 1
24-27 байт 22050 // эти значения как я понял должны стоять на оборот
28-31 байт 44100 // ведь это частота дискрет. значит должна быть на месте 24-27 байт, но тут они почему-то стоят наоборот. Пробовал другой файл там все ок.
32-33 байт 2
34-35 байт 16
36-39 байт           // здесь должна быть data, но её нет
40-43 байт 291939

Как извлечь данные не пойму. Пишу

  _FileWAVE.Read(subChunk2Size, SizeOf(subChunk2Size));
  Memo1.Lines.Add('40-43 байт ' + IntToStr(subChunk2Size));

  SetLength(data,subChunk2Size);
  _FileWAVE.Read(data, length(data));

в массив ничего не пишется, одни нули. Далее по размыслив, понял что нужно читать по размера сэмпла т.е. получается у меня семпл равен 16, значит читать по 16 бит.

Подскажите как правильно считать данные WAV файла в массив?

И когда мы получаем данные WAV, эти данные мы уже прогоняем через Фурье, правильно?

[Человек_Борща]

Подкину вам хорошую библиотеку обработки WAV.
ещё есть NewAC.

//Сам не в теме.. но мб Вы оттуда что-то да подберете.

[doktor255]

Глубоко сомневаюсь, что учитель примет работу со сторонней библиотекой. Ибо задачу решит не ученик, а чужой код.

[Человек_Борща]

doktor255, я не говорил про компонентокидательство, я про то что готовый код облегчает понимание его задачи. Можно написать свой код быстрее.

[doktor255]

Сообщение от Человек_Борща

я про то что готовый код облегчает понимание его задачи.

Код или DLL библиотека? Если второе, то вообще не облегчит.

[Человек_Борща]

Ну вы хоть глянули что я ему кидаю, прежде чем предполагать что-то с чем-то.

[doktor255]

Ну а оно, по сути мне надо? Я просто по опыту знаю, что большинство библиотек для работы со звуком идут в DLL, вот и высказался. Не DLL и слава богу. Значит задача сведётся к весёлому выдиранию функций из чужой библиотеки.