The purpose of this work is to study a class of singular elliptic system involving the p(x)p(x)-Laplace operator of the form
{−Δp(x)u=λa(x)|u|q(x)−2u(x)+1−α(x)2−α(x)−β(x)c(x)|u|−α(x)|v|1−β(x) in Ω,−Δp(x)v=μb(x)|v|q(x)−2v(x)+1−β(x)2−α(x)−β(x)c(x)|u|1−α(x)|v|−β(x) in Ω,u=v=0, on ∂Ω,
where Ω⊂ℝN,(N≥2) is a bounded domain with C2 boundary, λ,μ are two parameters, a,b,c∈C(¯Ω) are non-negative weight functions with compact support in Ω and p,q,α,β∈C(¯Ω) are assumed to satisfy the assumptions (A0)–(A2) in Sec. 1. We employ the Nehari manifold approach combined with some variational techniques in order to show the existence and the multiplicity of positive solutions.