检测二维地块中的倾斜
我需要自动检测2D绘图中的倾斜，如下图中用红色圆圈标记的区域。我只对“主”倾斜感兴趣，这意味着倾斜必须跨越x轴的最小长度。凹陷的数量是未知的，即不同的地块将包含不同的凹陷数量。有什么想法吗？ 更新： 按照要求，下面是样本数据，并尝试按照VINES的建议使用中值滤波对其进行平滑处理。 看起来我现在需要一种健壮的方法来近似每个点的导数，这样就可以忽略数据中保留的小点。是否有标准方法？
..
signal-processing相关内容
二阶巴特沃斯低通过滤系数的计算
使用 采样频率：10 KHz 截止频率：1 kHz 我如何实际计算下面差分方程的系数？ 我知道差分方程将是这种形式，但不知道如何实际计算并得出系数b0、b1、b2、a1、a2的数字 y(n) = b0.x(n) + b1.x(n-1) + b2.x(n-2) + a1.y(n-1) + a2.y(n-2) 我最终将用C++实现此LPF，但我需要先知道如何实际计算系
..
所需的软件锁相环示例代码
有人知道在哪里可以找到软件锁相环(SPLL)的实际代码示例吗？ 我需要一个能够跟踪1.1 kHz到1.3 kHz之间的PSK调制信号的SPLL。在谷歌上搜索会找到大量的学术论文和专利，但没有有用的东西。即使是去一趟大学图书馆，里面摆满了硬件锁相环的书籍，也只有一本关于SPLL的书中的一个章节，这更多的是理论而不是实践。 感谢您抽出时间。 Ian 推荐答案 这是一个交互式设
..
检测分段恒定信号中的步长
我有一个分段恒定信号，如下所示。我想检测步骤过渡的位置(用红色标记)。 我目前的做法： 使用移动平均过滤平滑信号(http://www.mathworks.com/help/signal/examples/signal-smoothing.html) 执行离散小波变换以获得间断 定位不连续点以获得阶跃过渡的位置 我目前正在实施检测不连续的最后一步。但是，我无法获得准确的位置，
..
学习数字信号处理
了解DSP有哪些好的资源(包括实际了解这些资源所需的数学和算法)？ 让我们假设我的数学技能也因为缺乏使用而生疏，那么下面是这样的路线图： 统计信息刷新器 微积分刷新器 FFT的扎实新手解释 (50步后.) 会很好，有望使DSP技能和知识接近“胜任”水平。 推荐答案 How Do I Learn DSP? A Beginner's Guide to Digi
..
如何检测(心电图)波的模式？
我正在试着从心电图上读取一幅图像，并检测其中的每一个主要波(P波、QRS波和T波)。我可以读取图像并获得向量(如(4.2; 4.4; 4.9; 4.7; ...))。我需要一个算法，它可以遍历这个矢量，并检测每个波的开始和结束时间。示例： 如果它们的大小总是相同的，或者如果我事先知道ECG有多少个波，那就很容易了。给定波形： 我提取向量： [0; 0; 20; 20; 20; 1
..
Python：实时或从WAV文件获取音量分贝级别
对于项目工作，我需要通过录制的音频文件或使用麦克风实时录制来测量音量。经过初步研究，我尝试使用声音文件库。使用soundfile读取该文件，并使用20*np.log10(np.sqrt(np.mean(np.absolute(a)**2)))计算了dB值。我得到了声音文件的负值。但是正常的声音可能在50-70分贝的范围内，我得到的是负值。有人能帮我解决这个问题吗？ 推荐答案 简答：db
..
Librosa基音跟踪-STFT
我正在使用this算法来检测 this音频文件。正如你所听到的，这是在吉他上演奏的E2音符，背景中有一点噪音。 我使用STFT生成了此频谱图： 我使用上面链接的算法，如下所示： y, sr = librosa.load(filename, sr=40000) pitches, magnitudes = librosa.core.piptrack(y=y, sr=sr, fmin=
..
一维数据中的步长检测
在Python中是否有用于检测一维数据中的步骤的现有实现？ 例如，检测到此数据中的一个步骤的内容： 关于算法的描述很多out there但我想知道在Python中是否有适合这项工作的内容？ 我不确定是否/如何提供该数据，但以下是： [ 594. 568.55555556 577.22222222 624.55555556 546.66666667 5
..
如何获得信号的高包络和低包络
我有一个相当嘈杂的数据，我正在试图计算出信号的高包络和低包络。它有点像MATLAB中的这个例子： http://uk.mathworks.com/help/signal/examples/signal-smoothing.html 在&Q；中提取峰值包络&Q；。Python中有没有类似的函数可以做到这一点呢？我的整个项目都是用Python编写的，最坏的情况是，我可以提取Numpy数组并将其
..
去除时间序列中的峰值和台阶之类的跳跃
我在野外有相当多的传感器可以测量水压。在过去，这些传感器的高度已经改变了好几次，造成了时间序列中的跳跃。因为这些时间序列是连续的，并且我有一个手动测量，所以从技术上讲我应该能够删除跳转(手动这很容易，但是测量太多了，所以我需要用python进行测量)。 我已尝试使用过滤中值删除跳转，但没有真正起作用。 我的代码： # filter out noise in signal
..
从CSD中寻找二维空间谱的正确方法
我尝试从上面的方程式实现空间谱(附后) 其中Kx，Ky是k空间中的网格点， C(w，r)-第i和第j个传感器之间的交叉谱密度(这里是大小为ns*ns>；no.传感器的数量)。 x，y是传感器之间的距离。(Kx，Ky的NK网格密度) 我寻找上述等式的合适的python实现。我有34个传感器，它们产生[row*column]=[n*34]大小的数据。首先，我在每个传感器的数据中找出了的交叉
..
基于numpy.fft.fft的功率谱
我通过下面的代码绘制的数字只是零附近的峰值，无论我如何更改数据。我的数据只是一列，它记录了某种信号的每个定时点。time_step是我应该根据数据中两个邻接点的间隔定义的值吗？ data=np.loadtxt("timesequence",delimiter=",",usecols=(0,),unpack=True) ps = np.abs(np.fft.fft(data))**2 tim
..
如何对R中的数据点运行高通或低通过滤？
我是R的初学者，我已尝试查找有关以下内容的信息，但未找到任何信息。 图片中的绿色图形由红色和黄色图形组成。但是假设我只有类似于绿色图表的数据点。如何使用low pass/high pass filter提取低/高频(即近似红色/黄色图形)？ 更新：图形是用 生成的 number_of_cycles = 2 max_y = 40 x = 1:500 a = number_of_cy
..
如何用scipy和lfilter来实时过滤？
免责声明：我可能不太擅长DSP，因此遇到的问题比让此代码正常工作更多。 我需要在传入信号发生时对其进行过滤发送。我试图让这段代码工作，但到目前为止我还没能做到。 引用scipy.signal.lfilter doc import numpy as np import scipy.signal import matplotlib.pyplot as plt from lib import
..
Goertzel算法在C中的实现
..
如何在时间序列中找到顶部和底部?
起初，这个问题听起来很愚蠢，但这不是根本问题.也许，它看起来似乎无法通过任何算法完全解决，但我假装是这样. 所以问题.我有图表，例如黄金.我需要找到时间轴上的顶部和底部在哪里.问题是我需要找到主要好转和主要低迷的开始.问题是有很多不相关的小涨和跌. 这是为了更好地理解的图片 - 红点是我想要找到的(不完全是，但在某种程度上像这样). 所以我可能需要过滤掉小的上调和下调，但不知道该
..
噪声正弦时间序列中的实时峰值检测
我一直在尝试实时检测正弦时间序列数据中的峰值，但是到目前为止我还没有成功.我似乎无法找到一种实时算法，可以以合理的准确度检测正弦信号中的峰值.我要么没有检测到峰值，要么沿着正弦波得到无数个点被检测为峰值. 对于类似正弦波且可能包含一些随机噪声的输入信号，什么是好的实时算法? 作为一个简单的测试用例，考虑一个频率和幅度始终相同的平稳正弦波.(确切的频率和幅度无关紧要；我任意选择了 60
..
使用 Swift 对增长的时间序列进行峰值检测
有没有人有一个好的算法来使用 Swift (v3) 测量增长的时间序列数据中的峰值?因此，在数据流入时检测峰值. 例如平滑 z 波算法的 Swift 版本.那个算法似乎很合适. 我需要检测如下所示的峰值.数据包含正数和负数.输出应该是峰值的计数器，和/或该特定样本的真/假. 样本数据集(上一个系列的总结): 让样本 = [0.01, -0.02, -0.02, 0.01, -0
..
估计两个时间序列之间的小时间偏移
我有两个时间序列，我怀疑它们之间存在时间偏移，我想估计这个时间偏移. 这个问题之前有人问过:求两个(非和谐)波之间的相位差和
..