摘要： The finite element method (FEM) has been one of the major numerical techniques for electromagnetics (EM) ever since the introduction of the spurious-free edge element. On the other hand, a so-called least-squares finite element method (LSFEM) is becoming increasingly popular due to its unique properties. For instance, the generated stiffness matrix is positive definite. Effort has been made to apply this node-element-based LSFEM to model EM problems, mostly in the theoretical aspect. In this paper, we will use normalized Maxwell’s equations incorporated with the constitutive relations of bianisotropic media, which represents the most general cases in the linear regime. Subsequently, the corresponding variational formulation is derived, followed by the generation of matrix systems. We then implement and validate various boundary conditions. Further we apply this scheme to simulate wave scattering in the presence of bianisotropic media that are categorized into 3 types. A good convergence rate and agreement between the LSFEM results and the referenced solutions has been achieved, which proves the Lagrange-element-based LSFEM as a feasible alternative to the edge-element-based FEM currently used in EM modeling.