数值分析重点公式
- 格式：docx
- 大小：327.46 KB
- 文档页数：10
第一章 非线性方程和方程组的数值解法
1）二分法的基本原理，误差：~12kbax
2）迭代法收敛阶：1lim0ipiic，若1p则要求01c
3）单点迭代收敛定理：
定理一：若当,xab时，(),xab且'()1xl，,xab，则迭代格式收敛于唯一的根；
定理二：设()x满足：①,xab时，(),xab，
②121212,,, ()(),01xxabxxlxxl有
则对任意初值0,xab迭代收敛，且：
定理三：设()x在的邻域内具有连续的一阶导数，且'()1，则迭代格式具有局部收敛性；
定理四：假设()x在根的邻域内充分可导，则迭代格式1()iixx是P阶收敛的()()()0,1,,1,()0jPjP（Taylor展开证明）
4）Newton迭代法：1'()()iiiifxxxfx，平方收敛
5）Newton迭代法收敛定理：
设()fx在有根区间,ab上有二阶导数，且满足：
①：()()0fafb；
②：'()0,,fxxab；
③：'',,fxab不变号
④：初值0,xab使得''()()0fxfx；
则Newton迭代法收敛于根。 6）多点迭代法：1111111()()()()()()()()()iiiiiiiiiiiiiiifxfxfxxxxxfxfxfxfxfxfxxx
收敛阶：152P
7）Newton迭代法求重根（收敛仍为线性收敛），对Newton法进行修改
①：已知根的重数r，1'()()iiiifxxxrfx（平方收敛）
②：未知根的重数：1''()(),()()()iiiiuxfxxxuxuxfx，为()fx的重根，则为()ux的单根。
8）迭代加速收敛方法：
2211211212()()iiiiiiiiiiixxxxxxxxxxx当不动点迭代函数()x在的某个邻域内具有二阶导数，'()1,0L平方收敛
9）确定根的重数：当Newton迭代法收敛较慢时，表明方程有重根
10）拟Newton法
其中11112222'1212()iiiniiiininnniiinfffxxxfffxxxAFxfffxxx
11）秩1拟Newton法：
Broyden秩1方法
第二章 线性代数方程组数值解法
1）向量范数：
①：非负性：0x，且0x的充要条件是0x；
②：齐次性：xx
③：三角不等式：xyxy
1范数：11niixx 2范数：12221()niixx
范数：1maxiinxx
p范数：11()nppipixx
2）矩阵范数：
①：非负性：0A，且0A的充要条件是0A；
②：齐次性：AA
③：三角不等式：ABAB
④：乘法不等式：ABAB
F范数：12211nnijFijAa
1范数：111maxnijjniAa，列和最大
范数：111maxnijinjAa，行和最大
2范数：2()HAAA，其中1()maxHiinAA，i为HAA的特征值，()AA
3）Gauss消元法（上三角阵）：313Mn；
Gauss-Jordan消元法（对角阵）：312Mn；
列选主元消元法：在消元之前进行行变换，将该列最大元素换置对角线主元位置；（可用于求逆矩阵）
全选主元消元法：全矩阵搜索矩阵最大元素进行行变换和列变换至其处于对角线主元位置；
4）三角分解法：
①：Doolittle分解法：A=LU，L单位下三角阵，U上三角阵
②：Crout分解法：A=LU，L下三角阵，U单位上三角阵
③：Cholesky分解法：A对称正定，TALL，L为单位下三角阵
④：改进的Cholesky分解法：A对称正定，TALDL，L为单位下三角阵，D为对角阵
⑤：追赶法：Crout分解法解三对角方程 5）矩阵的条件数1()1condAAA，谱条件数：1222()condAAA
6）如果1B，则IB为非奇异阵，且11()1IBB
7）迭代法基本原理：
①：迭代法：1iixBxK
②：()1B(lim0iiB，迭代格式收敛)
③：至少存在一种矩阵的从属范数，使1B
8）Jacobi迭代：ALDU
9）Gauss-Seidel迭代：111()()iixLDUxLDb
10）超松弛迭代法11iiixxr
11）二次函数的一维搜索：2111xxP
12）最速下降法：
选择方向0000()ZgradfxrbAx
进行一维搜索：1000xxr，其中00000(,)(,)rrArr
13）共轭梯度法：
第一步：最速下降法，00Pr，11rbAx，01(,)0rr
第二步：过1x选择0P的共轭方向110PrP，其中1000(,)(,)rAPPAP，过1x以1P为方向的共轭直线为11xxtP，进行二次函数的一维搜索211111111(,)(,)xxPrPAPP
14）一般的共轭梯度法：
第三章 插值法与数值逼近
1）Lagrange插值：0()()()nnjjjLxlxfx，
余项：(1)1()()()(1)!nnfExPxn 2）Newton插值：差商表
00100101010()()[ ]()[ ]()()[ ]()()nnnnfxfxfxxxxfxxxxxxxfxxxxxxxx余项(1)0101()()[ ]()()()(1)!nnnnfExfxxxxxxxxPxn
3）反插值
4）Hermite插值（待定系数法）'210()[()()()()]nnjjjjjHxxfxxfx
其中2'''1,1()()(),2(),12(),()njjjjjjjjjkkjjkxaxblxalxbxlxlxxx
余项：(22)21()()()(22)!nnfExPxn
5）分段线性插值：1111()()()jjjjjjjjjxxxxLxfxfxxxxx
插值基函数：0110101011110,,(),(),0,nnnnnnnnxxxxxxxxxxlxlxxxxxxxxxxx
余项：分段余项2(2)22,max()8MhMfx
6）有理逼近：反差商表
有理逼近函数式：000111122()()()()()nnnxxfxvxxxvxxxvxvx
7）正交多项式的计算：
定理：在[,]ab上带权函数()x的正交多项式序列0()nx，若最高项系数唯一，它便是唯一的，且由以下的递推公式确定
其中(,)()bijijaxdx
定理3.8
8）连续函数的最佳平方逼近：在2{1,,,,}nSpanxxx上，法方程为nHad， 其中1121(1)12131(2)1(1)1(2)1(21)nnnHnnn，10(,)()kkkdffxdx
均方误差：22**21(,)(,)niiiffPffad
最大误差：*01maxxfP
9）离散函数的最佳平方逼近（曲线的最小二乘拟合）：
法方程0(,)(,)njkjkjaf
其中00(,)()()(,)()()mjkijikiimkiikiixxffxx
第四章 数值积分
1）代数精度的概念及应用：对r次多项式的精确成立，以及代入法求解系数。
2）Lagrange插值代入
Lagrange插值基函数011011()()()()()()()()jjnjjjjjjjnxxxxxxxxlxxxxxxxx
0()()nbjjajfxdxHfx，其中()bjjaHlxdx
误差：(1)1()()()(1)!nbnafEfPxdxn
定理：数值积分公式具至少有n次代数精度其是差值型的
3）等距节点的Newton-Cotes公式
将拉格朗日差值积分公式中的差值节点ixaih即可，其中bahn；
00,(1)()!()!njnnjiijhHtidtjnj，令jjHCba（Cotes系数）则：
N-C公式的数值稳定性：当jC同号时是稳定的，否则不稳定，0()njjbaC（其中0maxjjn） N-C公式至少具有n次代数精度，若n为偶数，则其代数精度可提高到n+1次；
余项：
当n为偶数时，(2)1()()()(2)!nbnafEfxPxdxn
当n为奇数时，(1)1()()()(1)!nbnafEfPxdxn
4）复化的N-C公式
复化的梯形公式：将积分区间n等分，然后在每个区间上应用梯形公式
复化的Simpson公式：将积分区间n等分，然后在每个区间上应用Simpson公式
5）Romberg积分法
()mTh逼近()If的阶为2(1)mh
6）求积节点为n+1的机械求积公式的代数精度<=2n+1；
7）Gauss求积公式
1()nPx在[a，b]上与所有次数<=n的多项式带权1正交上式为Gauss求积公式、
8）Gauss-Legendre求积公式
给出1()nPx公式：0()1Px、1()Pxx、22(31)2xP······21()(1)2!nnnnndPxxndx
给出区间[1,-1]上的求积公式，取()nPx的零点为求积节点
① 取1()Px零点为0
② 取2P零点为33
对于区间[a,b]上的Gauss求积公式，令,[,]22abbaxttab，()()()22abbafxftgt，则：
余项：2(1)121101()()(),()()()2(22)!nnnnbagEfPtdtPtttttn
第五章 乘幂法
1）基本定理：
定理一：若12,,,n为A的特征值，()Px为某一多项式，则矩阵()PA的特征值是12(),(),,()nPPP。特别地，kA的特征值是12,,kkkn。
