<p style="padding:0 0 0 0; margin:0 0 0 0;">zdravim,</p>

<p style="padding:0 0 0 0; margin:0 0 0 0;">to zní dobře, byl byste prosím ochotrný se podělit o kód použitý pro níže uvedený test?</p>

<p style="padding:0 0 0 0; margin:0 0 0 0;">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á..</p>

<p style="padding:0 0 0 0; margin:0 0 0 0;">díky</p>

<p style="padding:0 0 0 0; margin:0 0 0 0;">v.</p>

<p style="padding:0 0 0 0; margin:0 0 0 0;"> </p>

<p style="padding:0 0 0 0; margin:0 0 0 0;"> </p>

<p style="padding:0 0 0 0; margin:0 0 0 0;"> </p>

<p style="padding:0 0 0 0; margin:0 0 0 0;">______________________________________________________________<br />
> Od: "Miroslav Mraz" <mrazik@volny.cz><br />
> Komu: hw-list@list.hw.cz<br />
> Datum: 01.06.2021 18:39<br />
> Předmět: Re: protonovy magnetometr, slaby signal<br />
></p>

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