[数字信号处理]IIR高通滤波器加速收敛算法实现

原创于 2019-01-03 22:19:52 发布 · 1.1k 阅读

4 ·

CC 4.0 BY-SA版权

文章标签：

#IIR #滤波器 #DSP #加速收敛

数字信号处理专栏收录该内容

9 篇文章

订阅专栏

本文介绍了一种在Matlab中实现高通滤波器的加速收敛方法，对比了传统Matlab方法与加速方法的效果。通过设计滤波器系数并应用到信号上，加速方法能在保持滤波效果的同时，提高处理速度。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

clc;clear all;close all;

% 原信号
sample_rate = 100;
n           = 5000;
x           = 1:1:n;
data        = 50000 + 1000*sin((x/sample_rate)*2*pi);

% matlab 慢收敛方法高通滤波
[b,a] = butter(2,0.33/50,'high');
data_fil_matlab = filter(b,a,data);

% 加速收敛方法
y_buff  = zeros(1,3);x_buff  = zeros(1,3);
data_fil_accelerate = zeros(n,1);

for data_index = 1:1:n
     %滤波器
    if(data_index==1)
        x_buff(1) = data(data_index); 
        x_buff(2) = data(data_index);
        x_buff(3) = data(data_index);
    end
   
     %IIR HP*******************************************************************
	y_buff(1) = y_buff(2);y_buff(2) = y_buff(3);
	x_buff(1) = x_buff(2);x_buff(2) = x_buff(3);
   
    %赋值
	x_buff(3)  = data(data_index);
   
    y_buff(3) = (-1)*a(2)*y_buff(2)+(-1)*a(3)*y_buff(1)+...
					b(1)*x_buff(3)+b(2)*x_buff(2)+b(3)*x_buff(1);
   
    data_fil_accelerate(data_index) = y_buff(3);
end

figure(1);
subplot(3,1,1);
plot(data);
title('Original Data');grid on;
subplot(3,1,2);
plot(data_fil_matlab);
title('Matlab Method Data');grid on;
subplot(3,1,3);
plot(data_fil_accelerate);
title('Accelerate Method Data');grid on;