Abstract
We use the sampling representations associated with Sturm-Liouville difference operators to derive generalized integral-valued trigonometric sums. This extends the known results where zeros of Chebyshev polynomials of the first kind are involved to the use of the eigenvalues of difference operators, which leads to new identities. In these identities Bernoulli's numbers play a role similar to that of Euler's in the old ones. Our technique differs from that of Byrne-Smith (1997) and Berndt-Yeap (2002).