matlab firpm函数
时间: 2023-12-12 19:35:08 浏览: 33
`firpm`函数是MATLAB中用于设计FIR滤波器的函数之一。它使用逼近法设计线性相位FIR滤波器,可以指定通带、阻带和过渡带的权重,从而实现对滤波器的精确控制。`firpm`函数的语法如下:
```matlab
b = firpm(N, F, A)
b = firpm(N, F, A, W)
b = firpm(N, F, A, W, V)
```
其中,`N`是滤波器的阶数,`F`是一个向量,包含了滤波器的频率取样点,`A`是一个向量,包含了这些取样点的幅度响应,`W`是一个向量,包含了这些取样点的权重,`V`是一个标志,用于指定滤波器类型。如果不指定`W`和`V`,则默认为等权重线性相位滤波器。
下面是一个使用`firpm`函数设计FIR滤波器的例子:
```matlab
% 设计一个50阶的FIR滤波器
N = 50;
% 设计一个通带频率为200Hz,截止频率为150Hz和350Hz的带阻滤波器
F = [0, 150, 200, 300, 350, 500];
A = [0, 0, 1, 1, 0, 0];
% 指定三个频段的权重
W = [3, 1, 100];
% 使用firpm函数计算滤波器系数
b = firpm(N, F/(500), A, W);
% 绘制滤波器的频率响应曲线
fvtool(b, 1);
```
上述代码将设计一个50阶的带阻滤波器,其中通带频率为200Hz,截止频率为150Hz和350Hz,三个频段的权重分别为3、1和100。最后,使用`fvtool`函数绘制滤波器的频率响应曲线。
相关问题
matlab firpm
firpm是MATLAB中用于设计最小最大无突变线性相位(Minimum Maximum Equiripple)FIR滤波器的函数。该函数通过指定频率点和对应的幅度来设计滤波器。具体用法如下所示:
b = firpm(n,f,a)
其中,n是滤波器的阶数,f是频率点矢量,a是与f对应的期望幅度。
f的取值范围是,其中0对应于直流频率,1对应于采样频率的一半。
f必须按升序排列,并且第一个元素必须为0,最后一个元素必须为1。
a是与f对应的期望得到的滤波器的幅度。
另外,调用firpm函数后,可以使用freqz函数绘制滤波器的频率响应曲线。
[h,w = freqz(b,1,512);
plot(f,a,w/pi,abs(h))
legend('Ideal','firpm Design')
xlabel 'Radian Frequency (\omega/\pi)', ylabel 'Magnitude'
上述代码中,h是滤波器的频率响应,w是频率点。
总结一下,firpm函数是用于设计最小最大无突变线性相位FIR滤波器的MATLAB函数,根据给定的频率点和对应的幅度来设计滤波器,并可通过freqz函数绘制频率响应曲线。<span class="em">1</span><span class="em">2</span><span class="em">3</span>
#### 引用[.reference_title]
- *1* [matlab使用firpm函数设计自定义响应的FIR数字滤波器](https://blog.csdn.net/leokingszx/article/details/80195577)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v92^chatsearchT0_1"}}] [.reference_item style="max-width: 33.333333333333336%"]
- *2* [信号处理常用matlab之数字滤波器及滤波函数](https://blog.csdn.net/weixin_37315722/article/details/117692169)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v92^chatsearchT0_1"}}] [.reference_item style="max-width: 33.333333333333336%"]
- *3* [Matlab中firpm函数设计具有非对称衰减的 FIR 带通滤波器](https://blog.csdn.net/D_handsome/article/details/126083500)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v92^chatsearchT0_1"}}] [.reference_item style="max-width: 33.333333333333336%"]
[ .reference_list ]
如何知道MATLAB中的firpm函数计算出的滤波器阶数
MATLAB中的firpm函数是根据一组设计参数(如通带频率、通带衰减、阻带频率、阻带衰减等)自动计算FIR滤波器的阶数。可以使用MATLAB中的length函数或size函数,来查看firpm函数输出的滤波器系数向量的长度,从而确定计算出的滤波器阶数。
例如,假设使用firpm函数设计一个通带截止频率为0.2*pi、通带最大衰减为0.1dB、阻带截止频率为0.3*pi、阻带最小衰减为60dB的低通滤波器,代码如下:
```matlab
fs = 1000; % 采样频率
fpass = 0.2*pi; % 通带截止频率
fstop = 0.3*pi; % 阻带截止频率
Ap = 0.1; % 通带最大衰减(dB)
Ast = 60; % 阻带最小衰减(dB)
[N, Fc, W, B] = firpmord([fpass, fstop]/(fs/2), [1, 0], [Ap, Ast]); % 计算设定参数
h = firpm(N, Fc, W, B); % 设计滤波器
```
运行上述代码后,可以使用length函数查看滤波器系数向量的长度,从而确定计算出的滤波器阶数:
```matlab
>> length(h)
ans =
101
```
可以看到,使用firpm函数计算出的滤波器阶数为101。