本文考虑如下$ p $-Laplacian 的非齐次拟线性 Kirchhoff-Schr\"odinger-Poisson~系统解的存在性:
\begin{equation*}
\begin{cases}
- \Big(a-b \displaystyle{\int_{\mathbb{R}^3}}|\nabla u|^p{\rm d}x \Big)\Delta_p u+|u|^{p-2}u+\lambda\phi_uu= |u|^{q-2}u+h(x), &\text{ }x\in \mathbb{R}^3,\\
-\Delta\phi=u^2, &\text{ }x\in \mathbb{R}^3,\\
\end{cases}
\end{equation*}
其中 $ a,b>0 $, $ \frac{4}{3} < p < \frac{12}{5} $, $ p < q < p^* =\frac{3p}{3-p} $, $ \lambda > 0 $. 当 $ h(x) $ 满足适当的条件时, 运用Ekeland变分原理和山路定理给出系统多重解的存在性.