  Rare $b$ decay processes are analyzed in the multi-Higgs doublet model.
Taking account of the constraint from the $b\to s \gamma$ process, the
branching ratio and the forward-backward asymmetry of the final leptons for the
$b\to s \ell^+ \ell^-$ process are calculated. It is shown that the branching
ratio can be a few times larger than the standard model prediction and the
asymmetry can be significantly different from that in the standard model.
Combining these observable quantities it is possible to determine complex
coupling constants associated with the charged Higgs mixing matrix. The CP
violating charge asymmetry in the $b\to s \ell^+ \ell^-$ process and the
branching ratio of the $b\to s \nu\bar{\nu}$ process are also calculated.
