There are different equivalent formulations of quantum mechanics, and how an interaction is represented in each formulation is different.
In the Schrodinger picture, operators are time independent and the state (or wavefunction) $\Psi$ depends on time. The Schrodinger equation governs the dynamics of the wavefunction:
\begin{equation}
i \hbar \frac{\partial \Psi}{\partial t} = H \Psi
\end{equation}
Note that the Schrodinger equation is always linear in the wave function $\Psi$, even in more complicated situations such as when you generalize to multi-particle states, or field theory. Interactions appear via terms in the Hamiltonian operator $H$ which are higher-than-quadratic-order in the position operator $x$ and momentum operator $p$.
For example, consider the following Hamiltonian:
\begin{equation}
H = \frac{p^2}{2 m} + \frac{m\omega^2}{2}x^2 + \lambda x^3
\end{equation}
- If $\lambda = 0$, then the Hamiltonian is the same as a harmonic oscillator. When we generalize this case to field theory, we would say there are no interactions. Note that the classical equations of motion following from this Hamiltonian (the equations of motion for a harmonic oscillator) are linear in $x$. We can solve the Schrodinger equation exactly in this case.
- If $\lambda \neq 0$, then when generalizing to field theory we would say there is a cubic interaction. Note that the classical equations of motion are not linear. We cannot usually solve the Schrodinger equation in cases like this. However note that the Schrodinger equation is still linear in the state $\Psi$.
In the Heisenberg picture, operators are time dependent and the state $\Psi$ is time-independent. For an operator $\mathcal{O}$ without explicit time dependence (such as the position or momentum operators), the equations of motion
\begin{equation}
\frac{{\rm d}\mathcal{O}}{{\rm d}t} = \frac{i}{\hbar}[H,\mathcal{O}]
\end{equation}
If we take our example above, and use the fact that $[x,p]=i \hbar$, then we would find the position and momentum have the Heisenberg equations of motion
\begin{eqnarray}
\frac{{\rm d} x}{{\rm d}t} &=& \frac{i}{\hbar}[H,x] = \frac{1}{m} p \\
\frac{{\rm d} p}{{\rm d}t} &=& \frac{i}{\hbar}[H,p] = - m \omega^2 x - 3 \lambda x^2
\end{eqnarray}
This system of equations probably more closely matches your intuition. When $\lambda=0$, the Heisenberg equations of motion are linear (in the position and momentum operator), and non-interacting. When $\lambda \neq 0$, the Heisenberg equations of motion are not linear and are interacting.
Finally I'll mention the Feynman path integral approach. There the basic quantity is the path integral
\begin{equation}
Z = \int \mathcal{D} x e^{\frac{i}{\hbar} S[x]}
\end{equation}
where $S$ is the classical action. The measure $\mathcal{D} x$ is constructed via a limiting process and represents an integral over all classical paths.
In our example, the action is
\begin{equation}
S = \int {\rm d} t \left( \frac{m}{2} \dot{x}^2 - \frac{m}{2} \omega^2 x^2 - \lambda x^3 \right)
\end{equation}
When $\lambda=0$ (free theory), the path integral is Gaussian and can be done exactly. When $\lambda \neq 0$ (interacting theory), the path integral is not Gaussian and usually can't be done exactly. Usually at this point a physicist would expand the exponential in powers of $\lambda$, assuming it is small, and use perturbation theory and Feynman diagrams.
This is not an exhaustive list of all possible formulations of quantum mechanics, but these are the most common, and at least should give you the idea that a physical idea (like whether a particle has interactions) can be mathematically represented in different ways. "A particular equation being linear or not linear" is a statement about the formalism and can depend on how you formulate the problem, whereas "a particle having interactions or not" is a statement about physics and cannot.
