首页 > 其他分享 >2023-06-28《计算方法》- 陈丽娟 - 向量和矩阵基础.md

2023-06-28《计算方法》- 陈丽娟 - 向量和矩阵基础.md

时间:2023-06-28 17:22:34浏览次数:46  
标签:md 06 导数 不等式 28 矩阵 范数 任意 向量

2023-06-28《计算方法》- 陈丽娟 - 向量和矩阵基础

Matlab计算方法矩阵范数导数条件数

本问补充向量和矩阵范数的相关知识,为下一章节的线性方程组的迭代法以及误差分析做准备。除了参考《计算方法》一书,还参考了华东师范大学数学学院的课程材料《迭代方法与预处理》以及陈新宇、伍元凯、赵熙乐和孙立君所著《面向时空交通数据修复及预测的低秩机器学习模型》, 有兴趣者可以去看他们在Github发布的文章 https://xinychen.github.io/books/spatiotemporal_low_rank_models.pdf

1. 向量范数

定义1:
, 实值函数的范数,如果

  1. 非负性;
  2. 齐次性:即, ;
  3. 三角不等式:对任意两个向量, 都有

向量范数所表示即中的一个度量,表示向量的长度,下面给出向量范数的几个例子:

  1. 1-范数:;
  2. 2-范数:;
  3. -范数:;
  4. p-范数:, ;

需要注意向量0-范数(计数)不满足向量范数的定义。另外向量和内积有以下关系:

推论1:
上的内积,则上的一个向量范数。

既然在中的向量范数有无穷多个定义,那么我们该如何选择使用哪个范数呢?下面的定理将告诉我们虽然向量范数有无穷多个,但是任意两个向量范数之间是等价的。

定理1:
上定义的任一向量范数都与向量1-范数等价,即存在正数, 对一切, 下述不等式成立:

按照书中的方式给出证明:
, 由三角不等式可知是连续函数,则在单位球面上有界。设其最大值为, 最小值为,则有

另设为任意非零向量,则


(即变为单位向量),然后即得

也即是

虽然定理1只说明任意向量范数都与范数等价,但是很容易就能得出对于任意都是等价的,特别的:

更广义的,
定理2:
有限维赋范线性空间上的所有范数都是等价的。

Cauchy-Schwartz 不等式
上的内积,则对任意, 有


且等号成立的充要条件是线性相关。

Holder 不等式
上的Euclidean内积, 则对任意, 有


其中.

备注
关于p-范数是范数,Cauchy-Schwartz 不等式,Holder 不等式的证明可以在《矩阵理论》-黄廷祝 上找到。

2. 矩阵范数

矩阵范数的定义和向量范数如出一辙:
定义2:
, 实值函数的范数,如果

  1. 非负性;
  2. 齐次性:即, ;
  3. 三角不等式:对任意两个向量, 都有
  4. 相容性:, for all .

后续叙述均在上讨论。

下面给出几个常见的矩阵范数:

  1. -范数
  2. -范数 (算子范数,诱导范数,导出范数)
  3. 1-范数(列范数):;
  4. -范数(行范数):;
  5. 2-范数(谱范数):,其中的谱半径,的特征值

简单计算即可验证-范数不是算子范数(找不到这样的

定理
上的所有范数都是等价的,特别地


引理
, 则, 且

矩阵范数的性质(对复数域成立)
1. 对任意相容范数有(由相容性定义可得);
2. 任意算子范数都是相容范数;
3. -范数是相容范数;
4. 是酉不变范数,即对任意酉矩阵 ,有


5. ,
6. 若是正规矩阵, 则,因此, 其中是任意算子范数。(,即.)

相应的结论可以退化到实数域中。

3. 病态矩阵

所谓病态矩阵即在线性方程组的求解过程中的舍入误差对求解会产生严重影响的一类系数矩阵。为了方便判断矩阵是否是病态的,给出条件数的定义:
定义
阶非奇异矩阵,则称数为矩阵的条件数,记为, 其具有性质
1. ;
2. , ;
3. 若,则.

由于条件数的计算涉及矩阵求逆,因此实际上难以实现,如何对条件数进行合理估计在《矩阵理论》一书中有所涉及。下面按照《计算方法》给出一些经验判据:
1. 用列主元高斯消元法时,出现绝对值很小的主元素;
2. 系数矩阵行列式的绝对值很小或很大;
3. 系数矩阵元素的量级差别很大;
4. 系数矩阵的行列式几乎线性相关。

定理
, 为非奇异矩阵,为方程组的精确解,为求得的近似解,其剩余向量为, 则有误差估计


证明: 由, 得到, 即, 则有



因此

4. 向量和矩阵(含范数)导数

向量函数的导数即是多元函数的导数,其中相关的梯度和黑塞矩阵相关知识相同。

写出的偏导数.

由于, 故


从而, .

写出的偏导数.


写出函数的导数。


其中
.

标签:md,06,导数,不等式,28,矩阵,范数,任意,向量
From: https://www.cnblogs.com/NEFPHYS/p/17511979.html

相关文章

  • 2023-06-28 小程序、h5、App各端的条件编译
    //表示代码仅在H5平台上面执行,其他平台不执行<!--#ifdefH5-->需条件编译的代码<!--#endif-->//表示代码在H5平台上面不执行,其他平台上面执行<!--#ifndefH5-->需条件编译的代码<!--#endif-->//表示代码在H5平台和App......
  • 20230628
    最近忙了三四个的准备去参展的项目被取消了,领导没有正面通知,听负责采购的同事说的,采购关于展会用的物料被拒批,忙了几个月,突然一下松了下来来,这几天基本上都是摸摸鱼,工作状态一下一下子,松懈了下来。今天看到一个博客园的哥们的在博客园上分享了近十年的经历,感觉挺有意思的,我关注了......
  • 20230628水题选做
    约束条件题意给定一些关系\(x=y或x\neqy\)。求是否能满足。分析显然并查集。我们考虑将约束条件排序,先使形如\(x_{i}=x_{j}\)的\(x_{i}\)和\(x_{j}\)合并。而后我们观察是否存在\(x_{i}\)和\(x_{j}\)已经合并但是关系是\(\neq\)。代码#include<bits/stdc++.h>usingname......
  • 错误 C2872 “IServiceProvider”: 不明确的符号
      <imgsrc="./%E6%A0%B8%E5%BC%B9/typora-user-images/image-20230628154332422.pngstyle="zoom:200%;"/><imgsrc="https://images.cnblogs.com/cnblogs_com/blogs/796407/galleries/2320826/o_230628075017_image-20230628154332422.png&q......
  • md常用快捷键
    一:菜单栏文件:alt+F编辑:alt+E段落:alt+P格式:alt+O视图:alt+V主题:alt+T帮助:alt+H二:文件新建:Ctrl+N新建窗口:Ctrl+Shift+N打开:Ctrl+O快速打开:Ctrl+P保存:Ctrl+S另存为:Ctrl+Shift+S偏好:Ctrl+,关闭:Ctrl+W三:编辑撤销:Ctrl+Z重做:Ctrl+Y剪切:Ctrl+X复制:Ctrl+C粘贴:Ctr......
  • curl -I 192.168.74.128 拒绝连接
     查看直接的linuxweb服务器信息 curlwww.baidu.com都是正常的,唯独curl本机的ip地址就一直报错 由于使用的yum安装的nginx,nginx的一些模块依赖一些lib库,所以在安装nginx,也必须先安装这些lib库,这些依赖库主要有g++、gcc、openssl-devel、pcre-devel和zlib-devel所以执......
  • 【寒假每日一题】AcWing 4728. 乘方
    目录一、题目1、原题链接2、题目描述二、解题报告1、思路分析2、时间复杂度3、代码详解 一、题目1、原题链接4728.乘方-AcWing题库2、题目描述小文同学刚刚接触了信息学竞赛,有一天她遇到了这样一个题:给定正整数 a 和 b,求 a^b 的值是多少。a^b 即 b 个 a 相乘的值,例如......
  • 2023.6.28 - vue项目打包内存堆栈溢出JS stacktrace
    vue项目打包时报错,JSstacktrace:ReachedheaplimitAllocationfailed-JavaScriptheapoutofmemory这是因为node打包时是有内存空间限制的,node能分配多少空间,默认是根据电脑内存占比来算的。在内存比较小的电脑里,默认分配给node的内存可能不足以支撑起项目运行或者打包......
  • C/C++自助点餐系统[2023-06-28]
    C/C++自助点餐系统[2023-06-28]面向对象程序课程设计任务书【题目】自助点餐系统【目的】通过设计一个小型的自助点餐系统,训练综合运用所学知识处理实际问题的能力,强化面向对象的程序设计理念,使自己的程序设计与调试水平有一个明显的提高。【要求】1、每个学生必须独立完成;......
  • JS中字符串28种常用API总结,substring、slice、JSON.stringify、match、split、search
    一、引言在前端开发中,处理字符串是一项常见的任务。JavaScript提供了一系列的字符串API,用于操作和处理字符串数据。字符串常用的API方法有很多,包括查找字符串、截取字符串、替换字符串、分割字符串、大小写转换、字符串拼接和字符串比较等等。本文将介绍一些常用的字符串API......