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 01:22 | Nová verze

    Byla vydána (𝕏) zářijová aktualizace aneb nová verze 1.105 editoru zdrojových kódů Visual Studio Code (Wikipedie). Přehled novinek i s náhledy a videi v poznámkách k vydání. Ve verzi 1.105 vyjde také VSCodium, tj. komunitní sestavení Visual Studia Code bez telemetrie a licenčních podmínek Microsoftu.

    Ladislav Hagara | Komentářů: 0
    včera 15:33 | Komunita

    Ve Firefoxu bude lepší správa profilů (oddělené nastavení domovské stránky, nastavení lišt, instalace rozšíření, uložení hesla, přidání záložky atd.). Nový grafický správce profilů bude postupně zaváděn od 14.října.

    Ladislav Hagara | Komentářů: 0
    včera 12:44 | Nová verze

    Canonical vydal (email) Ubuntu 25.10 Questing Quokka. Přehled novinek v poznámkách k vydání. Jedná se o průběžné vydání s podporou 9 měsíců, tj. do července 2026.

    Ladislav Hagara | Komentářů: 0
    včera 12:22 | Nová verze

    ClamAV (Wikipedie), tj. multiplatformní antivirový engine s otevřeným zdrojovým kódem pro detekci trojských koní, virů, malwaru a dalších škodlivých hrozeb, byl vydán ve verzi 1.5.0.

    Ladislav Hagara | Komentářů: 0
    včera 01:22 | Nová verze

    Byla vydána nová verze 1.12.0 dynamického programovacího jazyka Julia (Wikipedie) určeného zejména pro vědecké výpočty. Přehled novinek v příspěvku na blogu a v poznámkách k vydání. Aktualizována byla také dokumentace.

    Ladislav Hagara | Komentářů: 0
    8.10. 15:11 | Bezpečnostní upozornění

    V Redisu byla nalezena a v upstreamu již opravena kritická zranitelnost CVE-2025-49844 s CVSS 10.0 (RCE, vzdálené spouštění kódu).

    Ladislav Hagara | Komentářů: 5
    8.10. 14:00 | IT novinky

    Ministr a vicepremiér pro digitalizaci Marian Jurečka dnes oznámil, že přijme rezignaci ředitele Digitální a informační agentury Martina Mesršmída, a to k 23. říjnu 2025. Mesršmíd nabídl svou funkci během minulého víkendu, kdy se DIA potýkala s problémy eDokladů, které některým občanům znepříjemnily využití možnosti prokázat se digitální občankou u volebních komisí při volbách do Poslanecké sněmovny.

    Ladislav Hagara | Komentářů: 18
    8.10. 12:33 | Zajímavý software

    Společnost Meta představila OpenZL. Jedná se o open source framework pro kompresi dat s ohledem na jejich formát. Zdrojové kódy jsou k dispozici na GitHubu.

    Ladislav Hagara | Komentářů: 0
    8.10. 03:33 | IT novinky

    Google postupně zpřístupňuje českým uživatelům Režim AI (AI Mode), tj. nový režim vyhledávání založený na umělé inteligenci. Režim AI nabízí pokročilé uvažování, multimodalitu a možnost prozkoumat jakékoliv téma do hloubky pomocí dodatečných dotazů a užitečných odkazů na weby.

    Ladislav Hagara | Komentářů: 0
    7.10. 18:11 | Nová verze

    Programovací jazyk Python byl vydán v nové major verzi 3.14.0. Podrobný přehled novinek v aktualizované dokumentaci.

    Ladislav Hagara | Komentářů: 2
    Jaké řešení používáte k vývoji / práci?
     (38%)
     (46%)
     (15%)
     (17%)
     (21%)
     (15%)
     (17%)
     (15%)
     (15%)
    Celkem 203 hlasů
     Komentářů: 13, poslední 8.10. 07:41
    Rozcestník

    Dotaz: fast fourier v Lua

    25.4.2022 10:33 martin
    fast fourier v Lua
    Přečteno: 638×
    Ahoj, Potřebuji využít v jedné aplikaci FFT.
    Něco jsem našel na Rosetta code.
    FFT funguje ale inverzní FFT vrací něco jiného.
    help needed :(
    díky
    complex = {__mt={} }
    
    function complex.new(r, i)
      local new={real=r, imaginary=i or 0}
      setmetatable(new,complex.__mt)
      return new
    end
    
    function complex.__mt.__add(c1, c2)
      return complex.new(c1.real + c2.real, c1.imaginary + c2.imaginary)
    end
    
    function complex.__mt.__sub(c1, c2)
      return complex.new(c1.real - c2.real, c1.imaginary - c2.imaginary)
    end
    
    function complex.__mt.__mul(c1, c2)
      return complex.new(c1.real * c2.real - c1.imaginary * c2.imaginary,
    					 c1.real * c2.imaginary + c1.imaginary * c2.real)
    end
    
    function complex.expi(i)
      return complex.new(math.cos(i),math.sin(i))
    end
    
    function complex.__mt.__tostring(c)
      return "("..c.real..","..c.imaginary..")"
    end
    
    
    --[[---------------------------------------------------------------------
    Cooley–Tukey Fast Fourier Transformation                                -
    Zdroj: https://en.wikipedia.org/wiki/Cooley-Tukey_FFT_algorithm         -
    ]]-----------------------------------------------------------------------
    
    function FFT(vect)
      local n=#vect
      if n<=1 then return vect end
    
      local odd,even={},{}
      for i=1,n,2 do
        odd[#odd+1]=vect[i]
        even[#even+1]=vect[i+1]
      end
    
      FFT(even)
      FFT(odd)
    
      for k=1,n/2 do
        local t=even[k] * complex.expi(-2*math.pi*(k-1)/n)
        vect[k] = odd[k] + t
        vect[k+n/2] = odd[k] - t
      end
      return vect
    end
    
    --[[-----------------------------------------------------------------------
    Inverse Fast Fourier Transformation                                       -
    ]]-------------------------------------------------------------------------
    
    function IFFT(amplitudes)
    
    	local N = #amplitudes
    	local input = 1 / N
    
    	local i=nil
    	for i=1, N do
    		amplitudes[i].imaginary = -amplitudes[i].imaginary
    	end
    
    	FFT(amplitudes)
    
    	for i=1, N do
    		amplitudes[i].imaginary = -amplitudes[i].imaginary
    		amplitudes[i].real = amplitudes[i].real * input
    		amplitudes[i].imaginary = amplitudes[i].imaginary * input
    	end
    
    	return amplitudes
    end
    
    function toComplex(vector)
      vect={}
      for r,i in ipairs(vector) do
        vect[r]=complex.new(i)
      end
      return vect
    end
    
    -- test
    data = toComplex{1, 1, 1, 1, 0, 0, 0, 0}
    
    print("orig:", unpack(data))
    print("fft :", unpack(FFT(data)))
    print("ifft:", unpack(IFFT(data)))
    

    Řešení dotazu:


    Odpovědi

    Gréta avatar 25.4.2022 19:48 Gréta | skóre: 37 | blog: Grétin blogísek | 🇮🇱==❤️ , 🇵🇸==💩 , 🇪🇺==☭
    Rozbalit Rozbalit vše Re: fast fourier v Lua

    neumim dělat lulu :D :D :D :D ale ten tvuj lulu zdrojáček mi vrací

    orig:   (1,0)   (1,0)   (1,0)   (1,0)   (0,0)   (0,0)   (0,0)   (0,0)
    fft :   (4,0)   (1,-2.4142135623731)    (0,0)   (1,-0.41421356237309)   (0,0)   (1,0.41421356237309)    (0,0)   (1,2.4142135623731)
    ifft:   (1,-0)  (1,-5.5511151231258e-17)        (1,2.4894981252574e-17) (1,-5.5511151231258e-17)        (5.5511151231258e-17,0) (5.5511151231258e-17,5.5511151231258e-17) (0,-2.4894981252574e-17) (5.5511151231258e-17,5.5511151231258e-17)
    

    noa to é na mínus sedmnáctou (jako nakalkulačce hele :O ;D) sou tak děsně mrňavý čísla že je asi jako mužem považovat za nuly takže to máš asi jako dobře :D ;D

    Gréta avatar 25.4.2022 19:49 Gréta | skóre: 37 | blog: Grétin blogísek | 🇮🇱==❤️ , 🇵🇸==💩 , 🇪🇺==☭
    Rozbalit Rozbalit vše Re: fast fourier v Lua

    joa ty mrňavý čisilka se při výpočtech berou z nepřesnosti s děláním s floatama/číslama s plovoucí desetinou tečkou hele :O ;D

    25.4.2022 22:33 z_sk | skóre: 34 | blog: analyzy
    Rozbalit Rozbalit vše Re: fast fourier v Lua
    Čísla ináč zapísané. Viď: Vědecký zápis čísel.
    debian.plus@protonmail.com
    26.4.2022 02:56 .
    Rozbalit Rozbalit vše Re: fast fourier v Lua
    On nemluví o zápisu, ale o hodnotě, hňupe.
    Řešení 1× (Gréta)
    26.4.2022 05:44 martin
    Rozbalit Rozbalit vše Re: fast fourier v Lua
    Ahoj, nějak se mi tu potratil muj další příspěvek
    Nebo jsem dal jen náhled a pak ho neodeslal :/
    co jsem psal hned po tom.
    Bylo to ve spěchu takže jsem si nevšiml že číslo končí
    s e-017. Ve své slepotě jsem viděl jen -5,587543235696549
    To je celé.
    Díky za reakce.

    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.