taofly 发表于 2010-5-15 22:58

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

我的程序是这样的,但不知道为什么,每次得到的中值频率总有好多点在500Hz,然后其他好多点在50Hz附近,不知道是程序问题还是数据问题,我进行工频滤波了。
=stft(c);%semg列向量
% TF=TF';
TF=abs(TF);
=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;
页: [1]
查看完整版本: 请问中值频率的程序哪有问题