Clustering is to divide a dataset into multiple clusters based on the similarity between samples. Most existing clustering methods face two challenges. On the one hand, when defining the similarity between samples, the spatial distribution structure of the samples is often not considered, making it difficult to construct a stable similarity matrix. On the other hand, the sample graph structure constructed by graph clustering is too complex and has high computational costs. To solve these two problems, a Multi-order Nearest Neighbor Graph Clustering algorithm by fusing transition probability matrix (MNNGC) was proposed. Firstly, the nearest neighbor relationship and spatial distribution structure of samples were comprehensively considered, the similarity defined by shared nearest neighbor was weighted for densification, and the densification affinity matrix between nodes was obtained. Secondly, by utilizing multi-order probability transition between nodes, the correlation degrees of non-adjacent nodes were predicted, and a stable inter-node affinity matrix was obtained by fusing the multi-order transition probability matrix. Then, to further enhance the local structure of the graph, the multi-order nearest neighbor graph of nodes was reconstructed, and hierarchically clustered. Finally, the edge node allocation strategy was optimized. Positioning experimental results show that MNNGC achieves the highest Accuracy (Acc) among comparison clustering algorithms on all the synthetic datasets and 8 UCI datasets. The Acc, Adjusted Mutual Information (AMI), Adjusted Rand Index (ARI) and Fowlkes and Mallows Index (FMI) of MNNGC algorithm are improved by 38.6, 27.2, 45.4 and 35.1 percentage points, respectively, compared with Local Density Peaks-based Spectral Clustering (LDP-SC) algorithm.