Numpy求空间向量夹角余弦以及平面法向量

最新推荐文章于 2024-07-05 14:33:15 发布

2401_84009993

最新推荐文章于 2024-07-05 14:33:15 发布

阅读量437

点赞数 4

分类专栏：程序员文章标签： numpy 平面

本文链接：https://blog.csdn.net/2401_84009993/article/details/137473909

版权

这篇博客介绍了如何利用Numpy在3D空间中计算向量的夹角余弦以及平面的法向量。通过示例展示了如何创建3D图形并使用线性插值来求解特定点处的法向量值。

摘要由CSDN通过智能技术生成

xy_xz_dot=np.dot(xy_v,xz_v)#求点积
norm_xy=np.linalg.norm(xy_v)#求模
norm_xz=np.linalg.norm(xz_v)
cos_xy_xa=xy_xz_dot/(norm_xy*norm_xz)#cos(π/3)
angle_radians=np.arccos(cos_xy_xa)#π/3
angle_degree=np.degrees(angle_radians)#60度   
vector_P=np.cross(xy_v,xz_v)    #叉乘，计算法向量。1,1,1
print(vector_P)

calculate_p_vector()


由上可得余弦角以及法向量


接下来需要利用matplotlib的三角剖分库tri对一系列点进行三角剖分，


验证1：可以用自己定义的三角形网格进行3d绘图，ax.plot\_trisurf（）


验证2：可以利用tri库自带的mtri.LinearTriInterpolator()的\_\_call\_\_()方法对3d图形面上的点的z值进行计算


验证3：是否可以利用tri自带的calculate\_plane\_coefficients()函数对所有三角形面的平面方程系数进行计算

fig_4=plt.figure()
ax_4=fig_4.add_subplot(111,projection=‘3d’)
x=np.array([0,1,2])
y=np.array([0,1])
u,v=np.meshgrid(x,y)
u,v=u.flatten(),v.flatten()
t=np.array([[0,1,3],[1,3,4],[1,4,5],[1,2,5]])
z=np.array([1,0,