我做推广找不到我的网站,自学网站平面设计,天津建设工程,铁岭手机网站建设补偿 FIR 滤波器引入的延迟
对信号进行滤波会引入延迟。这意味着相对于输入#xff0c;输出信号在时间上有所偏移。此示例向您说明如何抵消这种影响。
有限冲激响应滤波器经常将所有频率分量延迟相同的时间量。这样#xff0c;我们就很容易通过对信号进行时移处理来针对延迟…补偿 FIR 滤波器引入的延迟
对信号进行滤波会引入延迟。这意味着相对于输入输出信号在时间上有所偏移。此示例向您说明如何抵消这种影响。
有限冲激响应滤波器经常将所有频率分量延迟相同的时间量。这样我们就很容易通过对信号进行时移处理来针对延迟进行校正。
以 500 Hz 的频率对心电图读数采样采样时间为 1 秒。添加随机噪声。重置随机数生成器以获得可再现性。
Fs 500;
N 500;
rng defaultxn ecg(N)0.25*randn([1 N]);
tn (0:N-1)/Fs;
使用滤波器阻挡 75 Hz 以上的频率以消除一部分噪声。使用 designfilt 设计一个阶数为 70 的滤波器。
nfilt 70;
Fst 75;d designfilt(lowpassfir,FilterOrder,nfilt, ...CutoffFrequency,Fst,SampleRate,Fs);
对信号进行滤波并绘图。与原始信号相比结果更平滑但存在滞后。
xf filter(d,xn);plot(tn,xn)
hold on, plot(tn,xf,-r,linewidth,1.5), hold off
title Electrocardiogram
xlabel Time (s), legend(Original Signal,Filtered Signal) 使用 grpdelay 检查滤波器造成的延迟是否等于滤波器阶数的一半。
grpdelay(d,N,Fs) delay mean(grpdelay(d))
delay 35对滤波后的信号进行时移以对齐数据。删除信号的前 delay 个采样。删除原始采样和时间向量的最后 delay 个采样。
tt tn(1:end-delay);
sn xn(1:end-delay);sf xf;
sf(1:delay) [];
对信号绘图并验证它们是否对齐。
plot(tt,sn)
hold on, plot(tt,sf,-r,linewidth,1.5), hold off
title Electrocardiogram
xlabel(Time (s)), legend(Original Signal,Filtered Shifted Signal) 另请参阅
designfilt | filter | filtfilt | grpdelay
相关主题
补偿 IIR 滤波器引入的延迟数字滤波实践介绍
