In this dissertation, the discovery potential of neutral Higgs bosons has been discussed respectively in the Standard Model (SM), minimal supersymmetric standard model (MSSM), and unified supersymmetric models. In the SM, we consider the Higgs boson being produced together with one jet, i.e., pp→jh0+X where j=g,q or q¯ (q=u,d,c,s); in supersymmetric theories, we focused on such channel as bg→bϕ (ϕ=h0,H0 or A0). The interesting decay mode is Higgs to tau pair followed with one tau decaying to e± or μ± and the other to mesons (π,a1 or ρ). By comparing with the relevant background processes, we found that the discovery of the SM Higgs boson is pessimistic at s=8 TeV, but it is very promising to search for neutral Higgs bosons in supersymmetric models owing to a large tan⁡β. In addition, we considered the experimental constraints from Br(b→sγ), Br(Bs0→μ+μ−), Δaμ and mh0=125 GeV. In the MSSM, the constraint from light Higgs mass favors the parameter space which has an intermediate tan⁡β and is explorable at the next run of LHC. The minimal supergravity unified model (mSUGRA) has large space allowed by mh0≈125 GeV, but it seems difficult to reconcile the flavor problems. In the model with minimal gauge mediated supersymmetry breaking (mGMSB), the situation is totally opposite.