General orthogonal polynomials are employed to perform the analysis, parameter identification and optimal control of linear time-varying systems. Applying the derived operational matrix of integration and operational matrix of time multiplication, the differential equation is simplified to a set of algebraic equations which is more suitable for digital computation. Illustrative examples are included for demonstration.