Bandpass fundamental and harmonics in Matlab
You have a problem with the comb filter design. According to Matlab documentation the following line creates the filter you want:
[b,a] = iircomb(round(25600/285), 2*285/25600/35,'peak');

Looking at the frequency response of the designed filter: 
Seems to be able to solve your problem...