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

版权声明:欢迎大佬们指点 https://blog.csdn.net/bloke_come/article/details/77994542

在前一段时间遇到了一个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;
展开阅读全文

没有更多推荐了,返回首页