C ++程序使用行列式计算三角形的面积
在本节中,我们将看到如何使用矩阵行列式在2D坐标空间中找到三角形的面积。在这种情况下,我们考虑的空间是2D。因此,我们将每个点放在矩阵中。将x值放在第一列,y放在第二列,并以1作为第三列。然后找到它们的决定因素。三角形的面积将是行列式值的一半。如果行列式为负,则只需取其绝对值即可。
$$Area\:=\:absolute\:of\begin{pmatrix}\frac{1}{2}\begin{vmatrix}x_1\:\:y_1\:\\:1\\x_2\:\:y_2\:\:1\\x_3\:\:y_3\:\:1\end{vmatrix}\end{pmatrix}$$
在这里我们假设这是3x3矩阵,所以行列式函数无法找到不是3x3的矩阵的行列式。
范例程式码
#include<iostream>
#include<cmath>
using namespace std;
double det(double M[3][3]) {
double t1 = (M[1][1] * M[2][2])-(M[1][2] * M[2][1]);
double t2 = (M[1][0] * M[2][2])-(M[1][2] * M[2][0]);
double t3 = (M[1][0] * M[2][1])-(M[1][1] * M[2][0]);
return (M[0][0]*t1) + (-M[0][1]*t2) + (M[0][2]*t3);
}
main() {
double M[3][3];
cout << "Enter Point p1 (x, y):";
cin >> M[0][0] >> M[0][1];
M[0][2] = 1;
cout << "Enter Point p2 (x, y):";
cin >> M[1][0] >> M[1][1];
M[1][2] = 1;
cout << "Enter Point p3 (x, y):";
cin >> M[2][0] >> M[2][1];
M[2][2] = 1;
int determinant = det(M);
cout << "The area is: " << fabs(determinant) * 0.5;
}输出结果
Enter Point p1 (x, y):3 4 Enter Point p2 (x, y):6 4 Enter Point p3 (x, y):3 9 The area is: 7.5