cosmetic.m

偏最小二乘算法在MATLAB中的实现

% cosmetic data. 17 samples, 8 inputs, and 11 outputs;
A=[-.319  1.246  1.9574  1.7677  1.7788  .2649  -.449  .1288  .1475  1.81  -.8346  1.7301  -1.0629  .7277  -.6592  .7324  -.6409  -.5196  1.4716;
1.7056  -1.5052  -2.3426  -1.9721  .2086  .2649  -1.6328  -1.475  .1475  -1.0344  -.8346  -.859  .8351  -.8685  -.6592  -1.1656  -.6409  -.5196  -.879;
-.319  -.6044  -1.4142  -.3173  .5125  -1.719  .1225  1.4773  -.9038  -1.0344  1.3362  .2506  .0759  -.8685  -.6592  .7324  1.9533  1.6887  -.879;
.1148  -1.0426  -1.5608  -1.7735  -1.387  .2649  -.5715  -.5638  .1475  .0756  -.602  -.859  1.2147  -.8685  -.6592  .7324  -.6409  -.5196  -.879;
1.4887  .4669  -.1926  .1129  .5125  .2649  -1.0613  -.5638  1.1988  -.8263  -1.3385  -.859  1.5943  1.5258  -.6592  -.4064  -.6409  -.5196  1.2198;
-1.0059  1.5382  .2716  1.3374  .0566  .6403  -1.6328  .8213  .1475  .0756  -.098  .9904  -1.0629  .7277  -.6592  1.8712  1.9533  -.5196  .8;
-2.5967  -2.1139  .125 .7417  .3605  -1.0756  .449  1.04  .6327  .4919  1.1036  .8794  -1.0629  .3286  .1272  -.1786  1.3048  -.5196  .3802;
.1148  .4669  .5648  -.5159  -1.0577  .6403  .6531  -.0899  1.1988  -.5835  .6385  1.2123  -1.0629  .7277  -.6592  .7324  -.6409  .9526  -.879;
.8018  -.2879  -.046  -.7145  .3605  .2649  -.7348  .8213  .1475  .9428  -.602  -.859  1.2147  1.5258  -.6592  -1.9248  -.6409  -.5196  .3802;  
-.5721  -.434  .7358  .3115  -1.387  -2.0943  .3674  -.7824  1.1988  -1.4506  .6385 .9904  -1.0629  .7277  -.6592  .7324  .3048  -.5196  .879;
-.319  -.434 .125  -.3173  -.9058  -1.0756  .9797  -1.0376  .6327  .2837  .1345 -.859  -1.0629  -.8685  2.29  -1.5452  -.6409  -.5196  -.879;
.5849  1.0756  .7358  -.3173  -1.5389  -.3785  -.8981  -1.0376  .1475  -.3753  -1.3385  .2506  .0759  1.5258  -.6592  .7324  -.6409  -.5196  -.879;
-.1021  -.2879  .4182  .7417  .3605  .6403  .8981  -.3451  .1475  .4919  .6385  -.859  .911  -.8685  1.3069  .7324  -.6409  -.5196  1.2198;
.5849  .6373  .2716  .7417  .6898  1.9808  -.1225  .5661  -2.3595  -1.2425  -1.5711  -.859  .911  -.8685  1.3069  -.786  -.1221  2.4248  1.2198;  
.8018  .9295  -.3392  -.5159  -.2726  -.0568  1.8369  -.5638  -.9038  .9428  1.6075  1.4268  -1.0629  -.8685  -.6592  -.0268  -.6409  1.6887  -.879;
-.3913  -.2879  .5648  -.2511  -.07  .533  .6531  -.5658  .1475  -.3772  .9873  -.859  .1518  -.8685  .9137  -.4064  -.6409  -.5196  -.879;
-.5721  .6373  .125  .9403  1.7788  .6403  1.143  2.1699  -1.8743  1.81  .1345  -.859  .4555  -.8685  1.3069  -.5582  .6562  -.5196  1.2198];
y=A(1:11,1:11);
ytest=A(12:17,1:11);
x=A(1:11,12:19);
xtest=A(12:17,12:19);
[ax,mx,sx]=auto(x);
sxtest=scale(xtest,mx,sx);
[ay,my,sy]=auto(y);
sytest=scale(ytest,my,sy);