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

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

%原始數(shù)據(jù)回歸方程的常數(shù)項(xiàng)
beta3（1:）=m

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