In this paper, we study a three-coupled variable-coefficient nonlinear Schrödinger system, crucial for picosecond pulse propagation in inhomogeneous optical fibers and modern communication advancements. Using the Riemann-Hilbert method, we analyze this system based on a \(4 \times 4\) matrix spectral problem from its Lax pair. After constructing relevant analytic functions and solving the Riemann-Hilbert problem, we obtain the N-soliton solutions under non-reflective boundary conditions. Subsequently, we analyze the dynamical behavior of multi-soliton solutions, revealing a series of physically meaningful phenomena like elastic collisions, soliton reflection, and picosecond pulse propagation patterns. Additionally, we investigate how the nonlinear coefficient \(\gamma (z)\) and the group velocity dispersion coefficient \(\beta (z)\) affect soliton solutions, clarifying their role in governing soliton trajectories and amplitude. This research provides essential theoretical support for precisely controlling picosecond pulses in heterogeneous optical fibers. This work enriches the soliton dynamics theory of the aforementioned system, provides critical guidance for multi-component optical fiber systems, and offers further theoretical support for advancing high-performance optical fiber communication technologies.