RTV forum PL | NewsGroups PL

Jak efektywnie obliczyć transformację FFT dla 15000 punktów bez całek?

FFT

NOWY TEMAT

elektroda NewsGroups Forum Index - Elektronika Polska - Jak efektywnie obliczyć transformację FFT dla 15000 punktów bez całek?

Goto page Previous  1, 2, 3  Next

J.F.
Guest

Fri Nov 24, 2006 12:53 pm   



On Fri, 24 Nov 2006 12:43:48 +0100, PAndy wrote:
Quote:
"J.F." <jfox_xnospamx@poczta.onet.pl> wrote in message
On Thu, 23 Nov 2006 22:23:10 +0100, seba wrote:
chcialbym obliczyc FFT dla np 15000 punktow

Raczej dla 16384, choc sa i dla niekraglych liczb.

Dostawia sie zera i po klopocie

Wbrew pozorom to jest jednak spory klopot ..

J.

PAndy
Guest

Fri Nov 24, 2006 3:22 pm   



"J.F." <jfox_xnospamx@poczta.onet.pl> wrote in message
news:i7ndm29q9kmo6deo1k7binkra54sg0ti9p@4ax.com...
Quote:
On Fri, 24 Nov 2006 12:43:48 +0100, PAndy wrote:
"J.F." <jfox_xnospamx@poczta.onet.pl> wrote in message
On Thu, 23 Nov 2006 22:23:10 +0100, seba wrote:
chcialbym obliczyc FFT dla np 15000 punktow

Raczej dla 16384, choc sa i dla niekraglych liczb.

Dostawia sie zera i po klopocie

Wbrew pozorom to jest jednak spory klopot ..

nie przesadzalbym...

J.F.
Guest

Fri Nov 24, 2006 3:36 pm   



On Fri, 24 Nov 2006 15:22:07 +0100, PAndy wrote:
Quote:
"J.F." <jfox_xnospamx@poczta.onet.pl> wrote in message
Raczej dla 16384, choc sa i dla niekraglych liczb.
Dostawia sie zera i po klopocie
Wbrew pozorom to jest jednak spory klopot ..

nie przesadzalbym...

Dostawic zaden, ale zdajesz sobie sprawe ze dalej liczysz
transformate czegos innego ?

J.

BartekK
Guest

Fri Nov 24, 2006 4:14 pm   



J.F. napisał(a):
Quote:
On Thu, 23 Nov 2006 23:36:30 +0100, BartekK wrote:
http://www.archelon.com/fft.html
A moze masz gdzies pod reka gotowca, ktory by poszedl bez floatow,
najlepiej bez liczenia cosinusow w locie (chocby z tablicy) - zalezy mi
glownie na predkosci, jakby sie udalo to na AVR bym chcial to robic :)

Przeciez tam wszystkiego jest jeden cosinus i to stalej ..
?

Moze ja cos zle widze, a z pewnoscia sie nie znam, ale widze tam
obliczenia na float (jak to odpale na avr to wyjdzie mi kolos z
predkoscia programatora do pralki)

--
| Bartlomiej Kuzniewski
| sibi@drut.org GG:23319 tel +48 696455098 http://drut.org/
| http://www.allegro.pl/show_user_auctions.php?uid=338173

J.F.
Guest

Fri Nov 24, 2006 5:07 pm   



On Fri, 24 Nov 2006 16:14:36 +0100, BartekK wrote:
Quote:
J.F. napisał(a):
On Thu, 23 Nov 2006 23:36:30 +0100, BartekK wrote:
http://www.archelon.com/fft.html
A moze masz gdzies pod reka gotowca, ktory by poszedl bez floatow,
najlepiej bez liczenia cosinusow w locie (chocby z tablicy) - zalezy mi
glownie na predkosci, jakby sie udalo to na AVR bym chcial to robic :)

Przeciez tam wszystkiego jest jeden cosinus i to stalej ..

?
Moze ja cos zle widze, a z pewnoscia sie nie znam, ale widze tam
obliczenia na float (jak to odpale na avr to wyjdzie mi kolos z
predkoscia programatora do pralki)

Obliczenia sa na float, bo taka natura tej transformaty - trzeba
mnozyc przez liczby z zakresu -1 ... +1, i to dosc dokladnie.

Mozesz probowac liczyc to staloprzecinkowo, ale trzeba sie bedzie
troche wysilic nad projektem, i wolno nadal bedzie.

Ale cosinusow tam praktycznie nie ma.

J.

PAndy
Guest

Fri Nov 24, 2006 5:24 pm   



"J.F." <jfox_xnospamx@poczta.onet.pl> wrote in message
news:5p0em2ts2h3hcmsig66cp1fa1fdb36n8m9@4ax.com...
Quote:
On Fri, 24 Nov 2006 15:22:07 +0100, PAndy wrote:
"J.F." <jfox_xnospamx@poczta.onet.pl> wrote in message
Raczej dla 16384, choc sa i dla niekraglych liczb.
Dostawia sie zera i po klopocie
Wbrew pozorom to jest jednak spory klopot ..

nie przesadzalbym...

Dostawic zaden, ale zdajesz sobie sprawe ze dalej liczysz
transformate czegos innego ?

Przepraszam - metode uzupelnienia zeremi przy FFT podaje sie praktycznie
w kazdym podreczniku DSP...

J.F.
Guest

Fri Nov 24, 2006 5:53 pm   



On Fri, 24 Nov 2006 17:24:01 +0100, PAndy wrote:
Quote:
"J.F." <jfox_xnospamx@poczta.onet.pl> wrote in message
Dostawia sie zera i po klopocie
Wbrew pozorom to jest jednak spory klopot ..
nie przesadzalbym...

Dostawic zaden, ale zdajesz sobie sprawe ze dalej liczysz
transformate czegos innego ?

Przepraszam - metode uzupelnienia zeremi przy FFT podaje sie praktycznie
w kazdym podreczniku DSP...

A pisza jak potem interpretowac wyniki ?

J.

PAndy
Guest

Fri Nov 24, 2006 6:47 pm   



"J.F." <jfox_xnospamx@poczta.onet.pl> wrote in message
news:er8em29d0c7okqn4efislrrks7k597pbe5@4ax.com...
Quote:
On Fri, 24 Nov 2006 17:24:01 +0100, PAndy wrote:
"J.F." <jfox_xnospamx@poczta.onet.pl> wrote in message
Dostawia sie zera i po klopocie
Wbrew pozorom to jest jednak spory klopot ..
nie przesadzalbym...

Dostawic zaden, ale zdajesz sobie sprawe ze dalej liczysz
transformate czegos innego ?

Przepraszam - metode uzupelnienia zeremi przy FFT podaje sie
praktycznie
w kazdym podreczniku DSP...

A pisza jak potem interpretowac wyniki ?

http://zone.ni.com/devzone/cda/tut/p/id/4880

J.F.
Guest

Fri Nov 24, 2006 6:53 pm   



On Fri, 24 Nov 2006 18:47:35 +0100, PAndy wrote:
Quote:
"J.F." <jfox_xnospamx@poczta.onet.pl> wrote in message
Przepraszam - metode uzupelnienia zeremi przy FFT podaje sie
praktycznie
w kazdym podreczniku DSP...

A pisza jak potem interpretowac wyniki ?

http://zone.ni.com/devzone/cda/tut/p/id/4880

"Unfortunately, it can lead to results that are easy to misinterpret"

J.

PAndy
Guest

Fri Nov 24, 2006 7:14 pm   



"J.F." <jfox_xnospamx@poczta.onet.pl> wrote in message
news:8ccem290dp6car7l0oet4q4mtlkk99e7rd@4ax.com...
Quote:
On Fri, 24 Nov 2006 18:47:35 +0100, PAndy wrote:
"J.F." <jfox_xnospamx@poczta.onet.pl> wrote in message
Przepraszam - metode uzupelnienia zeremi przy FFT podaje sie
praktycznie
w kazdym podreczniku DSP...

A pisza jak potem interpretowac wyniki ?

http://zone.ni.com/devzone/cda/tut/p/id/4880

"Unfortunately, it can lead to results that are easy to misinterpret"

Zgadza sie ale nie pisalem ze to metoda ktora jest panaceum na wszystko
i mozna ja stosowac bez zrozumienia...
natomiast jest alternatywa dla DFT w wielu sytuacjach gdy nie mamy
mozliwosci zrobienia zwyczjanie FFT bo mamy za malo probek...

Wojciech Bartczak
Guest

Fri Nov 24, 2006 9:50 pm   



Dnia Thu, 23 Nov 2006 22:23:10 +0100
"seba" <gennex@o2.pl> napisał(a):

Quote:
[...]

Hej !!

Proponuje zajrzec do ksiazki o matematyce dyskretnej. Ogolnie nie
jest to trudne. Pozatym calke mozna zastapic z pewna dokladnoscia
skonczona suma skladnikow Smile. Polcam jakas ksiazke o metodach
numerycznych np. Bjorck i Dalquist "Metody numeryczne" tam powinno
sie to znalezc. Wszystko bedzie jasne, a sam algorytm tez zrobi sie
duzo prostszy niezaleznie od jezyka programowania.

Pozdrawiam,

Wojciech Bartczak

Adam Dybkowski
Guest

Sat Nov 25, 2006 12:02 am   



BartekK napisał(a):

Quote:
http://www.archelon.com/fft.html

Moze ja cos zle widze, a z pewnoscia sie nie znam, ale widze tam
obliczenia na float (jak to odpale na avr to wyjdzie mi kolos z
predkoscia programatora do pralki)

Możesz przyjąć skończoną dokładność i robić obliczenia na liczbach
stałoprzecinkowych. Np. przeskalować wszystkie dane wejściowe do zakresu
<-1,1) i zapisać je jako liczy 16 bitowe (1 bit znaku i 15 bitów części
ułamkowej). Obliczenia na liczbach stałoprzecinkowych robi się szybko i
przyjemnie, tak pracuje większość obecnie produkowanych procesorów
sygnałowych (DSP), z czego Texas Instruments robi z 80% z nich.

--
Adam Dybkowski
http://www.amwaw.edu.pl/~adybkows/

Uwaga: przed wysłaniem do mnie maila usuń cyfry z adresu.

Sebastian Bialy
Guest

Sat Nov 25, 2006 12:53 am   



seba wrote:
Quote:
Drodzy Panstwo,chcialbym programowo obliczyc transformate FFT,znam troche
programowanie ale nie wiem jakiego zapisu uzyc,tak zeby obliczanie bylo w
miare szybkie

A w jakim jezyku? Jeśli w C to FFT najlepiej liczyć chyba fftw.
http://www.fftw.org/. Widziałem stosowaną na AVR i sam stosowałem na
różnych procesorach. Zazwyczaj 0 kłopotów, to naprawdę dość solidny
kawałek kodu i raczej sam nie wymyślisz szybciej.

Piotr Wyderski
Guest

Sat Nov 25, 2006 1:42 pm   



PAndy wrote:

Quote:
Raczej dla 16384, choc sa i dla niekraglych liczb.

Dostawia sie zera i po klopocie

Po klopocie, ale w sytuacjach wymagajacych zgrubnego
podzialu na podpasma, takich jak filtrowanie itp. Jesli
jednak chce sie policzyc _DFT_, to po dostawieniu zer
wyniki beda opisywaly cos innego, niz dana funkcja.

Quote:
A ja ze swej strony polecam sprzetowe FFT

A po co to liczyc sprzetowo, skoro algorytm jest latwy,
dobrze implementowalny i wydajny, zwlaszcza przy wsparciu
sprzetowym do adresowania w odwróconym porzadku bitowym?

Pozdrawiam
Piotr Wyderski

BartekK
Guest

Sat Nov 25, 2006 3:53 pm   



Adam Dybkowski napisał(a):
Quote:
BartekK napisał(a):
http://www.archelon.com/fft.html
Moze ja cos zle widze, a z pewnoscia sie nie znam, ale widze tam
obliczenia na float (jak to odpale na avr to wyjdzie mi kolos z
predkoscia programatora do pralki)
Możesz przyjąć skończoną dokładność i robić obliczenia na liczbach
stałoprzecinkowych. Np. przeskalować wszystkie dane wejściowe do zakresu
-1,1) i zapisać je jako liczy 16 bitowe (1 bit znaku i 15 bitów części
ułamkowej).
Myslalem o tym, ale dla danych wejsciowych o rozdzielczosci 12bit (a

docelowo raczej 16bit) skalowanie do 16bit nie byloby sensowne chyba
(przeciez i tak juz maja zakres 0-0xffff) wiec mozna po prostu przyjac
ze to co jest zapisane jako 0-65535 to jest -1 do +1. Tylko ze jak
nakarmilem ten algorytm (tak na piechote na pc) dosc losowymi danymi o
rozpietosci 16bitowej to floaty zaczely mi siegac od bardzo malych
(ktore nijak na 16bit -1 do +1 bym nie zapisal na 1bicie) do bardzo
duzych - ktore by mi sie w 32bit nie zmiescily... Sproboje jeszcze
przerobic to wszystko na int, i zobaczyc czy z poobcinania bitow
cokolwiek wyjdzie, czy wszystkie wspolczynniki powpadaja w nasycenie od
dolu lub gory.

--
| Bartlomiej Kuzniewski
| sibi@drut.org GG:23319 tel +48 696455098 http://drut.org/
| http://www.allegro.pl/show_user_auctions.php?uid=338173

Goto page Previous  1, 2, 3  Next

elektroda NewsGroups Forum Index - Elektronika Polska - Jak efektywnie obliczyć transformację FFT dla 15000 punktów bez całek?

NOWY TEMAT

Regulamin - Zasady uzytkowania Polityka prywatnosci Kontakt RTV map News map