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

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

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

Ниже приведен пример вычисления модуля спектра действительного массива чисел на основе реализации быстрого преобразования Фурье, написанный на 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; 
}
Я уже всю голову себе сломала... Что за модуль спектра действительного массива чисел?
и как проверить правильность работы этой ф-ции? Есть ли подходящая для проверки ф-ла?
Или вообще как находится этот "спектр массива"?

Возможно стоит обратиться в другую ветку форума
Но если дадите ответ, буду премного благодарна!
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
21.06.2016, 11:13     Проверка работы кода по быстрому преобразованию Фурье
Посмотрите здесь:

Проверка и доработка кода. C++
C++ проверка кода на правильность
C++ Проверка на правильность кода
C++ Проверка кода программы
Какой прогой можно по быстрому посмотреть таблицу, создаваемую SQL файлом? C++
Обнулить кусок памяти (выделенной под указатели) по быстрому C++
проверка кода, функции C++
C++ Подскажите, как сократить время работы кода? Проверка на простое число

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

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

Текущее время: 21:48. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru