振动数据分析性能比较多

fni=input('频域带通滤波-输入数据文件名：','s');
频域带通滤波-输入数据文件名：4444.txt
>> fid=fopen(fni,'r');
sf=fscanf(fid,'%f',1);
fmin=fscanf(fid,'%f',1);
fmax=fscanf(fid,'%f',1);
>> sx=fscanf(fid,'%s',1);
>> sy=fscanf(fid,'%s',1);
>> fno=fscanf(fid,'%s',1)
x=fscanf(fid,'%f',[1,inf]);
status=fclose(fid);
n=length(x);
>> t=(0:1/sf:(n-1)/sf)';
>> nfft=2^nextpow2(n);
ni=round(fmin*nfft/sf+1);
na=round(fmax*nfft/sf+1);
y=fft(x,nfft);
>> a=zeros(1,nfft);
>> a(ni:na)=y(ni:na);
>> a(nfft-na+1:nfft-ni+1)=y(nfft-na+1:nfft-ni+1);
>> y=ifft(a,nfft);
>> y=(real(y(1:n)))';
>> subplot(2,1,1);
>> plot(t,x);
>> xlabel('时间')
ylabel('加速度')
grid on
>> subplot(2,1,2);
plot(t,y);
xlabel('时间')