- android - 多次调用 OnPrimaryClipChangedListener
- android - 无法更新 RecyclerView 中的 TextView 字段
- android.database.CursorIndexOutOfBoundsException : Index 0 requested, 光标大小为 0
- android - 使用 AppCompat 时,我们是否需要明确指定其 UI 组件(Spinner、EditText)颜色
我正在尝试使用 buttord
函数设计一个模拟巴特沃斯滤波器(实际上,我正在移植一个从 MATLAB 调用此函数的程序到 Python)。
我的参数是:
通带频率 (Fp) = 10 Hz,给出 Wp = 2*pi*10 Hz
阻带频率 (Fs) = 100 Hz,给出 Ws = 2*pi*100 Hz
通带和阻带损耗/衰减(Rp、Rs)分别为 3 和 80 dB。
在 MATLAB 中我使用这段代码:
Wp = 2 * pi * 10
Ws = 2 * pi * 100
Rp = 3
Rs = 80
[N, Wn] = buttord(Wp, Ws, Rp, Rs, 's')
这给了我 N = 5
,Wn = 99.581776302
。
在 SciPy 中我尝试做同样的事情:
from numpy import pi
from scipy import signal
Wp = 2 * pi * 10
Ws = 2 * pi * 100
Rp = 3
Rs = 80
(N, Wn) = signal.buttord(Wp, Ws, Rp, Rs, analog=True)
我得到 N = 5
和 Wn = 62.861698649592753
。 Wn 与 MATLAB 给出的值不同,而且出奇地接近 Wp。这里有什么问题?
深入研究 SciPy 的来源和问题,我发现了 this pull request这可能解释了一些事情:原来 MATLAB 和 SciPy 有不同的设计目标(MATLAB 试图优化匹配阻带频率,而 SciPy 试图优化匹配通带频率)。
如果重要的话,我正在使用 MATLAB R2013a、Python 3.4.2 和 SciPy 0.15.0。
最佳答案
(我还在 scipy 邮件列表上发布了以下内容。)
当您使用 buttord 设计巴特沃斯滤波器时,没有足够的完全满足所有设计约束的自由度。所以那里是选择过渡区域的哪一端达到约束哪一端是“过度设计”的。 scipy 0.14.0 中所做的更改将该选择从阻带边缘切换到通带边缘。
一张图就清楚了。下面的脚本生成以下图。 (我将 Rp 从 3 更改为 1.5。-3 dB 与 Wn 处的增益一致,这就是您的 Wn 与 Wp 相同的原因。)使用旧约定或新约定生成的滤波器都满足设计约束。使用新约定,响应只会在通带末端碰到约束。
import numpy as np
from scipy.signal import buttord, butter, freqs
import matplotlib.pyplot as plt
# Design results for:
Wp = 2*np.pi*10
Ws = 2*np.pi*100
Rp = 1.5 # instead of 3
Rs = 80
n_old = 5
wn_old = 99.581776302787929
n_new, wn_new = buttord(Wp, Ws, Rp, Rs, analog=True)
b_old, a_old = butter(n_old, wn_old, analog=True)
w_old, h_old = freqs(b_old, a_old)
b_new, a_new = butter(n_new, wn_new, analog=True)
w_new, h_new = freqs(b_new, a_new)
db_old = 20*np.log10(np.abs(h_old))
db_new = 20*np.log10(np.abs(h_new))
plt.semilogx(w_old, db_old, 'b--', label='old')
plt.axvline(wn_old, color='b', alpha=0.25)
plt.semilogx(w_new, db_new, 'g', label='new')
plt.axvline(wn_new, color='g', alpha=0.25)
plt.axhline(-3, color='k', ls=':', alpha=0.5, label='-3 dB')
plt.xlim(40, 1000)
plt.ylim(-100, 5)
xbounds = plt.xlim()
ybounds = plt.ylim()
rect = plt.Rectangle((Wp, ybounds[0]), Ws - Wp, ybounds[1] - ybounds[0],
facecolor="#000000", edgecolor='none', alpha=0.1, hatch='//')
plt.gca().add_patch(rect)
rect = plt.Rectangle((xbounds[0], -Rp), Wp - xbounds[0], 2*Rp,
facecolor="#FF0000", edgecolor='none', alpha=0.25)
plt.gca().add_patch(rect)
rect = plt.Rectangle((Ws, ybounds[0]), xbounds[1] - Ws, -Rs - ybounds[0],
facecolor="#FF0000", edgecolor='none', alpha=0.25)
plt.gca().add_patch(rect)
plt.annotate("Pass", (0.5*(xbounds[0] + Wp), Rp+0.5), ha='center')
plt.annotate("Stop", (0.5*(Ws + xbounds[1]), -Rs+0.5), ha='center')
plt.annotate("Don't Care", (0.1*(8*Wp + 2*Ws), -Rs+10), ha='center')
plt.legend(loc='best')
plt.xlabel('Frequency [rad/s]')
plt.ylabel('Gain [dB]')
plt.show()
关于MATLAB 和 SciPy 对 'buttord' 函数给出了不同的结果,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/28056404/
我有 table 像这样 -------------------------------------------- id size title priority
我的应用在不同的 Activity (4 个 Activity )中仅包含横幅广告。所以我的疑问是, 我可以对所有横幅广告使用一个广告单元 ID 吗? 或者 每个 Activity 使用不同的广告单元
我有任意(但统一)数字列表的任意列表。 (它们是 n 空间中 bin 的边界坐标,我想绘制其角,但这并不重要。)我想生成所有可能组合的列表。所以:[[1,2], [3,4],[5,6]] 产生 [[1
我刚刚在学校开始学习 Java,正在尝试自定义控件和图形。我目前正在研究图案锁,一开始一切都很好,但突然间它绘制不正确。我确实更改了一些代码,但是当我看到错误时,我立即将其更改回来(撤消,ftw),但
在获取 Distinct 的 Count 时,我在使用 Group By With Rollup 时遇到了一个小问题。 问题是 Rollup 摘要只是所有分组中 Distinct 值的总数,而不是所有
这不起作用: select count(distinct colA, colB) from mytable 我知道我可以通过双选来简单地解决这个问题。 select count(*) from (
这个问题在这里已经有了答案: JavaScript regex whitespace characters (5 个回答) 2年前关闭。 你能解释一下为什么我会得到 false比较 text ===
这个问题已经有答案了: 奥 git _a (56 个回答) 已关闭 9 年前。 我被要求用 Javascript 编写一个函数 sortByFoo 来正确响应此测试: // Does not cras
所以,我不得不说,SQL 是迄今为止我作为开发人员最薄弱的一面。也许我想要完成的事情很简单。我有这样的东西(这不是真正的模型,但为了使其易于理解而不浪费太多时间解释它,我想出了一个完全模仿我必须使用的
这个问题在这里已经有了答案: How does the "this" keyword work? (22 个回答) 3年前关闭。 简而言之:为什么在使用 Objects 时,直接调用的函数和通过引用传
这个问题在这里已经有了答案: 关闭 12 年前。 Possible Duplicate: what is the difference between (.) dot operator and (-
我真的不明白这里发生了什么但是: 当我这样做时: colorIndex += len - stopPos; for(int m = 0; m < len - stopPos; m++) { c
思考 MySQL 中的 Group By 函数的最佳方式是什么? 我正在编写一个 MySQL 查询,通过 ODBC 连接在 Excel 的数据透视表中提取数据,以便用户可以轻松访问数据。 例如,我有:
我想要的SQL是这样的: SELECT week_no, type, SELECT count(distinct user_id) FROM group WHERE pts > 0 FROM bas
商店表: +--+-------+--------+ |id|name |date | +--+-------+--------+ |1 |x |Ma
对于 chrome 和 ff,当涉及到可怕的 ie 时,这个脚本工作完美。有问题 function getY(oElement) { var curtop = 0; if (oElem
我现在无法提供代码,因为我目前正在脑海中研究这个想法并在互联网上四处乱逛。 我了解了进程间通信和使用共享内存在进程之间共享数据(特别是结构)。 但是,在对保存在不同 .c 文件中的程序使用 fork(
我想在用户集合中使用不同的功能。在 mongo shell 中,我可以像下面这样使用: db.users.distinct("name"); 其中名称是用于区分的集合字段。 同样我想要,在 C
List nastava_izvjestaj = new List(); var data_context = new DataEvidencijaDataContext();
我的 Rails 应用程序中有 Ransack 搜索和 Foundation,本地 css 渲染正常,而生产中的同一个应用程序有一个怪癖: 应用程序中的其他内容完全相同。 我在 Chrome 和 Sa
我是一名优秀的程序员,十分优秀!