Форум программистов, компьютерный форум CyberForum.ru

Проверка работы кода по быстрому преобразованию Фурье - C++

Восстановить пароль Регистрация
Другие темы раздела
C++ Типы данных с плавающей точкой http://www.cyberforum.ru/cpp-beginners/thread1767866.html
Есть какие-нибудь типы данных фиксированной длины, вроде int16_t или __u16, но для работы с плавающей точкой? Вроде бы __fp32 встречается, но вообще, какова тут ситуация?...
C++ Универсальная программа тестирования Напишите универсальную программу тестирования. Тест, последовательность вопросов и варианты ответов должны находиться в текстовом файле. Имя файла теста программа должна получать из командной строки запуска программы. Количество вопросов теста неограничено. Вместе с тем, предлагается ввести следующее ограничение: текст вопроса и альтернативных ответов не должен занимать более одной строки экрана.... http://www.cyberforum.ru/cpp-beginners/thread1767864.html
Как в с++ реализуются обновления программ C++
Есть допустим у пользователя программка, я решил добавить туда кнопочку например, как обновить его программу не заменяя её? Допустим в играх которые весят 10+Гб выходит обнова на 100Мб которая меняет программу. Буду крайне благодарен если расскажите или посоветуете какую либо литературу. Вроде гуглом пользоваться умею, но адекватного ничего не нашёл. Пишу на Visual C++ 2008.
Перегрузка операторов класса Банковских вкладов C++
Класс deposit c полями данных number, investor, amount, percent, time для хранения данных о банковских вкладах (номер вклада, вкладчик, сумма, процентная ставка, длительность). Поле investor имеет тип «строка», поле time хранит данные о длительности вклада в месяцах, поле percent – процент годового дохода. Перегрузить: • операцию == таким образом, чтобы она проверяла на равенство сумму двух...
C++ Как определить наличие конструктора с заданными параметрами? http://www.cyberforum.ru/cpp-beginners/thread1767813.html
Есть шаблонный класс - обёртка вложенного объекта. Его задача: принимать объект формата JSON и передавать его в конструктор нашего объекта, если таковой имеется, иначе ничего не делать - наличие такого конструктора нужно определить на этапе компиляции. Также хочу сразу спросить, можно ли определять существование метода класса с заданным именем и нужными параметрами.
C++ Передать код под шаблонную функцию Ребят помоги передать код под одну шаблонную функцию У меня сделано через перезагрузку функции Сам код: #pragma hdrstop #pragma argsused #ifdef _WIN32 #include <tchar.h> #else typedef char _TCHAR; подробнее

Показать сообщение отдельно
Sarancha
0 / 0 / 0
Регистрация: 21.06.2016
Сообщений: 1
21.06.2016, 11:13     Проверка работы кода по быстрому преобразованию Фурье
Здравствуйте)
Я пытаюсь разобраться с БФП (курсовая работа у меня по нему) и набрела на вот этот код на википедии:

Ниже приведен пример вычисления модуля спектра действительного массива чисел на основе реализации быстрого преобразования Фурье, написанный на C++:
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
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
 
// AVal - массив анализируемых данных, Nvl - длина массива должна быть кратна степени 2. 
// FTvl - массив полученных значений, Nft - длина массива должна быть равна Nvl. 
 
const double TwoPi = 6.283185307179586; 
 
void FFTAnalysis(double *AVal, double *FTvl, int Nvl, int Nft) { 
int i, j, n, m, Mmax, Istp; 
double Tmpr, Tmpi, Wtmp, Theta; 
double Wpr, Wpi, Wr, Wi; 
double *Tmvl; 
 
n = Nvl * 2; Tmvl = new double[n]; 
 
for (i = 0; i < n; i+=2) { 
Tmvl[i] = 0; 
Tmvl[i+1] = AVal[i/2]; 
} 
 
i = 1; j = 1; 
while (i < n) { 
if (j > i) { 
Tmpr = Tmvl[i]; Tmvl[i] = Tmvl[j]; Tmvl[j] = Tmpr; 
Tmpr = Tmvl[i+1]; Tmvl[i+1] = Tmvl[j+1]; Tmvl[j+1] = Tmpr; 
} 
i = i + 2; m = Nvl; 
while ((m >= 2) && (j > m)) { 
j = j - m; m = m » 1; 
} 
j = j + m; 
} 
 
Mmax = 2; 
while (n > Mmax) { 
Theta = -TwoPi / Mmax; Wpi = sin(Theta); 
Wtmp = sin(Theta / 2); Wpr = Wtmp * Wtmp * 2; 
Istp = Mmax * 2; Wr = 1; Wi = 0; m = 1; 
 
while (m < Mmax) { 
i = m; m = m + 2; Tmpr = Wr; Tmpi = Wi; 
Wr = Wr - Tmpr * Wpr - Tmpi * Wpi; 
Wi = Wi + Tmpr * Wpi - Tmpi * Wpr; 
 
while (i < n) { 
j = i + Mmax; 
Tmpr = Wr * Tmvl[j] - Wi * Tmvl[j-1]; 
Tmpi = Wi * Tmvl[j] + Wr * Tmvl[j-1]; 
 
Tmvl[j] = Tmvl[i] - Tmpr; Tmvl[j-1] = Tmvl[i-1] - Tmpi; 
Tmvl[i] = Tmvl[i] + Tmpr; Tmvl[i-1] = Tmvl[i-1] + Tmpi; 
i = i + Istp; 
} 
} 
 
Mmax = Istp; 
} 
 
for (i = 0; i < Nft; i++) { 
j = i * 2; FTvl[i] = 2*sqrt(pow(Tmvl[j],2) + pow(Tmvl[j+1],2))/Nvl; 
} 
 
delete []Tmvl; 
}
Я уже всю голову себе сломала... Что за модуль спектра действительного массива чисел?
и как проверить правильность работы этой ф-ции? Есть ли подходящая для проверки ф-ла?
Или вообще как находится этот "спектр массива"?

Возможно стоит обратиться в другую ветку форума
Но если дадите ответ, буду премного благодарна!
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
 
Текущее время: 22:52. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru