Re: protonovy magnetometr, slaby signal

Vojtěch Petrucha petrucha na volny.cz
Čtvrtek Červen 3 08:24:39 CEST 2021


zdravim,
to zní dobře, byl byste prosím ochotrný se podělit o kód použitý pro níže uvedený test?
já jsem na tom mizerně s časem a je zcela evidentní, že matematice i programování rozumíte minimálně řádově lépe než já..
díky
v.
 
 
 
______________________________________________________________
> Od: "Miroslav Mraz" <mrazik na volny.cz>
> Komu: hw-list na list.hw.cz
> Datum: 01.06.2021 18:39
> Předmět: Re: protonovy magnetometr, slaby signal
>
Tak jsem si s tím po večerech trochu hrál, protože mě to docela zajímá. 
 Našel jsem nějaká stará data z roku 2018, zkusil jsem na to pustit tu 
 Hilbertovu transformaci ale moc mě to nenadchlo. Tak jsem si opět pohrál 
 s Fourierovou transformací. Podle obrázků na vašem webu (a moje 
 zkušenost to potvrzuje) je šířka hlavní spektrální čáry nepřímo úměrná 
 době měření, přitom celkem nezáleží na tom jaký je řád FFT. Jde tedy o 
 to jak najít co nejpřesněji její vrchol aniž by bylo nutné použít 
 extrémně velký řád FFT.
 Co se mi osvědčilo je použít FFT třeba jen 2^12 jako startovací metodu. 
 Dostanu tak hlavní frekvenci s nepřesností cca 10 Hz. Ten 20 Hz 
 (+/-10Hz) interval kolem ní pak rozsekám na 8. částí a pro každou z 
 těchto frekvencí vypočtu výkon signálu Goertzelovým algoritmem. Určím 
 maximum a mám novou hlavní frekvenci s nepřesností 20/8 Hz. To pak 
 opakuji dokud nedostanu vrchol s požadovanou přesností. Otázkou zůstává 
 jaká je skutečná nejistota měření. Ale z toho co mám k dispozici - asi 5 
 měření v jednom souboru těsně za sebou to vypadá, že jednotlivé vrcholy 
 se liší o 0.01 až 0.02 Hz a to není tak špatné.
 A co je na tom ještě dobré - není potřeba dělat FFT vysokého řádu, stačí 
 40 průchodů Goertzelovým algoritmem, ten je sice jednoduchý, ale pro 
 tuto požadovanou přesnost je nutné ho počítat v double. Asi by to šlo 
 počítat i v mcu, buffer pouze musí pobrat těch pár sekund vzorků - mohou 
 být 16.bit integer, na double se to převede až při výpočtu.
 
 Mrazík
 
 Dne 30. 05. 21 v 13:44 Vojtěch Petrucha napsal(a):
 > díky za komentář, na pc jsme to zkoušeli (teda konkrétně student), viz 
 > data na tom webu níže, fungovalo dobře, teď hledám jestli/jak to lze 
 > jednoduše udělat ve štole, kam chci magnetometr šoupnout k vyvíjenému 
 > nízkošumovému fluxgate, přičemž ideální řešení je levné, jednoduché, 
 > spolehlivé, mnoho proudu nežeroucí...  ..a nezabere příliš lidských 
 > zdrojů, kterých se v poslední době silně nedostává..
 > 
 > přinejhorším se bude frekvence měřit jen čítačem v stm32 určováním 
 > periody signálu z komparátoru...
 > 
 > v.
 > 
 _______________________________________________
 HW-list mailing list  -  sponsored by www.HW.cz
 Hw-list na list.hw.cz
 http://list.hw.cz/mailman/listinfo/hw-list <http://list.hw.cz/mailman/listinfo/hw-list>

------------- další část ---------------
HTML příloha byla odstraněna...
URL: <http://list.hw.cz/pipermail/hw-list/attachments/20210603/d96fe499/attachment.html>


Další informace o konferenci Hw-list