signal-processing

1熱度

1回答

我創建了一個python Gnuradio塊，現在我用C++重新編寫它。我注意到的是一件非常意外的事情 - 與Python版本（18％）相比，C++塊（python流程圖進程）消耗更多的CPU（〜125％）。我必須做一些錯誤的...所以 - 我創建的投入和產出比設定變量類型爲float沒有其他自定義代碼和編號的新塊1，我看到相同的行爲。我一定是做錯了什麼，但我不能告訴什麼... $ gnuradi

0熱度

1回答

巴特沃斯濾波器x平移

我有一個信號，想帶通濾波器是： 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

0熱度

1回答

在Matlab上使用FFT計算和繪製信號的頻譜

我在Matlab中使用FFT計算我的信號的頻譜時缺少一些東西。我的代碼：你會這麼好心來告訴我，我搞砸了？我試圖檢查是否algorythm工程，並使用相同的信號的這兩個取樣（相同的採樣頻率;在兩個不同的時間範圍0-10和0-100）： fs=1000; time10 = [0:1/fs:10]; time100 = [0:1/fs:100]; data10 = sin(2*pi*0

2熱度

1回答

Matlab fft函數交換索引

我寫了一個簡單和非常小的Matlab代碼，用於計算給定數組（或矢量）的離散傅里葉變換。我工作過，並得到了答案，我的Matlab代碼也給出了相同的答案。但fft通過交換索引給出了與此不同的答案。以下是我已經做了mannual算了一筆賬：這是第二圖像：這是第三個形象：從這些計算中可以明顯看出我的答案是{12, -3-3j, -2, -3+3j} 這裏是Matlab代碼，我用： clc; cl