Форум программистов, компьютерный форум, киберфорум
С++ для начинающих
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 4.83/6: Рейтинг темы: голосов - 6, средняя оценка - 4.83
1 / 1 / 0
Регистрация: 02.06.2013
Сообщений: 10

Ошибка в функции преобразования Фурье

16.06.2013, 01:02. Показов 1300. Ответов 1
Метки нет (Все метки)

Студворк — интернет-сервис помощи студентам
Здравствуйте!
Делаю гитарный тюнер, и застряла на моменте определения частоты уже снятого с микрофона звука. Звук с микрофона снимается в массив типа char, потом я загоняю его в функцию преобразование Фурье. Функцию честно писала с учебника Кинтцеля, немного подправив.
Допустим, при настройке шестой струны, частота должна быть в районе 82.41 герц. Мне же функция выдает массив из двух значений {-1701, 0}.
Подскажите, пожалуйста, в чем моя ошибка.

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
 
char * massiv1  = new char[400];
complex<double> b[1];
 
void ForwardDft(char samples[], int length, complex<double> *result)
{
   static const double
  twoPi = 2*3.1415926535897932384626;
 
   complex<double> samples2[400];
   for (int g = 0; g < 400; g++) {
    samples2[g].real((double)samples[g]);
}
 
 
   for(int f = 0; f<length; f++)
   {
      result[f] = 0.0;
      for(int t = 0; t < length; t++)   { 
      result[f] += samples2[t] * polar(1.0,-twoPi*f*t/length);
     }
   }
}
 
//Вызов функции
ForwardDft(massiv1, 400, b);
Добавлено через 8 часов 17 минут
Пришла к выводу, что многое делаю неправильно, переделала так:

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
double ForwardDft(char samples[], int length)
{
   static const double
twoPi = 2*3.1415926535897932384626;
 
complex<double> result[400];
complex<double> samples2[400];
double aa[400];
double max;
max=aa[0];
for (int g = 0; g < 400; g++) {
samples2[g].real((double)samples[g]);
}
 
 
   for(int f = 0; f<length; f++)
   {
      result[f] = 0.0;
      for(int t = 0; t < length; t++)
         {
         result[f] += samples2[t] * polar(1.0,-twoPi*f*t/length);
          }
   }
   for (int z = 0; z < 400; z++) {
    aa[z]=result[z].real();
      }
 
for( i = 0; i < 400; ++i )
{
 
max += aa[i];
}
max=max/400;
return max;
}
 
fr = ForwardDft(massiv1, 400);
Упорно выдает частоту 110 (частота пятой струны) на все струны. Подозреваю, что я неправильно работаю с полученными после преобразования Фурье значениями.
0
Programming
Эксперт
39485 / 9562 / 3019
Регистрация: 12.04.2006
Сообщений: 41,671
Блог
16.06.2013, 01:02
Ответы с готовыми решениями:

Ошибка в функции преобразования времени
Как можно или нужно написать программу, которая, при вводе с клавиатуры, делает преобразования времени при помощи двух функций. Одна...

Нужна рабочая функция быстрого/прямого и обратного преобразования Фурье
Нужна рабочая функция быстрого/прямого и обратного преобразования Фурье. В гугле много всего, рабочего не нашел.

Как из Wav файла получить комплексные значения для быстрого преобразования Фурье?
Пишу программу для вывода спектра Wav файла Беру данные из Wav файла, прочел шапку, дошел до самих значений. Как понимаю это амплитудные...

1
1 / 1 / 0
Регистрация: 02.06.2013
Сообщений: 10
18.06.2013, 21:03  [ТС]
Задачу решила, если вдруг кому-то тоже будет интересно, напишу, как:
Чтобы получить частоту, необходимо использовать индекс максимального элемента, а не сам максимальный элемент, т.к. его значение - это не частота, а максимальная энергия отрезка.
Дальше работаем по формуле:
частота = найденный индекс/размер массива * частота дискретизации

Плюс, как мне подсказали, индекс необходимо искать не по всему массиву, а только по его первой половине из-за зеркального эффекта.
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
inter-admin
Эксперт
29715 / 6470 / 2152
Регистрация: 06.03.2009
Сообщений: 28,500
Блог
18.06.2013, 21:03
Помогаю со студенческими работами здесь

Функции преобразования
Hi all Вот первый код #include &quot;stdafx.h&quot; #include &quot;conio.h&quot; using namespace std; class String{ private: char* ptr; ...

Функции преобразования
Подскажите, пожалуйста ф-ции преобразования в С++ из int в string и наоборот.

Вычисление функции разложением в ряд Фурье
Объясните для чайника как это решать и что вообще здесь нужно делать . sin x =...

Вычисление функции разложением в ряд Фурье
Надо доделать программу, в общем нужно сделать: 1) Чтобы значение s и y были равными, ну или хотя бы приближенными; 2) Чтобы в...

Вычислить значение функции разложением в ряд фурье
Вычислите s по формуле: s=1+cos(x)/1!+cos(2x)/2!+⋯cos(nx)/n!


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
2
Ответ Создать тему
Новые блоги и статьи
PhpStorm 2025.3: WSL Terminal всегда стартует в ~
and_y87 14.12.2025
PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .
Access
VikBal 11.12.2025
Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.
Новый ноутбук
volvo 07.12.2025
Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro
Музыка, написанная Искусственным Интеллектом
volvo 04.12.2025
Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .
От async/await к виртуальным потокам в Python
IndentationError 23.11.2025
Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .
Поиск "дружественных имён" СОМ портов
Argus19 22.11.2025
Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .
Сколько Государство потратило денег на меня, обеспечивая инсулином.
Programma_Boinc 20.11.2025
Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .
Ломающие изменения в C#.NStar Alpha
Etyuhibosecyu 20.11.2025
Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .
Мысли в слух
kumehtar 18.11.2025
Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .
Создание Single Page Application на фреймах
krapotkin 16.11.2025
Статья исключительно для начинающих. Подходы оригинальностью не блещут. В век Веб все очень привыкли к дизайну Single-Page-Application . Быстренько разберем подход "на фреймах". Мы делаем одну. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2025, CyberForum.ru