A fairly straightforward approach, that generalizes to $${\raise{-1ex}\mathop{\huge\mathrm{K}}_{n=0}^\infty}\frac{a_2 n^2+a_1 n+a_0}{b_1 n+b_0}$$ (using Gauss' notation), is to consider exponential generating functions of the convergents.
The value of the above CF (if it exists) is $\lim\limits_{n\to\infty}(P_n/Q_n)$, where the vectors $R_n=(P_n,Q_n)$ satisfy $$R_0=(1,0),\quad R_1=(0,1),\quad R_{n+2}=(b_1 n+b_0)R_{n+1}+(a_2 n^2+a_1 n+a_0)R_n,\quad(n\geqslant 0)$$ and we find that the EGF $R(z)=\sum_{n=0}^\infty R_n z^n/n!$ satisfies $$(1-b_1 z-a_2 z^2)R''(z)-\big(b_0+(a_1+a_2)z\big)R'(z)-a_0 R(z)=0,$$ an ODE of hypergeometric $_2F_1$ type (after a linear change of variable).
In our case, $b_0=2x$, $b_1=0$, $a_0=1$, and $a_1=a_2=4$; we suppose $x>0$ is real (for simplicity). Then $$R(z)=C_1\times{}_2F_1\left(\frac12,\frac12;1+\frac{x}{2};\frac12-z\right)+C_2\times\ldots,$$ where "$\ldots$" stands for the second solution of the ODE, which is unbounded at $z\to 1/2$, and $C_1,C_2$ are some constants (which of course depend on $x$, and determined by $R(0)=R_0$ and $R'(0)=R_1$).
Writing $R(z)=\big(P(z),Q(z)\big)$ [recall that $R_n$ are vectors] and assuming that our CF converges (in fact it's not hard to prove) to $f(x)$, we know that $$f(x)=\lim_{z\to 1/2^-}\frac{P(z)}{Q(z)}.$$
This is equal to $A/B$, where $(A,B)=C_2$ is our second "constant". Computing it, we find $$f(x)=\frac{F'(1/2)}{F(1/2)},\qquad F(z)={}_2F_1\left(\frac12,\frac12;1+\frac{x}{2};z\right).$$
The value of $F(1/2)=\Gamma\big((x+1)/4\big)/\Gamma\big((x+3)/4\big)$ can be computed using the integral representation of $_2F_1$ (the substitution $t=x(2-x)$ helps; there is also a known formula). The value of $F'(1/2)$ is obtained similarly, after integration by parts.
Finally, we obtain the CF $(1)$ in the OP as $4/\big(x+f(x)\big)$. This gives the expected result.
An equivalent (and a bit easier) approach is to use the corresponding recurrence for $_2F_1$ instead of the ODE.