I have noticed a certain similarity between Riccati 1st Order ODEs and linear 1st Order ODEs. Specifically, the general solution for each is given by any particular solution plus some function of the independent variable.
In the linear case, this function is a scalar multiple of any particular solution to the associated homogeneous ODE, so the answer takes the form $y = y_p + Cy_h$, where $y$ is the general solution, $y_p$ is any particular solution, $C$ is an arbitrary constant, and $y_h$ is any particular solution to the homogeneous ODE. This is much the same as the observation that the solution to a linear system of equations is given by a particular solution plus the nullspace.
The answer to a Riccati ODE takes the form $y = y_p + u$, where $u$ is some function of the independent variable. Given this similarity, can we think of $u$ as the zero set (the generalization of the nullspace into nonlinear contexts), or as in some way highly related to $Cy_h$?
If so, is there an analog to the linear ODE fact that $y_p + Cy_h = y_h \int \frac{f}{y_h} dx$, where $f$ is the non-homogenous term of the linear ODE, which could possibly give rise to another solution method for the Riccati?