函数y =溢出(b, x,重置)如果输入参数个数< 3,重置= true;结束持续的z p如果isempty (z) | |重置p = 0;z = 0(大小(b));结束(y, z, p) = fir_filter (b, x, z, p);结束函数(y, z, p) = fir_filter (b, x, z, p) y = 0(大小(x));nx =长度(x);nb =长度(b);为n = 1: nx p = p + 1;如果p > nb, p = 1;结束z (p) = x (n);acc = 0;k = p;为j = 1: nb acc = acc + b (j) * z (k);k = k - 1;如果k < 1, k =注;结束结束y (n) = acc;结束结束

函数overflow_test%计算滤波器系数%使用FIR1函数%信号处理工具箱。0.25 b % = fir1(11日);b = (-0.004465461051254 - -0.004324228005260 + 0.012676739550326 + 0.074351188907780 + 0.172173206073645 + 0.249588554524763 + 0.249588554524763 + 0.172173206073645 + 0.074351188907780 + 0.012676739550326 -0.004324228005260 -0.004465461051254) ';%输入信号nx = 256;t = linspace(0, 10 *π,nx) ';%的冲动x_impulse = 0 (nx, 1);x_impulse (1) = 1;%最大增益%时,将发生的最大增益滤波器%的输入线滤波器的迹象%脉冲响应。x_max_gain = (b)的迹象;x_max_gain = repmat (x_max_gain装天花板(nx /长度(b)), 1);x_max_gain = x_max_gain (1: nx);%的正弦f0 = 0.1;f1 = 2;x_sines =罪(2 *π* t * f0) + 0.1 *罪(2 *π* t * f1);%唧唧声f_chirp = 1/16;%的目标频率x_chirp =罪(π* f_chirp * t ^ 2);%线性啁啾x = [x_impulse x_max_gain、x_sines x_chirp);标题= {“冲动”,“马克斯增益”,正弦的总和,“唧唧喳喳”};y = 0(大小(x));为i = 1:尺寸(x, 2)重置= true;y (:, i) =溢出(b, x (:, i),重置);结束test_plot(1、标题、t, x, y)结束函数test_plot(图、标题、t, x, y₁)图(图)clf sub_plot = 1;font_size = 10;为i = 1:尺寸(x, 2)次要情节(4 1 sub_plot) sub_plot = sub_plot + 1;情节(t) x(:,我)“c”、t、y₁(:,我)“k”)轴(“紧”)包含(“t”,“字形大小”,font_size);标题(标题{},“字形大小”,font_size);甘氨胆酸ax =;斧子。字形大小= 10;结束图(gcf)结束