Phase-amplitude coupling (PAC) measures the relationship between the phase of low-frequency oscillations (LFO) and the amplitude of high-frequency oscillations (HFO).It plays an important functional role in neural information processing and cognition. Thus, we propose a novel method based on the Jensen-Shannon (JS) divergence and correlation matrix. The method takes the amplitude distributions of the HFO located in the corresponding phase bins of the LFO as multichannel inputs to construct a correlation matrix, where the elements are calculated based on the JS divergence between pairs of amplitude distributions. Then, the omega complexity extracted from the correlation matrix is used to estimate the PAC strength. The simulation results demonstrate that the proposed method can effectively reflect the PAC strength and slightly vary with the data length. Moreover, it outperforms five frequently used algorithms in the performance against additive white Gaussian noise and spike noise and the ability of detecting PAC in wide frequency ranges. To validate our proposed method with real data, it was applied to analyze the local field potential recorded from the dorsomedial striatum in a male Sprague-Dawley rat. It can replicate previous results obtained with other PAC metrics. In conclusion, these results suggest that our proposed method is a powerful tool for measuring the PAC between neural oscillations.