C ++ Угол между двумя векторами с acos - PullRequest
Новая
       33

C ++ Угол между двумя векторами с acos

0 голосов
/ 02 декабря 2018

Я пытаюсь заставить этот код работать, и он выдает углы, но не правильные, хотя они не достигают 180 °, и я не могу найти проблему.Я реализовал структуру Vector2, которая напоминает двумерный вектор, и я пытался определить для него оператор ^ как функцию, которая задает углы между двумя векторами.

Заранее спасибо за помощь! 

#include <iostream>
#include <cstdlib>
#include <cmath>

using namespace std;

struct Vector2 {
    double x, y;
};

ostream& operator<<(ostream& out, const Vector2& w){
    cout  << '[' << w.x << " , " << w.y << ']';
    return out;
}

double operator*(const Vector2& a, const Vector2& b){
    double w = a.x*b.x + a.y*b.y;
    return w;
}

double absc(const Vector2& d) {
    return sqrt(d*d);  
}

constexpr Vector2 operator+(const Vector2& a, const Vector2& b){
    Vector2 y{a.x+b.x, a.y+a.y};
    return y;
}

constexpr Vector2 operator-(const Vector2& a,const Vector2& b){
    Vector2 z{a.x-b.x, a.y-b.y};
    return z;
}

double operator^(const Vector2& a,const Vector2& b){
    double n = (absc(a))*(absc(b));
    double c = acos((a*b)/n);
    return c;
}

int main(){
    Vector2 a = {1.0, 1.0};
    Vector2 b = {4.0, 7.0};
    Vector2 c = {-2.0, 5.0};
    double d =c^b;
    double e = a^c;
    double f = a^b;
    auto grad = [](double rad) { return rad * (45./atan(1.0)); };

    cout<<"Angle Alpha: "<<grad(c^b)<<"rad "<<d<<"\n";
    cout<<"Angle Beta: "<<grad(a^c)<<"rad "<<e<<"\n";
    cout<<"Angle Gamma: "<<grad(a^b)<<"rad "<<f<<"\n";

    return EXIT_SUCCESS;
}

1 Ответ

0 голосов
/ 02 декабря 2018

Вы не сравниваете правильные углы, попробуйте это вместо этого (создайте векторы из точек, которые у вас есть): 

double d = (a-c)^(a-b);
double e = (b-a)^(b-c);
double f = (c-b)^(c-a);
auto grad = [](double rad) { return rad * (45./atan(1.0)); };

cout<<"Angle Alpha: "<<grad((a-c)^(a-b))<<"rad "<<d;
cout<<"\nAngle Beta: "<<grad((b-a)^(b-c))<<"rad "<<e;
cout<<"\nAngle Gamma: "<<grad((c-a)^(c-b))<<"rad "<<f;

И теперь у вас есть 180 ° или пи.

...