如何用matlab代码实现窄带滤波器被一系列脉冲激励后的响应
时间: 2023-05-25 14:07:10 浏览: 51
下面是一个实现窄带滤波器被一系列脉冲激励后的响应的 Matlab 代码示例:
```matlab
% 定义滤波器参数
fc = 1000; % 中心频率
bw = 500; % 带宽
fs = 8000; % 采样率
order = 6; % 阶数
% 计算滤波器系数
[b,a] = butter(order,(2/bw)*(fc/fs));
% 生成脉冲激励信号
t = 0:(1/fs):0.005; % 生成时域范围为 5ms 的采样时间点
x = [1,zeros(1,length(t)-1)]; % 生成一个单个脉冲
% 滤波器频率响应
f = 0:10:5000;
H = freqz(b,a,f,fs);
% 对每个脉冲进行滤波,并绘制响应结果
y = filter(b,a,x);
plot(t,y);
xlabel('时间 (s)');
ylabel('幅度');
title('滤波器响应');
% 绘制滤波器频率响应曲线
subplot(2,1,1)
plot(f,abs(H));
xlabel('频率 (Hz)');
ylabel('幅度响应');
title('滤波器频率响应');
subplot(2,1,2);
plot(f,angle(H));
xlabel('频率 (Hz)');
ylabel('相位响应');
title('滤波器相位响应');
```
这段代码首先定义了一个窄带滤波器的参数:中心频率 fc 和带宽 bw、采样率 fs 和阶数 order,然后使用 butter 函数计算出滤波器系数 b 和 a。接着,生成了一个单个脉冲 x 和采样时间点 t。然后,使用 filter 函数对每个脉冲进行滤波,并绘制响应结果。最后,使用 freqz 函数计算并绘制出滤波器的频率响应曲线。
运行以上 Matlab 代码即可实现窄带滤波器被一系列脉冲激励后的响应。