NumPy叉乘

1、叉乘-数学解释

向量积,数学中又称外积、叉积,物理中称矢积、叉乘,是一种在向量空间中向量的二元运算。与点积不同,它的运算结果是一个向量而不是一个标量。并且两个向量的叉积与这两个向量和垂直。其应用也十分广泛,通常应用于物理学光学和计算机图形学中。

叉乘是三维设计和游戏开发的底层最简单逻辑,计算机从业者都应用理解叉乘。

2维空间中的叉乘是

看起来像个标量,事实上叉乘的结果是个向量,方向在z轴上。我们来做一个推导,我们可以来尝试计算一下两个向量围成的平行四边形的面积,其实可以求向量围成的三角形的面积。我们知道面积是底*高,所有推导如下:

设:


我找到向量A',使得A⟂A',



也就出现了



也就是

就有如下推导:


这个推导不是很容易直接看懂,我来解释一下,首先是A’是如何求得的,这非常简单,但是我们程序员可能忘了:

a,b是两个向量:a=(a1,a2) b=(b1,b2)

a垂直b:a1b1+a2b2=0

由此可以得到A’,然后是|A|=|A`|,因为是垂直关系,向量的膜肯定是相等的啊,然后就是


这是其实更加简单直接查询点积公式即:


同样的后面的利用点积代数定义:

很容易就得出结论,这个说的这样细致,就是因为我们程序员可能很久不碰数学了,这个推导对于高中生是非常简单的。好了,展示:

继续说叉乘:

叉乘不是面积,叉乘的膜才是面积,严格来说,叉乘是有向面积,正负代表了面积的方向。好像这样没什么意思啊,计算面积没有简化都少啊,来看三维向量。

3维空间中的叉乘是:


这里涉及到行列式的计算,后面在说明,我们来看一下结合意义


不难发现,在3D图像学中,叉乘的概念非常有用,可以通过两个向量的叉乘,生成第三个垂直于a,b的法向量,从而构建X、Y、Z坐标系,所以上面描述的有方向的面积现在应该明白了,因为二维向量的法向量其实就是Z轴,我们把二维向量补齐成三维向量计算,a=(a1,a2,0) b=(b1,b2,0) 发现结果为(0,0,a1b2-a2b1)的一个平行于Z轴的向量。

2、叉乘-程序应用

P1(-1,0,1)、P2(0,2,2)、P3(0,-1,2)三个点正好可以围成一个三角形,这三角形的面积是多少?

import numpy as np

# 构建点
P1 = np.array([-1, 0, 1])
P2 = np.array([0, 2, 2])
P3 = np.array([0, -1, 2])

# A和B两个向量尾部相连
A = P3 - P1
B = P3 - P1
# 计算叉乘
A_B = np.cross(A, B)
# 计算叉乘的膜
AB_mo = np.linalg.norm(A_B)
# 计算面积
Area = AB_mo / 2
print("三角形的面积为:", Area)

三角形的面积为: 2.1213203435596424


已知四个点p1(2,0,0),p2(1,0,2),p3(-1,0,0),p4(2,5,0)可以构成平行四面体,求它的体积?

这个题目我们利用小学的知识一样可以解决,但是如何利用程序的方式去解决?我们不可能让计算机自己画图,然后低乘高吧,而且如果四面体的所有面全部平行四边形(没有90度的角),我们即使绘图也是不容易求解的。所以,展示:



空间向量A,B,C组成平行6面体的三条邻边,H是高,我们知道向量的叉乘的模是面积,所以底面积就是|B×C|,关键看H如何表示,我们设i为H方向的单位向量,则:

这里参考点积公式



即可明白一个向量点积一个单位向量,结果为这个单位向量方向上的投影。继续展示:


得到:


接着展示:

import numpy as np

# 构建点
P1 = np.array([2, 0, 0])
P2 = np.array([1, 0, 2])
P3 = np.array([-1, 0, 0])
P4 = np.array([2, 5, 0])

# 构建向量,以P1为原点
A = P2 - P1
B = P4 - P1
C = P3 - P1

# 计算叉乘
A_B = np.cross(B, C)
# 计算点积
volume = np.inner(A, A_B)

print("四点围成的体积:", volume)

四点围成的体积: 30

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 201,552评论 5 474
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 84,666评论 2 377
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 148,519评论 0 334
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,180评论 1 272
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,205评论 5 363
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,344评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,781评论 3 393
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,449评论 0 256
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,635评论 1 295
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,467评论 2 317
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,515评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,217评论 3 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,775评论 3 303
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,851评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,084评论 1 258
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,637评论 2 348
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,204评论 2 341