功率谱c语言 微盘,Cython的功率谱

我正在尝试使用Cython优化我的代码.它正在做一个功率谱,而不是使用FFT,因为这是我们在课堂上被告知要做的事情.我曾尝试在Cython中编写代码,但没有看到任何区别.这是我的代码

#! /usr/bin/env python

# -*- coding: utf8 -*-

from __future__ import division

cimport numpy as np

import numpy as np

cimport cython

@cython.boundscheck(False)

def power_spectrum(time, data, double f_min, double f_max, double df,w=1 ):

cdef double com,f

cdef double s,c,sc,cc,ss

cdef np.ndarray[double, ndim=1] power

cdef np.ndarray[double, ndim=1] freq

alfa, beta = [],[]

m = np.mean(data)

data -= m

freq = np.arange( f_min,f_max,df )

for f in freq:

sft = np.sin(2*np.pi*f*time)

cft = np.cos(2*np.pi*f*time)

s = np.sum( w*data*sft )

c = np.sum( w*data*cft )

ss = np.sum( w*sft**2 )

cc = np.sum( w*cft**2 )

sc = np.sum( w*sft*cft )

alfa.append( ( s*cc-c*sc )/( ss*cc-sc**2 ))

beta.append( ( c*ss-s*sc )/( ss*cc-sc**2 ))

com = -(f-f_min)/(f_min-f_max)*100

print "%0.3f%% complete" %com

power = np.array(alfa)**2 + np.array(beta)**2

return freq,power,alfa,beta

时间和数据通过numpy.loadtxt加载并发送到此函数.当我做

cython -a power_spectrum.pyx

.html文件很黄,所以效率不高.特别是整个for-loop和电源的计算并返回一切.

我曾尝试阅读Cython的官方指南,但由于我从未用C语言编写,因此有点难以理解.

所有帮助都非常精确:)