Стас писал(а): Если правильно понял коло-во отсчетов влияет на разрешающую способность, правильно ли понял?

Да, совершенно правильно.
И конечно же количество отсчётов на частоту основной гармоники сигнала никак не влияет.

8745 ...книгу по ЦОС Лайонс не читал, но осуждаю! :D Что значит 8745 сделать 16 точечное БПФ ? Ты берёшь вектор данных размерностью 16. Делаешь БПФ. Что получаем на выходе? Комплексный вектор размерностью 16. Что в этом векторе? Комплексный спектр на частотах: 0, 31.25, 62.5, 93.75, 125.0, 156.25, ...

Теперь с помощью функции взаимной корреляции найдём координаты фрагмента изображения. В середине - исходная картинка 132х132. Справа вверху - фрагмент картинки 56х56 (дополнен нулями до 132х132) Внизу - функция взаимной корреляции. 103.jpg Полезная ссылка: http://forum.algolist.ru/algorithm-image-gr...

Обрежем исходную картинку до размера 132х132 и выполним прямое и обратное БПФ Рис.1 - исходная картинка.
Рис.2 -логарифм модуля двухмерного преобразования Фурье
Рис.3 - восстановленная картинка

Итак исходная картинка: 01.JPG Картинка, переведена массив в 256 градаций серого и выведена в Surfer8: 101.jpg Качество, конечно низкое, т.к. в Surfer8 не удалось подобрать соответствующую шкалу. Но на этом этапе это не критично, так как файл A(x,y) мы получили и теперь можем его подвергнуть 2D прео...

Прежде чем приступить к обработке, необходимо "вытащить" информацию из файла картинки. Проще всего использовать bmp- файлы. Структура файла приведена здесь: http://c-site.h1.ru/infa/bmp_struct.htm Как превратить RGB в градации серого написано здесь: http://www.cyberforum.ru/csharp-beginners/thread90...

Heavy писал(а): ...И еще: если на входе алгоритма присутствует шум, или просто звук, который не содержит искомой компоненты, то модуль растет со временем.

Может на этот эффект Вам намекали здесь?
http://electronix.ru/forum/index.php?sh ... 33098&st=0

Прямое 2D Преобразование Фурье: subroutine FFTF_2D(N,M,U,UFFT) use IMSL integer*4 i,j,n double complex , dimension (256,256) :: U,UFFT double complex , dimension (256) :: XK,YK do j=1,M do i=1,N XK(i)=U(i,j) ! по строкам end do call DFFTCF(n,XK,YK) do i=1,N UFFT(i,j)=YK(i) end do end do ! UX(i,j) do...

Рассмотрим сначала одномерный случай. Пусть x(t), y(t) две действительные функции. Можно представить представить комплексную функцию f(t)=x(t)+i y(t) Обозначив F(W) преобразование Фурье функции f(t),получим: X(W)=[F(W)+F*(N-W)]/2 Y(W)=[F(W)-F*(N-W)]/(2i) А можно ли аналогично поступить в двухмерном ...

trofim2 писал(а): В последней формуле при к-те а_2 Y[n-2]?

Да. Исправил.

Да, на частотах среза близких к 0 и частоте Найквиста - представление коэффициентов - float, желательно...