In this paper, a suboptimal state feedback integral decentralized tracking control synthesis for interconnected linear time-variant systems is proposed by using orthogonal polynomials. Particularly, the use of operational matrices allows, by expanding the subsystem input states and outputs over a shifted Legendre polynomial basis, the conversion of time-varying parameter differential state equations to a set of time-independent algebraic ones. Hence, optimal open-loop state and control input coefficients are forwardly determined. These data are used to formulate a least-square problem, allowing the synthesis of decentralized state feedback integral control gains. Closed-loop asymptotic stability LMI conditions are given. The proposed approach effectiveness is proved by solving a nonconstant reference tracking problem for coupled inverted pendulums.