Re: ricerca radici di polinomio di grado qualsiasi (anche non intersecante, ma solo tangente, l'asse X)

From: Soviet_Mario <Soviet.Mario_at_CCCP.MIR>
Date: Wed, 19 Dec 2012 17:14:06 +0100

Il 17/12/2012 20:51, pippo ha scritto:
> T'interessano le radici reali o quelle complesse? O entrambe? Ti
> interessa la molteplicit�?

soltanto le radici reali, anzi positive, e ristrette in un
dominio ben preciso, come ho scritto a Tommaso Russo
nell'altra branca, dove la "genesi" dei polinomi �
formalizzata con precisione

>
> In buona sostanza il metodo delle secanti e delle tangenti sono dei
> Newton modificati, nei quali sostituisci la derivata della funzione con
> una sua approssimazione.

Cmq le funzioni di derivazione le ho gi� implementate,
quindi anche usare le forme esatte non mi cruccia.

> Trattandosi di polinomi, ovvero di funzioni 'lisce', piene di buone
> qualita', e curandoti troppo di velocit� di convergenza, anche se
> consideri gradi elevati, tutto sommato un metodo vale l'altro, visto che
> difficilmente sarai in presenza di oscillazioni tali da costituire un
> problema per la convergenza.

ah, senza i trucchi per isolare gli zeri che mi ha suggerito
T. Russo, sono incappato pi� volte in 1) overflow di memoria
con la bisezione, minimi con le tangenti, non convergenza o
convergenza lentissima per polinomi che erano solo tangenti
l'asse.
Cmq la situazione sta migliorando, con Horner, anche se non
ho ancora finito, il primo codice fa schifo, � pieno di
gestioni ad hoc per le derivate finali (rettilinee), quelle
prefinali (di II grado) e le precedenti, da computare a loro
volta come il polinomio originale
ciao
CCCP

>
> In ogni caso esistono metodi specializzati per determinare zeri di
> polinomi, ma sono praticamente esercizi di stile...
>
>
> Il 10/12/2012 18:01, Soviet_Mario ha scritto:>
> > Vorrei un parere comparativo, relativo ai due metodi, delle tangenti


--
1) Resistere, resistere, resistere.
2) Se tutti pagano le tasse, le tasse le pagano tutti
Soviet_Mario - (aka Gatto_Vizzato)
Received on Wed Dec 19 2012 - 17:14:06 CET