This paper explores the application of Probabilistic Neural Network (PNN), SupportVector Machine (SVM) and k-means clustering as tools for automated classification of massive stellar spectra. The data set consists of a set of stellar spectra associated with the Sloan Digital Sky Survey (SDSS) SEGUE-1 and SEGUE-2, which consists of 400,000 data from 3850 to 8900Å with 3646 data points each. We investigate the application of principal components analysis (PCA) to reducing the dimensionality of data set to 280, 400 and 700 components.We show that PNN can give fairly accurate spectral type classifications σ RM S = 1.752, σ RM S = 1.538 and σ RM S = 1.391 and K-means can classify these spectra with an accuracy of σ RM S = 1.812, σ RM S = 1.731 and σ RM S = 1.654 and SVM with the accuracy of σ RM S = 1.795, σ RM S = 1.674 and σ RM S = 1.529 across the 280, 400 and 700 components, respectively. By using K-means the classification of the spectra renders 38 major classes. Furthermore, by comparing the results we noticed that PNN is more successful than K-means and arXiv:1609.03147v1 [astro-ph.IM] 11 Sep 2016 2 SVM in automated classification.