5.2 FIR滤波器的卷积(脉冲响应、频率响应,幅频响应,暂态-搞不清楚的有救了啊)

-5.2 FIR滤波器的卷积
如果将抽头系数的个数用M表示:卷积公式如下:
在这里插入图片描述
上述困扰大多大学生的卷积公式如何理解呢?
h(k)就是滤波器的信道系数,x(n-k)是所谓经过延迟k个器件后的数据。
这句话不理解的话,可以结合上一节,我们所化的FIR框图来理解。

单位脉冲响应
一个滤波器的单位脉冲响应就如同它的名字一样:它是滤波器对一个单位脉冲序列的响应,也就是当滤波器的输入是一个单位脉冲时的时域输出序列。

单位脉冲序列:就是在某一个时刻有一个有用信息进入,以表示为1,其余时刻全部没有输入信号,表示为0.

从某种意义上讲:当某人问你一个FIR滤波器的脉冲响应时,他等于是在问你FIR滤波器的系数。因为FIR滤波器的系数是个数有限的,所以它的脉冲响应在时域上也是有限的,即FIR滤波器有限脉冲响应。
在这里插入图片描述
脉冲序列与信道系数卷积获得脉冲响应。
如上图所示,我们一直把信道卷积系数定义为1/5,也就是每次都是计算的过去5分钟,平均每分钟的车流量。
很明显,我们可以通过更改滤波器的系数,可以使它具有更好的表现。

如何评价一个滤波器是否表现良好
我们可以通过定义滤波器的频率响应来判断滤波器是否“”表现“”良好
频率响应:y(n)=h(k)*x(n) 其中h(k)代表的是**脉冲响应(FIR滤波器的系数)**也就是过去h(k)分钟,每分钟平均车流量。

利用时域相卷等于频域相乘,滤出掉我们不想要的频率的部分,我们可以得到想要的频率部分,然后再对信号做IDFT。

那么可以对h(k)做FFT,那么就可以得到滤波器的幅频响应|H(m)|和相频特性

下图是五个点1/5,我们再扩展到64个点,再做FFT,可以得到幅频响应。b图是归一化后的幅频特性
在这里插入图片描述

  • 暂态

    可以观察到输出信号和输入信号频率相同,这是LTI系统的本质,输入某个频率的正弦信号,输出同样是该频率的正弦信号,只不过幅度和相位发生了改变。

    其次,可以注意到输出信号的前面4个点不是正弦形状的,这些点只是滤波器的暂态过渡响应。在这种抽头延迟线FIR滤波器中,过渡响应点数与滤波器的延拓单元数量D相同,在过渡响应之后才是滤波器的稳态时域响应

    它意味着对于拥有D个延迟器的抽头延迟线FIR滤波器而言前面D个系统的输出并不能代表系统对该序列的正确响应。直到D+1个输出才是系统对该序列的准确响应。

    在这里插入图片描述

    如果换一个更高频率的的输入信号,输出信号也是同频率的正弦波,但是幅度却降为了原来的0.69.
    这是由低通滤波器的本质决定的,低通滤波器对高频信号的衰减大于低频信号的衰减幅度。

  • 总结
    FIR滤波器的时域卷积是通过对移位的输入信号与信道系数的乘积求和得到的。

FIR滤波器的输出序列等于输入序列与信道脉冲响应(信道抽头系数)的卷积。

FIR滤波器的幅频响应等于该滤波器脉冲响应的DFT

FIR滤波器的输出信号幅频响应等于输入信号幅频响应和滤波器自身幅频响应的乘积

时域相卷等于频域相乘

  • 9
    点赞
  • 111
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
好的,您的问题是关于FIR滤波器设计和仿真分析的。下面是我的回答: 等波纹逼近法(Equiripple Approximation Method)是一种常用的FIR滤波器设计方法,其主要思想是在通带和阻带上均匀分布若干个波纹,使得滤波器幅频响应在通带和阻带的最大值和最小值之间波动最小。其设计过程一般包括以下几个步骤: 1. 确定滤波器的通带和阻带边界频率以及通带和阻带的最大衰减量。 2. 根据通带和阻带的频率响应要求,选择适当的滤波器类型和阶数。 3. 使用等波纹逼近法设计滤波器的单位脉冲响应,得到滤波器的系数序列。 4. 利用MATLAB等工具对设计的FIR滤波器进行仿真验证,包括单位脉冲响应频率响应的分析。 具体的设计和仿真过程可以参考以下步骤: 1. 确定滤波器的通带和阻带频率以及通带和阻带的最大衰减量。例如,假设我们需要设计一个抗混叠的低通滤波器,通带频率为1kHz,阻带频率为5kHz,通带最大衰减量为1dB,阻带最小衰减量为40dB。 2. 选择适当的滤波器类型和阶数。根据通带和阻带的频率响应要求,选择适当的滤波器类型和阶数。常用的滤波器类型包括矩形窗、汉明窗、布莱克曼窗等。一般来说,滤波器的阶数越高,频率响应的拟合精度就越高,但是计算量也会增加。例如,假设我们选择了一个10阶的矩形窗滤波器。 3. 设计滤波器的单位脉冲响应。利用等波纹逼近法,设计滤波器的单位脉冲响应。在MATLAB中,可以使用firpm函数进行设计。例如,下面的代码可以生成一个10阶的矩形窗滤波器的系数序列: ``` fs = 10e3; % 采样率 f = [0 1e3 5e3 fs/2]/(fs/2); % 频率向量 m = [1 1 0 0]; % 通带和阻带的增益值 dev = [0.05 0.01 0.05]; % 通带和阻带的最大波纹 b = firpm(10,f,m,dev); % FIR滤波器系数 ``` 4. 对设计的滤波器进行仿真验证。利用MATLAB等工具,对设计的滤波器进行仿真验证,包括单位脉冲响应频率响应的分析。例如,下面的代码可以对设计的滤波器进行单位脉冲响应频率响应的仿真分析: ``` % 单位脉冲响应分析 n = 0:99; % 采样点数 h = filter(b,1,[1 zeros(1,99)]); % FIR滤波器的单位脉冲响应 subplot(2,1,1); stem(n,h); title('单位脉冲响应'); xlabel('采样点数'); ylabel('幅度'); % 频率响应分析 [H,f] = freqz(b,1,4096,fs); % FIR滤波器频率响应 subplot(2,1,2); plot(f/1e3,20*log10(abs(H))); title('频率响应'); xlabel('频率/kHz'); ylabel('增益/dB'); ``` 以上就是使用等波纹逼近法设计FIR滤波器单位脉冲响应频率响应仿真分析的具体步骤和代码示例。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值