Let $\gamma, \omega, c$ be positive constants, let $\mathbf{Q}_{a}$ and $\mathbf{Q}_{b}$ be three-dimensional vectors, and let $\mathbf{B}(\mathbf{r})=\mathbf{B}(x,y,z)$ be a vector field. Let $\mathbf{B}_{\circ}= \mathbf{B}(\mathbf{r}_{\circ})$ for some position $\mathbf{r}_{\circ}$ and write $\mathbf{B}=\mathbf{B}_{\circ}+\delta\mathbf{B}$. Consider the following coupled system for the unknown vectors $\mathbf{a}$ and $\mathbf{b}$: $$ -\frac{\omega^2}{\gamma}\mathbf{a}=c\mathbf{B}_{\circ}\times\mathbf{Q}_{a}+\omega\mathbf{b}\times\mathbf{B}_{\circ} \tag{1} $$ $$ -\frac{\omega^2}{\gamma}\mathbf{b}=c\mathbf{B}_{\circ}\times\mathbf{Q}_{b}-\omega\mathbf{a}\times\mathbf{B}_{\circ} \tag{2} $$ If it helps, $\mathbf{Q}_{a}=\int_{0}^{2\pi/\omega}\mathbf{B}\cos(\omega t)\,dt$ and $\mathbf{Q}_{b}=\int_{0}^{2\pi/\omega}\mathbf{B}\sin(\omega t)\,dt$. All the $x, y, z$ are functions of $t$ but $\mathbf{B}$ is (basically) abitrary so those integrals can't be done.
Is it possible to solve $(1)$ and $(2)$ for $\mathbf{a}$ and $\mathbf{b}$ while still keeping everything a vector? I believe that one could simply write out all the components of $(1)$ and $(2)$ to get two linear relationships between $\mathbf{a}$ and $\mathbf{b}$, which is solvable in principle, but I'd really rather not do that. I've tried just subbing one equation into the other, but you end up with the thing to solve for being buried in a double cross product, so I'm not sure.
Hint: as double cross product can be computationally exhausting (a bit), why don't you have a preliminary look at the vector components? If we forget about constants for a second putting $\omega=\gamma=c=1$, then we obtain
$$a=(Q_a-b)\times B, $$ $$b=(Q_b+a)\times B. $$
As $a_i=\epsilon_{ilr}(Q_a-b)_{lr}B_r$, then $b_j=\epsilon_{jik}(Q_b)_iB_k+ \epsilon_{jik}\epsilon_{ilr}(Q_a-b)_{lr}B_r,$ which can be manipulated with
$$\epsilon_{jik}\epsilon_{ilr}=-\delta_{jl}\delta_{kr}+\delta_{jr}\delta_{lk}. $$