MATLAB带通滤波器开始端和结尾端数据异常（解决的小技巧）

在前一段时间遇到了一个MATLAB带通滤波方面的问题，滤波后首尾段数据都出现异常状况，然后就开始疯狂的查资料，但是这方面的资料说实话确实很少，尤其是这种滤波器性能导致的问题，后来偶然的巧合下，在MATLAB中文论坛里面看到了一片因为时延导致数据异常的情况，就借鉴他那个方法成功的解决了我这个问题，总体滤完波后没有什么明显的数据异常问题
先看一下我以前的异常数据情况：

添加了异常处理部分后：

很明显异常部分处理掉了，我的异常处理思路就是，既然出错部分是由于滤波器性能决定的，那么计算出滤波器的延迟是多少，那么在数据的开始部分和结尾部分添加一部分正常数据，正好是延迟部分就可以，虽然会有一点影响，但是影响不会很大，然后滤完波后，再将异常数据去掉，就可以了
下面是我的部分异常处理代码：

fs = 1000;
fp1 = 1; fs1 = 40;%过渡带范围，单位KHz
Wp1 = fp1*pi/(fs/2);Ws1 = fs1*pi/(fs/2);
Wdelta = Ws1-Wp1;%计算过渡带宽度
N = ceil(6.2*pi/Wdelta);
if(rem(N,2) == 1)
    N = N-1;
end
Wn = (Ws1+Wp1)/2/pi;
window = hanning(N+1);
b = fir1(N,Wn,window);
delay=length(b)/2-1/2;%计算出延迟
addData = data1(length(data1)-delay + 1:end);
input=[addData input addData];
[wfdata_fd delay1]=filter(b,1,input);
output=wfdata_fd(length(delay1)/2+1:end);%去掉前半部分异常数据
fzdata=fliplr(output);%数据翻转
output2=fzdata(length(delay1)/2+1:end);%去掉后半部分异常数据
fzdata1=fliplr(output2);%再次翻转成正常数据
EEGData = fzdata1;