python correlate_Python signal.correlate方法代码示例

本文整理汇总了Python中scipy.signal.correlate方法的典型用法代码示例。如果您正苦于以下问题：Python signal.correlate方法的具体用法？Python signal.correlate怎么用？Python signal.correlate使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在模块scipy.signal的用法示例。

在下文中一共展示了signal.correlate方法的27个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于我们的系统推荐出更棒的Python代码示例。

示例1: test_pdos_1d

​点赞 6

​

# 需要导入模块: from scipy import signal [as 别名]

# 或者: from scipy.signal import correlate [as 别名]

def test_pdos_1d():

pad=lambda x: pad_zeros(x, nadd=len(x)-1)

n=500; w=welch(n)

# 1 second signal

t=np.linspace(0,1,n); dt=t[1]-t[0]

# sum of sin()s with random freq and phase shift, 10 frequencies from

# f=0...100 Hz

v=np.array([np.sin(2*np.pi*f*t + rand()*2*np.pi) for f in rand(10)*100]).sum(0)

f=np.fft.fftfreq(2*n-1, dt)[:n]

c1=mirror(ifft(abs(fft(pad(v)))**2.0)[:n].real)

c2=correlate(v,v,'full')

c3=mirror(acorr(v,norm=False))

assert np.allclose(c1, c2)

assert np.allclose(c1, c3)

p1=(abs(fft(pad(v)))**2.0)[:n]

p2=(abs(fft(mirror(acorr(v,norm=False)))))[:n]

assert np.allclose(p1, p2)

p1=(abs(fft(pad(v*w)))**2.0)[:n]

p2=(abs(fft(mirror(acorr(v*w,norm=False)))))[:n]

assert np.allclose(p1, p2)

开发者ID:elcorto，项目名称:pwtools，代码行数:25，

示例2: test_rank0

​点赞 6

​

# 需要导入模块: from scipy import signal [as 别名]

# 或者: from scipy.signal import correlate [as 别名]

def test_rank0(self, dt):

a = np.array(np.random.randn()).astype(dt)

a += 1j * np.array(np.random.randn()).astype(dt)

b = np.array(np.random.randn()).astype(dt)

b += 1j * np.array(np.random.randn()).astype(dt)

y_r = (correlate(a.real, b.real)

+ correlate(a.imag, b.imag)).astype(dt)

y_r += 1j * (-correlate(a.real, b.imag) + correlate(a.imag, b.real))

y = correlate(a, b, 'full')

assert_array_almost_equal(y, y_r, decimal=self.decimal(dt) - 1)

assert_equal(y.dtype, dt)

assert_equal(correlate([1], [2j]), correlate(1, 2j))

assert_equal(correlate([2j], [3j]), correlate(2j, 3j))

assert_equal(correlate([3j], [4]), correlate(3j, 4))

开发者ID:Relph1119，项目名称:GraphicDesignPatternByPython，代码行数:19，

示例3: test_consistency_correlate_funcs

​点赞 6

​

# 需要导入模块: from scipy import signal [as 别名]

# 或者: from scipy.signal import correlate [as 别名]

def test_consistency_correlate_funcs(self):

# Compare np.correlate, signal.correlate, signal.correlate2d

a = np.arange(5)

b = np.array([3.2, 1.4, 3])

for mode in ['full', 'valid', 'same']:

assert_almost_equal(np.correlate(a, b, mode=mode),

signal.correlate(a, b, mode=mode))

assert_almost_equal(np.squeeze(signal.correlate2d([a], [b],

mode=mode)),

signal.correlate(a, b, mode=mode))

# See gh-5897

if mode == 'valid':

assert_almost_equal(np.correlate(b, a, mode=mode),

signal.correlate(b, a, mode=mode))

assert_almost_equal(np.squeeze(signal.correlate2d([b], [a],

mode=mode)),

signal.correlate(b, a, mode=mode))

开发者ID:Relph1119，项目名称:GraphicDesignPatternByPython，代码行数:20，

示例4: calculateBackgroundSignal

​点赞 6

​

# 需要导入模块: from scipy import signal [as 别名]

# 或者: from scipy.signal import correlate [as 别名]

def calculateBackgroundSignal(self, mat, vmat, nuc_cov):

offset=self.start-mat.start-vmat.w

if offset<0:

raise Exception("Insufficient flanking region on \

mat to calculate signal")

self.vmat = vmat

self.bias_mat = mat

self.cov = CoverageTrack(self.chrom, self.start, self.end)