load hdata.txt
h=hdata;
[m,n]=size(h)
subplot(3,1,1)
plot(-h')
title('Profiles')
xlabel('x')
a=h'*h/(m*n)
[e,lambda]=eig(a)
cnk=h*e
subplot(3,1,2)
plot(-e)
xlabel('x')
title('Eigenvectors')
subplot(3,1,3)
plot(cnk)
xlabel('n')
ylabel('C_{nk}');
title('Weights');