问题条件数(Conditioning of a problem)

该博客是学习《Numerical Linear Algebra with Applications Using MATLAB》的一些总结，仅供学习使用。

通常即使使用一个稳定的算法来解决一个问题，该问题对数据中小的改变或扰动仍然是敏感的。这些扰动可能来自舍入误差(roundoff error)，收集实验数据时的小的测量误差，没有被信号过滤的噪音，近似无穷级数时的截断误差等。如果该问题落入此类，则称该问题为ill-conditioned（病态的）。

Def 1: 不管用什么算法来解决一个问题，如果它的数据中的一个小的相对变化会导致它的计算解的一个大的相对误差， 则称这个问题是病态的（ill-conditioned），如果问题数据中的小扰动导致计算解的小相对误差，则称问题是良态的(well-conditioned)。

需要搞清楚两个概念：stability 和 conditioning。

1. Stable or unstable 是对于算法而言的；
2. Well or ill-conditioned 是指特定的问题而非使用的算法。

显然，将不稳定的算法应用于病态的问题是在自找麻烦。（asking for disaster）.

但是，对于求解某个问题而言，一个稳定的算法也会产生差的结果，如果输入数据时病态的。 我们假定$f(x)$表示一个算法，$x$表示数据。 我们可以利用$f(x)$来定义well 和ill-conditioning.

Def 2. 定义$x$和$\bar{x}$分别是初始数据和扰动数据， $f(x)$和$f(\bar{x})$是对应的解，则有
如果$|x-\bar{x}|$很小，$|f(x)-f(\bar{x})|$也很小，则该问题关于数据$x$是well-conditioned.
如果$|x-\bar{x}|$很小，$|f(x)-f(\bar{x})|$很大，则该问题关于数据$x$是ill-conditioned. 如下图*

问题对数据的敏感性可以通过定义条件数(condition number)来测量。 条件数越大，则该问题对数据的改变的敏感性越高。

对于一个特定的数据$x$,假定数据中有很小的误差，问题的输入为$\bar{x}=x+\Delta x$，则计算的值为$f(\bar{x})$而非$f(x)$, 结果的相对误差除以输入数据的相对误差为：
$$\frac{\frac{|f(x)-f(\bar{x})|}{|f(x)|}}{\frac{|x-\bar{x}|}{|x|}}$$

这个比例测量了一个函数对于输入数据的误差或者改变的敏感性程度。 这导致一个特定问题$f$的条件数的数学定义。

Def 3. 一个问题$f$关于输入数据$x$的条件数定义为
$$C_f(x)=\underset{ϵ\to 0^{+}}{\lim }\underset{\Vert \delta x\Vert \le ϵ}{\sup }\frac{\frac{\Vert f(x+\delta x)-f(x)\Vert }{\Vert f(x)\Vert }}{\frac{\Vert \delta x\Vert }{\Vert x\Vert }}.$$

例1. 寻求下面三个函数的条件数
a. $f(x)=x^2,C_f(x)=|x||2x|/\left|x^2\right|=2$, 因此$f$ 是 well-conditioned.
b. $f(x)=1/(1-x),C_f(x)=|x|\left|1/(1-x{)}^2\right|/|1/(1-x)|=|x/(1-x)|.f(x)$ 是ill-conditioned 在 $x=1$ 附近，并且是 well-conditioned在其他任何地方.
c. $f(x)=e^x,C_f(x)=|x|\left|e^x/e^x\right|=|x|.f(x)$对于大的$x$是 ill-conditioned.
(b) 和 (c) 表明在$x$的特定范围内，$x$ 相对小的改变将引起$f(x)$的计算的巨大误差.

定义4. 矩阵$A$的条件数定义为$\Vert A\Vert \Vert A^{-1}\Vert$，表示为$\kappa (A)$. $\Vert \cdot \Vert$表示任意的矩阵范数.

定义5. 矩阵$A$是ill-conditioned如果$A$的条件数很大(large)；否则$A$是well-conditioned.
当然，这里的large是模糊的，条件数如果是在10^4或者更大的范围内则肯定是ill-conditioned,但对于一些矩阵而言，一个较小的条件数也可能表示ill-conditioned,决定矩阵是否是ill-conditioned并不是一门准确的科学。

cond(A); % 2-norm
cond(A,1); % 1-norm
cond(A,'inf'); % infty-norm
cond(A,'fro'); % Frobenius-norm