Harmonic amplification phenomena could appear at the point of common connection (PCC) of the cable line terminal. However, the distributed parameter model of the cable line contains hyperbolic functions with plural variables, which makes it challenging to obtain the harmonic amplification factor (HAF). Hence, a time-domain method combining the Kalman filter and convolution inversion (KFCI) methods is proposed to address this problem. First, the Kalman filter method optimizes the square wave pulse response (SWPR) with measurement error. Then, the optimized SWPR data are used to get the HAF by the convolution inversion method. Next, the harmonic amplification characteristics of cable lines are explored. Finally, an experimental simulation model is built on the PSCAD software, verifying the optimization effectiveness of Kalman filter for the SWPR with error and the accuracy of the HAF calculated by the proposed method. The analysis rationality of harmonic amplification properties is also demonstrated.