這里提供的matlab函數(shù)有四個(gè)，分別是移動(dòng)窗口最小二乘多項(xiàng)式平滑(sgsmooth.m)，移動(dòng)窗口中位數(shù)平滑(median_smooth.m)，粗糙懲罰平滑(rpsmooth.m)以及kernel平滑(kernel.m)。
      移動(dòng)窗口最小二乘多項(xiàng)式平滑采用的格式是：
       y=sgsmooth(x,width,order)
       x是數(shù)據(jù)：光譜或者是色譜，而且輸入形式必須為一個(gè)向量，width是平滑時(shí)選擇的窗口寬度，order是多項(xiàng)式擬和的次數(shù)。需要注意的是多項(xiàng)式擬和的次數(shù)order必須小于width，而且width必須為一個(gè)奇數(shù)，否則提示error。

       移動(dòng)窗口中位數(shù)平滑的格式：
        y=median_smooth(x,width)
        x是光譜或者色譜數(shù)據(jù)，同樣輸入必須為一個(gè)向量，width是移動(dòng)平滑窗口寬度。

        粗糙懲罰平滑格式：
        f=rpsmooth(x,nambda)，x為光譜或者色譜數(shù)據(jù)，也必須為一個(gè)向量。nambda為懲罰參數(shù)，可以為一個(gè)數(shù)，也可以為一個(gè)含有兩個(gè)數(shù)的向量，nambda向量中的元素但是不能超過2，因?yàn)槿绻麑Υ植趹土P有所認(rèn)識，可以直接輸入一個(gè)數(shù)進(jìn)行計(jì)算，而如果對參數(shù)選擇沒有把握，則可以以[up low]的向量形式輸入兩個(gè)值進(jìn)行計(jì)算。nambda不能為0。如果nambda為0，無平滑效果。如果沒有選擇nambda值，默認(rèn)的是從1：100進(jìn)行循環(huán)計(jì)算，步長為10。
        需要注意的是，在matlab中產(chǎn)生一個(gè)矩陣Q和R(詳見算法文檔或者分析化學(xué)手冊第十分冊)時(shí)間花費(fèi)比較大，如果數(shù)據(jù)點(diǎn)大于1000，建議分段計(jì)算，我計(jì)算了一下，用1000個(gè)數(shù)據(jù)點(diǎn)進(jìn)行計(jì)算，時(shí)間花費(fèi)是70~80s。

         kernel平滑格式：
        smooth_d=kernel(t,x,h,fun)
         x同樣是一個(gè)光譜或者色譜數(shù)據(jù)向量。t為光譜中的波長點(diǎn)或者色譜中的保留時(shí)間點(diǎn)，如果沒有這些數(shù)據(jù)的，可以用點(diǎn)數(shù)來代替，比如x數(shù)據(jù)點(diǎn)有500個(gè)，可用其索引1,2,3,…等來組成一個(gè)向量t。h為帶寬，必須為正數(shù)。fun為kern函數(shù)，在matlab文件中提供了四個(gè)函數(shù)：均勻函數(shù)，二次函數(shù)，一階導(dǎo)數(shù)估計(jì)和二階導(dǎo)數(shù)估計(jì)，詳細(xì)請見算法文檔，但是一階導(dǎo)數(shù)和二階導(dǎo)數(shù)估計(jì)在平滑中沒有用到，因此沒有給予介紹。分別用'uni'，'qua'，'fde'，'sde'表示。默認(rèn)的是采用二次函數(shù)平滑。