<div dir="ltr">Zdravím,<div><br></div><div>zajímavá otázka. Je pravda, že ten jmenovatel 2B z hlediska výpočtu nevypadá moc dobře, ale co, je to konstanta a po rozumné normalizaci to půjde. Zajímavý, praktičtější polynom pro výpočet je zde: <a href="http://www.analog.com/en/design-center/reference-designs/hardware-reference-design/circuits-from-the-lab/CN0381.html#rd-description">http://www.analog.com/en/design-center/reference-designs/hardware-reference-design/circuits-from-the-lab/CN0381.html#rd-description</a> , část  4-Wire RTD Configuration, vzorec 6. </div><div><br></div><div>Dnešní procesory by asi výpočet zvládly, ale mnohem rychlejší bude setříděná tabulka odporů (nebo přímo výstupu A/D převodníku) a teplot, vyhledání intervalu v tabulce a mezi tím lineární aproximace. </div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">S pozdravem<br>ing. Pavel Poucha<br>jednatel<br><a href="mailto:pavel.poucha@papouch.com" target="_blank">pavel.poucha@papouch.com</a><br>Tel. +420 777 232 485<br><br>Papouch s.r.o. - vývoj<br>Papouch store s.r.o. - obchod</div><div>Papouch production s.r.o. - výroba</div><div dir="ltr"><br></div><div dir="ltr">Máte-li chuť, navštivte naše stránky <a href="http://www.papouch.com/" target="_blank">http://www.papouch.com/</a></div></div></div></div></div></div></div>
<br><div class="gmail_quote">2018-02-11 15:22 GMT+01:00 balu <span dir="ltr"><<a href="mailto:balu@k-net.fr" target="_blank">balu@k-net.fr</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">presne tak... vsetky funkcie som tak, ako vsetci tu uz videl milion krat. Ale az dnes to prvy krat prakticky implementujem, takze ma zaujimaju detaily. Presne ako hovorite, zarazilo ma to delenie velmi malym cislom.<br>
<br>
Pri hladani nejakych relevantnych zdrojov som az teraz zistil, ze vsade su len tie zbytocne rovnice, ale chyba pouzitelne riesenie.<br>
<br>
Vzhladom na to, ze ide o najpouzivanejsi priemyselny standard som predpokladal, ze uzitocne riesenie bude znacne preflaknute :-) Pritom to nie je az tak zlozity problem, ktory vyzaduje tazko ziskane know-how.<span class="HOEnZb"><font color="#888888"><br>
<br>
b.</font></span><div class="HOEnZb"><div class="h5"><br>
<br>
<br>
<br>
On 11/02/2018 15:13, Miroslav Mraz wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
1. To řešení je správně.<br>
2. Nelámal bych si s tím hlavu. Pro numeriku bude to řešení stejně asi nepoužitelné - budete dělit skoro nulou zřejmě něco co bude také velmi malé protože vznikne odečtením dvou skoro stejných čísel.<br>
Takže vzít původní data a nafitovat inverzní funkci na polynom požadovaného stupně a zkontrolovat odchylky.<br>
Je zajímavé, že v datasheetu bývají tyto aproximace tak jak uvádíte, ale pro inverzní funkce většinou nikoli. Ačkoli je to prakticky o hodně potřebnější.<br>
<br>
Mrazík<br>
<br>
Dne 11.2.2018 v 14:31 balu napsal(a):<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
skusim menej kontroverznu temu, aj ked urcite aj na tomto poli urobili ruski matematici velke pokroky uz v 15-tom storoci :-)<br>
<br>
Implementujem vypocet teploty z RTD senzora (napriklad PT100) a zaujala ma jedna vec.<br>
<br>
Najjednoduchia verzia rovnice zavislosti odporu elementu od teploty je linearna:<br>
<br>
Rt = R0 + A R0 t<br>
<br>
ktora ma riesenie:<br>
<br>
t = (-R0 + Rt)/(A R0)<br>
<br>
kde R0 je odpor pri danej teplote (napr. 0 stupnov C), Rt je namerany odpor a A je linearny clen teplotneho koeficientu odporu.<br>
<br>
Rozsirena rovnica so zvysenou presnostou obsahuje aj kvadraticky clen<br>
<br>
Rt = R0*(1 + A*t + B*t^2)<br>
<br>
jej riesenie je<br>
<br>
t = (-A R0 - Sqrt[R0] Sqrt[A^2 R0 - 4 B R0 + 4 B Rt])/(2 B R0)<br>
(alebo v grafickej forme vo wikipedii <a href="https://wikimedia.org/api/rest_v1/media/math/render/svg/b9446f92a09e799c680dd99e23ff005f6f1adc4e" rel="noreferrer" target="_blank">https://wikimedia.org/api/rest<wbr>_v1/media/math/render/svg/b944<wbr>6f92a09e799c680dd99e23ff005f6f<wbr>1adc4e</a> <br>
)<br>
<br>
Idem si zlomit mozog, ale nedokazem pochopit, preco z riesenia vypadol clen linearnej zavislosti, vsetko z jeho okolia sa zmenilo na konstantu. Rovnica teraz funguje len na kvadratickom clene B*Rt.<br>
<br>
Velmi naivne by som ocakaval, ze rovnica t = f(Rt) bude mat viacej clenov, ktore budu postupne spresnovat jej vysledok. T.j. ze bude fungovat aj ked pouzijem len A a B=0. V tomto pripade, ale kompletnejsie riesenie pre B=0 exploduje.<br>
<br>
?? Preco ??<br>
<br>
b.<br>
<br>
</blockquote>
______________________________<wbr>_________________<br>
HW-list mailing list  -  sponsored by <a href="http://www.HW.cz" rel="noreferrer" target="_blank">www.HW.cz</a><br>
<a href="mailto:Hw-list@list.hw.cz" target="_blank">Hw-list@list.hw.cz</a><br>
<a href="http://list.hw.cz/mailman/listinfo/hw-list" rel="noreferrer" target="_blank">http://list.hw.cz/mailman/list<wbr>info/hw-list</a><br>
</blockquote>
______________________________<wbr>_________________<br>
HW-list mailing list  -  sponsored by <a href="http://www.HW.cz" rel="noreferrer" target="_blank">www.HW.cz</a><br>
<a href="mailto:Hw-list@list.hw.cz" target="_blank">Hw-list@list.hw.cz</a><br>
<a href="http://list.hw.cz/mailman/listinfo/hw-list" rel="noreferrer" target="_blank">http://list.hw.cz/mailman/list<wbr>info/hw-list</a><br>
</div></div></blockquote></div><br></div>