diff --git a/oldMcode/compare.m b/oldMcode/compare.m deleted file mode 100644 index 3383e9b..0000000 --- a/oldMcode/compare.m +++ /dev/null @@ -1,26 +0,0 @@ -maxSize=6; - -mseErrorAutoenc=zeros([maxSize,1]); -for i=1:maxSize - ae=trainAutoencoder(trainingData,i,'ShowProgressWindow',false); - predicted=predict(ae,trainingData); - mseErrorAutoenc(i)=mse(trainingData-predicted); -end - -fig=figure(); -subplot(1,3,1); -plot(mseErrorAutoenc) - -%PCA -[COEFF,SCORE,latent] = princomp(trainingData); -subplot(1,3,2); -plot(cumsum(latent)./sum(latent)) - -%NNMF -D=cell(maxSize,1); -for i=1:maxSize - [W,H,D{i}]=nnmf(trainingData,i); -end - -subplot(1,3,3); -plot(diff(cell2mat(D))) \ No newline at end of file diff --git a/oldMcode/noOfSynergies.m b/oldMcode/noOfSynergies.m new file mode 100644 index 0000000..6dea77d --- /dev/null +++ b/oldMcode/noOfSynergies.m @@ -0,0 +1,42 @@ +maxSize=6; + +mseAutoenc=zeros([maxSize,1]); +for i=1:maxSize + ae=trainAutoencoder(trainingData,i,'ShowProgressWindow',false); + predicted=predict(ae,trainingData); + mseAutoenc(i)=mse(trainingData-predicted); +end + +fig=figure(); +subplot(1,3,1); +plot(mseAutoenc) +title('Autoencoder') +xlabel('number of syergies') +ylabel('mse') + +%PCA +[COEFF,SCORE,latent] = pca(trainingData,'Centered',false); +msePCA=zeros([maxSize,1]); +for i=1:maxSize + predicted=SCORE(:,1:i)*COEFF(:,1:i)'; + msePCA(i)=mse(trainingData-predicted); +end +subplot(1,3,2); +plot(msePCA) +title('PCA') +xlabel('number of syergies') +ylabel('mse') + +%NNMF +mseNNMF=zeros([maxSize,1]); +for i=1:maxSize + [W,H]=nnmf(trainingData,i); + predicted=W*H; + mseNNMF(i)=mse(trainingData-predicted); +end + +subplot(1,3,3); +plot(mseNNMF) +title('NNMF') +xlabel('number of syergies') +ylabel('mse') \ No newline at end of file