请问中值频率的程序哪有问题

taofly

我的程序是这样的，但不知道为什么，每次得到的中值频率总有好多点在500Hz，然后其他好多点在50Hz附近，不知道是程序问题还是数据问题，我进行工频滤波了。
[F,T,TF]=stft(c);%semg列向量
% TF=TF';
TF=abs(TF);
[TFrow,TFcol]=size(TF);
for (i=1:TFrow)
    medianf(i)=fft_medianf(TF(i,(TFcol/2+1):TFcol));
    end；
其中的调用函数fft_medianf是
function medianf=fft_medianf(x)
fs=1000;
N=length(x);
d=sum(x)/4;
sum1=x(1);
for (t=1:1:N)
if(sum1<d)
    t=t+1;
    sum1=sum1+x(t);
    if(sum1>d)  break;
    end;
end;
end;
% medianf=(t*fs)/N;
medianf=t*fs/512;