<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 />