clc
clear
close?all
load?leafname.mat;???%載入光譜數據，一行代表一個樣本
matrix?=?leafspectral_matrixZT;???????%?將載入的光譜數據賦給matrix???
leafmsc?=?msc（matrix）;????????????????%?對光譜數據進行msc多元散射校正
%%???????????????????????????????繪制圖像進行對比?msc前后的對比
plot（leafmsc‘）
figure（2）
plot（matrix‘）
leafmsc?=?leafmsc‘;
%%???????????????????????????????進行?sgolya平滑?和SG-1階平滑
N?=?input（‘請輸入擬合次數‘）;????%?設置SG濾波器中多項式的階數
F?=?input（‘設置窗口參數（奇數）‘）;???%?設置SG濾波器中框架的邊長
%%
[bg]=sgolay（NF）;?????????????%?sgolay濾波器函數調用
???????????????????????????????%?b的維度是F*F的.b的每一行代表一個濾波器的系數。中間那行用于穩態數據的平滑，
???????????????????????????????%b中間那行的前面所有行代表信號的末端數據，
???????????????????????????????%后面所有行代表信號起始階段的數據。
[m?n]?=?size（leafmsc）
%%???使用b中數據進行SG平滑
figure（3）
for?i?=?1:n
????ycenter?=?conv（leafmsc（:i）b（（F+1）/2:）‘valid‘）;????%?求穩態數據的平滑結果
????ybegin?=?b（end:-1:（F+3）/2:）?*?leafmsc（F:-1:1i）;???%?求初始階段數據的平滑結果（b和intensity的數據是矩陣想乘的關系，但要注意二者數據的順序進行了顛倒）
????yend?=?b（（F-1）/2:-1:1:）?*?leafmsc（end:-1:end-（F-1）i）;??%?求末階段段數據的平滑結果（b和intensity的數據是矩陣想乘的關系，但要注意二者數據的順序進行了顛倒）

????Y?=?[ybegin;ycenter;yend];
????plot（Y）
????hold?on;????
end
%%?使用g中數據進行SG平滑。g中數據代表不同階的SG平滑效果。第一列代表0階，第二列代表1階.....
figure（4）?????%?繪制SG_2階的數據
for?j?=?1:n
???dy?=?zeros（m4）;
???for?p?=?0:3
???????dy（:p+1）?=?conv（leafmsc（:j）?factorial（p）/（-4.44）^p?*?g（:p+1）?‘same‘）;??%?g（:1）:平滑系數，g（:2）:一階導數系數。利用求導的方法求解系數，式中/（4.44）^p
???end
????
????plot（dy（:2））;????%?改變??dy（:value）中的value可以畫出不同階數的SG平滑，1代表0階，2代表1階
????hold?on
end
?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1859??2018-09-01?10:08??spectralmatrix_msc_sgolay.m

-----------?---------??----------?-----??----

?????????????????1859????????????????????1