orthogonalize.m

来自「动态模糊神经网络的例程」· M 代码 · 共 16 行

function [w,a]=orthogonalize(p)
[u,v]=size(p);
w(:,1)=p(:,1);
a=eye(v);
for k=2:v
b=zeros(u,1);
   for i=1:k-1
      a(i,k)=w(:,i)'*p(:,k)/(w(:,i)'*w(:,i));
      b=b+a(i,k)*w(:,i);
   end
   
   w(:,k)=p(:,k)-b;
end