Due to the complex vertical distribution of aerosols in the atmosphere, the commonly used gradient method and wavelet covariance transform (WCT) method still have large uncertainties in automatically and continuously identifying the planetary boundary layer height (PBLH). For this reason, we applied the 2-D matrix method to the retrieval of PBLH in two cases observed by a dual-wavelength lidar over Jinhua city. The 2-D matrix method optimizes PBLH retrieval results from the distance and time dimensions. The results show that in the case of profiles with multi-layer aerosols, the 2-D matrix method is superior to the WCT method in terms of PBLH retrieval. The correlation coefficient between the PBLH retrieved by the 2-D matrix method at 532 nm and that by the WCT method at 1064 nm is 0.87. However, the correlation coefficient is only 0.37 if the 2-D matrix method is replaced by the WCT method at 532 nm. Additionally, the consistent variation trend in PBLH and surface temperature indicates reliable PBLH retrieval results. The 2-D matrix method provides a reference for the automatic and continuous identification of the PBLH by lidar in the absence of low clouds (no clouds appearing below 3 km).