  We present a hybridized discontinuous Galerkin (HDG) method for stationary
linearized incompressible magnetohydrodynamics (MHD) equations. At the heart of
the paper is the introduction of an HDG flux of the dual saddle-point form of
the MHD equations that facilitates the hybridization of discontinuous Galerkin
(DG) method. We carry out the $\textit{a priori}$ error estimates for the
proposed HDG method on simplicial meshes in both two- and three-dimensions. The
analysis provides optimal convergence for the fluid velocity and the magnetic
variables, and quasi-optimal convergence for the remaining quantities.
Numerical examples are presented to verify the theoretical findings.
