clear;clc;
close all
sf=50000*1/0.0006;%采样频率
fmin=80000;%最小截止频率
fmax=300000;%最大截止频率
load('D:\课程学习\教研室\盲源分离\DEMO_BSS\50noise\s1.dat');%信号数据
x=s1;
n=length(x);
t=(0:1/sf:(n-1)/sf)';
nfft=2^nextpow2(n);%取大于最接近n的2的幂次方为FFT长度
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('时间');title('滤波前');
subplot(2,1,2);plot(t,y);xlabel('时间');title('滤波后');