线性方程组的方法-数值分析-王兵团-北京交通大学


[Class]数值分析.王兵团.北京交通大学.全128讲[48:35:32]_哔哩哔哩_bilibili

注解:

1.线性代数中线性方程组的方法:克拉默法则。

线性方程组:

Ax=b

解:

xi=Di/D

如果A可逆,还可以写成:x=A-1/b

方程组的解是:系数行列式某一项换成等式右端常数项/系数行列式。

既然可以有这么好的公式,那为何还要学习其它解法呢?答:好多数学的公式一旦用到计算机里面,就不行了。

有人实验过,100万/s的计算量,解算40阶的线性方程组的解,要算一年。

天气预报的话有几百万几千万的方程组,怎样快速解出来?

大规模集成电路也需要解大规模方程组的。

人们需要快速求解大规模的线性方程组,这样,理论解就不行了。

数值分析讲的是怎样用计算机快速求出数学问题的解。

注解:

1.如果方程组的数量>未知数个数,没有解,或者有最小二乘解。

2.如果如果方程组的数量<未知数个数,没有解,有无穷多解。

3.如果如果方程组的数量=未知数个数,有唯一解。计算机做的最多的是这种情况:即方(阵)的情况。

4.线性方程组怎么得来的?答:每次实验得来的。

5.一个系统,有n个元器件,xi代表第i个元器件,每个元器件相当于一个变量xi,它们之间的变化有一定的关系。每实验一次,得到一个它们之间相互关系的方程。实验n次,得到n个方程。通过方程组,求出n个元器件情况。

注解:

1.计算机求的解都是近似解。

2.学一个东西,怎样学好?答:通过类比去看。

注解:

1.简单迭代法是怎样做的?

答:

数值分析6(2简单迭代法)苏州大学_哔哩哔哩_bilibili

注解:

1.初值可以给成:[0,0,0,...].

2.x是一个向量。

注解:

1.构造迭代格式所用的等价形式一定是有的。

2.未必都收敛的意思:比如,如果c给的合适,就收敛,如果不合适,就不收敛。

3.前文的引例就是example9.

 注解:

1.xk代表第k步的迭代值。xk是一个向量,所以ε(k)也是一个向量,是指第k步的迭代误差。

2.红色部分的等式是一个递推式子。

3.。。。完全取决于迭代矩阵B,跟初值怎样选择是没有关系的。那看一下什么样的矩阵迭代B,它的(k+1)次幂作用到一个非零向量上,当k趋近于无穷大时能让误差逐步收敛到0?

注解:

1.或者后面的等式:相当于范数收敛于0数值。

2.||B||,这个范数未必是2范,1范,2范,3范,4范。。。无穷范都行的。只要找到B的某范数小于1就行。

3.不同的范数,值不一样,有可能2范小于1,但无穷范数可能>1,也有可能1、2、3范数>1,但这并不能说明所有范数都是>1的,有可能是有<1的范数,没有找到。

4.范数有无穷多个,不可能逐一检验。

注解:

1.所有的范数,都是>=谱半径的。

2.这样就不用着1范,2范,3范,或者几范了,只要谱半径<1,迭代就会收敛。只要谱半径>=1,就不收敛。

 注解:

1.等价变形等式的两端都有未知向量,也就是都有未知参数。 

注解:

1.线性方程组迭代求解,做的比较早的是Jocobi,

 2.如果是方程两边只有一个变量,那么可以直接求出这个参数变量的值,但是线性方程组的变量有很多,比如有30个,那么可以:利用第一个方程组把第一个变量x1解出来,利用第二个方程把第二个参数变量x2表示出来,。。。

注解:

1.其实(3.5)式可以表示成矩阵的的形式,但是Jacobi时代矩阵的概念可能还没出来。

2.那我可以第一个方程解第二个变量,第二个方程我解第一个变量,这样可以吗?答:可以,但是没有规律。(3.5)式有唯一性。

3.那也可以第一个方程解第n个变量,第n个变量解第1个变量,这样也可以,但是这样的做法不是最简单的做法。

注解:

1.这个公式就是著名的Jacobi迭代公式。

2.使用Jacobi迭代公式计算线性方程组的解就叫做Jacobi迭代法。

3.在欣赏别人做的东西的同时,看还能不能改进了,不要迷信、神话专家。做不了大的突破,小的突破也可以,这种意识一定要有。----王兵团

4.Jacobi迭代过程中,看方程组里面的第一个方程,当使用第k次的参数向量的值带入第1个式子解出等式左边的x1的时候,x1的值是第k+1次的值了。此时,解第2个方程的时候,第2个方程右边的x1具备了新的值,那可不可以把更新后值带入呢?答:可以。应该带入更新后的值,按照这个思路,当利用第n个方程对第n个参数进行迭代计算的时候,方程右边所有的参数变量都更新了新的值了。这样的话,就得到了如下的Seidel公式了:

注解:

1.有人在Jacobi迭代公式的基础上做了一次小小的加工,做的工作很少,但是留名了。他就是Seidel.

2.数学公式构造都要有一个特点,要准确,要快速。

3.只看第n个方程,方程左边的参数(也是参数向量的最后一个参数)的值已经是第k+2次的值了。

4.对比Jacobi迭代来说,Seidel从第k次迭代到第k+1次后参数向量的第n个分量,是Jacobi迭代法从第k+1次迭代到k+2次后参数向量的第n个分量。

5.都收敛的情况下,Seidel迭代比Jacobi迭代能更快的收敛。

1是Jacobi迭代是收敛范围,3是Seidel迭代的收敛范围,2是他们共同的收敛范围。在区域2中,Seidel迭代比Jacobi迭代收敛的更快,也就是前者比后者有更快的收敛速度。

6.随着电子计算机的发展,现在有并行计算,可以很多计算机同时做一个任务。Jacobi迭代更适合于并行计算。在Jacobi迭代中,给出参数向量xk,可以并行算出参数向量xk+1。而Seidel迭代是一环套一环的,只能依次计算,这样的话,Jacobi迭代又比Seidel迭代快了。

7.有一些旧的知识点随着新技术尤其是计算机技术和算法的发展,可能又会焕发出青春。

注解:

1.把变量的差作为迭代的补偿,为了使得补偿变得更加丰富,加上一个松弛因子,通过放大或者缩小Δx,以调控新的值,最终目的是使得迭代更快的收敛。

注解:

1.不动点的理解:在方程的右边带入参数向量(x1,x2,x3),得到的点还是(x1,x2,x3),但是值是新值。

然后,写出线性方程组的Jacobi和Seidel迭代格式:

注解:

1.考试题目:写出一个线性方程组的Jacobi迭代格式,Seidel迭代格式。

注解:

1.这是Sor(超松弛)迭代格式。

下面的内容是用矩阵写线性方程组。那能否通过矩阵描述线性方程组的迭代计算呢?

 注解:

1.迭代矩阵B是一个常数矩阵。

注解:

1.一个数列收敛的定义是:

lim xk=x*(k趋近于无穷)<=> lim |xk-x*|(k趋近于无穷)=0

绝对值代表距离,两个数的距离趋近于0,说明它们靠的越来越近。

2.向量收敛怎样定义?

答:向量可以看成是空间的一个点。

向量的收敛还是从距离的角度进行定义。

3.范数就是距离,是空间的距离。

注解:

1.还有个p范数,可以包含这3中范数。

||x||p=  

2.以平面点(3, 4)为例,1范是7,2范是5,无穷范数是4. 1范是走直角后的距离,2范是走斜边的距离,无穷范数是先像y轴投影,再走较长边的距离。

3.研究这么多范数的目的,是为了研究距离,研究距离的目的是为了定义收敛。距离为趋近于0意味着无限靠近,这意味着收敛。

注解:

1.人们把矩阵看成是特殊的向量。