|
|
Я на предыдущей странице Вам писал, этого вопроса еще не было.
Во-первых я глюк у себя обнаружил. Не в fft, при вычислении resul, т.е. модуля спектра. Там у меня написано умножение на v (x*=v,y*=v),
нужно умножение на w (x*=w,y*w). Извините. main вчера на ходу сочинил, описался. Вроде больше не должно быть глюков. Но вообще я на них горазд, осторожнее :-).
Теперь по вопросу. DAT - комплексный массив. 64 комплексных значения.
Четные real,нечетные -imidg. Для действительного сигнала imidg=0, т.е. нечетные=0.
По частотам. Частота дискретизации к размеру Фурье никак не относится.
Исходная формула FFT:
S(k)=SUM(s(i)*exp(-j*2pi*(k*Fd/N)*i*Td)=SUM(s(i)*exp(-j*2pi*(k*/N)*i), поскольку Fd=1/Td.
N - размер Фурье, длина массива. В моем варианте N=64. Шаг между частотами(отсчетами) в спектре=Fd/N. k- номер частоты. Fd параметр свободный. Любым можно выбрать.
А синус нужной частоты генерируется просто: sin=(2*pi*F0*i/Fd). F0 -требуемая частота.
Если Fd/F0 делится нацело, то в спектре будет палка на номере k=Fd/F0.
И комплексно-сопряженная (для действ сигнала) ей на номере l=N-k.
Если Вы выбираете Fd=128 Гц и F0=16 Гц, то палка в спектре будет на 8-м отсчете, сопряженная ей на 64-8=56 отсчете.
Если Fd/F0 не делится нацело, то палки в спектре не будет, она размажется по закону 1/x. С максимумом на ближайшем целом к Fd/F0. Аналогично для сопряженной.
Последнее. В спектре (в DAT) после fft возможны маленькие значения на частотах, где их быть вроде и не должно. Это целочисленные эффекты.
E-mail: 
info@telesys.ru
