发布网友 发布时间:2022-04-25 20:05
共1个回答
热心网友 时间:2023-10-14 19:38
这段代码是级联积分梳状滤波器中梳状滤波器的部分。
[h1,f1]=freqz(ones(1,D),1,1000,fs)在matlab help中查freqz有详细解释及举例。该句用于分析长度为5的FIR滤波器的频率响应,返回值h1为频率响应,f1为频率轴,已经是以Hz为单位了(0~fs/2)
因此下句也就用错了
plot(f1/(fs/2),20*log10(abs(h1))-max(20*log10(abs(h1)))应改为
plot(f1,20*log10(abs(h1))-max(20*log10(abs(h1)))
此时画出的是幅频图,幅度进行了归一化,幅度单位是dB。