许多问题关于DSP /钢管/巴特沃斯滤波器(模拟离散)与脉冲不变的信息检索方法

6视图(30天)
我工作在一个小的学校项目。我需要创建一个6阶巴特沃斯模拟滤波器转换成离散的版本和脉冲不变的方法。最终产品将IIR 6阶巴特沃斯滤波器和找到它的两极。我知道我缺乏一些基本的理解和我试图赶上。无论如何,这里是我的代码到目前为止,我的问题。
清晰的所有;
clc;
关闭所有;
%的模拟滤波器
num = 1;
穴= [1 3.864 7.464 9.142 7.464 3.864 1);
filterTF =特遣部队(num穴)% 6阶巴特沃斯滤波器传递函数
[b] = lp2lp (num,窝,12566.4);%植入截止频率在2000赫兹
f = 0:1:4000;% 4 k赫兹频率0赫兹之间
w = 2 *π* f;%转换为w
BUTT6 =频率(b、w);%模拟滤波器的频率响应
(z,磷、钾)= tf2zp (b);%找到模拟滤波器的极点
p
图(1)
情节(f, abs (BUTT6));%幅度响应的模拟滤波器
图(2)
pzmap (filterTF)%对模拟滤波器极点
%数字滤波器
(bz, az) = impinvar (b, a, 10000)%转换为数字脉冲不变的方法,采样率为10000 hz
[BUTT6Z, wd] = freqz (bz, az);%频率resposne数字滤波器
阿兹BUTT6digitalTF = filt (bz)
(z1, p1, k1) = tf2zp (bz, az);%找到数字滤波器的极点
p1
图(3)
情节(wd /π/ 2 * 10000,abs (BUTT6Z));%级响应数字滤波器
图(4)
p1 zplane (0)%为数字滤波器极点
% 0用于零因为它去e + 14 ? ? ?
1。谁能确认我的波兰人在模拟和数字滤波器?(我知道波兰人对IIR将RHS zplane。
2。有可能硬编码这个不使用impinvar函数
3所示。我需要重复使用6阶巴特沃斯高通滤波。我怎么能去呢?在s域,分子(s ^ 6) /(我们)^ 6。我该如何改变我的分子?我必须改变分母系数?
4所示。离散传递函数是正确的使用BUTT6digitalTF = filt (bz, az) ?我不很regonize。
5。对离散zplane,为什么我的零e + 14 . .我必须用0 z的地方为了看到两极。
谢谢你!

答案(0)

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!