In order to make the distribution of lateral physical properties and layer parameters between adjacent measuring points smoother and more continuous and reduce the limitations of a single geophysical inversion method，a pseudo-two-dimensional lateral constrained joint inversion study of controlled source audio frequency magnetotelluric method (CSAMT) and micro-motion spectrum ratio method was carried out. The microtremor data is numerically simulated using the spectral ratio method，combined with the CSAMT-based limited memory BFGS (L-BFGS) inversion algorithm，introducing the lateral constraint theory，and adding the cross-gradient function to achieve the mutual coupling of two different physical parameters. A set of quasi-two-dimensional lateral constraint joint inversion algorithms was developed，and the accuracy and effectiveness of the algorithms were verified through two sets of theoretical models. Meanwhile，the inversion algorithm is used to invert the measured data in Yanqing，Beijing. The results show that there is a good correspondence between the abrupt interface morphology of resistivity and shear wave velocity，which proves the practical value of the laterally constrained joint inversion algorithm.