Higher-order Hamiltonians for the trigonometric Gaudin model