Ich möchte ein sin cos für Optimierungstabelle nachschlagen erstellen, indem ein Array-Index 0-UCHAR_MAX
verwenden, so dass 0 Radiant Index 0, pi/2
Radiant ist UCHAR_MAX/4
:Warum ist diese sin cos Tabelle falsch, wenn das Radiant groß ist?
sincos.h
#include <limits.h>
#include <math.h>
int sini[UCHAR_MAX];
int cosi[UCHAR_MAX];
#define MAGNIFICATION 256
#define SIN(i) sini[i]/MAGNIFICATION
#define COS(i) cosi[i]/MAGNIFICATION
void initTable(){
for(int i=0;i<UCHAR_MAX;i++){
sini[i]=sinf(i*2*M_PI/UCHAR_MAX)*MAGNIFICATION;
cosi[i]=cosf(i*2*M_PI/UCHAR_MAX)*MAGNIFICATION;
}
}
der Grund der Verwendung von UCHAR_MAX
als max ich die guten Verwendung von unsigned char Überlauf machen will den Radianten das ist von 0 bis 2*pi
nur variiert simuliert: zum Beispiel, wenn der Wert von Radianten 2*pi
ist, wird der Index des Array UCHAR_MAX
weil es überläuft, es automatisch ly wird 0 und es ist keine Mod erforderlich (wenn ich 0 bis 360 als Domäne verwende, muss ich jedes Mal index%360
berechnen). Dann teste ich es mit einigen Radiant Werte:
float rad[]={2.0f,4.0f,6.0f,8.0f,10.0f,-2.0f,-4.0f,-6.0f,-8.0f,-10.0f};
wie folgt aus:
#include "sincos.h"
#include <stdio.h>
int main(){
initTable();
unsigned char radToIndex;
float rad[]={2.0f,4.0f,6.0f,8.0f,10.0f,-2.0f,-4.0f,-6.0f,-8.0f,-10.0f};
int scalar=123;
printf("scalar=%d\n",scalar);
for(int i=0;i<sizeof(rad)/sizeof(float);i++){
radToIndex=rad[i]*UCHAR_MAX/2/M_PI;
printf("%d*sin(%f) : %f , %d\n",scalar,rad[i],scalar*sinf(rad[i]),scalar*SIN(radToIndex));
}
return 0;
}
ich den Tisch testen mit 123*sin(radian)
, fand gehen die Ergebnisse beginnt über den eigentlichen man, wenn die Größe von Radiant zunimmt (wenn radian beträgt 10 oder -10):
scalar=123
123*sin(2.000000) : 111.843582 , 111
123*sin(4.000000) : -93.086708 , -92
123*sin(6.000000) : -34.368107 , -35
123*sin(8.000000) : 121.691063 , 122
123*sin(10.000000) : -66.914597 , -61
123*sin(-2.000000) : -111.843582 , -112
123*sin(-4.000000) : 93.086708 , 90
123*sin(-6.000000) : 34.368107 , 38
123*sin(-8.000000) : -121.691063 , -122
123*sin(-10.000000) : 66.914597 , 59
und Prüfung mit einem anderen Datum:
float rad[]={0.01f,0.1f,1.0f,10.0f,100.0f,1000.0f,-0.01f,-0.1f,-1.0f,-10.0f,-100.0f,-1000.0f};
Ausgabe:
scalar=123
123*sin(0.010000) : 1.229980 , 0
123*sin(0.100000) : 12.279510 , 12
123*sin(1.000000) : 103.500931 , 102
123*sin(10.000000) : -66.914597 , -61
123*sin(100.000000) : -62.282974 , -97
123*sin(1000.000000) : 101.706184 , -25
123*sin(-0.010000) : -1.229980 , 0
123*sin(-0.100000) : -12.279510 , -8
123*sin(-1.000000) : -103.500931 , -100
123*sin(-10.000000) : 66.914597 , 59
123*sin(-100.000000) : 62.282974 , 98
123*sin(-1000.000000) : -101.706184 , 22
Der Fehler erhöhen, wenn Größe zunimmt, so bin ich ganz sicher, dass die Tabelle ungenau wird, wenn Radiant groß ist. In sincos.h gibt es einen Wert MAGNIFICATION die Genauigkeit zu steuern, ich habe es 256 bis 4096 geändert, aber es scheint nicht viel Verbesserung:
scalar=123
123*sin(0.010000) : 1.229980 , 0
123*sin(0.100000) : 12.279510 , 12
123*sin(1.000000) : 103.500931 , 102
123*sin(10.000000) : -66.914597 , -62
123*sin(100.000000) : -62.282974 , -97
123*sin(1000.000000) : 101.706184 , -25
123*sin(-0.010000) : -1.229980 , 0
123*sin(-0.100000) : -12.279510 , -9
123*sin(-1.000000) : -103.500931 , -100
123*sin(-10.000000) : 66.914597 , 59
123*sin(-100.000000) : 62.282974 , 99
123*sin(-1000.000000) : -101.706184 , 22
warum das passieren? Gibt es einen logischen Fehler der Tabelle?
Wie Sie die Genauigkeit messure? – Zich