clcclear
light=importdata（‘data0.txt‘）;%讀取數據
py=light.data（:2）;%碳含量
px=light.data（:[3:end]）;%光譜強度

pz=[pxpy];?%廣義矩陣????
mu=mean（pz）;?sig=std（pz）;????%求均值和標準差
rr=corrcoef（pz）;????%求相關系數矩陣
%plot（rr（3649[1:3648]）‘k‘）;
%legend（‘光譜強度與碳含量相關系數‘1）
ab=zscore（pz）;??%數據標準化
a=ab（:?[1:end-1]）;b=ab（:end）;???%提出標準化后的自變量和因變量數據
[XLYLXSYSBETAPCTVARMSEstats]=plsregress（ab）;
xw=a\XS??;??%求自變量提出成分系數，每列對應一個成分，這里xw等stacts.W
yw=b\YS??;??%求因變量提出成分的系數
%ncomp=input（‘請根據PCTVAR的確定提出成分的個數ncomp=‘）;
ncomp=10;
[XL2YL2XS2YS2BETA2PCTVAR2MSE2stats2]=plsregress（abncomp）;
n=size（a2）;?m=size（b2）;???%n是自變量的個數，m是因變量的個數

%原始數據回歸方程的常數項
beta3（1:）=m

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1533??2015-05-18?17:47??Plsregress.m