漫步线性代数一——引言
我们以线性代数的中心问题开启我们的航程:解决线性方程。最重要并且也是最简单的情况就是位置函数的数目等于方程的数目。现在我们有包含nn个未知变量的nn个方程,首先从n=2n=2开始:
\begin{equation} 1x+2y=3\\ 4x+5y=6\tag1 \end{equation}未知变量是 x,yx,y。我打算用消元和行列式两种方式来求解方程。当然 x,yx,y由 1,2,3,4,5,61,2,3,4,5,6这些数确定,问题就是如何利用这六个数求解出方程。
消元 第二个方程减去第一个方程的四倍。从而消去第二个方程中的xx,只留下了关于yy的方程:
(方程2)−4(方程1)−3y=−6(2)\begin{equation}(\text{方程2})-4(\text{方程1})\qquad -3y=-6\tag2\end{equation}这样就得到y=2y=2。然后xx从第一个方程1x+2y=31x+2y=3计算出来:
1x+2(2)=3x=−1(3)\begin{equation}1x+2(2)=3\qquad x=-1\tag3\end{equation}计算出来后,x,yx,y也应该满足第二个方程。代入得:4×(−1)+5×2=64\times (-1)+5\times 2=6
行列式 y=2y=2完全依赖于方程中的六个数。对于yy存在一个公式(当然xx也有),是两个行列式的比值:
y=∣∣∣1436∣∣∣∣∣∣1425∣∣∣=1⋅6−3⋅41⋅5−2⋅4=−6−3=2(4)\begin{equation}y=\frac{\begin{vmatrix}1&3\\4&6\end{vmatrix}}{\begin{vmatrix}1&2\\4&5\end{vmatrix}}=\frac{1\cdot 6-3\cdot 4}{1\cdot 5-2\cdot 4}=\frac{-6}{-3}=2\tag4\end{equation}如果你知道2×22\times 2方阵的行列式,那么它就没有那么神秘了。它同样给出了解y=2y=2。同样利用行列式,我们可以求出xx:
x=∣∣∣3625∣∣∣∣∣∣1425∣∣∣=3⋅5−2⋅61⋅5−2⋅4=3−3=−1(5)\begin{equation}x=\frac{\begin{vmatrix}3&2\\6&5\end{vmatrix}}{\begin{vmatrix}1&2\\4&5\end{vmatrix}}=\frac{3\cdot 5-2\cdot 6}{1\cdot 5-2\cdot 4}=\frac{3}{-3}=-1\tag5\end{equation}
我们来比较这两种方法,考虑nn非常大(在科学计算中n=1000n=1000是非常适中的大小)。事实就是直接对1000个方程使用行列式将是一个大灾难。左边将会有百万级别的数目,既然是正确的,但是效率很低。之后会提到该公式如何得到(克莱姆法则),但是目前我们需要一个很好的办法来解决这1000个方程。
最好的办法是高斯消元法。这个算法一直被用于解决大型的方程组。以后的大部分例子都是n=3n=3,这是看不出太大区别。方程(2)(4)(2)(4)基本使用相同的步骤得到y=2y=2。之后通过回带到方程(3)中很快得出xx。对于更大的nn值,依然有效。消去法比计算行列式要好。
消去法的想法看起来很简单,通过几个例子就能掌握它。它是非常基础的内容,通过简化矩阵我们就能理解它。在此我想讲四点更深层次的内容:
线性代数带来了平面几何。在十维空间中可视化九维平面不太容易。而理解相交于十个方程解的那些平面更难。但是不见得是不可能的。在图1 中有两条直线,相交于点(x,y)=(−1,2)(x,y)=(-1,2)。线性代数将图像放到十维空间里,在这个空间里,我们的直觉不得不去想象其几何形状。
图1:左边是例子的单个解,中间和右边是奇异情况,分别是无解和多个解现在考虑矩阵符号,将nn个未知量表示成向量xx,nn个方程表示成Ax=bAx=b。我们用消去矩阵乘以AA得到上三角矩阵UU。这些步骤将矩阵AA分解为L×UL\times U,其中LL是下三角矩阵:
A=∣∣∣1425∣∣∣=∣∣∣1401∣∣∣∣∣∣102−3∣∣∣=L×U(6)\begin{equation}A=\begin{vmatrix}1&2\\4&5\end{vmatrix}=\begin{vmatrix}1&0\\4&1\end{vmatrix}\begin{vmatrix}1&2\\0&-3\end{vmatrix}=L\times U\tag6\end{equation}首先我们需要介绍矩阵和向量以及乘法规则。每个矩阵都存在转置ATA^{\mathrm{T}}。这个矩阵还存在逆矩阵A−1A^{-1}
大部分情况下,消去法不会存在问题。矩阵可逆的话,方程Ax=bAx=b还有一个解。但是在特殊的情况下这种方法就被打破了,既可能是方程组的顺序出错(通过交换一笑就能产生),也可能是方程没有唯一解。如果将例子中的5换成8,就出现了奇异的情况:
1x+2y=34x+8y=6(7)\begin{equation}1x+2y=3\\4x+8y=6\tag7\end{equation}消去法依然用第二个减去第一个的四倍,那么结果是:
(方程2)−4(方程1)0=−6(\text{方程2})-4(\text{方程1})\qquad 0=-6
对于nn个方程组,我们希望粗略算出需要多少步消去运算。计算代价经常决定着模型的精度。一百个方程需要一百万步(乘法和减法)的三分之一。计算机可以很快地计算出来。在一百步之后,舍入误差就已经很明显了。(有些问题敏感,而有些不敏感)在不知道全部细节的情况下,我们想明白实际中出现的大型系统以及他们是如何被解决的。
之后我会尽可能有效的介绍消去算法。这个算法用于各种各样的应用中,同时,用矩阵的方式(系数矩阵AA,消去矩阵EE,行交换矩阵PP以及因子L,UL,U)理解它也是非常重要的。我希望接下来的一系列文章能够让大家感到轻松。
漫步线性代数一——引言相关推荐
- 漫步线性代数二——线性方程的几何形状
漫步线性代数二--线性方程的几何形状 2016年08月15日 23:10:10 会敲键盘的猩猩 阅读数:1818 几何形状 理解这个主题的方法是举例说明.我们以两个极其简单的方程开始,可以说大家在没有 ...
- 漫步线性代数二十一——行列式引言
在一百年前,行列式不是线性代数的中心,但是数学的方向一直在变换!毕竟,仅仅一个数就能告诉我们许多矩阵的信息. 对行列式的一种理解是:它对A−1,A−1bA^{-1},A^{-1}b的每一项给出了明确的 ...
- 漫步线性代数二十四——行列式应用
本篇文章介绍四个应用:AA的逆,求解Ax=bAx=b,盒子的体积以及主元.他们都是线性代数里面非常关键的计算,而行列式给出了这些答案的公式. 1.计算A−1A^{-1}.2×22\times 2矩阵展 ...
- 漫步线性代数十九——快速傅里叶变换(上)
傅里叶级数是无限维度上的线性代数,它里面的向量实际上是函数f(x)f(x):他们分别投影到正弦和余弦上:然后乘以傅里叶系数ak,bka_k,b_k.用ak,bka_k,b_k乘以无限的正弦和余弦序列就 ...
- 漫步线性代数十七——正交基和格拉姆-施密特正交化(上)
对于一个正交基,每个向量和其他所有向量垂直,坐标轴就是互相正交的.我们还可以进一步改善:每个向量除以它的长度得到单位向量,这样的话正交基变成了标准正交基: 16.如果 qTiqj={01i≠j,给出正 ...
- 漫步线性代数十一—— 四个基本子空间
上篇文章处理了定义而不是,我们知道基是什么,但不知道如何找到他们.现在,从一个明确描述的子空间开始,我们开始计算显式的基. 子空间能用两种方式描述.第一,我们可以给出一个生成空间的向量集合.(例如:列 ...
- 漫步线性代数十——线性无关,基和维数
m,nm,n没有给出线性系统实际大小的真实信息,在我们上文的例子中有三行和四列,但是第三行仅仅是前两行的组合,在消元后得到了零行,它对奇次问题Ax=0Ax=0 没有影响.第四列同样是相关的,列空间减到 ...
- 漫步线性代数七——特殊矩阵和应用
本篇文章有两个目标.第一是解释实际问题中大型线性方程组Ax=bAx=b的一种解法,事实是,工程或经济学中大型和现实的问题能够引导我们更深入理解这些知识.但是有一个很重要应用却不需要大量的准备工作. 另 ...
- 漫步线性代数五——三角分解和行交换
我们继续消去法,看看它对矩阵意味着什么.我们从Ax=bAx=b开始: Ax=⎡⎣⎢24−21−67102⎤⎦⎥⎡⎣⎢uvw⎤⎦⎥=⎡⎣⎢5−29⎤⎦⎥=b(1) \begin{equation} Ax ...
最新文章
- CVPR 2019 | 百度17篇论文入选,AI巨头都在关注什么?(附论文下载)
- 【华为云技术分享】探索软件复杂性简洁之道
- linux挂载磁盘没有权限,linux肿么知道哪个盘没有被挂载
- 华三配置telnet
- 最好免费的 HTML5 JS 网站视频播放器收集
- PourOver – 快速筛选和排序大的数据集合
- 局部变量AND全局变量
- 三星app图标不见了怎么恢复_手机便签的内容不见了怎么恢复?试试这款便签软件的时间轴功能...
- python 图形库介绍_数据分析Python手绘图形库有哪些
- MFC 通用对话框之颜色对话框
- matlab中测交流电压,间接测量交流电压有效值的方法与应用
- Adobe Photoshop CC已停止工作重装显卡驱动解决方法
- 《国际学术论文写作与发表》课后题
- ensp中ap获取不到ip_[网络求助]华为ap无法获取到ip
- 【kali技巧】查看宿主机windows的ip地址
- 央视CCTV摄制组进驻领存深入报道
- 配置MyBatis Plus 的乐观锁功能
- 关于python中面向对象(一)
- 如何将ios键盘返回由英文设置成中文
- c 语言形式参数,高级语言C的形式参数.doc
热门文章
- 给大一师弟师妹的一些建议
- 页面的访问速度如何提高?(我所知道的,哪位高手有其它的高见,敬请提醒....)...
- 百度启用Baidu.co.jp域名,有利于其在日本推广
- 面试官系统精讲Java源码及大厂真题 - 42 常用的 Lambda 表达式使用场景解析和应用
- Jeecg-Boot 一款基于SpringBoot 快速开发平台
- 容器编排技术 -- Kubernetes kubectl create secret tls 命令详解
- Chrome谷歌插件开发-01
- Python第十二章-多进程和多线程01-多进程
- ‘packaging‘ with value ‘jar‘ is invalid. Aggregator projects require ‘pom‘ as packaging.
- js 弹出框 背景不滑动 方案