【正文】 )附錄10：求白葡萄與白葡萄酒理化指標(biāo)之間聯(lián)系程序（）clc,clear ,format long gpz=。)。,sol,39。xlswrite(39。sheet539。39。 %計(jì)算原始數(shù)據(jù)回歸方程的系數(shù) end sol=[ch0。 %計(jì)算原始數(shù)據(jù)回歸方程的常數(shù)項(xiàng) for i=1:m xish(:,i)=xishu(:,i)./sig_x39。sig_y=sig(n+1:end)。mu_y=mu(n+1:end)。 %求 Y*關(guān)于 t 的回歸系數(shù) xishu=w_star*beta_z。,i)。 end if Q_h2(i) fprintf(39。 Q_h2(1)=1。 %求殘差向量 press_i(j)=sum(cancha.^2)。 %刪除第 j 個(gè)觀測值 beta1=[t1,ones(num1,1)]\f1。 %把舍去的第 j個(gè)樣本點(diǎn)保存起來 t1(j,:)=[]。 she_t=t1(j,:)。 %求誤差平方和 %以下計(jì)算 press(i) for j=1:num t1=t(:,1:i)。 %求回歸方程的系數(shù)，數(shù)據(jù)標(biāo)準(zhǔn)化，沒有常數(shù)項(xiàng) cancha=f0t*beta。 %計(jì)算殘差矩陣 e0=e。)。*t(:,i))。 %計(jì)算成分 ti 的得分 alpha=e039。 %提出最大特征值對應(yīng)的特征向量 w_star(:,i)=chg*w(:,i)。)。 %提出對角線元素，即提出特征值 [val,ind]=sort(val,39。 [vec,val]=eig(matrix)。*f0*f039。%求樣本點(diǎn)的個(gè)數(shù) chg=eye(n)。f0=data(:,n+1:end)。y0=pz(:,n+1:end)。m=3。 %求相關(guān)系數(shù)矩陣 data=zscore(pz)。sig=std(pz)。)。,39。附錄9：求紅葡萄與紅葡萄酒理化指標(biāo)之間聯(lián)系程序（）clc,clear ,format long gpz=xlsread(39。sheet539。39。)。,stf,39。xlswrite(39。sheet339。39。)。,y,39。xlswrite(39。sheet139。39。, ind=ind39。)。 %計(jì)算綜合得分 [stf,ind]=sort(tf,39。 %num 為選取的主成分的個(gè)數(shù) df=jiubiao*x(:,1:num)。 %構(gòu)造與x 同維數(shù)的元素為177。 %數(shù)據(jù)標(biāo)準(zhǔn)化 r=corrcoef(jiubiao)。)。,jiubiao14,39。xlswrite(39。bjfxw=bjfxw39。 bjfxw(:,i)=t1.*fenzim。bjfxw=zeros(n,28)。)。,39。,39。fenzim=xlsread(39。)。,39。,39。%下面考慮的是釀酒白葡萄的芳香物質(zhì)指標(biāo)fxbj=xlsread(39。)/4。)。,39。,39。seze=xlsread(39。)/3。)。,39。,39。dpph=xlsread(39。)/2。)。,39。,39。chun=xlsread(39。huangtong=huangtong39。huangtong=sum(huangtong39。P33:R6039。jiu39。39。zongfen=zongfen.*。)/3。)。,39。,39。zongfen=xlsread(39。danning=danning39。danning=sum(danning39。H33:J6039。jiu39。39。)。,ind,39。xlswrite(39。sheet439。39。)。,z,39。xlswrite(39。sheet239。39。)。,x,39。 xlswrite(39。 %把得分按照從高到低的次序排列 stf=stf39。descend39。 %計(jì)算各個(gè)主成分的得分 tf=df*z(1:num)/100。1 的矩陣 x=x.*f %修改特征向量的正負(fù)號，每個(gè)特征向量乘以所有分量和的符號函數(shù)值 num=3。 %計(jì)算相關(guān)系數(shù)矩陣 %下面利用相關(guān)系數(shù)矩陣進(jìn)行主成分分析，x 的列為 r 的特征向量，即主成分的系數(shù) [x,y,z]=pcacov(r) %y為 r 的特征值，z 為各個(gè)主成分的貢獻(xiàn)率 f=repmat(sign(sum(x)),size(x,1),1)。jiubiao=zscore(jiubiao)。sheet139。39。jiubiao=[hjfxw huagan danning zongfen huangtong chun dpph seze]jiubiao13=jiubiao(:,1:3)。endhjfxw=sum(hjfxw)。for i=1:27 t1=fxhj(:,i)。n=length(fenzim)。C2:C7439。hongjiu39。39。fxhj(isnan(fxhj))=0。E2:AE7439。hongjiu39。39。seze=seze39。seze=sum(seze39。AD3:AH2939。jiu39。39。dpph=dpph39。dpph=sum(dpph39。Z3:AB2939。jiu39。39。chun=chun39。chun=sum(chun39。T3:X2939。jiu39。39。huangtong=huangtong.*。)/3。)。,39。,39。huangtong=xlsread(39。zongfen=zongfen39。zongfen=sum(zongfen39。L3:N2939。jiu39。39。danning=danning.*。)/3。)。,39。,39。danning=xlsread(39。)/3。)。,39。,39。附錄7：求紅葡萄酒的分級程序（）clc,clear %紅葡萄酒huagan=xlsread(39。sheet539。39。)。,stf,39。xlswrite(39。sheet339。39。)。,y,39。xlswrite(39。sheet139。39。, ind=ind39。)。 %計(jì)算綜合得分 [stf,ind]=sort(tf,39。 %num 為選取的主成分的個(gè)數(shù) df=putao*x(:,1:num)。 %構(gòu)造與x 同維數(shù)的元素為177。 %數(shù)據(jù)標(biāo)準(zhǔn)化 r=corrcoef(putao)。)。,putao110,39。xlswrite(39。putao=[putao1 putao2 putao3]。putao2=[tongchun,zongtang,yuantang,guowu,ph,dingsuan,huoli,hebiandu,ziyouji]。bjfxw=bjfxw39。 bjfxw(:,i)=t1.*fenzim。bjfxw=zeros(n,28)。)。,39。,39。fenzim=xlsread(39。)。,39。,39。%下面考慮的是釀酒白葡萄的芳香物質(zhì)指標(biāo)fxbj=xlsread(39。EC34:EC6139。putao39。39。)。,39。,39。axing=xlsread(39。DU34:DU6139。putao39。39。)。,39。,39。guopi=xlsread(39。DM34:DM6139。putao39。39。)。,39。,39。guogengbi=xlsread(39。DE34:DE6139。putao39。39。)。,39。,39。guosui=xlsread(39。CW34:CW6139。putao39。39。)。,39。,39。gusuan=xlsread(39。CO34:CO6139。putao39。39。)。,39。,39。ph=xlsread(39。CG34:CG6139。putao39。39。)。,39。,39。yuantang=xlsread(39。BY34:BY6139。putao39。39。)。,39。,39。tongchun=xlsread(39。BO34:BO6139。putao39。39。)。,39。,39。huangtong=xlsread(39。BE34:BE6139。putao39。39。)。,39。,39。zongfen=xlsread(39。AW34:AW6139。putao39。39。)。,39。,39。hebiandu=xlsread(39。)/3。)。,39。,39。huoli=xlsread(39。AJ34:AJ6139。putao39。39。)。,39。,39。%shisuan=shisuanguosuan=xlsread(39。AF34:AF6139。putao39。39。huagan=huagan39。huagan=sum(huagan39。AB34:AD6139。putao39。39。vc=vc39。vc=sum(vc39。Y34:Z6139。putao39。39。baizhi=baizhi39。baizhi=sum(baizhi39。V34:W6139。putao39。39。suan=suan39。suan=sum(suan39。C34:T6139。putao39。39。)。,ind,39。xlswrite(39。sheet439。39。)。,z,39。xlswrite(39。sheet239。39。)。,x,39。 xlswrite(39。 %把得分按照從高到低的次序排列 stf=stf39。descend39。 %計(jì)算各個(gè)主成分的得分 tf=df*z(1:num)/100。1 的矩陣 x=x.*f %修改特征向量的正負(fù)號，每個(gè)特征向量乘以所有分量和的符號函數(shù)值 num=9。 %計(jì)算相關(guān)系數(shù)矩陣 %下面利用相關(guān)系數(shù)矩陣進(jìn)行主成分分析，x 的列為 r 的特征向量，即主成分的系數(shù) [x,y,z]=pcacov(r) %y為 r 的特征值，z 為各個(gè)主成分的貢獻(xiàn)率 f=repmat(sign(sum(x)),size(x,1),1)。putao=zscore(putao)。sheet339。39。putao19=putao(:,1:9)。putao3=[vc,gusuan,ganwu,guosui,baili,guogengbi,chuzhilv,guopi,lxing,axing,bxing]。putao1=[hjfxw,suan,danning,zongfen,huagan,baizhi,shisuan,guosuan,mengsuan]。endhjfxw=sum(hjfxw)。for i=1:27 t1=fxhj(:,i)。n=length(fenzim)。C3:C5739。hongpu39。39。fxhj(isnan(fxhj))=0。E3:AE5739。hongpu39。39。)。,39。,39。bxing=xlsread(39。DY3:DY2939。putao39。39。)。,39。,39。lxing=xlsread(39。DQ3:DQ2939。putao39。39。)。,39。,39。chuzhilv=xlsread(39。DI3:DI2939。putao39。39。)。,39。,39。baili=xlsread(39。DA3:DA2939。putao39。39。)。,39。,39。ganwu=xlsread(39。CS3:CS2939。putao39。39。)。,39。,39。dingsuan=xlsread(39。CK3:CK2939。putao39。39。)。,39。,39。guowu=xlsread(39。CC3:CC2939。putao39。39。)。,39。,39。zongtang=xlsread(39。BU3:BU2939。putao39。39。)。,39。,39。luchun=xlsread(39。BI3:BI2939。putao39。39。)。,39。,39。danning=xlsread(39。BA3:BA2939。putao39。39。)。,39。,39。ziyouji=xlsread(39。AS3:AS2939。putao39。39。huoli=huoli39。huoli=sum(huoli39。AL3:AN2939。putao39。39。)。,39。,39。mengsuan=xlsread(39。AH3:AH2939。putao39。39。)。,39。,39。shisuan=xlsread(39。)/3。)。,39。,39。huagan=xlsread(39。)/2。)。,39。,39。vc=xlsread(39。)/2。)。,39。,39。baizhi=xlsread(39。)/2。)。,39。,39。subplot(2,2,1),hist(bai1)subplot(2,2,2),hist(bai2)subplot(2,2,3),normplot(bai1)subplot(2,2,4),normplot(bai2)[muhat1,sigmahat1,muci1,sigmaci1]=normfit(bai1)[muhat2,sigmahat2,muci2,sigmaci2]=normfit(bai2) 附錄5：求紅葡萄的分級程序（）clc,clear %紅葡萄suan=xlsread(39。bai1=bai1(:)。M3:V3039。sheet539。39。)。,39。,39。subplot(2,2,1),hist(hong1)subplot(2,2,2),hist(hong2)subplot(2,2,3),normplot(hong1)subplot(2,2,4),normplot(hong2)[muhat1,sigmahat1,muci1,sigmaci1]=normfit(hong1)[muhat2,sigmahat2,muci2,sigmaci2]=normfit(hong2)附錄4：求兩組紅葡萄酒樣品的總體分布圖程序（）clc,clearbai1=xlsread(3