yinzi.txt

matlab中的移植算法

pn=size(x);
n=pn(1);
m=pn(2);
x
for j=1:m
   xx(j)=sum(x(:,j))/n;
end 	
for j=1:m
  s1(j)=sum((x(:,j)-xx(j)).^2);
  s(j)=sqrt(s1(j)/n);
end
xx;
s1;
for i=1:n
   for j=1:m
     x(i,j)=(x(i,j)-xx(j))/s(j); 
   end     
end  
x;    
for j=1:m
   xx(j)=sum(x(:,j))/n;
end 
for j=1:m
  s1(j)=sum((x(:,j)-xx(j)).^2);
  s(j)=sqrt(s1(j)/n);
end
xx;
s1;
for t=1:m
   for s=1:m
      r(t,s)=sum((x(:,t)-xx(t)).*(x(:,s)-xx(s)));
      r(t,s)=r(t,s)/sqrt(s1(t)*s1(s));
   end
end
r
[v1,d]=eig(r)
aa=v1*sqrt(d)
nmd1=eig(r);
nmd2=nmd1;
nmd1=sort(nmd1);
for i=1:m
   nmd(i)=nmd1(m+1-i);
end
nmd
for i=1:m
   d(i,i)=nmd(i);
end
for i=1:m
   for j=1:m
   if nmd(i)==nmd2(j)
     xihao(i)=j;
   end
   end
end
for i=1:m
   v(:,i)=v1(:,xihao(i));
end
v
aa=v*sqrt(d)
sum2=0;
sum1=nmd/sum(nmd)
for i=1:m
   sum2=sum2+sum1(i);
   if sum2>0.85
      i
      break
    end  
end
k=i;
A=aa(:,1:k)
for j=1:n
   y(:,j)=A'*inv(r)*(x(j,:))';
end 
b=y'
for i=1:m
   h(i)=sum(A(i,:).^2);
end	
for j=1:k
   g(j)=sum(A(:,j).^2);
end   
hold off
hold on
for i=1:n
   plot(b(i,1),b(i,2),'.');
%text(b(i,1)+0.02,b(i,2)+0.02,num2str(i));
end
h
g