aa.m

MATLAB科学计算与工程分析源代码源程序3

 X=[191 36 50;
189 37 52;
193 38 58;
162 35 62;
189 35 46;
182 36 56;
211 38 56;
167 34 60;
176 31 74;
154 33 56;
169 34 50;
166 33 52;
154 34 64;
247 46 50;
193 36 46;
202 37 62;
176 37 54;
157 32 52;
156 33 54;
138 33 68
];
Y=[5 162 60;
2 110 60;
12 101 101;
12 105 37;
13 155 58;
4 101 42;
8 101 38;
6 125 40;
15 200 40;
17 251 250;
17 120 38;
13 210 115;
14 215 105;
1 50 50 ;
6 70 31;
12 210 120;
4 60 25;
11 230 80;
15 225 73;
2 110 43];
[n,p]=size(X);
[n,q]=size(Y);
meanX=mean(X);%junzhi
varX=var(X);%fangcha
meanY=mean(Y);%junzhi
varY=var(Y);%fangcha
for i=1:p      
    for j=1:n
    X0(j,i)=(X(j,i)-meanX(i))/((varX(i))^0.5);
    end
end
for i=1:q
    for j=1:n
    Y0(j,i)=(Y(j,i)-meanY(i))/((varY(i))^0.5);
    end
end
V1=inv(X0'*X0)*X0'*Y0*inv(Y0'*Y0)*Y0'*X0;
V2=inv(Y0'*Y0)*Y0'*X0*inv(X0'*X0)*X0'*Y0;
[VV1,lamda1]=eig(V1);
[VV2,lamda2]=eig(V2);

lamda11=sum(lamda1);
lamda21=sum(lamda2);
k=1;
while(k<=(length(lamda1)))
    [maxlamda1(k),I]=max(lamda11);   
    maxVV1(:,k)=VV1(:,I);
    lamda11(I)=[];
    VV1(:,I)=[];
     [maxlamda2(k),I]=max(lamda21);
    maxVV2(:,k)=VV2(:,I);
     lamda21(I)=[];
    VV2(:,I)=[];
    k=k+1;
end


F=X0*maxVV1
G=Y0*maxVV2