Abstract: To improve and perfect the previous reliability analysis for aggregated Markov repairable systems with multi-working levels, the further study has been done in this research. By using the ion-channel modeling theory, Laplace transformation and matrix method, the probability density functions within sojourning time at each working lever, the mean sojourning time and the mean numbers that the system visits distinct working levels could be obtained. A numerical example was presented to show the results obtained in this research. It could be found that the obtained results are more realistic to reflect the evolution of the multi-working levels systems, compared to the results from previous work.