logaritm cu PIC

[Blacksmith]

Are cineva o rutina care sa calculeze logaritm dintr-un numar ?

Pt PIC-uri pe 8 biti.

 

Am gasit ceva pe net dar are o instructiune care nu exista, "SKPNC" :

http://www.dattalo.com/technical/software/pic/piclog.html

[sofian]

mplab stie skpnc sau skpz si sunt niste macrouri

vezi aici

http://www.mv.com/ipusers/ke1g/pic_std_macros.html

[Liviu M]

mplab stie...

MPASM stie, nu mplab.

[Blacksmith]

I-auzi ce bine !Dar tocmai acum cand gasii functie logaritmica mi-am dat seama ca nu mai imi trebuie. Acum imi trebuie ridicare la o putere fractionala si functia radical ... :(

[cirip]

Acum imi trebuie ridicare la o putere fractionala si functia radical

y=x^(a/b)y=x^a + x^(-b)log(y)=a*log(x)-b*log(x)=(a-b)*log(x)=c10^log(y)=10^cy=10^c, unde c=(a-b)*log(x)am notat x^a = "x la puterea a" . Sper ca nu am gresit pe undeva ...Dupa cum vezi, nu e chiar asa de rau ca stii sa calculezi logaritmul. Daca mai afli si cum se face ridicare la putere, poti sa calculezi rezultatul oricarei ridicari la o putere fractionara, inclusiv radical. Radicalul este un caz particular in care a=1 si b=2.Am folosit baza 10 in exemplu, dar se aplica la logaritmul in orice baza. Daca logaritmul pe care il ai este cel natural, atunci va trebui sa afli cum se calculeaza e^x.

[cirip]

y=x^a + x^(-b)

am facut-o de oaie si nu mai pot sa corectez.y=x^a + x^(-b) se va inlocui cuy=x^(a * b^(-1))

[cirip]

Pffff ... va rog sa ignorati cele doua posturi de mai sus. Inca o data, cu exemplu. y=x^(a/b)log(y)=(a/b)*log(x), notam c=(a/b)*log(x)10^(log(y))=10^((a/b)*log(x))=10^cy=10^cexemplu in excelul atasat.Asa e cand nu le faci pe hartie. Mintea mai joaca feste. Dar am impresia ca Dudi si-a pierdut interesul. [attachment=0]radical_log.xls[/attachment]

[10vid]

De obicei mai toate functiile (log, ln, e, sqrt, x^y, sin, cos, tan etc.....) se aproximeaza prin serii Taylor, care constau din cateva adunari si inmultiri.Cirip, foarte interesant si adevarat ceea ce spui, dar logaritmul ala nu vine moca, impartirea e groaznica si mai e si o inmultire. Toate astea pentru a calcula un biet radical?! Cand de fapt ar fi mai simplu sa calculezi direct radicalul. Sunt rutine pe net cam in jur de 100 ciclii. Nu e mult, avand in vedere ca o inmultire 16x16 biti e vreo 140 ciclii pe PIC16.

[cirip]

Intru totul de acord, dar Dudi vrea si ridicarea la putere fractionara, nu numai cazul particular al radicalului. Si m-am gandit ca daca tot are logaritmul ala acolo, de ce sa nu-l foloseasca?

Oricum, toata chestia e flexare la muschiul neuronului. Un pic de artimetica asa de intretinerea tonusului mental. Ca daca e vorba sa optimizam, depinde de context.

 

Cat despre aproximari, intr-adevar seria Taylor este destul de populara, dar are dezavantajul convergentei lente. Relativ recent am folosit pentru e^x aproximarea Padé, care este de tip rational-polinomiala, care converge mai repede si are precizie foarte buna pe un domeniu mai larg de valori ale exponentului. Detalii aici: http://en.wikipedia.org/wiki/Pad%C3%A9_approximant

[Blacksmith]

Am gasit niste Application Notes de la Microchip cu functiile astea, sa vedem ce e de capu lor :

 

- PDF-ul

 

- si softul !

[Blacksmith]

Le-am testat, merg foarte bine !