signal-processing

    0热度

    2回答

    我正在使用PsychToolbox在fMRI扫描仪上运行实验。我想收集扫描仪发出的脉冲时间点。如果您对此主题不熟悉:扫描仪发出一个信号 - 相当于按定义的时间间隔按键盘数字“5”。第一个这样的信号启动代码,使得扫描器的测量和代码的开始时间同步。我的代码是顺序的,如下图所示。我怎么能实现一个循环,在任何时候检查这些输入的“5”,即使我的主代码运行在“for”循环中? 我的代码: % here I w

    4热度

    4回答

    我有两个不同长度的Python列表。人们可能会认为其中一个列表比另一个列表大数倍。 这两个列表都包含相同的物理数据,但采集的采样率不同。 我的目标是对较大的信号进行降采样,使其具有与较小信号一样多的数据点。 我想出了下面的代码,基本上没有工作,但既不是非常Python化,也不能够在一个高性能的方式处理非常大的列表: import math a = [1,2,3,4,5,6,7,8,9,10]

    0热度

    1回答

    我有两个滤波器效果(低通和高通类型)。但是,在处理完audiosample缓冲区后,两种效果在运行时都不会听到。要声明的第二个过滤器只能听到。我已经通过交换调用共振频率位的顺序进行了很多实验,第二个要声明的过滤器类型只是使它在运行时输出。 我得到了一个建议,描述了一个名为“音频链接”的过程。有人对此有什么想法吗?或者不能在音频系统中使用两种不同滤波器类型的超级DSP处理缓冲器?请任何帮助都将非常感

    0热度

    1回答

    我有一个包含4500个长度为1800的向量的矩阵,为此我需要计算矩阵中每两个向量之间的DTW(动态时间翘曲)距离。 我已经使用一个嵌套循环填满一个4500x4500矩阵(其看起来像一个三角形)的一半: matr = zeros(4500,4500); %initializing empty matrix to fill dtw distance for i=1:4500 x = new

    0热度

    1回答

    我正在使用来自Surge库的FFT功能。返回数组中的每个元素与什么相关?我怎么知道哪些频率对应哪些值? 例如,我试图在某个频率范围内找到峰值。

    2热度

    1回答

    我的测量值和相应的权重一大熊猫数据帧的一系列加权值的: df = pd.DataFrame({'x': np.random.randn(1000), 'w': np.random.rand(1000)}) 我要平滑的测量值(x)同时服用逐元素 权重( w)。这与滑动窗户的重量无关,我还希望应用其中的 (例如,三角形窗口或更有用的东西)。因此,为了计算每个窗口内的平滑值,该函数不仅应该通过窗函数

    -2热度

    1回答

    我正在编写从加速计获取信号的代码,并使用TVD算法来消除信号。我们可以在下面的图表看,我没做到这一点: 但我要检测的高峰,我想,我可以用MATLAB中的findchangepts功能,但是当我运行它,我具备以下条件: 我们可以看到它不能检测到峰值。有没有人知道该功能中的任何设置,或者在峰内去噪的方法?

    -2热度

    2回答

    我想确定估计模型对未来新数据的适用程度。为此,经常使用预测误差图。基本上,我想比较测量的输出和模型输出。我使用最小均方算法作为均衡技术。有人可以请帮助什么是适当的方式来绘制模型和测量数据之间的比较?如果估计值接近真值,那么曲线应该非常接近。以下是代码。 u是均衡器的输入,x是有噪声的接收信号,y是均衡器的输出,w是均衡器权重。应该使用x和y*w来绘制图表吗?但是x是嘈杂。我很困惑,因为测量的输出x

    0热度

    1回答

    得到传递函数F(1/Z)我定义的传递函数如下: I =[1,zeros(1,65800)]; f = filter(num, den, I); F = fftshift(fft(f)); 其中F为F(Z)。使用F,我怎样才能得到F(1/z)?

    3热度

    1回答

    我有一些音频数据(浮点数组),我用它来绘制一个简单的 波形。 绘制时,波形不会在边缘出现最大值。 没问题 - 数据只需要标准化。我迭代一次找到最大值,然后再次迭代,再除以最大值。再画一次,一切都很棒! 但是等待有一个大声的介绍或大爆炸的视频导致其余的波形仍然很小。 经过一番研究,我遇到了应该解决这个问题的RMS。我遍历样本并计算RMS,并再次将每个样本除以RMS值。这导致相当大的 “裁剪”: 什么