通过二分法求任意sin正弦值
首先在一个水平坐标图的第一象限内,画一1/4圆,其半径为1,圆心为(0,0),该圆与x轴y轴相交点分别为R和L,也就是所求J角的范围。
通过不断二分缩小LR之间的差距,缩小查找的误差,在误差小到一定范围就能得到非常近似的正弦值
#include <iostream>
#include <cmath>
using namespace std;
struct node{//LR点所包含的信息不仅有两点的坐标,还有该点与x轴的夹角jiao
double x,y;
double jiao;
};
node l,r,mid;
int main()
{
l.x = 0;//L点夹角90
l.y = 1;
l.jiao = 90;
r.x = 1;//R点夹角0
r.y = 0;
r.jiao = 0;
double j,k;
cin >> j;//所求正弦的角
while(l.jiao-r.jiao>1e-7)//LR两夹角之间的角度没达到10的负7次方,
{//即精度不足小数点后7位就进行二分
mid.x = (l.x+r.x)/2;//取两点中点
mid.y = (l.y+r.y)/2;
k = 1/sqrt(mid.x*mid.x+mid.y*mid.y);
mid.x *=k;//中点与圆点的射线在圆上的交点坐标(mid.x,mid.y)
mid.y *=k;
mid.jiao = (l.jiao+r.jiao)/2;//该射线夹角为(L+R)/2
if(j<mid.jiao) l = mid;//若mid角大于j则在下半部分查找
else r = mid; //反之在上半部分查找
}
cout << mid.y/1;//该交点的纵坐标除以斜边长1就是正弦值
}