Reseni rovnic embedded

Ales Prochaska, Divesoft prochaska na divesoft.cz
Neděle Leden 14 13:40:08 CET 2018


A proč není možné použít přímo vzorce pro výpočet kořenů? Přiznávám,
že když se brala stabilita numerických algoritmů tak jsem asi chyběl,
ale čistě prakticky mi podobné věci vycházejí tak, že když úloha
vypadá jako špatně podmíněná, tak to znamená, že jsem zvolil špatné
rovnice, protože praktické úlohy obvykle vycházejí "dobře". Podle mě
by stačilo kontrolovat, jestli se v průběhu výpočtu nějaký dělitel
příliš nepřiblížil nule a v takovém případě hlásit, že jsou divná
vstupní data.

Aleš Procháska

> Ahoj,
> Tak něco konkrétního. Já bych si snížil řád polynomu na kubický pomocí
> derivace. Ta je v tomhle případě triviální.
> Pak bych pomocí Vietových vzorců spočítal  kořeny toho kubického plynomu
> a tím mám dost dobrý vstup do metody Regula falsi.
> Pokud do toho vhodně promítnete specifičnost požadavků (reálné kořeny,
> omezený rozsah, stanovená přesnost), mohlo by to být i docela svižné.

> JP



> Dne 13.01.2018 v 20:46 Jan Matuska napsal(a):
>> Dovolím si mírné OT, když je zmínka o STM32H7 na 400 MHz.....
>>
>> No ve škole jsem používal tohle:
>>
>> https://en.wikipedia.org/wiki/TI-89_series
>>
>> Dodnes mě uchvacuje, tak dobře to s tím HW fungovalo - vlastně funguje 
>> dodnes...
>>
>> Honza
>>
>>
>>
>> Dne 12.01.2018 v 21:12 Jaroslav Buchta napsal(a):
>>> Ono asi zalezi na tom, jaky MCU.  S treba STM32H7 na 400 MHz s FPU uz 
>>> se IMHO daji pouzit vsechny bezne postupy.
>>>
>>> Dne 12.01.2018 v 19:27 Ondrej napsal(a):
>>>> Řešení rovnic pomocí umělé inteligence - to mě zajímá. Kde začít 
>>>> hledat?
>>>>
>>>> Jinak víc než jsem napsal bohužel není k dispozici. Docela mě to 
>>>> překvapuje, že není rozumný "solver" pro ARM/DSP atd.
>>>>
>>>> O.
>>>>
>>>>
>>>> Dne 12.1.2018 v 1:00 Jiří Nesvacil napsal(a):
>>>>> Zdravim, nevidim spatny vysledek pokud Vam vyjde hodnota blizici se 
>>>>> nekonecnu. Pokud ma rovnice vice reseni nebo ji nezjednodusite, tak 
>>>>> to je horsi. Musite trosku vic popsat co chcete, jinak Vam asi 
>>>>> nikdo nepomuze. Jinak rovnice a derivace pro reseni extremu, ciste 
>>>>> numericky vypocet, optimalizace jako skokove metody nebo az 
>>>>> algoritmy oznacovane jako umele inteligence.
>>>>>
>>>>> Jirka
>>>>>
>>>>> Dne 11.01.2018 v 19:37 Ondrej napsal(a):
>>>>>> Teoreticky ano, jenomže když vezmeme obyčejnou kvadratickou 
>>>>>> rovnici (-b +- sqrt(D)) / 2a, tak pokud bude a relativně malé, tak 
>>>>>> bude řešení numericky mimo. A to je problém.
>>>>>>
>>>>>> Tohle je samo takový základní trik, pro větší stupně se to celé 
>>>>>> komplikuje a tam už se ztrácím a navíc nechci znovu objevovat kolo.
>>>>>>
>>>>>> Ondřej
>>>>>>
>>>>>>
>>>>>> Dne 10.1.2018 v 21:30 Slavomir Skopalik napsal(a):
>>>>>>> Zkuste zacit zde:
>>>>>>>
>>>>>>> https://en.wikipedia.org/wiki/Root-finding_algorithm
>>>>>>>
>>>>>>> treba Vam to pomuze.
>>>>>>>
>>>>>>> Slavek
>>>>>>>
>>>>>>> Ing. Slavomir Skopalik
>>>>>>> Executive Head
>>>>>>> Elekt Labs s.r.o.
>>>>>>> Collection and evaluation of data from machines and laboratories
>>>>>>> by means of system MASA (http://www.elektlabs.cz/m2demo)
>>>>>>> -----------------------------------------------------------------
>>>>>>> Address:
>>>>>>> Elekt Labs s.r.o.
>>>>>>> Chaloupky 158
>>>>>>> 783 72 Velky Tynec
>>>>>>> Czech Republic
>>>>>>> ---------------------------------------------------------------
>>>>>>> Mobile: +420 724 207 851
>>>>>>> icq:199 118 333
>>>>>>> skype:skopaliks
>>>>>>> e-mail:skopalik na elektlabs.cz
>>>>>>> http://www.elektlabs.cz
>>>>>>>
>>>>>>> On 10.1.2018 20:44, Ondrej wrote:
>>>>>>>> Zdravím konferu
>>>>>>>>
>>>>>>>>
>>>>>>>> potřeboval bych v rámci embedded systému řešit rovnici ax^4 + 
>>>>>>>> bx^3 + cx^2 + dx + e = 0 (nebo-li najít kořeny polynomu čtvrtého 
>>>>>>>> řádu). Teoreticky by to mělo jít přímo vzorečky, ale obávám se 
>>>>>>>> numerické nestability (oprávněně?) Použít něco jako 
>>>>>>>> Jenkins–Traub mi zase přijde trochu overkill a navíc pomalé.
>>>>>>>>
>>>>>>>> Rád bych využil toho, že úloha není zcela obecná, protože mě 
>>>>>>>> zajímají jen reálné kořeny, a-e jsou reálná čísla, navíc všechna 
>>>>>>>> čísla mají rozumný rozsah (řekněme +- 1e5 s rezervou) a kořeny 
>>>>>>>> mě zajímají jen v tomto rozsahu. Nepožaduji ani extra velkou 
>>>>>>>> přesnost.
>>>>>>>>
>>>>>>>>
>>>>>>>> Máte s tím někdo zkušenosti? Existuje nějaká knihovna, která jde 
>>>>>>>> zkompilovat pro embedded systém (ARM, gcc)? Případně doporučení 
>>>>>>>> čeho se vyvarovat.
>>>>>>>>
>>>>>>>>
>>>>>>>> Díky za rady
>>>>>>>>
>>>>>>>> Ondřej
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> _______________________________________________
>>>>>>>> HW-list mailing list  -  sponsored by www.HW.cz
>>>>>>>> Hw-list na list.hw.cz
>>>>>>>> http://list.hw.cz/mailman/listinfo/hw-list
>>>>>>>
>>>>>>>
>>>>>>> _______________________________________________
>>>>>>> HW-list mailing list  -  sponsored by www.HW.cz
>>>>>>> Hw-list na list.hw.cz
>>>>>>> http://list.hw.cz/mailman/listinfo/hw-list
>>>>>>
>>>>>> _______________________________________________
>>>>>> HW-list mailing list  -  sponsored by www.HW.cz
>>>>>> Hw-list na list.hw.cz
>>>>>> http://list.hw.cz/mailman/listinfo/hw-list
>>>>>
>>>>> _______________________________________________
>>>>> HW-list mailing list  -  sponsored by www.HW.cz
>>>>> Hw-list na list.hw.cz
>>>>> http://list.hw.cz/mailman/listinfo/hw-list
>>>>
>>>> _______________________________________________
>>>> HW-list mailing list  -  sponsored by www.HW.cz
>>>> Hw-list na list.hw.cz
>>>> http://list.hw.cz/mailman/listinfo/hw-list
>>>
>>>
>>> _______________________________________________
>>> HW-list mailing list  -  sponsored by www.HW.cz
>>> Hw-list na list.hw.cz
>>> http://list.hw.cz/mailman/listinfo/hw-list
>>
>> _______________________________________________
>> HW-list mailing list  -  sponsored by www.HW.cz
>> Hw-list na list.hw.cz
>> http://list.hw.cz/mailman/listinfo/hw-list

> _______________________________________________
> HW-list mailing list  -  sponsored by www.HW.cz
> Hw-list na list.hw.cz
> http://list.hw.cz/mailman/listinfo/hw-list




Další informace o konferenci Hw-list