王济那本书上的程序，感觉程序中的加速度的有效值不对呀，数值的级数完全不对，对带通滤波的数据随然进行了傅里叶逆变换，但是没有用，这个对吗？
a=fft(x,nfft);
for j=1:nc
fl=fc(j)/oc6;
fu=fc(j)*oc6;
nl=round(fl*nfft/sf+1);
nu=round(fu*nfft/sf+1);
if fu>sf/2
m=j-1;break
end
b=zeros(1,nfft);
b(nl:nu)=a(nl:nu);
b(nfft-nu+1:nfft-nl+1)=a(nfft-nu+1:nfft-nl+1);
c=ifft(b,nfft);
yc(j)=sqrt(var(real(b(1:n))));
end
哪位大侠用过，给点介绍呗。
[复制链接]
发表于 2010-12-22 20:39
提升卡
变色卡
显身卡
发表于 2010-12-24 14:32
