注意FIR滤波器获得结果的采样点数偏移

eric2013

14个点的采样点偏移，高通：

1. fs=1000;                  %设置采样频率 1K
2. N=320;                    %采样点数      
3. n=0:N-1;
4. t=n/fs;                     %时间序列
5. f=n*fs/N;                  %频率序列
7. x1=sin(2*pi*50*t);
8. x2=sin(2*pi*200*t);
9. x=sin(2*pi*50*t)+sin(2*pi*200*t);  %50Hz和200Hz正弦波混合               
11. plot(n, x2, 'b');   %绘制信号x的波形  
12. xlabel('时间');
13. ylabel('幅值');
14. title('原始信号');
15. hold on;
16.   
17. b=fir1(28, 125/500,  'high');     %获得滤波器系数，截止频率125Hz，高通滤波。
18. y=filter(b, 1, x);
19. plot(n(1:305), y(15:319), 'r');
20. legend('原始信号','滤波后信号');
21. grid on;

复制代码

eric2013

14个采样点的偏移，低通

1. fs=1000;                  %设置采样频率 1K
2. N=320;                    %采样点数      
3. n=0:N-1;
4. t=n/fs;                     %时间序列
5. f=n*fs/N;                  %频率序列
7. x1=sin(2*pi*50*t);
8. x2=sin(2*pi*200*t);
9. x=sin(2*pi*50*t)+sin(2*pi*200*t);  %50Hz和200Hz正弦波混合               
11. plot(n, x1, 'b');   %绘制信号x的波形  
12. xlabel('时间');
13. ylabel('幅值');
14. title('原始信号');
15. hold on;
16.   
17. b=fir1(28, 125/500);     %获得滤波器系数，截止频率125Hz.
18. y=filter(b, 1, x);
19. plot(n(1:305), y(15:319), 'r');
20. legend('原始信号','滤波后信号');
21. grid on;

复制代码

xiaosir

学习....

caicaptain2

eric2013:14个采样点的偏移，低通
fs=1000;                  %设置采样频率 1K
N=320;                & .. (2014-12-21 14:12) 

什么意思？

是说经过数字滤波后的波形发生了相移吗？

eric2013

caicaptain2:什么意思？

是说经过数字滤波后的波形发生了相移吗？
 (2015-03-12 16:45) 

是的，线性相位FIR的群延迟。

fcdxsdxsza

管理員你好:

看到這篇   忽然讓我想到  示波器程式中  在   DSO_FirFilter_Step1KHz.c 裡面  有

    /* 80階FIR濾波 */
    arm_fir_f32(&S, FirDataInput, FirDataOutput, BLOCK_SIZE);
   
    /* 輸出濾波器數值 */
    for(i = 0; i < DSOSCREEN_LENGTH; i++)
    {
        g_usWaveBuf =FirDataOutput[i+50];       /*  這裡的  +50 是因為 " 群延迟 "嗎? */
    }

  又請問  是如何算出+50的?  謝謝