firls
最小二乘线性相位数字滤波器的设计
描述
例子
设计一个低通滤波器的过渡带
下面的代码展示了如何设计一个低通滤波器过渡带的订购225。
创建频率和振幅矢量,f和一个。
0.25 - 0.3 f = [0 1]
f =1×40 0.2500 0.3000 1.0000
一个= [1 1 0 0]
一个=1×41 1 0 0
使用firls获取n + 1订单的系数n低通滤波器。
b = firls (255 f);
滤波器的脉冲响应
fvtool (b,“冲动”)
设计一个反对称与分段线性滤波器通频带
下面显示了如何设计一个24阶反对称与分段线性滤波器通频带,和情节预期的和实际的振幅响应。
创建频率和振幅矢量,f和一个。
f = [0 0.3 0.4 0.6 0.7 0.9];一个= [0 1 0 0 0.5 - 0.5];
使用firls获得25滤波器的系数。
f b = firls(24日,,希尔伯特的);
画出理想的振幅响应随着过渡地区。
情节(f。*π,“o”,“markerfacecolor”,(1 0 0));持有在;情节(f。*π,“r——”,“线宽”2);
使用freqz获得设计滤波器的频率响应和策划滤波器的幅度响应。
[H F] = freqz (b, 1);情节(F, abs (H));集(gca),“xlim”π[0])传说(“过滤规范”,“过渡区”,级响应的)
输入参数
输出参数
更多关于
算法
firls设计一个线性相位冷杉过滤器。这个过滤器最小化加权、集成平方误差之间的理想的分段线性函数和滤波器的幅度响应所需的频段。
参考[1]描述背后的理论方法firls。函数解决了涉及一个内积矩阵的线性方程组的大小差不多n / 2使用MATLAB®\运营商。
这个函数设计类型I, II, III, IV线性相位滤波器。I和II型分别是n偶数和奇数的违约。的希尔伯特的和“区别”旗帜生产类型III (n)和第四(n奇怪)过滤器。不同滤波器类型有不同的对称性和限制他们的频率响应(见[2]详情)。
| 线性相位滤波器类型 | 过滤器订单 | 对称系数 | 反应H (f), f = 0 | 反应H (f), f = 1(奈奎斯特) |
|---|---|---|---|---|
I型 |
甚至 |
(k) = b (n + 2 k), k = 1,……n + 1 |
没有限制 |
没有限制 |
II型 |
奇怪的 |
(k) = b (n + 2 k), k = 1,……n + 1 |
没有限制 |
H(1)= 0 |
类型III |
甚至 |
b (k) = - b (n + 2 k), k = 1,……n + 1 |
H(0)= 0 |
H(1)= 0 |
IV型 |
奇怪的 |
b (k) = - b (n + 2 k), k = 1,……n + 1 |
H(0)= 0 |
没有限制 |
引用
[1]公园卢铁荣,和C.S. Burrus,数字滤波器设计,1987岁的约翰·威利& Sons 54 - 83页。
[2]奥本海姆前任所长A.V.以R.W.谢弗,离散时间信号处理新世纪,1989年,页256 - 266。
版本历史
介绍了R2011a
