· Text 2261 · w8 · inlagg
Finns rutin för detta i banken?
M.v.h. Micke
Retro Preservation Workspace
Monitorn
Browse archived Monitorn möten and the ABC80 mailing list, with full-text and text-number search.
SQR med ASCII-aritmetik
11 posts in thread
-
-
· Text 2265 · w8 · kommentarMan kan räkna ut SQR ur tal så här: NYARE TIPS=(X/TIPS+TIPS)/2 X är talet att dra roten ur, X är ett tips på vad SQR(X) är (det går bra med ett första tips på 2). Upprepa beräkningen men använd då NYARE TIPS istället för det gamla. Upprepa igen osv ända till tipset inte ändrar sig mer, då har du fått fram SQR(X). Jag tror att man kan få ännu bättre fart på det hela om man som första tips väljer X/2. Finns det möjligen någon som vet ännu bättre första tips?
-
· Text 2266 · w8 · kommentarTar inte det TID i ASCII? Fast det kanske inte finns något bättre sätt... M.v.h. Micke
-
· Text 2273 · w8 · kommentarDen metoden är snabb, d.v.s. den antar snabbt ett närmevärde till roten du söker. Dessutom ger den ett fel som du kan bortse ifrån. Det är nämligen så att när NYARE TIPS = TIPS har du roten med maximalt antal RÄTTA siffror.
-
· Text 2277 · w8 · kommentarOm det nu handlar om räkning med ADD¤, DIV¤ osv, så föreslår jag att man i de fall talområdet räcker till, beräknar det första tipset med rotutdragning med vanliga flyttal!
-
· Text 2278 · w8 · kommentarSå här: A¤=det som ska tas roten ur T¤=första tips T¤=num¤(sqr(val(a¤))? Det förutsätter tal<100.100, eller har jag gjort fel? M.v.h. Micke
-
· Text 2279 · w8 · kommentarDet var så jag tänkte, ja, fast gränsen för A¤ blir väl att 10-exponenten får hålla sig mellan -127 och +128, eller vad det nu är på ABC80. (På ABC800: +/-38). Men det bör väl inte vara något besvär, om den ger error i fall talet inte går in, då kan man ju starta med att dela med två i DIV¤. .... Eller så kan man i det läget halvera exponenten. Då har man ju nästan redan dragit roten ur talet i alla fall.
-
· Text 2280 · w8 · kommentarASCII-aritmetiken i ABC80 kan inte hantera exponenter. För att NUM¤(SQR(VAL(A¤))) inte ska ge exponentformat måste SQR(VAL(A¤)) ligga i intervallet 1E-6 ... 1E+5 eller 0.000001 ... 99999.9 Det kommer det att göra om A¤ (som man drar roten ur) ligger i intervallet 1E-12 ... 1E+10 När man använder NUM¤(SQR( )) som första försök så verkar det som att redan en omgång i NYTEST=(A/TEST+TEST)/2 ger minst 10 decimalers nogranhet i NYTEST. Ett problem med ABC80 (i alla fall på checksumma 11273) är att vid division med DIV¤ så får man inte ha speciellt många decimaler. Om jag försöker ha mer än 10 eller 11 decimaler så får jag error.
-
· Text 2285 · w8 · kommentarGunnar Forsell visade mig en metod som hette något krångligt, men som enligt honom var den absolut bästa. Om intresse finns så kan jag lägga den i INLÅDA/EJKLARA. Ska jag det? (Den hette något med 'Newman...'...) M.v.h. Micke
-
· Text 2290 · w8 · kommentarNewton-Raphsons metod. Den absolut bästa metoden, ge inga fel med given noggrannhet, konvergerar sanbbt mot nollstället för funktionen ifråga. Dessutom är den lätt. Ny Tips = Tips - f(tips)/f'(tips) Så ser den ut. f'(x) är derivatan av f(x).
-
· Text 2294 · w8 · kommentarOm Newton-Raphsons metod tillämpas på SQR så ser det ut så här: Nytips = Tips - f(Tips)/f'(Tips) där f(x)=xÜ2-A och f'(x)=2*x A = tal till vilket kvadratroten sökes Nytips = Tips - (TipsÜ2-A)/(2*Tips) förenkla Nytips = Tips - Tips/2 + A/2/Tips Nytips = Tips/2 + A/2/Tips Nytips = (A/Tips + Tips)/2 Den sista algoritmen visar att metoden i inlägg 2265 utgår från Newton-Raphsons metod.
-
-
-
-
-
-
-
-
-
-