主要内容

mslowess

用非参数法平滑峰值信号

语法

= mslowess (X强度
mslowess(…“秩序”,OrderValue,……)
mslowess(…“跨越”,SpanValue,……)
mslowess(…“内核”,KernelValue,……)
mslowess(…,'RobustIterations',RobustIterationsValue,……)
mslowess(…,'ShowPlot',ShowPlotValue,……)

参数

X 一组具有峰值的信号的分离单位值的向量。向量中的元素数等于矩阵中的行数强度.分离单元可以量化波长、频率、距离、时间或m/z,这取决于产生信号数据的仪器。
强度 具有相同分离单位范围的一组峰值的强度值矩阵。每一行对应一个分离单元值,每一列对应一组有峰值的信号或保留时间。行数等于向量中的元素数X

描述

提示

使用以下语法处理来自产生信号数据的任何分离技术的数据,如光谱学、核磁共振、电泳、色谱或质谱。

= mslowess (X强度平滑原始噪声信号数据,强度,使用默认跨度为的局部加权线性回归(Lowess)方法10样本。

请注意

mslowess假设输入向量,X,可能没有均匀间距的分离单元。因此,用于平滑的滑动窗口的居中使用最接近的样本X价值而非以之而论X索引。

请注意

当输入向量,X,没有重复值或NaN值,算法的速度大约是原来的两倍。

mslowess (X强度,……”PropertyName',PropertyValue,……)调用mslowess使用属性名/属性值对的可选属性。可以以任意顺序指定一个或多个属性。每一个PropertyName必须用单引号括起来,不区分大小写。这些属性名/属性值对如下:

mslowess(…“秩序”,OrderValue,……)指定顺序(OrderValue)。输入1(线性多项式拟合或Lowess),2(二次多项式拟合或黄土),或0(相当于加权局部均值估计器,而且可能更快,因为只执行平均值计算,而不是最小二乘回归)。默认值为1

请注意

曲线拟合工具箱™软件也指低阶平滑2黄土平滑。

mslowess(…“跨越”,SpanValue,……)指定平滑内核的窗口大小。如果SpanValue大于1,窗口等于SpanValue独立于分离单位向量的样本数,X.默认值为10样本。较高的值将使信号更加平滑,但要牺牲计算时间。如果SpanValue小于1,窗口大小被认为是数据中点数的一个分数。例如,当SpanValue0.005,窗口大小等于0.50中点数的%X

mslowess(…“内核”,KernelValue,……)所指定的函数KernelValue用于对观测到的强度进行加权。靠近被平滑的分离单元位置的样本在确定估计值时具有最大的权重。KernelValue可以是以下字符向量(或字符串)中的任何一个:

  • “tricubic”(默认)(1 - (dist /距离)。^ 3)。^ 3

  • “高斯”- - - - - -exp (- (2 * dist /距离)^ 2)。

  • “线性”- - - - - -1-dist /距离

mslowess(…,'RobustIterations',RobustIterationsValue,……)指定迭代次数(RobustValue)适合强健的身体。如果RobustIterationsValue0(默认),不执行健壮拟合。对于鲁棒平滑,在每个跨度上的小残差值被忽略,以改进新的估计。12健壮的迭代通常就足够了,而较大的值在计算上可能会比较昂贵。

请注意

对于一个X具有均匀间距分离单元的向量,具有非鲁棒平滑OrderValue等于0相当于用核向量对信号进行滤波。

mslowess(…,'ShowPlot',ShowPlotValue,……)在原始信号上绘制平滑信号。当你打电话mslowess没有输出参数,信号被绘制,除非ShowPlotValue.当ShowPlotValue真正的,只有第一个信号强度是绘制。ShowPlotValue也可以包含到?中的一个信号的索引强度

例子

  1. 加载一个mat文件,包括生物信息学工具箱™软件,其中包含一些示例数据。

    负载sample_lo_res
  2. 平滑频谱,画出原始信号和平滑信号的第一个频谱图。

    y = mslowess (MZ_lo_res Y_lo_res,“Showplot”,真正的);

  3. 放大图中的一个区域,看看原始信号和平滑信号的区别。

    轴([7350 7550 0.1 1.0])

之前介绍过的R2006a