MATLAB代码:给定FIR滤波器H(Z)=0.2*(1 + 1z^-1 + 5z^-2+3z^-3++z^4) 请画出其级联型、频率抽样型结构
时间: 2023-10-15 20:04:33 浏览: 117
以下是MATLAB代码:
```matlab
% FIR滤波器系数
h = 0.2 * [1, 1, 5, 3, 1];
% 级联型结构
cascadedH = dfilt.dffir(h);
fvtool(cascadedH); % 绘制幅频响应图
% 频率抽样型结构
upsampleFactor = 2;
downsampleFactor = 3;
upsampledH = upfirdn(h, ones(1, upsampleFactor), upsampleFactor);
downsampledH = downsample(upsampledH, downsampleFactor);
freqSampledH = dfilt.dffir(downsampledH);
fvtool(freqSampledH); % 绘制幅频响应图
```
其中,`dfilt.dffir`函数用于创建FIR滤波器对象,`fvtool`函数用于绘制幅频响应图,`upfirdn`函数用于进行上采样和下采样操作。运行以上代码可以得到级联型结构和频率抽样型结构的幅频响应图。
相关问题
MATLAB代码:采用dir2fs(h)函数实现给定FIR滤波器H(Z)=0.2*(1 + 1z^-1 + 5z^-2+3z^-3++z^4) 请画出其频率抽样型结构级联型
这里是MATLAB代码,使用dir2fs(h)函数实现给定FIR滤波器H(Z)=0.2*(1 + 1z^-1 + 5z^-2+3z^-3++z^4) 的频率抽样型结构级联型:
```
h = [0.2 0.2 1 5 3 1];
fs = dir2fs(h);
fvtool(fs);
```
解释一下:首先定义给定FIR滤波器的系数 $h$,然后使用 dir2fs 函数将其转换为频率抽样型结构级联型,最后使用 fvtool 函数画出其频率响应。
MATLAB代码:采用dir2fs(h)函数实现给定FIR滤波器H(Z)=0.2*(1 + 1z^-1 + 5z^-2+3z^-3++z^4) 请画出其级联型、频率抽样型结构
以下是MATLAB代码实现:
```matlab
% FIR滤波器系数
h = 0.2 * [1, 1, 5, 3, 1];
% 级联型结构
sos = dir2fs(h);
fvtool(sos, 'Analysis', 'freq')
% 频率抽样型结构
z = exp(1i * pi * (0:4) / 5);
b = poly(z);
H = fft(h, 5);
B = fft(b, 5);
G = H ./ B;
fvtool(G, 'Analysis', 'freq')
```
其中,`sos = dir2fs(h)`将FIR滤波器系数转换为级联型结构的系数,`z = exp(1i * pi * (0:4) / 5)`生成5个等间距的频率点,`b = poly(z)`得到插值多项式的系数,`H = fft(h, 5)`和`B = fft(b, 5)`分别计算FIR滤波器系数和插值多项式在5个频率点处的频率响应,`G = H ./ B`得到频率抽样型结构的系数。最后,`fvtool`函数可以用于绘制滤波器的频率响应图。
阅读全文