f0=200;
x=load('200_0.txt');
N=length(x);
N2=round(N/2);
%subplot 211; plot(x); grid;
x=x-mean(x);
xf=fft(x);
xf=xf(1:N2)/N*2;
A=abs(xf);
=max(A);
indsecL=A(index-1)>A(index+1);
df=indsecL.*A(index-1)./(Amax+A(index-1))-(1-indsecL).*A(index+1)./(Amax+A(index+1));
ddf=f0/(index-1-df);
fs=N*ddf
其它文件只要改变f0和文件名就可以了,这几个数据文件求出的采样频率为 64710Hz左右。 谢谢大家的热心帮助,今天真是获知良多阿:handshake :handshake :handshake
页:
1
[2]