The parametrizations of experimental yields of K ± , 0 , ϕ and Λ + Σ 0 are proposed as function of available energy, s NN , and number of participants, ⟨ A part ⟩ b , for s NN from 2.15 to 3 GeV. For all the dataset the ⟨ A part ⟩ b was extracted using the Glauber Monte Carlo method. The α exponent of yield dependency on ⟨ A part ⟩ b appears not to change with beam energy and is found to be 1.30 ± 0.02. Our parametrization and the predictions of public versions of RQMD.RMF, SMASH and UrQMD transport models are compared to the HADES experimental data for Ar+KCl at s NN of 2.61 GeV. The phenomenological parametrization currently offers the best overall description of these yields. Predictions are given for yields from Ag + Ag collisions at available energies of 2.41 and 2.55 GeV, analysed by HADES, Au + Au experiment at 2.16 and 2.24 GeV planned by this collaboration, some yields for STAR’s Au + Au collisions at 3 GeV, and for Au + Au collisions planned by CBM, up to 3.85 GeV.