signal-processing

0热度

1回答

我有一个信号，想带通滤波器是： def butter_bandpass_prep(lowcut, highcut, fs, order=5): """Butterworth bandpass auxilliary function.""" nyq = 0.5 * fs # Minimal frequency (nyquist criterion) low = lowc

-2热度

1回答

如何计算频带内的功率谱？

我有离散的脑电信号，我试图从每个通道中提取绝对功率。采样率为30Hz。我感到困惑的是我不知道如何准确编程。首先我做： y = fft(eegdata); 这给了我一个复杂的两倍。我必须在Hamming窗口+ 50％重叠的频带（alpha +）内进行计算，但我不知道如何去做。任何人都可以帮我解决这个问题吗？

-2热度

1回答

随着时间改变导入信号的音高对数/指数

如何随时间改变导入信号的音高对数/指数关系？请注意，使用的导入信号不是单一频率，所以简单的扫描或啁啾命令将不起作用，因为我将导入人声音频文件，我只是创建了下面的示例，可以测试/显示我遇到的问题。我可以随时间线性地改变信号的音高，这很好用，请参阅下面的测试代码和频率图的第1部分。由于Sheljohn的代码 %Sweep question part 1 clear all,clf reset,t

1热度

2回答

比较峰的算法：它们是否同相？

我正在开发一个比较两个数字列表的算法。这些列表代表使用robust peak detection method在信号中发现的峰值。我想提出一些方法来确定峰值是否处于同相，异相或两者之间（无法确定）。例如：这些阵列将在阶段考虑： [ 94 185 278 373 469], [ 89 180 277 369 466] 但这些阵列会出的阶段： [51 146 242 349], [99 200 30

0热度

1回答

查找条件y值的x轴位置

我有大约180000个数据点，看起来像这样 - Figure1。它的0和5以不规则的间隔。 data<- 0,0,0,5,5,5,0,0,0,5,5,5,0,0.. and so on. 我想找到开始和结束的y值是5（用蓝色标记）的索引。图中附有图。我正在使用R进行分析。我尝试使用changepoint软件包中的cpt.meanvar函数，但它并未提供所有更改点位置。有另一种方法吗？

1热度

1回答

如何从麦克风输入读取实时频率？

我想实时获取从麦克风获取的语音输入的频率。我搜索了这个，并了解了FFT和另外的2,3种算法，但是实现这些算法似乎非常复杂。我正在寻找一个C＃库，使我能够将频率简单地插入到数组中而无需实现它。

1热度

1回答

如何在matlab中同时绘制图表来播放音频文件？

我想播放音频文件同时绘制和更新图表。但是，我的音频文件被中断。我想在后台顺利播放音频文件并同时更新图形。 for i=1:10 player = audioplayer(audio, Fs); play(player); scatter(x(i),y(i),'r.') end

-4热度

1回答

C编程中的晶格滤波器？

我想在C中实现晶格fir过滤器。我的问题是如何使延迟时间（z-1）。我试着计算每次迭代的前一个样本，但是第一个样本的第一个样本有负值。您可以找到算法here

0热度

1回答

窗口大小和重叠PSD的评价matlab pwelch函数

任何人都可以请建议理想的窗口大小和重叠的样本在Matlab中的pwelch函数。我有几个200 ms的EEG信号，采样率为1000（信号长度或采样数量= 200），以评估频谱功率。默认情况下，pwelch使用汉明窗口并将数据分成8段，重叠度为50％。对于只有200个样本的信号，这些默认值是否正确？默认值工作正常，它给我一个PSD图。我想确保我所做的是在概念上是正确的，如果有人可以提出任何更好的方法

6热度

1回答

使用倍频程/ matlab代码随时间递增/逐渐改变信号音调

我可以使用resample对整个信号进行音高移位，并尝试了相位声码器代码here。我也试过repmat and interpolation，我看着fft and interp1 我怎样才能逐步/逐步改变信号的音高随着时间的推移？我已经包含了一个Original Signal的示例，我试图让Processed Signal听起来像（我使用Audacity创建了处理后的信号并使用其效果Sliding