The near-field resonance response of a plasmonic metamaterial lattice under an oblique incident field is theoretically investigated via applying the first-proposed M-K dynamic matrix method. By considering the electric, magnetic and field-dipole interactions, we construct a dissipative many-body Lagrange model for a reference lattice. A collective forced vibration equation, with the degree of freedom equals to the number of nanoparticles in a cell, is introduced to describe the lattice resonance under a polarized field. The resonance frequencies can be conveniently obtained from the poles of transfer function matrix. Based on this elegant matrix differential equation, one can calculate the dynamic response of plasmonic lattice and analysis the normal modes from dispersion relations. The feasibility of this method is shown in details from three examples: simple square lattice, binary chain and chessboard lattice. Reasonable and significant results show that M-K matrix method is highly universal in a large frequency band and perfectly matched with numerical simulations. It proves to be a powerful tool in treating resonance problem of metallic nanoparticle lattice.