Форум программистов «Весельчак У»
  *
Добро пожаловать, Гость. Пожалуйста, войдите или зарегистрируйтесь.
Вам не пришло письмо с кодом активации?

  • Рекомендуем проверить настройки временной зоны в вашем профиле (страница "Внешний вид форума", пункт "Часовой пояс:").
  • У нас больше нет рассылок. Если вам приходят письма от наших бывших рассылок mail.ru и subscribe.ru, то знайте, что это не мы рассылаем.
   Начало  
Наши сайты
Помощь Поиск Календарь Почта Войти Регистрация  
 
Страниц: [1]   Вниз
  Печать  
Автор Тема: интегралы через быстрое преобразование Фурье  (Прочитано 17886 раз)
0 Пользователей и 1 Гость смотрят эту тему.
Sommer
Молодой специалист

us
Offline Offline

« : 13-08-2003 11:06 » 

кто - нибудь занимался реализацией этого?...
Записан

когда-нибудь, я верю, ты будешь ехать по этому городу и поймёшь, что хочешь увидеть меня за рулём мчащейся по соседней полосе машины.
но тогда меня уже не будет
в этом городе
forever yours.
RXL
Технический
Администратор

ru
Offline Offline
Пол: Мужской

WWW
« Ответ #1 : 13-08-2003 13:50 » 

Сейчас разбераюсь...
Какой иммено вопрос?
Записан

... мы преодолеваем эту трудность без синтеза распределенных прототипов. (с) Жуков М.С.
Sommer
Молодой специалист

us
Offline Offline

« Ответ #2 : 13-08-2003 19:40 » 

RXL,
вопрос в дискретной свертке
в алгоритме ее реализации посредством FFT
если занимаешься этим - очень бы хотелось пообщаться
я не совсем понимаю как реализовать сам алгоритм дискретной сертки..
пользую FFT с www.fftw.org
Записан

когда-нибудь, я верю, ты будешь ехать по этому городу и поймёшь, что хочешь увидеть меня за рулём мчащейся по соседней полосе машины.
но тогда меня уже не будет
в этом городе
forever yours.
RXL
Технический
Администратор

ru
Offline Offline
Пол: Мужской

WWW
« Ответ #3 : 14-08-2003 07:12 » 

Свертка (convolution) и FFT - вещи разные.
Свертка делает умножение двух сигналов, а FFT разложение на состовляющие синусоиды.

Что именно ты хочешь проинтегрировать? Сам сигнал, или некоторую составляющую?
Записан

... мы преодолеваем эту трудность без синтеза распределенных прототипов. (с) Жуков М.С.
Sommer
Молодой специалист

us
Offline Offline

« Ответ #4 : 14-08-2003 17:04 » 

задача состоит в реализации дискретной свертки посредством FFT - об этом можно почитать у Баскакова в РЦИС

дело в том, что коэффициенты преобразования Фурье свертки являются произведениями коэффициентов ДПФ сввертываемых сигналов

моя проблема пока состоит в понимании как посчитать отсчеты свертываемых сигналов, если помнишь там у одного из них 2 параметра( обычно t и тау)
Записан

когда-нибудь, я верю, ты будешь ехать по этому городу и поймёшь, что хочешь увидеть меня за рулём мчащейся по соседней полосе машины.
но тогда меня уже не будет
в этом городе
forever yours.
RXL
Технический
Администратор

ru
Offline Offline
Пол: Мужской

WWW
« Ответ #5 : 15-08-2003 09:20 » 

Чесно говоря, я почти ни чего не понял... Видимо, мы разную литературу читаем  Отлично

Мне кажется, что я бы быстрее и правильнее понял, если бы ты сказал своими словами, без математики, тау и прочих китов.  Ага
Ты хочешь сделать свою реализацию FFT?
Записан

... мы преодолеваем эту трудность без синтеза распределенных прототипов. (с) Жуков М.С.
Sommer
Молодой специалист

us
Offline Offline

« Ответ #6 : 15-08-2003 10:10 » 

RXL, есть какая-то программа, которая вычисляет кучу вложенных интегнралов месяцы

моя задача используя FFT (его писать ИМХО самому бессмысленно - легче и умнее взять с www.fftw.org) решать те же задачи во много раз быстрее

используется свойство дискретной сверки(интеграл произведения 2 - х функций) - ее можно представиь через произведение коэф - в дискретного преобразования Фурье тех сигнналов, интеграл проризведения которых мы берем

ты сам чем конкретно занимаешься?
меня интересует алгоритм программирования этой дискретной свертки
Записан

когда-нибудь, я верю, ты будешь ехать по этому городу и поймёшь, что хочешь увидеть меня за рулём мчащейся по соседней полосе машины.
но тогда меня уже не будет
в этом городе
forever yours.
Sommer
Молодой специалист

us
Offline Offline

« Ответ #7 : 15-08-2003 10:14 » 

RXL,
нужен алгоритм дискретной свертки через FFT

FFT уже есть
Записан

когда-нибудь, я верю, ты будешь ехать по этому городу и поймёшь, что хочешь увидеть меня за рулём мчащейся по соседней полосе машины.
но тогда меня уже не будет
в этом городе
forever yours.
RXL
Технический
Администратор

ru
Offline Offline
Пол: Мужской

WWW
« Ответ #8 : 15-08-2003 10:26 » 

Я читаю литературу отсюда: http://www.dspguide.com/
Точнее: http://www.dspguide.com/pdfbook.htm
А DFT детально описывается тут:
http://www.dspguide.com/ch8.htm
http://www.spectrumsdi.com/ch8.pdf
На нем я пока и остановился.
-------
Мне ранее приходилось сталкиваться с вопросами обнаружения сигналов в телефонии, а сейчас решил изучить вопрос детальнее.
Записан

... мы преодолеваем эту трудность без синтеза распределенных прототипов. (с) Жуков М.С.
RXL
Технический
Администратор

ru
Offline Offline
Пол: Мужской

WWW
« Ответ #9 : 15-08-2003 10:59 » 

Честно говоря, даже походив по инету, так и не понял что есть "свертка через FFT" - я себе пока это не могу представить. в общем, когда пойму, обязательно напишу.

Сама свертка представляет собой простую формулу: y[n]=sum(h*x[n-i]), где x - входной сигнал, h - "импульсный отклик" (impulse responce), y - выходной сигнал, n - число отсчетов в выходном сигнале, i - длина h[], длина входного сигнала n-i+1.
Записан

... мы преодолеваем эту трудность без синтеза распределенных прототипов. (с) Жуков М.С.
Sommer
Молодой специалист

us
Offline Offline

« Ответ #10 : 15-08-2003 11:36 » 

RXL,
дискретная свертка через FFT :
C[f[k]] = C[x[k]] * C[y[k]]
литературу, по преобразованиям Фурье я уже указывал - есть в любой библиотеке

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

я вроде разобрался в алгоритме...
сделаю - выложу - он простой
проблема была не втом как пограммировать, а как вообще его себе представить\осознать...

--
а как ты обнаруживал сигналы в телефонии?....
Записан

когда-нибудь, я верю, ты будешь ехать по этому городу и поймёшь, что хочешь увидеть меня за рулём мчащейся по соседней полосе машины.
но тогда меня уже не будет
в этом городе
forever yours.
RXL
Технический
Администратор

ru
Offline Offline
Пол: Мужской

WWW
« Ответ #11 : 15-08-2003 12:13 » 

Несколько дальше нашел перемещение свертки из time domain в frequency domain - похоже ты от этом говоришь. Действительно простая вещь - оба сигнала (одной длины - самый короткий дополняется нулями) проходят FFT, а потом просто перемножаются коэффициенты с одним индексом. Упомянается, что это полезно при большой длине импульсного отклика (рекомендуют для длины больше 32).
 Ты об этом говорил?
-----
Есть алгоритм Goertzel-я для обнаружения заранее известного набора частот. Значительно проще и быстрее FFT. Кроме того, в АОН-е (всмысле, тел. аппарате с запросом АОН-а) на Z80 был реализован тот же алгоритм, только для однобитного (!) сигнала. Время вычисления для 7 частот составляло около 15мс на процессоре в 4МГц.
Записан

... мы преодолеваем эту трудность без синтеза распределенных прототипов. (с) Жуков М.С.
Diletant
Помогающий

de
Offline Offline

« Ответ #12 : 15-08-2003 13:15 » 

Мне тяжело вмешиваться в спор столь квалифицированной компании, но как мне помнится из вузовской математики, интеграл Фурье от свертки двух функций равен произведению интегралов Фурье от каждой функции. На этом основан метод быстрого вычисления свертки: Сначала вычисляется прямое FFT каждой функции, коэффициенты перемножаются, и вычисляется обратное FFT. Это имелось ввиду?
При этом не совсем понятно при чем тут форум по С?
Записан
Sommer
Молодой специалист

us
Offline Offline

« Ответ #13 : 15-08-2003 14:50 » 

Diletant,
да именно это имелось в виду
я внес тему  именно в этот форум т.к. нужен алгоритм на С++

правда пока полемика вокруг понимания алгоритма впринципе Улыбаюсь

RXL, ясно, спасибо
Записан

когда-нибудь, я верю, ты будешь ехать по этому городу и поймёшь, что хочешь увидеть меня за рулём мчащейся по соседней полосе машины.
но тогда меня уже не будет
в этом городе
forever yours.
Diletant
Помогающий

de
Offline Offline

« Ответ #14 : 15-08-2003 17:24 » 

Цитата: mandela

 нужен алгоритм на С++



Алгоритм на языке программирования, IMHO, нонсенс. Или вопрос в том как комплексную арифметику на С реализовать? А так берется "алгоритм" на Фортране и аккуратно переписывается. Тем более вся программка строк этак в двести. Если же речь идет об FFT, оптимизированном под конкретный процессор, то тема обсуждалась на  www.shelek.com/club, там еще ссылки на специализированные сайты приводились
Записан
Sommer
Молодой специалист

us
Offline Offline

« Ответ #15 : 15-08-2003 17:57 » 

Diletant,
я наверное не правильно выразился
в действительности же нужен был алгоритм и его реализация на С++
что касается предыдущего форума, тот там тему создавал опять же я
но тогда вопорс стоял не о свертке, а о реализации FFT, как ты уже отметил.
Записан

когда-нибудь, я верю, ты будешь ехать по этому городу и поймёшь, что хочешь увидеть меня за рулём мчащейся по соседней полосе машины.
но тогда меня уже не будет
в этом городе
forever yours.
Страниц: [1]   Вверх
  Печать  
 

Powered by SMF 1.1.21 | SMF © 2015, Simple Machines