基于频率方法的高通滤波故障诊断程序实现

%%%%%%%%%%%%%%%%%%%%adaptive filter1
clear
load O100W2D5D0D3A.csv
X=O100W2D5D0D3A(:,1);
x=X(512*100+1:512*200);
%%highpass filter with frequency method 3000 hz
fs=200000;%frequency sampling
freqmin=0;
freqmax=3000;
n=length(x);
t=(0:1/fs:(n-1)/fs)';
nfft=2^nextpow2(n);
ni=round(freqmin*nfft/fs+1);
na=round(freqmax*nfft/fs+1);
y=fft(x,nfft);
a=zeros(1,nfft);
y(ni:na)=a(ni:na);
y(nfft-na+1:nfft-ni+1)=a(nfft-na+1:nfft-ni+1);
a=ifft(y,nfft);
y=real(a(1:n));
figure(1)
subplot(3,1,1);plot(t,x);title('原始信号');
subplot(3,1,2);plot(t,y);title('高通滤波信号');
%%
F1=abs(hilbert(y));
k=2^nextpow2(length(F1));% get the value :power of two more than length(F1)
F1=F1-(ones(1,length(F1)))'*mean(F1);
fs=200000;%sampling frequency
f=fs*(0:((k/2)-1))/k;
Y2=fft(F1,k);

Pyy2=Y2.*conj(Y2)/k;
kk=1024;f1=f(1:kk);
subplot(3,1,3);plot(f1,Pyy2(1:kk));title('求包络后所得频谱' );xlabel('frequency (Hz)');ylabel('Signal Value');

--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
%%%  程序――平均求绝对值――功率谱法

clear
load O150W5D0D0D5A.csv
X=O150W5D0D0D5A(1:200000,1);
%%highpass filter with frequency method 3000 hz
fs=200000;%frequency sampling
freqmin=0;
freqmax=3000;
n1=length(X);
t1=(0:1/fs:(n1-1)/fs)';
nfft=2^nextpow2(n1);
ni=round(freqmin*nfft/fs+1);
na=round(freqmax*nfft/fs+1);
y=fft(X,nfft);
a=zeros(1,nfft);
y(ni:na)=a(ni:na);
y(nfft-na+1:nfft-ni+1)=a(nfft-na+1:nfft-ni+1);
a=ifft(y,nfft);
y=real(a(1:n1));
figure(1)
subplot(3,1,1);plot(t1,X);title('原始信号');
subplot(3,1,2);plot(t1,y);title('高通滤波信号');