tf2sos iirnotch给予不同的结果

4视图(30天)

显示旧的评论

Daoyuan 2014年12月11日

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/166350-tf2sos-of-iirnotch-giving-different-results

评论道: Daoyuan2014年12月15日

答:接受明星黾

你好,

我试图用iirnotch过滤(),我注意到,如果我从iirnotch跑系数(二阶)tf2sos(仍然是相同的数字),(b)的频率响应和(sos)显示在freqz()或fvtool()给非常不同的结果。

如:[b] = iirnotch(我们,bw);[sos1] = tf2sos (b);freqz (b, a, Fs);freqz (sos, Fs);

我想知道为什么有区别吗?

0评论
显示1年长的评论藏1年长的评论

登录置评。

在回答这个问题。

接受的答案

明星黾 2014年12月11日

0
链接

直接链接到这个答案

//www.tatmou.com/matlabcentral/answers/166350-tf2sos-of-iirnotch-giving-different-results answer_162090

你没有邮政代码的细节,所以我只能提供一个通用的回复。文档状态,second-order-section (sos)实现滤波器的典型生产最稳定的版本的过滤器。传递函数(tf)通常是不稳定的。因此,绘制过滤器passband-stopband特征可能是不同的两个实现。(我将使用“sos”实现当我设计一个过滤器)。

6个评论
显示5年长的评论隐藏5年长的评论

Daoyuan 2014年12月12日

谢谢你的回答,我的代码类似于如我给了,举个例子,我是测试电力线路干扰的陷波滤波器。

% 60 Hz噪音,采样频率为300赫兹

我们= 60 / (300/2);bw =我们/ 35;

[b] = iirnotch(我们,bw);

fvtool (b);

sos = tf2sos (b);

fvtool (sos);

从这段代码的结果,第一和第二图看起来完全不同,第一个图有大约-80 db切口平面带其他地方,和第二大切口的只有-40分贝,和后non-flatband等级。但是我很困惑,因为如果你看中的值b &和SOS值,它们是相同的价值观:

sos = 1 0.98236277 -0.607133581 0.98236277 -0.607133581 - 0.96472554

b = 0.98236277 -0.607133581 0.98236277

-0.607133581 - 0.96472554 = 1

然而,fvtool或freqz中的图显示了非常不同的结果。我想知道自切口提供了一个二阶iir滤波器,sos是二阶iir滤波器,相同的价值观,他们是如何产生不同的结果。据我所知,我认为sos应该产生相同的结果与原始iir滤波器如果都是二阶,自从tf2sos将高阶转换成级联二阶滤波器。

明星黾 2014年12月12日

我图你的数据时得到相同的结果 freqz 。我没有DSP系统工具箱所以我不能实验。

我得到好的结果,这段代码为您的过滤器在两种特遣部队和紧急求救信号实现(信号处理工具箱):

                                  Fs = 300;%采样频率
                                 
                                  Fn = f / 2;%奈奎斯特频率
                                 
                                  Ws = [59 61] / Fn;%阻带范围
                                 
                                  Wp = [57 63] / Fn;%通频带的限制
                                 
                                  [n, n] = buttord (Wp Ws 1 10);
                                 
                                  [b] =黄油(n, n,“停止”);%设计阻带滤波器
                                 
                                  (sos, g) = tf2sos (b);
                                 
                                  图(3)
                                 
                                  freqz (b)
                                 
                                  图(4)
                                 
                                  freqz (sos)

这里的滤波器系数:

                                  b = [949.9142 e + e - 003 -1.1746 000 2.2629 -1.1746 e + e + 000 000 949.9142 e - 003);
                                 
                                  a = [1.0000 e e + -1.2047 + 000 000 902.3386 2.2604 -1.1444 e + e + 000 000 e - 003);
                                 
                                  sos = [1.0000 e + e - 003 000 -618.2445 1.0000 1.0000 000 e + e + e - 003 -554.5812 000 949.5072 e - 003;
                                 
                                  1.0000 e + e - 003 000 -618.2445 1.0000 1.0000 000 e + e + e - 003 -650.1666 000 950.3231 e - 003);
                                 
                                  g = 949.9142 e - 003;

Daoyuan 2014年12月15日

好,谢谢你的建议,我将尝试发送错误报告。

= D

登录置评。

类别

信号处理 DSP系统工具箱滤波器的实现单频过滤器

找到更多的在单频过滤器在帮助中心和文件交换

社区寻宝

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

开始狩猎!