全通滤波器 相位校正 matlab,一文读懂滤波器的线性相位,全通滤波器,群延迟...

延迟和

全通滤波器,来解释这些概念。

先说结论:

线性相位能保证信号中各频率成分的相对相位关系不改变。

通俗解释是:

信号经过线性相位滤波器后,各个频率分量的延时时间是一样的。

1. 延迟

举一个最简单的FIR的例子,延迟。假设16kHz的采样频率,一个采样周期的延迟,可以用FIR来表示。利用Matlab来观看这个滤波器的频率响应,代码如下。采样频率为Fs = 16kHz, 采样周期为Ts,Ts = 1/Fs。

num = [0,1]

den = [1,0]

fvtool(num,den)

下图中,蓝色的实线表示的幅频响应,为0dB。红色实线表示相频响应,主要看相频响应。图中对1k,2k,4kHz频点的横纵坐标有截图,16kHz采样率下,1kHz的正弦信号一个完整的周期(这里说的周期指的是2 π 2\pi2π)内会得到16个采样值。一个采样周期的延迟,带来的相位变化是− 2 π / 16 = − π / 8 = − 22.5 ° -2\pi/16=-\pi/8=-22.5\degree−2π/16=−π/8=−22.5°。而2kHz信号一个完整周期(2 π 2\pi2π)内会得到8个采样值,那么一个采样周期的延迟带来的相位变化是− 2 π / 8 = − π / 4 = − 45 ° -2\pi/8=-\pi/4=-45\degree−2π/8=−π/4=−45°,同理,对4kHz的信号,相位变化是− 2 π / 4 = − π / 2 = − 90 ° -2

  • 2
    点赞
  • 19
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
要实现fir非线性相位全通滤波器,并处理信号,可以按照以下步骤进行: 1.设计fir非线性相位全通滤波器 可以使用matlab中的firls函数或firpm函数进行设计,具体方法如下: firls函数: ```matlab % 设计非线性相位全通滤波器 fs = 1000; % 采样率 f1 = 20; % 通带截止频率 f2 = 200; % 阻带截止频率 M = 100; % 滤波器阶数 % 生成频率响应 h1 = [ones(1, f1/(fs/2)), zeros(1, (f2-f1)/(fs/2)), ones(1, (fs/2-f2)/(fs/2))]; h2 = fliplr(h1); % 设计非线性相位全通滤波器 h_allpass = firls(M, [0 f1 f2 fs/2]/(fs/2), [h1 h2 h1 h2]); ``` firpm函数: ```matlab % 设计非线性相位全通滤波器 fs = 1000; % 采样率 f1 = 20; % 通带截止频率 f2 = 200; % 阻带截止频率 M = 100; % 滤波器阶数 % 生成频率响应 h1 = [ones(1, f1/(fs/2)), zeros(1, (f2-f1)/(fs/2)), ones(1, (fs/2-f2)/(fs/2))]; h2 = fliplr(h1); % 设计非线性相位全通滤波器 h_allpass = firpm(M, [0 f1 f2 fs/2]/(fs/2), [h1 h2 h1 h2], 'hilbert'); ``` 2.读入信号并进行滤波 读入信号可以使用matlab中的audioread函数或wavread函数,具体方法如下: ```matlab % 读入信号 [x, fs] = audioread('filename.wav'); ``` 然后进行滤波处理,可以使用matlab中的filter函数,具体方法如下: ```matlab % 进行滤波处理 y = filter(h_allpass, 1, x); ``` 3.播放或保存处理后的信号 播放处理后的信号可以使用matlab中的sound函数,具体方法如下: ```matlab % 播放处理后的信号 sound(y, fs); ``` 保存处理后的信号可以使用matlab中的audiowrite函数或wavwrite函数,具体方法如下: ```matlab % 保存处理后的信号 audiowrite('filename_processed.wav', y, fs); ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值