各位大虾，怎么用matlab做典型相关分析??THX

r=[X,Y];
a=zscore(r);
R=cov(a)
p=1;q=1;n=1144;%p,q为X,Y的列数，n为行数
R11=R(1:p,1:p);R12=R(1:p,p+1:p+q);R22=R(p+1:p+q,p+1:p+q);R21=R(p+1:p+q,1:p);
[v1,d1]=eig(R11),
[v2,d2]=eig(R22)
p1=inv(v1*sqrt(d1)*v1');
p2=inv(v2*sqrt(d2)*v2');
T1=p1*R12*inv(R22)*R21*p1;
T2=p2*R21*inv(R11)*R12*p2;
[va,da]=eig(T1),[vb,db]=eig(T2),
A1=p1*va,
B1=p2*vb,
r=sqrt(sum(da))