  The parametrizations of experimental yields of K$^{\pm,0}$, $\phi$ and
$\Lambda+\Sigma^0$ are proposed as function of available energy,
$\sqrt{s_\mathrm{NN}}$, and number of participants, $\langle A_\mathrm{part}
\rangle_\mathrm{b}$, for $\sqrt{s_\mathrm{NN}}$ from 2.15 to 3 GeV. For all the
dataset the $\langle A_\mathrm{part} \rangle_\mathrm{b}$ was extracted using
the Glauber Monte Carlo method. The $\alpha$ exponent of yield dependency on
$\langle A_\mathrm{part} \rangle_\mathrm{b}$ appears not to change with beam
energy and is found to be 1.30 $\pm$ 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
$\sqrt{s_\mathrm{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 unpublished yields for STAR's Au+Au collisions at 3 GeV,
and for Au+Au collisions planned by CBM, up to 3.85 GeV.
