Diskuze - danyk.cz

Vlákno z kategorie: Hlavní diskuze
Celkem 21 odpovědí.


25.04.2021 (20:10:06) K # IP X
loučák:
Zase budu otravovat s matlabem. Kdo jste toctu říkali, že jste dělali fft v matlabu? Dá se u té fce nějak nastavit rozsah a krok? Nepotřebuju analyzovat kmitočty v GHz ale za to potřebuju trochu jemnější krok než po 14Hz. Díky
26.04.2021 (11:06:35)  # IP X
VroutekB:
Dyk velikost binu máš snad jasně danou vzorkovacím kmitočtem a velikostí FFT. A je jedno, jestli matlab nebo jinde FFT...
26.04.2021 (13:43:18)  # IP X
loučák:
No to mě právě zarazilo, že ať mám na 70ms 14kSa nebo 7MSa tak furt vypadá začátek .csv nějak takle
0,14.286,28.571,42.857,57.143,71.429,85.714,100,114.29, prostě brutálně velikej krok a nevím jak je definovat jemněji
26.04.2021 (13:45:32)  # IP X
loučák:
a pak logicky když člověk chce analyzovat proud síťového zařízení, tam mu tam těch 50Hz a násobky dost chybějí. Popř. dá se zjišťovat amplituda jen na určených frekvencích, tedy nosná plus násobky?
26.04.2021 (13:48:30)  # IP X
loučák:
napsal jsem to kokotksy. Mám prostě časovou základnu 70ms a vemu si z ní 14k vzroků, pak 140k a pak 7M a to fft vypadá pořád stejně, hlavně stejný krok, těch cca 14Hz
26.04.2021 (17:51:35)  # IP X
VroutekB:
Tak ti znova říkám, že rozlišení binu FFT je přesně vzorkovací_f děleno velikost FFT.
Nevim čemu řikáš časová základna, chci vědět, jaká je vzorkovací frekvence a jaká velikost FFT Nic z toho co píšeš, totiž nedává smysl.
26.04.2021 (18:56:39)  # IP X
loučák:
Jo blbě jsem to napsal, vzorek je dlouhý 70ms. A vše mělo mělo směřovat k vzorkovací frekvenci 200MHz, tedy 140k vzorků. A čím je definovaná velikost fft?
26.04.2021 (19:41:47)  # IP X
StandaM (web) :
Boha jeho. :) Kdyz mas fs = 200e6 a N = 140kSa, tak bude krok:
f_bin = fs/N = 1.428kHz.
Vime? Asi nebudes analyzovat 50Hz se vzorkovanim 200e6. Tedy muzes, ale bude zaznam dlouhej jako prase a bude to i trvat a stejne z nej pouzijes jen malej kousek spektra ze zacatku.
27.04.2021 (13:45:42)  # IP X
loučák:
Dobře, chápu, dojebal jsem to. Tak teda další otázka. Nechce se mi jet extra do Ostravy abych to navzorkoval znovu. Je v matlabu nějaká fce, aby se z toho souboru vybrala každý 1000cí vzorek?
27.04.2021 (14:43:18)  # IP X
VroutekB:
Samozřejmě že je, řiká se tomu decimace.
Ale taky před její aplikací musíš použít nějaký loupás filtrování **02
27.04.2021 (17:10:25)  # IP X
loučák:
Tak na vyzobání 1000ciny dat z matice už jsem došel, průběh v čase vypadá ok, stejně jako v plne palbě, fft jako přimka přímé úměrnosti...tedy jak se prosím dělá to filtrování?
27.04.2021 (19:07:32)  # IP X
VroutekB:
Už si někdy slyšel o aliasingu? **02
27.04.2021 (19:14:29)  # IP X
loučák:
Ano, slyšel, nejvyšší analyzovaná f musí být nejvýše polivou f vzorkovací. Ale nevím jak to provést sw, chápu že osciloskop si tohle řeší sám a počítá s nastavenou vzorkovací f, nikoho nenapadne, že by se měla dodatečně 1000x zmenšit.
27.04.2021 (19:20:12)  # IP X
loučák:
A pokud to chápu dobře, tak je nutné prohnat dolní propustí signál ještě před změnou vzorkovací f. Což hádám bude fest oříšek v matlabu a míň času mi zabere otočka do Ostravy. Jenže se obávám, že ten škopek ani tak nízké vzorkovací f neumí, takže mě čeká návrat k matlabu. Oser toto, nejvíc mě sere, že je to práce pro práci, tenhle vyjebanej způsob mám přiřazenej jen proto, abych měl té práce co nejvíc, jinak bych to prej měl moc jednoduchý. Takle mi to bylo otevřeně řečeno.
27.04.2021 (19:24:29)  # IP X
StandaM (web) :
28.04.2021 (08:53:20)  # IP X
loučák:
Díky za tip, zkusím si poradit.
Prej pojď na energetiku, říkali. Vyhneš se programování, říkali...
29.04.2021 (16:12:06)  # IP X
VroutekB:
A matlab je programování od kedy? **02
29.04.2021 (16:20:11)  # IP X
RayeR (web) :
Tak asi od doby co basic, python, JS je taky programovani... :)
30.04.2021 (17:19:13)  # IP X
DR. DOKTOR:
Důkladně jsem to nečetl a možná budu plácat blbosti, ale nebylo by řešením prostě vzít přesně 3 celé periody, tedy 60 ms, pak rozlišení FFT bude 50/3 Hz a zachytí tedy přesně 50 Hz se všemi harmonickými
1.05.2021 (22:00:55)  # IP X
loučák:
No vlastně s délkou vzorku mě nenapadlo si pohrát. Díky za tip, prostě se asi tomu novému měření neubráním
3.05.2021 (15:51:23)  # IP X
loučák:
Díky za rady, pomohli mi ke snad přijatelnému výsledku. A mám poslední blbou otázku. Proč je matice výstupu fft komplexní? Tipnu si, imaginární hodnota značí fázové spektrum?
3.05.2021 (19:51:12)  # IP X
StandaM (web) :
Ona je komplexni i ta vstupni matice. Jen se to moc nepouziva. Vystupem jsou pochopitelne re-im, protoze jak rikas amplituda-faze. Prepocet z kartezskych na polarni:
Y_mag = abs(Y);
Y_arg = angle(Y);
A taky nezapomen, ze pul spektra jsou negativni kmitocty - zahodit. Pro real-y vstup jsou obe pulky komplexne sdruzene (a vzajemne v opacnem poradi).
Přezdívka:
Heslo:
Text: