MIT18.065 数据分析、信号处理和机器学习中的矩阵方法 学习笔记 -- Lecture 11 Minimizing_xSubject to Ax=b

本文探讨了在二维和三维空间中不同范数(如L1、L2和Lp)的最小化问题,以及如何通过Gram-Schmidt方法和Krylov-Arnoldi算法处理非正交矩阵的正交化,特别关注了稀疏矩阵在大维度情况下的处理策略。
摘要由CSDN通过智能技术生成

 1.minimized different norms 寻找最小化的范数

二维平面中的最小范数

找出在直线 3x_{1}+4x_{2}=1最小的范数1,范数2,范数3。

min ||x||_{1},||x||_{2},||x||_{\infty } with 3x_{1}+4x_{2}=1

附:L1范数 L1=|x_{1}|+|x_{2}|

        L2范数L2=\sqrt{\sum_{}^{}|x_{i}|^{2}}

        L\infty范数LL\infty =max|x_{i}|

p=1 即L1范数 :假设x1与x2均为正,则用y=x1+x2;当x1为负时,存在y=-x1+x2;当x2为负时,存在y=x1+-x2,y=-x1+-x2等情况,因此L1范数图像呈钻石状。

p=2 即L2范数 :算术平方根形式,即该直线与原点的距离,因此L2范数图像呈现圆形。

p=\infty即L\infty范数:与坐标轴距离最大,L\infty图像呈现矩形。

根据题目要求寻找最小范数p,图像求解如下:

随着范数p的增大,最小范数点在逐渐向下移动,直到(\frac{1}{7},\frac{1}{7})

推广到多维

In 3D,直线变成平面或两个方程一条直线,L1扩展为立体钻石,L2扩展为球形,L3扩展为立方体。在三维方向上扩展。

2.Gram-Schmidt

 板块内容包括:
- Standard way A->QR

- Column exchanges

- "Krylov - Arnoldi"

(1)Srandard way A->QR

i. what Gram-Schmidt about?

即从A获得Q的过程

因此Gram-Schmidt,从一个非正交非独立的矩阵A,通过构造获得一个标准正交矩阵Q与其他矩阵的乘积R。

由于Q是标准正交矩阵,则Q^{-1}=Q^{T}(前面视频讲过),因此具体可计算出R的值,如下:

ii. review:what Gram-Schmidt about? 如何从A获得Q?

  A=\left [ a_{1} a_{2} ... a_{n}\right ]  为非正交矩阵且非独立。因此每一个向量ai直接也not atrhogonal.

需要通过Gram-Schmidt算法将每一个向量ai调整为单位正交。

计算方法如下:

Gram-Schmidt 两个步骤:1.减去原有向量上的分量,获得A;2.除以长度 unit vector q

(2)Column exchanges

i. What is the risk?标准Gram-Schmidt顺序求解的问题在哪?

如果a2和a1非常接近,因此在与a1正交的方向上,只存在非常小的向量。

通过exchange向量改变问题。

解决方法:在剩余的所有向量中,均减去a1方向上的分量,后进行比较,选择剩余向量较大的做第二个基向量q2。即在选择基向量时,选择最大的。

在标准方式中,当计算到column i 时,才会减掉之前i-1方向上的分量;

在修改后的方式中,在每一步中,将所有列都减去上一步已知的方向分量,后选择最大的向量。

just like picking the biggest number to go into the pivot

- Are we permuting every time? 答:每次都需要进行排序,选择最大

(3)"Krylov - Arnoldi"

Ax=b

A是一个很大的矩阵,并且为一个稀疏矩阵,无法进行转置。

So what we can do?

答:matrix times a vector.

b, Ab, A(Ab), ..., A^{j-1}b (空间的一组基向量)组成krylov space K_{j}X_{j}=best vector in K_{j}.

当前基向量非独立,非正交。

正交化基向量:Arnoldi && Lanczos

因此Arnoldi -- Lanczos算法是Krylov基向量通过改进后的Gram-Schmidt方法变成正交基向量。

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值