王济那本书上的程序,感觉程序中的加速度的有效值不对呀,数值的级数完全不对,对带通滤波的数据随然进行了傅里叶逆变换,但是没有用,这个对吗?
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
哪位大侠用过,给点介绍呗。