偏最小二乘建模的全过程MATLAB程序与结果

一•利用160组数据建PLS回归模型。>>clear>>loadysj>>X=ysj(:,1:8);>>Y=ysj(:,9:11);>>E0=stand(X);>>F0=stand(Y);>>A=rank(E0);>>[W,C,T,U,P,R]=bykpcr(E0,F0);W：自变量轴权重；C：因变量轴权重；T：自变量系统主成分得分；U：因变量系统主成分得分；P：模型效应载荷量；R：因变量载荷量。.确定主成分个数三种方法:1)复测定系数：R2hR2h弋(灯刈：|2)复测定系数表示所提取的主成分的可解释变异信息占总变异的百分比。当h=m，复测定系数的值足够大时，可再第m步终止主成分的提取计算。通常R2>0.85m即可。>>RA=plsra(T,R,F0,A)RA=0.33900.48310.57310.63580.64880.65220.65310.33900.48310.57310.63580.64880.65220.65310.6537结论：利用这个方法，无法确定。(2)类似典型相关分析中的精度分析方法：>>[Rdx,RdX,RdXt,Rdy,RdY,RdYt]=plsrd(E0,F0,T,A)Rdx=0.30340.43480.05390.13260.00820.01320.03310.02080.26610.19180.05490.19320.18520.00010.04160.06710.04000.10100.32810.01910.45570.05290.00020.00300.02060.08130.48680.04920.04690.30260.00210.01040.00160.04720.58690.09210.01260.19550.01010.05400.26670.22290.25170.00020.04470.06380.06340.0866

0.27460.18590.01120.00410.00060.04340.45690.02330.54670.44300.00180.00010.00720.00030.00080.0001RdX=0.21500.21350.22190.06130.09510.08400.07610.0332RdXt=1.0000Rdy=0.00920.00020.13250.04380.01950.00190.00020.00030.07610.06130.01120.05680.00010.00250.00010.00060.45910.16970.00090.00000.00130.00100.00110.0001RdY=0.18140.07710.04820.03360.00700.00180.00050.0003RdYt=0.3498>>[V]=LJRdX(RdX)V=0.21500.42840.65040.71170.80680.89080.96681.0000(3)累计贡献率：>>[U]=LJGXL(X,T,A)U=0.17560.38460.57910.63080.71980.79810.87110.9043（4）交叉有效性由于不会编交叉有效性的MATLAB程序，因此，没再验证。>>h=6;•计算PLS回归方程的系数首先，求标准化变量F,F,F关于E,E,E,E,E,E,E,E的回归系数。0102030102030405060708>>SCOEFF=pls(h,8,W,P,R)SCOEFF=-0.0103-0.01190.00730.01140.0131-0.0081-0.0072-0.00820.00510.03980.0458-0.0284-0.0282-0.03240.0201-0.0221-0.02540.0157-0.0048-0.00550.00340.00070.0008-0.0005其次，求原始变量yl,y2,y3,关于xl,x2,x3,x4,x5,x6,x7,x8的回归系数。>>[COEFF,INTERCEP]=plsiscoeff(X,Y,SCOEFF)COEFF=-0.0000-0.00000.00390.00010.0000-0.0331-0.1944-0.019858.11930.06950.0071-20.7671-0.0192-0.00195.7253-0.0011-0.00010.3146-0.0013-0.00010.37680.00000.0000-0.0004INTERCEP=0.43710.0292309.2128（三）.建立PLS回归方程xl-风量;x2-热风温度;x3-热风压;x4-透气性;x5-富氧率;x6-煤粉量;x7-料速（料批）；x8-铁量差;yl-Si;y2-S;y3-铁元素还原速率；yl二0.4371—0xl+O.OOOlx2-0.1944x3+0.0695x4-0.0192x5-O.OOllx6-0.0013x7+0x8y2二0.0292-0xl+0x2-0.0198x3+0.0071x4-0.0019x5-O.OOOlx6-O.OOOlx7+0x8y3=309.2128+0.0039xl-0.0331x2+55.1193x3-20.7671x4+5.7253x5+0.3146x6+0.3768x7-0.0004x8二.利用40组数据来计算命中率。Si—0.1;S--0.01;铁元素还原速率—50。[M]=mzl(XX,YY)M=0.52500.32500.7250若改成：Si—0.15;S--0.015;铁元素还原速率—60。[M]=mzl(XX,YY)M=0.80000.82500.8750三.给出预测值与实际值的折线图

四.李志玲建的模型风量、风温、风压、透气性、富氧、喷煤、料速、铁量差；铁水硅含量［Si］、铁水硫含量［S］、炉渣碱度［R］、铁还原率。其中y=\y^yM/=凶乂2，耳九乂5乂