已知平面内三点,求其平面的法向量

三点平面法向量

设三点坐标为A(x1,y1,z1),B(x2,y2,z2),C(x3,y3,z3)

向量AB=(x2-x1,y2-y1,z2-z1),AC=(x3-x1,y3-y1,z3-z1)

AB、AC所在平面的法向量即AB×AC=(a,b,c),其中:

a=(y2-y1)(z3-z1)-(z2-z1)(y3-y1)

b=(z2-z1)(x3-x1)-(z3-z1)(x2-x1)

c=(x2-x1)(y3-y1)-(x3-x1)(y2-y1)

设 a=(ax,ay,az), b=(bx,by,bz)。i,j,k分别是X,Y,Z轴方向的单位向量,则:

a×b=(aybz-azby)i+(azbx-axbz)j+(axby-aybx)k

a·b=(axbx+ayby+az*bz)

Python代码

python 复制代码
def normal_vector(p1, p2, p3):
    x1, y1, z1 = p1
    x2, y2, z2 = p2
    x3, y3, z3 = p3

    a = (y2 - y1) * (z3 - z1) - (z2 - z1) * (y3 - y1)
    b = (z2 - z1) * (x3 - x1) - (x2 - x1) * (z3 - z1)
    c = (x2 - x1) * (y3 - y1) - (y2 - y1) * (x3 - x1)

    return (a, b, c)

if __name__ == '__main__':
    p1 = 1.0, 5.2, 0.0
    p2 = 2.8, 3.9, 1.0
    p3 = 7.6, 8.4, 2.0
    p4 = normal_vector(p1, p2, p3)
    print(p4)

    # output (-5.800000000000001, 3.0, 14.340000000000002)

C++代码

cpp 复制代码
#include<iostream>

using namespace std;

//三维double矢量
struct Vec3d
{
	double x, y, z;

	Vec3d()
	{
		x = 0.0;
		y = 0.0;
		z = 0.0;
	}
	Vec3d(double dx, double dy, double dz)
	{
		x = dx;
		y = dy;
		z = dz;
	}
	void Set(double dx, double dy, double dz)
	{
		x = dx;
		y = dy;
		z = dz;
	}
};

//计算三点成面的法向量
void Cal_Normal_3D(const Vec3d& v1, const Vec3d& v2, const Vec3d& v3, Vec3d &vn)
{
	//v1(n1,n2,n3);
	//平面方程: na * (x -- n1) + nb * (y -- n2) + nc * (z -- n3) = 0 ;
	double na = (v2.y - v1.y)*(v3.z - v1.z) - (v2.z - v1.z)*(v3.y - v1.y);
	double nb = (v2.z - v1.z)*(v3.x - v1.x) - (v2.x - v1.x)*(v3.z - v1.z);
	double nc = (v2.x - v1.x)*(v3.y - v1.y) - (v2.y - v1.y)*(v3.x - v1.x);

	//平面法向量
	vn.Set(na, nb, nc);
}

int main()
{	
	// Vec3d v1(1.0, 5.2, 3.7);
	// Vec3d v2(2.8, 3.9, 4.5);
	// Vec3d v3(7.6, 8.4, 6.2);   
	// 法向量为:-5.81 0.78 14.34;

	Vec3d v1(1.0, 5.2, 0.0);
	Vec3d v2(2.8, 3.9, 1.0);
	Vec3d v3(7.6, 8.4, 2.0);
	// 法向量为:-5.8 3 14.34
	
	Vec3d vn;
	Cal_Normal_3D(v1, v2, v3, vn);
	// cout <<"法向量为:"<< vn.x << '\t' << vn.y << '\t' << vn.z << '\n';
    cout <<"法向量为:"<< vn.x << " " << vn.y << " " << vn.z << '\n';

	return 0;
}
相关推荐
想跑步的小弱鸡4 小时前
Leetcode hot 100(day 3)
算法·leetcode·职场和发展
xyliiiiiL6 小时前
ZGC初步了解
java·jvm·算法
爱的叹息6 小时前
RedisTemplate 的 6 个可配置序列化器属性对比
算法·哈希算法
独好紫罗兰7 小时前
洛谷题单2-P5713 【深基3.例5】洛谷团队系统-python-流程图重构
开发语言·python·算法
每次的天空7 小时前
Android学习总结之算法篇四(字符串)
android·学习·算法
请来次降维打击!!!8 小时前
优选算法系列(5.位运算)
java·前端·c++·算法
qystca8 小时前
蓝桥云客 刷题统计
算法·模拟
别NULL8 小时前
机试题——统计最少媒体包发送源个数
c++·算法·媒体
weisian1519 小时前
Java常用工具算法-3--加密算法2--非对称加密算法(RSA常用,ECC,DSA)
java·开发语言·算法
程序员黄同学10 小时前
贪心算法,其优缺点是什么?
算法·贪心算法