函数插值
type
status
date
slug
summary
tags
category
icon
password
Property
 
已知平面上若干互不相同的点 ,其关系属于某未知对应关系,如何构造一函数,使得该函数可以穿过这些点并刻画
这种问题被称为插值问题,上述的函数被称为插值函数
 

Lagrange多项式法

Lagrange多项式法使用一个 次多项式来拟合上述的 个点,且可以得出唯一的结果
 
证明:
对于点集 ,找到一 次多项式 使得等价于方程组:
展开之,可使其变为一个系数矩阵是范德蒙矩阵的关于多项式系数向量的线性方程组,且范德蒙矩阵在互不相同时满秩,因此多项式系数向量存在且唯一。
 
函数逼近与曲线拟合
type
status
date
slug
summary
tags
category
icon
password
Property
 
 
函数逼近问题给出一目标函数 ,希望提出一(容易处理的)近似函数 来代替 ,以近似刻画并简化对 的诸多运算.
而插值问题是给出一系列点 ,希望提出一拟合函数 以近似刻画这些点所隐含的函数关系 .
函数逼近法的总体过程
  1. 提出近似函数 的形式并预留一些待定系数
  1. 使用某种误差估计手段写出误差关于上述待定系数 的关系
  1. 通过最小化误差 求出
由于其有着便于计算的性质,通常使用多项式函数 来进行拟合,同时使用最小二乘法来描述误差
在某些情况下,可能会将自变量 映射为其它的函数,如设置原型 逼近,则可以对逼近点列进行变换如
 

内积空间

数值微分
type
status
date
slug
summary
tags
category
icon
password
Property
 
 
数值微分问题求一函数 在其解析域内的某点 处的近似一阶导数 ,常用的方法是用一简单插值函数来代替 求导

插值法数值微分

在原函数 上找到 个点,据其构造一插值函数 ,以 来代替
误差分析:在点集上的 次插值余项可写作:
其中
对上式求导并将近似求导点 代入其中:
 

两点公式

数值积分
type
status
date
slug
summary
tags
category
icon
password
Property
 
 
数值积分求一函数在其连续域 内的积分近似值.

矩形近似

将不会求解的问题抽象为简单的问题——即使这种抽象是不精确的,因此想到积分中值定理:
可以在 内自由选择 的值以近似估计,这一式子实际上是用一矩形来近似代替积分,因此称为矩形公式。特殊地,若 时,上式被称为左矩形公式、右矩形公式和中矩形公式。
矩形近似本质上是用零次插值函数 来代替
 
 

插值近似

使用对 上的若干点进行次插值,得到插值函数 ,以 来代替 。易知 次插值公式的截断误差为
常微分方程的数值解
type
status
date
slug
summary
tags
category
icon
password
Property
 
 
对于一个附带了初值条件(否则将无法定义的截距)的一阶常微分方程:
求函数 在若干点 的函数值
假定方程右侧的多项式 是连续的,且关于 满足比连续更为高级的条件Lipschitz连续条件(收缩映射):
只要满足Lipschitz条件,存在单解
 

一阶近似法

既然只需要在离散的点集上进行求算,故上述方程便允许一定程度上的信息损失,而这种近似操作的首选目标便是微分或积分算子:

导数的差商近似

恰定方程组的解
type
status
date
slug
summary
tags
category
icon
password
Property
 
 
恰定方程组指系数矩阵秩与其阶数相同,且与增广矩阵秩相同的方程组,它有且仅有一个解向量。
它的一般形式是
我们需要通过某种手段解出符合上述条件的
 
 

Cramer法则

将上述列向量 替换 中的第 个列向量,把产生的新矩阵记作 ,那么解向量的第 项:
它要求系数矩阵是满秩的
矩阵的误差分析
type
status
date
slug
summary
tags
category
icon
password
Property
 
直观地,误差分析首先要求我们能够确定且统一地描述误差(即测定值与实际值之间的"距离"),而单靠矩阵显然做不到这一点,因此借用了范数这一工具.

度量方法:范数

范数实际上是实分析中的一条重要概念,是长度概念的推广
假设有一线性空间 ,存在某一函数 映射至非负实数空间 ,且这个函数满足:
  1. 当且仅当 (正定性)
  1. (绝对齐次性)
  1. (三角不等式)
则这一函数被称为 上的一个范数,二元体 被称为一个赋范线性空间或Banach空间.
一句话描述:范数是某种距离的刻画.
 
具有如下性质的范数是等价的:
超定线性方程组
type
status
date
slug
summary
tags
category
icon
password
Property
 

超定线性方程组的最小二乘解

看一个引例:对物理量 执行三次测量,结果分别是 ,则这个真实量最有可能是? 记函数 ,称其为方差函数。我们可以后验地认为,当方差函数取最小值时的 是最有可能取到的。这样的思想被称为“极大似然”,这种方法被称为“最小二乘”。
将上述的情况转化为矩阵形式:
上式看作是一个线性方程组,而若 不完全相同,则上述方程组的增广矩阵之秩比系数矩阵秩大,故这样的方程组属"无解"方程组,称为超定线性方程组。 但仍可以通过某一优化目标来将这个问题转化为一个最优化问题. 当这一优化目标是"均方误差最小化",或称为"2-范数最小化"时,这样的解法称为"超定线性方程组的最小二乘解法".
 

解的推导:几何理解

对于线性方程组 ,将 列向量记为 ,其中,考虑优化条件:
它的意义是: 的距离最小。其中可以由矩阵乘法法则展开为:
线性方程组迭代法与数值解
type
status
date
slug
summary
tags
category
icon
password
Property
 
向量和矩阵的极限
上的向量序列,若:
就称向量序列 收敛至
向量收敛的充分必要条件是向量列表中的每一个元素均收敛至,矩阵同理
 

迭代法的一般形式

对恰定线性方程组 ,构造与其同解的线性方程组
,代入上述方程:
矩阵特征值计算方式
type
status
date
slug
summary
tags
category
icon
password
Property
 

幂法

设置方阵 的各特征值 中仅有唯一的模长最大的值,且各特征矢 线性无关,即:
是任意不为零的 维向量:
由迭代公式:
可以证明,在 足够大时,矩阵特征值 是向量 与向量 的任意两个对应分量之商
 
证明
由于 维向量且线性无关,故它是 空间中的一组基底,故无论的取值,它总能被表示为 ,其中, 不全为0, 即是各个特征矢的线性组合
 
非线性方程组的数值解
type
status
date
slug
summary
tags
category
icon
password
Property
 

二分查找法

对于方程 ,其奇数重根 有如下的性质:存在 的某个邻域 ,使得邻域内分别在零点 左和右的两点 ,所以只需保证 ,每次将 替换为,直到误差达到要求即可。
二分查找法一定收敛,但收敛速度较慢。但二分查找法较为依赖初始值的选取,且无法处理偶数重根。
 
 

迭代法

不动点迭代

如果一函数 满足在 ,那么就可以尝试使用迭代格式
来找到这个 。但这个迭代格式在某区间 上不一定能够收敛,于是寻找使它收敛的条件:
首先,在这个区间 内方程 应有根,即满足:
集合约束和无约束优化
type
status
date
slug
summary
tags
category
icon
password
Property
 

优化问题

该优化问题的含义:寻找合适的 ,使得函数达到最小。
函数 称为目标函数或价值函数,是一个实值函数。
维向量,表示为相互独立,通常称为决策变量
集合维实数空间 的一个子集,称为约束集或可行集。如果 ,则该问题是无约束优化问题。
通常约束集可表示为 ,其中, 表示由函数组成的向量,这种形式的约束称为函数约束。
 
 

极小点

一个 元实值函数 ,定义域为 。对于定义域 中的点 如果存在 ,对于所有满足 的向量 ,不等式 都成立,则称 是函数 在定义域 上的一个局部极小点。如果对于所有 ,不等式 都成立,则称 是函数 在定义域 中的一个全局极小点。