matlab代码设计FIR高通滤波器,使其满足wp=0.8pi,Rp=0.1dB,ws=0.7pi,As=60dB
时间: 2023-08-21 12:01:50 浏览: 111
设计FIR高通滤波器的MATLAB代码如下:
```matlab
% 指定滤波器参数
wp = 0.8*pi; % 通带截止频率
ws = 0.7*pi; % 阻带截止频率
Rp = 0.1; % 通带最大衰减量
As = 60; % 阻带最小衰减量
% 计算滤波器阶数和截止频率
delta_p = (10^(Rp/20)-1)/(10^(Rp/20)+1);
delta_s = 10^(-As/20);
A = -20*log10(min(delta_p,delta_s));
if A<21
n = ceil((ws-wp)/(2*0.92*pi*(ws+wp)));
else
n = ceil(3.1*(ws-wp)/(2*pi*(ws+wp)));
end
wc = (ws+wp)/2;
% 设计FIR滤波器
b = fir1(n,wc/pi,'high',kaiser(n+1,A));
% 绘制滤波器幅频响应曲线
[H,w] = freqz(b,1,1024);
mag = 20*log10(abs(H));
plot(w/pi,mag)
grid on
xlabel('Normalized Frequency (\times\pi rad/sample)')
ylabel('Magnitude (dB)')
title('FIR Highpass Filter Frequency Response')
```
其中,`fir1`函数用于设计FIR滤波器,其参数说明如下:
- `n`:滤波器阶数;
- `wc`:滤波器截止频率;
- `'high'`:指定高通滤波器;
- `kaiser(n+1,A)`:使用Kaiser窗函数进行滤波器设计,其中`A`为滤波器的最小衰减量。
运行上述代码后,可以得到FIR高通滤波器的幅频响应曲线。
阅读全文