function [Ex, Ey, Ez, normE] = get_ele_vertices_bele(mesh, solver, phy)
%GET_ELE_VERTICES_BELE  重构散射场并叠加背景场 E_total = E_s + E_b

[Ex, Ey, Ez, ~] = get_ele_vertices(mesh, solver);

V = mesh.Vertex;
for iv = 1:mesh.NbrVertex
    eb = phy.bele.Eb(V(iv, 1), V(iv, 2), V(iv, 3));
    Ex(iv) = Ex(iv) + eb(1);
    Ey(iv) = Ey(iv) + eb(2);
    Ez(iv) = Ez(iv) + eb(3);
end

normE = sqrt(abs(Ex .* conj(Ex) + Ey .* conj(Ey) + Ez .* conj(Ez)));

end