Tenemos
$$ F(x,y)= \frac{Q}{x} \frac{\partial P}{\partial x} - \frac{P}{x} \frac{\partial Q}{\partial x} - \frac{Q}{y} \frac{\partial P}{\partial y} + \frac{P}{y} \frac{\partial Q}{\partial y} $$
donde
$$P = \sum_{i=1}^{N}(a_i x + b_i y)^2 $$ $$Q = \sum_{i=1}^{N}(c_{i} x + d_{i} y)^2$$
y $a_i, b_i, c_{i}, d_{i}$ son constantes los parámetros definidos para $1\leq i,j \leq N$.
Queremos encontrar las soluciones de la ecuación de $F(x,y)=0$. Para un caso más sencillo, donde $P$ se define como el anterior y $Q=1$, mediante un cambio de variable $m=\frac{y}{x}$, obtenemos una ecuación cuadrática para $m$ que podría resolverse simplemente y dar las dos soluciones para $m$. Me pregunto si con el mismo cambio de variable, se podría obtener las soluciones para el caso más general (con $Q$ de la forma indicada más arriba o similares).
Algunos antecedentes: $x,y$ están las coordenadas de un punto en un espacio de dos dimensiones y definir una línea que pasa a través de $(x,y)$ y el origen $(0,0)$. $P$ y $Q$ son derivados de la proyección de algunos otros de los puntos en esta línea. El uso de $m=y/x$ tiene sentido, ya que las proyecciones sólo depende de la pendiente de la línea y no los valores reales de a $(x,y)$. Nuestro objetivo es encontrar la línea (que se define por su pendiente $m$) que satisface $F(x,y)=0$.
EDIT: Este problema es el caso especial de un problema más general para encontrar la asignación de un $N$ espacio tridimensional a una $M$ espacio tridimensional que minimice la suma de las distancias de los pares de puntos similares ($P$) normalizado a la suma de las distancias de los pares de distintos puntos de ($Q$). La solución que se dio aquí lo soluciona para la asignación de 2D 1D. También he publicado una nueva pregunta para el caso de asignación de 3D a 1D aquí, que espero que se pueda resolver de una manera similar.