Вот это было бы замечательно, от вашей картинки толку мало 

Хмм, и что это даст ?

 Точность чего повысится ?

 В данном случае вы предлагаете интерполяцию но зачем? Можно взять частоту дискретизации 44100 Гц и запас будет выше крыши. Pavia, ход вашей мысли не уловил. Если можно разъясните по подробней.

Pavia, можете объяснить каким образом это может нам помочь ? Мы ведь просто увеличиваем ЧД , и взяв дпф мы просто получаем тот же сигнал только растянутый на большую частоту .
Этот способ сможет помочь при использовании АКФ , но при определении частоты с помощью дпф - только увеличит вычислительную нагрузку.

vedun, 
Kpeved, 
С утра времени подумать не хватило. 
Малость ошибся. Интерполировать надо в частотной области. Во временной выражается добавлением нулей к хвосту. Хотя я предпочту вместо нулей сделать так. БПФ как есть, а потом уже поиск максимума и после уже его уточнение через интерполяцию (кубическую).

Это сообщение отредактировал(а) Pavia - 1.3.2012, 15:18

 Собственно я делал тоже самое, только вместо кубической квадратичная интерполяция.

vedun, не знаю , подойдет ли тебе вот в таком вот формате - PCM 16 бит ,по 2 байта на значение .Чд 44100 и 8000.  Каждую струну я записывал отдельно по 10с по 2 удара (имеется ввиду извлекал звук из струны )). Они подписаны с 6(самая толстая ) по 1 (самая тонкая) , и ещё все вместе .Просто скинул считанные массивы в .txt файлы , так что никакой кодировки .
 Инструмент- дешёвая электруха , микрофон в 15ти см от струн .

http://dl.dropbox.com/u/29529555/44k-1.rar

 Думаю что подойдёт, на выходных попробую, потом отпишу что получилось.

Да, и совсем забыл сказать . Всего 6 струн , настроены по тюнеру кроме двух , кот-е я немного расстроил ) Если определишь какие - думаю задача решена ) . Строй стандартный , в инете все частоты есть . 

 Напишите ещё порядок байт, первый байт в файле это младший в слове, второй старший или наоборот? А то по спектрам сложно понять, они на глазок выглядят одинаково.  И ещё, я правильно догадался что числа со знаком?

 Мне это вообще ни о чём не говорит. Я посчитаю частоты а вы решайте что там и куда отстроено. 

Kpeved, 
А wav не судьба было выложить?

Pavia, я думал напрямую с массивом проще работать . Ок , позже wav выложу .

vedun, числа обычные signed short . Да , со знаком - амплитуда ведь и вверх и вниз идет ) 
 Порядок такой - младший , потом старший .

Померял на файлах с ЧД = 8000 Гц, везде второй импульс.
Там на спектрах стоит лес гармоник, мерял частоту 1-й (она везде почти самого низкого уровня, возможно это завал АЧХ микрофона на низких частотах её давит):

• 1 - 329,25 Гц
• 2 - 252,00 Гц
• 3 - 195,59 Гц
• 4 - 147,49 Гц
• 5 - 110,31 Гц
• 6 - 81,69 Гц

Алгоритм простой, как угол дома. Выделяем импульс, я делал это вручную, можно сделать амплитудный детектор. Берём ДПФ от этого импульса и находим на нём первый максимум. Берём его и 2 окружающих отсчёта. По 3-м точкам строим параболу и находим её максимум. Теперь зная ЧД, номер отсчёта ДПФ и дробную поправку находим положение пика частоты.
В данныхрасчётах я брал длину ДПФ 8192. Каков результат судить вам.

Это сообщение отредактировал(а) vedun - 2.3.2012, 21:14

vedun, 
Пока грубо посчитал с точностью до 0.5Гц через  БПФ N=16*1024 отсчёта fd=8000

1 - 328.6 Гц 
2 - 251 Гц 
3 - 194.8 Гц 
4 - 146 Гц
5 -  109.4  Гц
6 - 80 Гц

Что то у нас результаты разнятся на 1 Гц.

Пересчитал более точно через ПФ с шагом 0.1 Гц N=4*16*1024  fd=44100  
1 - 328.5
2 - 251.2
3 - 196.1
4 - 146.5
5 - 109.5
6 - 80.8

Пересчитал более точно через ПФ тс шагом 0.01 Гц N=4*16*1024  fd=44100
1 - 328.53
2 - 251.24
3 - 196.09
4 - 146.55
5 - 109.53
6 - 80.81

Видимо над точностью надо поработать.

Это сообщение отредактировал(а) Pavia - 3.3.2012, 12:11

А какой вы фильтр использовали?

vedun, 
Прямоугольный, по пробовал Ланцоша и Хэмминга разница 0.02 Гц от прямоугольного.
Но видимо из-за упругости струны часта плавает в пределах 1 Гц. Постепенно уменьшаясь. Но у вас всё равно на 1-2 ГЦ больше.

Само ПФ.
И детектирование частоты



Это сообщение отредактировал(а) Pavia - 3.3.2012, 12:37

 Если я правильно понял то ПФ - это полосовой фильтр. Меня интересовало какая у него полоса пропускания и центральная частота. Это FIR или IIR?

 А что от этого частота может менятся ? Мне казалось что параметры струны влияют только на скорость затухания звука а не на его частоту.