abclinuxu.cz AbcLinuxu.cz itbiz.cz ITBiz.cz HDmag.cz HDmag.cz abcprace.cz AbcPráce.cz
AbcLinuxu hledá autory!
Inzerujte na AbcPráce.cz od 950 Kč
Rozšířené hledání
×
dnes 15:22 | Nová verze

Byl vydán Mozilla Firefox 64.0. Přehled novinek v poznámkách k vydání a na stránce věnované vývojářům. Nejnovější verze tohoto webového prohlížeče přináší například ovládání více panelů, nebo správce úloh, který lze otevřít v nabídce Firefoxu > Více > Správce úloh, nebo napsáním about:performance do adresního řádku.

Ladislav Hagara | Komentářů: 0
dnes 13:00 | Zajímavý článek Ladislav Hagara | Komentářů: 0
včera 22:33 | Nová verze

Po 3 měsících vývoje od vydání verze 14 byla vydána nová stabilní verze 15 open source systému Nextcloud, forku ownCloudu, umožňujícího provoz vlastního cloudového úložiště. Přehled novinek i s náhledy v příspěvku na blogu. Pro vyzkoušení Nextcloudu je k dispozici demo.

Ladislav Hagara | Komentářů: 4
včera 18:00 | IT novinky

Počítačová hra Doom slaví 25 let. Společností id Software ji vydala 10. prosince 1993. Zahrát si ji lze například na Internet Archive.

Ladislav Hagara | Komentářů: 13
9.12. 23:55 | Zajímavý článek

Nakladatelství Raspberry Pi vydalo 244 stránkového průvodce pro úplné začátečníky s jednodeskovým počítačem Raspberry Pi The Official Raspberry Pi Beginner’s Guide (pdf). Programování ve visuálním programovacím jazyce Scratch je věnována nová příručka Code Club Book of Scratch Volume 1 (pdf). Vydáno bylo také třetí číslo časopisu věnovaného počítačovým hrám Wireframe (pdf).

Ladislav Hagara | Komentářů: 0
9.12. 23:44 | Nová verze

U příležitosti oslav jednoho roku prací na debianím balíčku, vyšlo GPXSee 7.0. Nová verze přináší zejména podporu vektorových map (Mapbox PBF) pomocí nově vzniklého Qt pluginu.

Martin Tůma | Komentářů: 8
7.12. 23:33 | Nová verze

UBports, nadace a komunita kolem Ubuntu pro telefony a tablety Ubuntu Touch (seznam podporovaných zařízení), oznámila vydání nové verze Ubuntu Touch OTA-6. Zdůrazněna jsou například vylepšení webového prohlížeče Morph.

Ladislav Hagara | Komentářů: 0
7.12. 21:33 | Komunita

Richard Hughes na svém blogu informuje, že se společnost AMI (American Megatrends Incorporated), největší dodavatel firmwarů a BIOSů (AMIBIOS, Aptio, ...), zapojila do projektu LVFS (Linux Vendor Firmware Service). Seznam projektem podporovaných zařízení by se měl tedy časem rozšířit.

Ladislav Hagara | Komentářů: 0
7.12. 20:11 | Komunita

Nadace pro svobodný software (FSF) informuje, že seznam doporučených svobodných distribucí GNU/Linuxu byl rozšířen o distribuci Hyperbola GNU/Linux-libre. Hyperbola vychází z Arch Linuxu a na rozdíl od něj splňuje požadovaná kritéria.

Ladislav Hagara | Komentářů: 7
7.12. 19:55 | IT novinky

Společnost Valve oznámila, že počítačová hra Counter-Strike: Global Offensive (Steam) je nově free to play. Hra Company of Heroes 2 (Steam) je do pondělí 10. prosince zdarma.

Ladislav Hagara | Komentářů: 0
Chystáte se přejít na Wayland na „desktopu“?
 (26%)
 (7%)
 (11%)
 (31%)
 (26%)
Celkem 104 hlasů
 Komentářů: 14, poslední včera 12:19
Rozcestník

Dotaz: Vyhladenie dát pomocou FFT

marulinko avatar 5.4.2006 18:30 marulinko | skóre: 12 | blog: Technohead
Vyhladenie dát pomocou FFT
Přečteno: 452×
Potreboval by som získať z rozhádzaných nameraných údajov vyhladenú krivku. Používam na spracovanie údajov Octave a na vykreslenie Gnuplot. Na vyhladenie by sa dalo použiť vyhladenie pomocou FFT (Fast Fourier Transformation). Máte niekto skúsenosti ako na to ? Už 2 deň hladám na nete, našiel som rôzne softy napr. BruteFIR ale nešlo mi to skompilovať. Neviem či by bol na to vhodný ale zo screenshotov by to možno fungovalo. A keďže neštudujem ani matematiku, presne ani neviem ako to vyhladzovanie vlastne funguje :(. Možno by sa to dalo aj naprogramovať. Potrebujem nakopnúť správnym smerom. Alebo doporučiť nejaký program. Díky za každú radu.
BruteFIR
FFTW
Octave
Gnuplot

Odpovědi

5.4.2006 20:54 fikusek | skóre: 1
Rozbalit Rozbalit vše Re: Vyhladenie dát pomocou FFT
jde o to jak ty data vypadaji?!?!
treba jsou to jednoducha data, ktera by slo prolozit polynomem.
nebo to bude nejaky signal, ktery se bude muset
vyhlazovat sloziteji.
nechej na sebe mejla a je to.
marulinko avatar 6.4.2006 10:06 marulinko | skóre: 12 | blog: Technohead
Rozbalit Rozbalit vše Re: Vyhladenie dát pomocou FFT
môj mail je majo(at)marulinko(dot)com
5.4.2006 22:39 Honza
Rozbalit Rozbalit vše Re: Vyhladenie dát pomocou FFT

A co přesně si pod pojmem vyhlazení dat představuješ? Pokud tím myslíš jejich vykreslení do grafu ve formě spojité křivky, tak se dá použít proložení nějakým polynomem jak už radil kolega výše.

FFT je pouze algoritmus výpočtu diskrétní Fourierovy transformace, která převádí řadu z časové oblasti do spektrální. Pokud jsou tvoje data signál a pod pojmem vyhlazení si představuješ například odstranění vyšších spektrálních složek, je možné FFT využít. Říká se tomu filtrace ve spektrální oblasti. Provedeš jednoduše FFT daného signálu, odstraníš část spektra, kterou nechceš (doporučuji zachovat spektrum pokud možno spojité) a modifikované spektrum pomocí IFFT transformuješ zpět do časové oblasti.

Jinak si jsem na 99% jistý, že FFT i IFFT budou v Octave implementované. Pravděpodobně to bude v nějaké knihovně DSP apod. Blíže neporadím, protože používám Scilab nebo Matlab.

5.4.2006 22:53 Jan Martinek | skóre: 43 | blog: johny | Brno
Rozbalit Rozbalit vše Re: Vyhladenie dát pomocou FFT
Těžko říct, jakou matematickou metodu bys měl pro svá data použít. Ale kdyby ses rozhodl pro FFT, tak existuje knihovna FFTW (kterou jsi zmínil) a používá ji spousta programů - i ten BruteFIR, který jsi uvedl. Základní funkce, která vypočítá jednorozměrnou FFT za použití FFTW může vypadat třeba takhle:
#include <math.h>
#include <stdlib.h>
#include <fftw3.h>

int calculate_fft(int N, double *in_dbl, double *out_dbl){
  int i;
  double re, im;


  fftw_plan p;
  double *in, *out;

  in = fftw_malloc(sizeof(double) * N);
  out = fftw_malloc(sizeof(double) * N);

  for(i=0; i<N; i++){
    in[i] = in_dbl[i]; /* Zero complex part */
  }


  p = fftw_plan_r2r_1d(N, in, out, FFTW_R2HC, FFTW_ESTIMATE);
  fftw_execute(p);
  fftw_destroy_plan(p);


  out_dbl[0] = out[0];

  for(i = 1; i < (N + 1) / 2; ++i){
    re = out[i];
    im = out[N-i];
    out_dbl[i] = sqrt(re*re+im*im);
  }

  fftw_free(in);
  fftw_free(out);

  return 0;
}
Nejdřív vypočte reálnou a imaginární část a pak z toho vezme absolutní hodnotu. Měl bych asi zmínit ještě spoustu dalších věcí, ale to už by mělo s linuxem pramálo společného.

Při překladu to musíš slinkovat s -lfftw3.
marulinko avatar 6.4.2006 10:02 marulinko | skóre: 12 | blog: Technohead
Rozbalit Rozbalit vše Re: Vyhladenie dát pomocou FFT
Díky za kód, len by som sa chcel spýtať či to zlinkovanie znamená ak dám gcc tento_kod -L/usr/include/fftw3.h ? Podobný problém som mal pri kompilácií brutefire. Tam mi chýba ešte súbor: -lfftw3f. -lfftw3 som vďaka tvojmu skriptu už našiel :)
6.4.2006 10:39 Jan Martinek | skóre: 43 | blog: johny | Brno
Rozbalit Rozbalit vše Re: Vyhladenie dát pomocou FFT
To, co jsem ti poslal, je pouze funkce, kterou když zavoláš, tak ti spočítá FFT. Bude se ti to hodit, až budeš mít napsaný celý svůj funkční kód v céčku a posledním chybějícím střípkem bude právě tato funkce. Jediné, co s tím můžeš v této chvíli udělat je zkompilovat do object souboru pomocí
gcc -c neco.c
Ten můžeš následně spojit se svým hotovým programem a teprve potom vznikne spusitelná binárka.

Linkují se knihovny, includují se hlavičkové soubory. Ne naopak ;-)

Upřímně řečeno - zkus raději tu fft funkci z octave. Na učení céčka si nech nějaký čas.
6.4.2006 07:39 LS_999
Rozbalit Rozbalit vše Re: Vyhladenie dát pomocou FFT
Dle meho nazoru je FFT na vyhlazeni dat je jako jit na mouchu s kulometem. A krome toho si myslim, ze diky efectum jako "spectral leaking" to bude vyhlazovat i tam, kde by nemelo (na "kraji"). Pokud nestudujete matematiku, usetrete si to utrpeni zapasit s fft. Doporucuju k vyhlazeni pouzit gnuplot, ktery je z octave vicemene pristupny pres prikaz plot atd. (viz jakysi serial na abclinuxu), ale mnohem lepe je pouzit zvlast. Lze k vyhlazeni pouzit spliny anebo pripadne fitovat nejakou vhodnou funkci, napr. polynom, jak uz nekdo psal anebo e^-x atd. Priklad priblizne v gnuplot:

Parabola:

fit a*x**2+b*x+c "mojedata.xy" using ($1):($2) via a,b,c

plot a*x**2+b*x+c

Vyhlazeni splinem: plot "mojedata.xy" using ($1):($2) smooth csplines

Vice v dokumentaci ke gnuplot (help fit, help plot)
marulinko avatar 6.4.2006 09:51 marulinko | skóre: 12 | blog: Technohead
Rozbalit Rozbalit vše Re: Vyhladenie dát pomocou FFT
Ďakujem za všetky odpovede. Tie dáta vyzerajú takto: tu je časť z nich
čas Mk F P
0.000 -383.789 200.195 0.000
0.020 -383.789 202.637 0.005
0.040 -351.807 205.078 -0.001
0.060 -351.807 204.773 -0.002
0.080 -199.890 -312.500 -0.001
0.100 -343.811 204.468 0.000
0.120 -679.626 199.585 -0.014
0.140 -351.807 202.484 -0.002
0.160 -343.811 183.563 0.009
0.180 -343.811 201.111 -0.004
A vykreslené to vyzerá takto: Graf (7.82kB)
Je to vstup zo snímačov do meracej karty a následne do PC.
Vyhladiť to potrebujem kôli eliminácií šumu a jednoduchšiemu odčítaniu hodnôt. Idem skúsiť tú funkciu fit z gnuplot.
môj mail je majo(at)marulinko(dot)com
6.4.2006 14:01 Jura
Rozbalit Rozbalit vše Re: Vyhladenie dát pomocou FFT
Hergot, bych tam proste od ruky(od mysi) nacmaral caru, aby to jako kopirovalo ten graf, napisu k temu FFT, a kdo pozna, ze ne? Se s tem moc seres :)))
6.4.2006 17:30 lada
Rozbalit Rozbalit vše Re: Vyhladenie dát pomocou FFT
ja bych na to sel takto:
udelal bych nekolik mereni takze bych dostal nekolik takovych krivek jak je na tvem obrazku (rekneme 10 nebo i 100 zalezi na narocnosti)
srovnal je tak, aby sedel zacatek toho skoku
udelal z nich prumer
vychutnaval si pohled na hladkou krivku
mozna nic moc ale rozhodne to funguje
27.7.2006 10:17 Petr "Glubo" Sýkora | skóre: 21 | blog: Glubnik
Rozbalit Rozbalit vše Re: Vyhladenie dát pomocou FFT
Rozhodně to funguje při "dostatečně náhodném" šumu ;)
„O mrtvých jen v dobrém." „Pojďme se bavit o Stalinovi."

Založit nové vláknoNahoru

Tiskni Sdílej: Linkuj Jaggni to Vybrali.sme.sk Google Del.icio.us Facebook

ISSN 1214-1267   www.czech-server.cz
© 1999-2015 Nitemedia s. r. o. Všechna práva vyhrazena.