草庐IT

【数理知识】向量数乘,内积,外积,matlab代码实现

Zhao-Jichao 2024-07-20 原文
1【数理知识】向量数乘,内积,外积,matlab代码实现
2【数理知识】矩阵普通乘积,哈达玛积,克罗内克积,点乘,点积,叉乘,matlab代码实现

文章目录

1. 向量基本形式

形如 ( a 1 a 2 ⋮ a n ) \left(\begin{matrix} a_1 \\ a_2 \\ \vdots \\ a_n \\ \end{matrix}\right) a1a2an 的形式称之为向量。


2. 向量的数乘

指用一个数乘以向量中的每个元素

b ∗ ( a 1 a 2 ⋮ a n ) = ( a 1 a 2 ⋮ a n ) ∗ b = ( a 1 ∗ b a 2 ∗ b ⋮ a n ∗ b ) \begin{aligned} b * \left(\begin{matrix} a_1 \\ a_2 \\ \vdots \\ a_n \\ \end{matrix}\right) &= \left(\begin{matrix} a_1 \\ a_2 \\ \vdots \\ a_n \\ \end{matrix}\right) * b &= \left(\begin{matrix} a_1 * b \\ a_2 * b \\ \vdots \\ a_n * b \\ \end{matrix}\right) \end{aligned} b a1a2an = a1a2an b= a1ba2banb


3. 向量的内积

等于对应位置相乘再相加,两个向量的内积的结果是变成一个标量(也叫点乘)

a ⋅ b = ( a 1 a 2 ⋮ a n ) ⋅ ( b 1 b 2 ⋮ b n ) = ( a 1 a 2 ⋯ a n ) ( b 1 b 2 ⋮ b n ) = a 1 b 1 + a 2 b 2 + ⋯ + a n b n \begin{aligned} a \cdot b = \left(\begin{matrix} a_1 \\ a_2 \\ \vdots \\ a_n \\ \end{matrix}\right) \cdot \left(\begin{matrix} b_1 \\ b_2 \\ \vdots \\ b_n \\ \end{matrix}\right) =\left(\begin{matrix} a_1 & a_2 & \cdots & a_n & \end{matrix}\right) \left(\begin{matrix} b_1 \\ b_2 \\ \vdots \\ b_n \\ \end{matrix}\right) &= a_1 b_1 + a_2 b_2 + \cdots + a_n b_n \end{aligned} ab= a1a2an b1b2bn =(a1a2an) b1b2bn =a1b1+a2b2++anbn

a = [1;
	 2;
     3];
b = [4;
	 5;
     6];

Matlab语法:dot(a, b)

>> dot(a, b)
ans =
    32

>> a .* b
ans =
     4
    10
    18

>> a * b
错误使用  * 
用于矩阵乘法的维度不正确。请检查并确保第一个矩阵中的列数与第
二个矩阵中的行数匹配。要单独对矩阵的每个元素进行运算,请使用
TIMES (.*)执行按元素相乘。

点乘的几何意义是可以用来表征或计算两个向量之间的夹角,以及在 b b b 向量在 a a a 向量方向上的投影,有公式:

a ⋅ b = ∣ a ∣ ∣ b ∣ cos ⁡ θ a \cdot b = |a| |b| \cos \theta ab=a∣∣bcosθ


详细推导过程请参考:向量点乘(内积)和叉乘(外积、向量积)概念及几何意义解读,这里直接放结论

θ = arccos ⁡ ( a ⋅ b ∣ a ∣ ∣ b ∣ ) \theta = \arccos (\frac{a \cdot b}{|a| |b|}) θ=arccos(a∣∣bab)

根据这个公式就可以计算向量 a a a 和向量 b b b 之间的夹角。从而就可以进一步判断这两个向量是否是同一方向,是否正交(也就是垂直)等方向关系,具体对应关系为:

  • a ⋅ b > 0 a \cdot b > 0 ab>0 方向基本相同,夹角在0°到90°之间
  • a ⋅ b = 0 a \cdot b = 0 ab=0 正交,相互垂直
  • a ⋅ b < 0 a \cdot b < 0 ab<0 方向基本相反,夹角在90°到180°之间

4. 向量的外积

叉乘的运算结果是一个向量而不是一个标量。更为熟知的叫法是法向量,该向量垂直于 a a a b b b 向量构成的平面。两个向量的叉积与这两个向量组成的坐标平面垂直(也叫向量积、叉乘、叉积)

a × b = ( x 1 y 1 z 1 ) × ( x 2 y 2 z 2 ) = ∣ X Y Z x 1 y 1 z 1 x 2 y 2 z 2 ∣ = ( y 1 z 2 − y 2 z 1 ) X − ( x 1 z 2 − x 2 z 1 ) X + ( x 1 y 2 − x 2 y 1 ) Z \begin{aligned}a \times b= \left(\begin{matrix} x_1 \\ y_1 \\ z_1 \\ \end{matrix}\right) \times \left(\begin{matrix} x_2 \\ y_2 \\ z_2 \\ \end{matrix}\right) &= \left|\begin{matrix} X & Y & Z \\ x_1 & y_1 & z_1 \\ x_2 & y_2 & z_2 \\ \end{matrix}\right| = (y_1 z_2 - y_2 z_1)X - (x_1 z_2 - x_2 z_1)X + (x_1 y_2 - x_2 y_1)Z \end{aligned} a×b= x1y1z1 × x2y2z2 = Xx1x2Yy1y2Zz1z2 =(y1z2y2z1)X(x1z2x2z1)X+(x1y2x2y1)Z

a = [1;
	 2;
     3];
b = [4;
	 5;
     6];

Matlab语法:dot(a, b)

>> cross(a, b)
ans =
    -3
     6
    -3

>> dot(a, cross(a, b))
ans =
     0

>> dot(b, cross(a, b))
ans =
     0

关于向量外积的几何含义,我们假设在存在一个三维空间。

在下图可以看到,向量 V V V U U U 的外积是指在三维空间中,由 V , U V, U V,U 所组成平面的法向量。

下面一幅图阐述了向量外积的交换律。

在二维空间中,叉乘还有另外一个几何意义就是:aXb等于由向量a和向量b构成的平行四边形的面积。

Ref

  1. 矩阵和向量的乘法—点乘、叉乘、内积、外积、数乘、哈达玛积、克罗内克积
  2. 向量点乘(内积)和叉乘(外积、向量积)概念及几何意义解读
  3. 向量积 - 百度百科

有关【数理知识】向量数乘,内积,外积,matlab代码实现的更多相关文章

  1. ruby - 如何在 buildr 项目中使用 Ruby 代码? - 2

    如何在buildr项目中使用Ruby?我在很多不同的项目中使用过Ruby、JRuby、Java和Clojure。我目前正在使用我的标准Ruby开发一个模拟应用程序,我想尝试使用Clojure后端(我确实喜欢功能代码)以及JRubygui和测试套件。我还可以看到在未来的不同项目中使用Scala作为后端。我想我要为我的项目尝试一下buildr(http://buildr.apache.org/),但我注意到buildr似乎没有设置为在项目中使用JRuby代码本身!这看起来有点傻,因为该工具旨在统一通用的JVM语言并且是在ruby中构建的。除了将输出的jar包含在一个独特的、仅限ruby​​

  2. ruby-on-rails - Rails 源代码 : initialize hash in a weird way? - 2

    在rails源中:https://github.com/rails/rails/blob/master/activesupport/lib/active_support/lazy_load_hooks.rb可以看到以下内容@load_hooks=Hash.new{|h,k|h[k]=[]}在IRB中,它只是初始化一个空哈希。和做有什么区别@load_hooks=Hash.new 最佳答案 查看rubydocumentationforHashnew→new_hashclicktotogglesourcenew(obj)→new_has

  3. ruby - 如何根据特征实现 FactoryGirl 的条件行为 - 2

    我有一个用户工厂。我希望默认情况下确认用户。但是鉴于unconfirmed特征,我不希望它们被确认。虽然我有一个基于实现细节而不是抽象的工作实现,但我想知道如何正确地做到这一点。factory:userdoafter(:create)do|user,evaluator|#unwantedimplementationdetailshereunlessFactoryGirl.factories[:user].defined_traits.map(&:name).include?(:unconfirmed)user.confirm!endendtrait:unconfirmeddoenden

  4. ruby-on-rails - 浏览 Ruby 源代码 - 2

    我的主要目标是能够完全理解我正在使用的库/gem。我尝试在Github上从头到尾阅读源代码,但这真的很难。我认为更有趣、更温和的踏脚石就是在使用时阅读每个库/gem方法的源代码。例如,我想知道RubyonRails中的redirect_to方法是如何工作的:如何查找redirect_to方法的源代码?我知道在pry中我可以执行类似show-methodmethod的操作,但我如何才能对Rails框架中的方法执行此操作?您对我如何更好地理解Gem及其API有什么建议吗?仅仅阅读源代码似乎真的很难,尤其是对于框架。谢谢! 最佳答案 Ru

  5. ruby - 模块嵌套代码风格偏好 - 2

    我的假设是moduleAmoduleBendend和moduleA::Bend是一样的。我能够从thisblog找到解决方案,thisSOthread和andthisSOthread.为什么以及什么时候应该更喜欢紧凑语法A::B而不是另一个,因为它显然有一个缺点?我有一种直觉,它可能与性能有关,因为在更多命名空间中查找常量需要更多计算。但是我无法通过对普通类进行基准测试来验证这一点。 最佳答案 这两种写作方法经常被混淆。首先要说的是,据我所知,没有可衡量的性能差异。(在下面的书面示例中不断查找)最明显的区别,可能也是最著名的,是你的

  6. ruby - 寻找通过阅读代码确定编程语言的ruby gem? - 2

    几个月前,我读了一篇关于ruby​​gem的博客文章,它可以通过阅读代码本身来确定编程语言。对于我的生活,我不记得博客或gem的名称。谷歌搜索“ruby编程语言猜测”及其变体也无济于事。有人碰巧知道相关gem的名称吗? 最佳答案 是这个吗:http://github.com/chrislo/sourceclassifier/tree/master 关于ruby-寻找通过阅读代码确定编程语言的rubygem?,我们在StackOverflow上找到一个类似的问题:

  7. ruby - Net::HTTP 获取源代码和状态 - 2

    我目前正在使用以下方法获取页面的源代码:Net::HTTP.get(URI.parse(page.url))我还想获取HTTP状态,而无需发出第二个请求。有没有办法用另一种方法做到这一点?我一直在查看文档,但似乎找不到我要找的东西。 最佳答案 在我看来,除非您需要一些真正的低级访问或控制,否则最好使用Ruby的内置Open::URI模块:require'open-uri'io=open('http://www.example.org/')#=>#body=io.read[0,50]#=>"["200","OK"]io.base_ur

  8. 程序员如何提高代码能力? - 2

    前言作为一名程序员,自己的本质工作就是做程序开发,那么程序开发的时候最直接的体现就是代码,检验一个程序员技术水平的一个核心环节就是开发时候的代码能力。众所周知,程序开发的水平提升是一个循序渐进的过程,每一位程序员都是从“菜鸟”变成“大神”的,所以程序员在程序开发过程中的代码能力也是根据平时开发中的业务实践来积累和提升的。提高代码能力核心要素程序员要想提高自身代码能力,尤其是新晋程序员的代码能力有很大的提升空间的时候,需要针对性的去提高自己的代码能力。提高代码能力其实有几个比较关键的点,只要把握住这些方面,就能很好的、快速的提高自己的一部分代码能力。1、多去阅读开源项目,如有机会可以亲自参与开源

  9. Matlab imread()读到了什么 (浅显 当复习文档了) - 2

    matlab打开matlab,用最简单的imread方法读取一个图像clcclearimg_h=imread('hua.jpg');返回一个数组(矩阵),往往是a*b*cunit8类型解释一下这个三维数组的意思,行数、数和层数,unit8:指数据类型,无符号八位整形,可理解为0~2^8的数三个层数分别代表RGB三个通道图像rgb最常用的是24-位实现方法,即RGB每个通道有256色阶(2^8)。基于这样的24-位RGB模型的色彩空间可以表现256×256×256≈1670万色当imshow传入了一个二维数组,它将以灰度方式绘制;可以把图像拆分为rgb三层,可以以灰度的方式观察它figure(1

  10. 华为OD机试用Python实现 -【明明的随机数】 2023Q1A - 2

    华为OD机试题本篇题目:明明的随机数题目输入描述输出描述:示例1输入输出说明代码编写思路最近更新的博客华为od2023|什么是华为od,od薪资待遇,od机试题清单华为OD机试真题大全,用Python解华为机试题|机试宝典【华为OD机试】全流程解析+经验分享,题型分享,防作弊指南华为o

随机推荐