Chebyshev matrix product states with canonical orthogonalization for spectral functions of many-body systems