I hate to post answers so soon after asking, but I think I solved the above problem, so I don't want to waste anyone's time.
Polynomial approximation of $H$. Let $x_0$ be a root of $x \to \det H(x)$ of multiplicity $1 \leq \nu$. Since $H(x)$ is bounded, we thus see that the smallest eigenvalue of $H(x)$ must be at least $\epsilon|x-x_0|^\nu$ for some $\epsilon>0$. Consider the power series of $H(x)$, and the Taylor polynomial $P(x)$ of order $m = \max\{\nu,d\}$: $$ H(x) = \sum_{k=0}^{\infty} H_k (x-x_0)^k \text{ and } P(x) = \sum_{k=0}^{m} H_k (x-x_0)^k. $$ Hopefully, $x \to \det P(x)$ also has a root of order $\nu$ at $x_0$ (this hopefully follows from the fact that $\det X$ is a polynomial of degree $d$ in the entries of the matrix $X$, so the first $m\geq \nu$ terms of $\det P(x)$ coincide with the first $m$ terms of the Taylor series of $\det H(x)$ near $x=x_0$).
$P$ is an "optimal preconditioner" for $H$. We also have the estimate $$ \|H-P\| \leq C|x-x_0|^{m+1}. $$ in a neighborhood of $x_0$. Thus, in a punctured neighborhood of $x_0$ and nonzero vector $v$, $$ \frac{v^TPv}{v^THv} = 1-\frac{v^T(H-P)v}{v^THv} = 1-R, $$ $$ |R| \leq {C \over \epsilon}|x-x_0|^{1+m-\nu} \leq 0.5. $$ This leads to our key estimate $$ 0.5v^TPv \leq v^THv \leq 1.5v^TPv, $$ in some neighborhood of $x_0$, say $x \in (a,b)$ and for arbitrary vector $v$.
The main result. Let $q(x)$ be an arbitrary polynomial of degree $n$. Note that $q^TPq$ is a polynomial of degree $m+2n$. From Rahman and Schmeisser, (13.3.2), which apparently is usable for the "$p$-norm" with $p=1/2$, we find that \begin{align} \sup_{(a,b)} q^THq & \leq 1.5\sup_{(a,b)} q^TPq \\\\ & \leq 1.5(2n+m+1)^4 \left({1 \over b-a}\int_a^b \sqrt{q^TPq}\right)^2 \\\\ & \leq 3(2n+m+1)^4 \left({1 \over b-a}\int_a^b \sqrt{q^THq}\right)^2. \end{align} The result follows by partitioning $[0,1]$ into finitely many intervals, and each interval $(a_k,b_k)$ is either a suitably small neighborhood of the finitely many roots of $\det H(x)$, or $H(x)$ is uniformly bounded below and above on $(a_k,b_k)$.
Commentary. The above was a simplified version of my actual problem so I'm not sure my solution will be usable. It turns out that my function $H(x)$ is also parametrized by some other variable, i.e. $H(x,y)$, but the supremum and integral are w.r.t. $x$. For that, it feels like some matrix version of the Weierstrass preparation theorem would be needed, but I'm currently stuck. I guess I've made the additional assumption that $\det H(x)$ is not uniformly $0$.
Edit. Hm. I just discovered this paper on matrix-valued preparation theorems. MaybyMaybe that is what I need.