cost 371 ms
fftpack 的分段錯誤

[英]Segmentation fault with fftpack

我正在嘗試使用從 https 下載的 fftpack(特別是 rffti.f、rffti1.f、rfftf.f、rfftf1.f、radf4.f、radf2.f、radf3.f、radf5.f、radfg.f) : //www.netlib.org/fftpack/index.html使用此測試代碼 ...

對於小數據量,MKL FFTW 比 FFTPACK 慢

[英]MKL FFTW is Slower than FFTPACK for small data size

20 年前我寫了一個矩陣計算 C++ 庫,我願意使用英特爾 MKL 庫來提高它的性能。 對於復值向量/矩陣,我的庫使用兩個拆分 arrays:一個用於實部,一個用於虛部。 以下是計時結果: N=65536, fftw 時間 = 0.005(s), fftpack 時間 = 0.001(s) N ...

如何將兩個 2D RFFT arrays (FFTPACK) 相乘以兼容 NumPy 的 FFT?

[英]How to multiply two 2D RFFT arrays (FFTPACK) to be compatible with NumPy's FFT?

我正在嘗試將兩個 2D arrays 相乘,這些 2D arrays 用fftpack_rfft2d() (SciPy 的 FFTPACK RFFT)轉換,結果與我從scipy_rfft2d() (SciPy 的 FFT RFFT)得到的結果不兼容。 下圖共享腳本的output,顯示: 兩個輸入a ...

以 FFT 的頻率計算 dFT

[英]computing dFT at the frequencies of the FFT

我正在計算在 x_i, i=0,1,...,N(具有已知 dx)處以頻率 u_j, j=0,1,...,N 采樣的函數 f(x) 的 dFT,其中u_j 是 np.fft.fftfreq(N, dx) 生成的頻率,並將其與 np.fft.fft(f(x)) 的結果進行比較。 我發現兩人不同意... ...

scipy.fftpack的FFT凍結

[英]scipy.fftpack's FFT freezing

當計算大小約為150萬個項目的FFT時: FFT計算永遠不會完成,並且程序被凍結。 如果我將1492828更改為1492827 ,它似乎可以工作。 但是,如果我將1492828更改為1492826 ,它仍然凍結,這有點奇怪。 這是一個已知的錯誤? 注意: CPU ...

如何對不同周期的信號進行FFT

[英]How to do FFT on signal of different period

我需要對非嚴格周期性的步態信號執行FFT。 下面是我為雙面創建的代碼。 但是,由於信號不是嚴格周期性的,我不知道這是否正確-人們的每一步走的方式都不一樣。 如果我每個周期采樣相同數量的樣本,那么我不知道如何賦予fftpack.fftfreq(n,d)參數d。 現在,我將d設置為樣本之間的 ...

為什么運行此腳本會凍結我的計算機?

[英]Why does running this script freeze my computer?

我用Python用SciPy編寫了一個腳本,對信號執行短時傅立葉變換。 當我在具有一千個時間點的信號上運行它時,它運行良好。 當我在具有一百萬個時間點的信號上運行它時,它凍結了我的計算機(計算機沒有響應,並且如果正在播放音頻,計算機將輸出跳過和循環的嗡嗡聲); 這是我嘗試過的所有3次一致的 ...

Python 的 MATLAB fftfilt 等價物

[英]MATLAB fftfilt equivalent for Python

我正在嘗試將在 MATLAB 中創建的以下函數轉換為 Python, 所以,我試圖在這里將 MATLAB 中定義的這個函數實現到 Python 中 我在想函數fftfilt是 MATLAB 中 fftfilt 的克隆,當我運行時出現以下錯誤 所以,我的問題是:Python 中是否有與 MATLAB ...

使用fftpack5.1進行傅立葉變換時出現問題

[英]Trouble with Fourier Transform using fftpack5.1

我在Fortran 90中使用FFTPACK5.1遇到問題,該問題包含用於計算離散傅里葉變換的子例程。 我設法安裝它並使用例程,但是當我檢查頻率為A的簡單正弦波是否一切正常時,得到的非零系數不是A(在頻率空間,頻譜),而是2A 。 頻譜發生了變化,我不明白為什么。 我幾乎可以肯定(但我有疑 ...

使用fftpack在對數軸上的實際數據功率譜

[英]Power spectrum of real data with fftpack on log axis

我已經閱讀了很多關於此主題的討論( 在lomb-scargle和fft之間進行比較 , 在python中繪制功率譜 , Scipy / Numpy FFT頻率分析等),但是仍然無法解決 ,因此我需要一些技巧。 我有一個光子事件列表(檢測與時間),數據可在此處獲取 。 這些列是time , c ...

Java fftpack的基本用法

[英]Basic use of Java fftpack

我在MATLAB中編寫了一些簡單的音頻處理算法,我將用於Android應用程序。 我在將FFT實現轉換為Java時遇到了一些麻煩。 我正在嘗試使用fftpack的Java版本,據說效率很高,但對我來說也有點神秘,也許是因為我對FFT缺乏一貫的把握。 我只是有一個雙數組,我想從中生成 ...

如何在 Android Studio 中導入 java 包

[英]How to import a java package in Android Studio

我是 android 新手,我想將 jfftpack 導入到我在 android studio 中的項目中,但我不知道如何導入它。 原始代碼是 import ca.uol.aig.realdoublefft ,我不知道將 java 文件放在哪里。 我試圖將它們放到 libs 文件夾中。 這是 j ...

以不同的分辨率獲取輸出ifft

[英]getting output ifft at a different resolution

我正在嘗試使用scipy.fftp平滑和插值python中的一些定期數據。 我設法獲取了數據的fft,刪除了wn以上的高階頻率(通過執行myfft [wn:-wn] = 0),然后使用ifft(myfft)重構了數據的“平滑”版本。 ifft創建的數組的點數與原始數據相同。 我如何使用該f ...

Scipy的fftpack dct和idct

[英]Scipy's fftpack dct and idct

假設您使用dct函數,然后不對數據進行操作並使用反轉變換; 倒置數據不會與預轉換數據相同嗎? 為什么浮點問題? 是報告的問題還是正常行為? 任何人都有解釋為什么? 當然,一個簡單的c*10**-1可以解決問題,但是如果你重復調用函數在幾個維度上使用它,那么錯誤就會變大: ...


 
粵ICP備18138465號  © 2020-2024 STACKOOM.COM