| 查看: 3021 | 回复: 14 | ||
| 【奖励】 本帖被评价5次，作者dnp增加金币 3.75 个 | ||
| 当前主题已经存档。 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[资源]
数据平滑方法简介（一）
|
||
|
我和lhdcsu分别发了一些帖子，并介绍了化学计量学里面最基础的两门学科：线性代数以及统计学，前者为我们的体系建立了一个模型，并提供了非常丰富而又灵活的解决方法，后者为我们提供了评价我们这个体系的方法以及如何最好地解决这个体系的思路。在化学中，数学方法解决了化学家非常多的问题，然而仍然有非常多的体系用目前的化学计量学方法无法解决，因而必须发展新的方法以解决化学中所存在的黑箱子已经非常迫切了。 数据平滑方法是色谱，波谱解析中必不可少的一个预处理步骤。因为对于我们仪器中所产生的数据，噪声存在是不可避免的，而如何有效的去除我们数据中存在的噪声以提高色谱或者波谱的信噪比，一直以来都是化学家研究的重大难题之一，并到现在已经提出了非常多的方法，从最简单的移动窗口平均法到移动窗口拟和多项式平滑，再到小波变换，傅利叶变换，以及各种各样的回归平滑方法，如分量回归平滑等。虽然平滑必不可少，然而选择一个好的平滑方法却非常重要，因为一个参数选择不好的话就可能使你的数据处理结果因为失真太严重而没有说服力，或者会达不到平滑的效果(尤其是在有异常点或者异常峰存在的时候)。 由于方法太多，因而无法穷尽所有的平滑方法，而且有很多方法我也是不知道的，所以在这里发表了一下我所知道的一些平滑方法，以为大家讨论之资，也请大家能够提出一些看法或者新的方法来互相讨论。这个帖子中我想讨论的是移动窗口多项式最小二乘拟和平滑方法，粗糙惩罚方法，以及kernel平滑方法。 移动窗口拟和多项式平滑方法(也称为Savitzky-Golay平滑)最初是由Savitzky & Golay在1964年提出来，后来经过Steinier (1973)纠正，并由Gorry给出了更一般的计算公式。我试了一下，由Gorry给出的计算公式中，用正交多项式对SG平滑的边界点进行计算，我觉得不是很好，因此只用了Steinier给出的纠正算法来平滑，在matlab里面也有一个命令 y = sgolayfilt(x,k,f)，x是平滑数据，k是平滑多项式窗口宽度，f是多项式次数。 然而，移动窗口多项式是采用的最小二乘拟和，而最小二乘拟和有一个缺陷就是在求逆的时候如果矩阵奇异的话，求逆就会得到非常大的误差。同时，如果你的数据存在异常点(outlier)的话，用最小二乘所得到的模型就会有崩溃的危险，因而在这种情况下需要用一个比较稳健的方法来进行平滑，最简单的就是移动窗口中位数平滑方法，在程序附件程序中已经给出，在图outlier中curve3就是经过中位数平滑后的结果，虽然有点失真，但蓝框框中的奇异点却被平滑。而且，上面的移动窗口平滑，对于最初的几个点和最后的几个点都无法估计，就像我在前面说的，Gorry给出了一般的SG平滑公式，边界点由正交多项式拟和，这个有时候得到的误差却很大。因此，粗糙惩罚就由统计中引入了化学计量学中，其最初由Silverman提出，他是在最小二乘判别后面加上一项惩罚项，这样既兼顾到了拟和的作用(这个是最小二乘的优势，不能丢)，同时又对数据进行了惩罚，越粗糙的数据惩罚越大，即使其对整个平滑过程的作用就越小，从而相对加大了好点的作用。 kernel方法在各个领域应用得非常广泛，其平滑方法更是在信息领域作用非凡。其也是一种移动窗口平滑方法，利用加权方法提取周围点的信息，而中心点的权重最大，这样就对窗口内的各个点进行拟和(本质上和移动窗口多项式平滑方法相同)。因此，kernel平滑方法就是利用不同的权重函数对数据进行平滑，当初最流行的是Nadaraya-Waston估计器，即权函数，后来由Gasser & Muller进行了改进，我的程序就是用改进的平滑方法。 matlab程序和算法文件都在附件的压缩文件(smooth.rar)中，程序都是我自己编的，请大家指正一下。特别是经过这么段时间的观察，论坛里面高手多多，请多多指教。 上传图片：两个图分别为有奇异点(蓝框框中，人为设置)(outlier.jpg)和无奇异点(normal.jpg)的色谱图，图中 curve 1:模拟的色谱图加高斯噪声； curve 2: 经过SG平滑，平滑窗口大小为7，多项式次数为3； curve 3: 经过中位数平滑，窗口大小为9(请注意蓝框框中的奇异峰位置)； curve 4: 用粗糙惩罚进行平滑(注意边界点)； curve 5: 经过平方kernel平滑。 最后祝大家新年快乐！  [ Last edited by 佳怡 on 2009-4-23 at 21:35 ] |
回复此楼
3楼2007-12-30 20:05:40
2楼2007-12-30 17:23:24
4楼2007-12-30 21:03:58
5楼2007-12-30 21:36:22
