Optimal design for nonlinear estimation of the hemodynamic response function