Python fftpack.diff函数代码示例

OGeek|极客世界-中国程序员成长平台 › 门户 › 编程› Python›Python编程经验

原作者: [db:作者] 来自: [db:来源] 收藏邀请

本文整理汇总了Python中scipy.fftpack.diff函数的典型用法代码示例。如果您正苦于以下问题：Python diff函数的具体用法？Python diff怎么用？Python diff使用的例子？那么恭喜您, 这里精选的函数代码示例或许可以为您提供帮助。

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

示例1: bench_random

 def bench_random(self):
     print()
     print('Differentiation of periodic functions')
     print('=====================================')
     print(' size  |  convolve |    naive')
     print('-------------------------------------')
     for size,repeat in [(100,1500),(1000,300),
                         (256,1500),
                         (512,1000),
                         (1024,500),
                         (2048,200),
                         (2048*2,100),
                         (2048*4,50),
                         ]:
         print('%6s' % size, end=' ')
         sys.stdout.flush()
         x = arange(size)*2*pi/size
         if size < 2000:
             f = sin(x)*cos(4*x)+exp(sin(3*x))
         else:
             f = sin(x)*cos(4*x)
         assert_array_almost_equal(diff(f,1),direct_diff(f,1))
         assert_array_almost_equal(diff(f,2),direct_diff(f,2))
         print('| %9.2f' % measure('diff(f,3)',repeat), end=' ')
         sys.stdout.flush()
         print('| %9.2f' % measure('direct_diff(f,3)',repeat), end=' ')
         sys.stdout.flush()
         print(' (secs for %s calls)' % (repeat))

开发者ID:ChadFulton，项目名称:scipy，代码行数:28，代码来源:bench_pseudo_diffs.py

示例2: test_period

 def test_period(self):
     for n in [17,64]:
         x = arange(n)/float(n)
         assert_array_almost_equal(diff(sin(2*pi*x),period=1),
                                   2*pi*cos(2*pi*x))
         assert_array_almost_equal(diff(sin(2*pi*x),3,period=1),
                                   -(2*pi)**3*cos(2*pi*x))

开发者ID:BranYang，项目名称:scipy，代码行数:7，代码来源:test_pseudo_diffs.py

示例3: test_int

 def test_int(self):
     n = 64
     x = arange(n)*2*pi/n
     assert_array_almost_equal(diff(sin(x),-1),-cos(x))
     assert_array_almost_equal(diff(sin(x),-2),-sin(x))
     assert_array_almost_equal(diff(sin(x),-4),sin(x))
     assert_array_almost_equal(diff(2*cos(2*x),-1),sin(2*x))

开发者ID:BranYang，项目名称:scipy，代码行数:7，代码来源:test_pseudo_diffs.py

示例4: test_random_odd

 def test_random_odd(self):
     for k in [0,1,2,3,4,5,6]:
         for n in [33,65,55]:
             f = random((n,))
             af = sum(f,axis=0)/n
             f = f-af
             assert_almost_equal(sum(f,axis=0),0.0)
             assert_array_almost_equal(diff(diff(f,k),-k),f)
             assert_array_almost_equal(diff(diff(f,-k),k),f)

开发者ID:BranYang，项目名称:scipy，代码行数:9，代码来源:test_pseudo_diffs.py

示例5: test_random_even

 def test_random_even(self):
     for n in [32,64,56]:
         f = random((n,))
         af = sum(f,axis=0)/n
         f = f-af
         # zeroing Nyquist mode:
         f = diff(diff(f,1),-1)
         assert_almost_equal(sum(f,axis=0),0.0)
         assert_array_almost_equal(direct_hilbert(direct_ihilbert(f)),f)
         assert_array_almost_equal(hilbert(ihilbert(f)),f)

开发者ID:BranYang，项目名称:scipy，代码行数:10，代码来源:test_pseudo_diffs.py

示例6: test_expr_large

 def test_expr_large(self):
     for n in [2048,4096]:
         x = arange(n)*2*pi/n
         f = sin(x)*cos(4*x)+exp(sin(3*x))
         df = cos(x)*cos(4*x)-4*sin(x)*sin(4*x)+3*cos(3*x)*exp(sin(3*x))
         ddf = -17*sin(x)*cos(4*x)-8*cos(x)*sin(4*x)\
              - 9*sin(3*x)*exp(sin(3*x))+9*cos(3*x)**2*exp(sin(3*x))
         assert_array_almost_equal(diff(f),df)
         assert_array_almost_equal(diff(df),ddf)
         assert_array_almost_equal(diff(ddf,-1),df)
         assert_array_almost_equal(diff(f,2),ddf)

开发者ID:BranYang，项目名称:scipy，代码行数:11，代码来源:test_pseudo_diffs.py

示例7: test_expr

 def test_expr(self):
     for n in [64,77,100,128,256,512,1024,2048,4096,8192][:5]:
         x = arange(n)*2*pi/n
         f = sin(x)*cos(4*x)+exp(sin(3*x))
         df = cos(x)*cos(4*x)-4*sin(x)*sin(4*x)+3*cos(3*x)*exp(sin(3*x))
         ddf = -17*sin(x)*cos(4*x)-8*cos(x)*sin(4*x)\
              - 9*sin(3*x)*exp(sin(3*x))+9*cos(3*x)**2*exp(sin(3*x))
         d1 = diff(f)
         assert_array_almost_equal(d1,df)
         assert_array_almost_equal(diff(df),ddf)
         assert_array_almost_equal(diff(f,2),ddf)
         assert_array_almost_equal(diff(ddf,-1),df)

开发者ID:BranYang，项目名称:scipy，代码行数:12，代码来源:test_pseudo_diffs.py

示例8: find_peaks

def find_peaks(signal, v_mph, v_mpd, v_cutMin, v_cutMax):
    #signal[-100:100] = 0.0
    
    
    foierSignal = filterFourier(signal, v_cutMin, v_cutMax)
    diffFoierSignal = diff(foierSignal,1)
    
    dev_flag=True
    if (dev_flag == True):
      import matplotlib.pyplot as plt
      from random import randint
      
      plt.figure()
      plt.plot(np.linspace(100, len(signal)-100, num=len(signal.real[100:-100])), 40*signal.real[100:-100], '-r')
      plt.plot(np.linspace(100, len(diffFoierSignal)-100, num=len(diffFoierSignal.real[100:-100])), -diffFoierSignal.real[100:-100], '-g')
      plt.savefig('test/%sdifStress.png'%( randint(0,99) ), dpi=150)

    from detect_peaks import detect_peaks
    ind = detect_peaks(-diffFoierSignal.real[150:-150], mph=v_mph, mpd=v_mpd, show=False)
    
    vmax = []
    vmin = []
    id_vmax = []
    id_vmin = []
    i = 0
    d_find = 50
    while i < len(ind):
        s = signal
        
        if (len(ind) != 0 ):
          if ( np.fabs(len(s)-ind[i]) > d_find and ind[i] > d_find ):
            vmax.append(max( s[ ind[i] - d_find : ind[i] + d_find] ) )     
            id_vmax.append([k for k, j in enumerate(s[ind[i]-d_find:ind[i]+d_find]) if j == vmax[i]][0] + ind[i] - d_find)  
            
          else:
            vmax.append(max( s[ ind[i] - d_find : ind[i]] ) )
            id_vmax.append([k for k, j in enumerate(s[ind[i]-d_find:ind[i]]) if j == vmax[i]][0] + ind[i]-d_find)
            
            
          if ( np.fabs(len(s)-ind[i]) > d_find and ind[i] > d_find ):
            vmin.append(min( s[ ind[i] - d_find : ind[i] + d_find] ) )
            id_vmin.append([k for k, j in enumerate(s[ind[i] - d_find : ind[i] + d_find]) if j == vmin[i]][0]+ind[i]-d_find)
            
          else:
            vmin.append(min( s[ ind[i] : ind[i] + d_find] ) )
            id_vmin.append([k for k, j in enumerate(s[ ind[i] : ind[i] + d_find ] ) if j == vmin[i]][0]+ind[i]-d_find)
          
        i += 1
    return [id_vmax,vmax,id_vmin,vmin]

开发者ID:zhekan，项目名称:MDrun，代码行数:49，代码来源:stress_analis.py

示例9: test_sin

 def test_sin(self):
     for n in [32,64,77]:
         x = arange(n)*2*pi/n
         assert_array_almost_equal(diff(sin(x)),cos(x))
         assert_array_almost_equal(diff(cos(x)),-sin(x))
         assert_array_almost_equal(diff(sin(x),2),-sin(x))
         assert_array_almost_equal(diff(sin(x),4),sin(x))
         assert_array_almost_equal(diff(sin(4*x)),4*cos(4*x))
         assert_array_almost_equal(diff(sin(sin(x))),cos(x)*cos(sin(x)))

开发者ID:BranYang，项目名称:scipy，代码行数:9，代码来源:test_pseudo_diffs.py

示例10: test_zero_nyquist

 def test_zero_nyquist(self):
     for k in [0,1,2,3,4,5,6]:
         for n in [32,33,64,56,55]:
             f = random((n,))
             af = sum(f,axis=0)/n
             f = f-af
             # zeroing Nyquist mode:
             f = diff(diff(f,1),-1)
             assert_almost_equal(sum(f,axis=0),0.0)
             assert_array_almost_equal(diff(diff(f,k),-k),f)
             assert_array_almost_equal(diff(diff(f,-k),k),f)

开发者ID:BranYang，项目名称:scipy，代码行数:11，代码来源:test_pseudo_diffs.py

示例11: time_diff

 def time_diff(self, size, soltype):
     if soltype == 'fft':
         diff(self.f, 3)
     else:
         direct_diff(self.f, 3)

开发者ID:1641731459，项目名称:scipy，代码行数:5，代码来源:fftpack_pseudo_diffs.py

示例12: test_definition

 def test_definition(self):
     for n in [16,17,64,127,32]:
         x = arange(n)*2*pi/n
         assert_array_almost_equal(diff(sin(x)),direct_diff(sin(x)))
         assert_array_almost_equal(diff(sin(x),2),direct_diff(sin(x),2))
         assert_array_almost_equal(diff(sin(x),3),direct_diff(sin(x),3))
         assert_array_almost_equal(diff(sin(x),4),direct_diff(sin(x),4))
         assert_array_almost_equal(diff(sin(x),5),direct_diff(sin(x),5))
         assert_array_almost_equal(diff(sin(2*x),3),direct_diff(sin(2*x),3))
         assert_array_almost_equal(diff(sin(2*x),4),direct_diff(sin(2*x),4))
         assert_array_almost_equal(diff(cos(x)),direct_diff(cos(x)))
         assert_array_almost_equal(diff(cos(x),2),direct_diff(cos(x),2))
         assert_array_almost_equal(diff(cos(x),3),direct_diff(cos(x),3))
         assert_array_almost_equal(diff(cos(x),4),direct_diff(cos(x),4))
         assert_array_almost_equal(diff(cos(2*x)),direct_diff(cos(2*x)))
         assert_array_almost_equal(diff(sin(x*n/8)),direct_diff(sin(x*n/8)))
         assert_array_almost_equal(diff(cos(x*n/8)),direct_diff(cos(x*n/8)))
         for k in range(5):
             assert_array_almost_equal(diff(sin(4*x),k),direct_diff(sin(4*x),k))
             assert_array_almost_equal(diff(cos(4*x),k),direct_diff(cos(4*x),k))

开发者ID:BranYang，项目名称:scipy，代码行数:20，代码来源:test_pseudo_diffs.py

示例13: print

# CsPbI3 - distorted
PB_X=0.469972*NGX
PB_Y=0.530081*NGY
PB_Z=0.468559*NGZ

# CsPbI3 - perfect cubic
#PB_X=0.5*NGX
#PB_Y=0.5*NGY
#PB_Z=0.5*NGZ

# MAPBI3 - pseudo cubic distorted
#PB_X=0.476171*NGX
#PB_Y=0.500031*NGY
#PB_Z=0.475647*NGZ

# Read out massive grad table, in {x,y,z} components
print(grad_x[PB_X][PB_Y][PB_Z],grad_y[PB_X][PB_Y][PB_Z],grad_z[PB_X][PB_Y][PB_Z])
# Norm of electric field at this point
print(np.linalg.norm([ grad_x[PB_X][PB_Y][PB_Z],grad_y[PB_X][PB_Y][PB_Z],grad_z[PB_X][PB_Y][PB_Z] ]))

# OK, let's try this with a Spectral method (FFT)
# JMF - Not currently working; unsure of data formats, need worked example
from scipy import fftpack
V_FFT=fftpack.fftn(grid_pot[:,:,:])
V_deriv=fftpack.diff(grid_pot[:,:,:]) #V_FFT,order=1)

# Standard catch all to drop into ipython at end of script for variable inspection etc.
from IPython import embed; embed() # End on an interactive ipython console to inspect variables etc.

开发者ID:Yvetteoo，项目名称:MacroDensity，代码行数:28，代码来源:FieldAtPoint.py

注：本文中的scipy.fftpack.diff函数示例由纯净天空整理自Github/MSDocs等源码及文档管理平台，相关代码片段筛选自各路编程大神贡献的开源项目，源码版权归原作者所有，传播和使用请参考对应项目的License；未经允许，请勿转载。

鲜花

握手

雷人

路过

鸡蛋

该文章已有0人参与评论

请发表评论

全部评论

专题导读

More+

10-27 六六分期app的软件客服如何联系？(六六分期

11-06 可心卡盟:win10系统火狐flash插件崩溃怎么

11-06 亲亲特价:怎么删除回收站图标

11-06 济南大学虚拟社区:鲁大师节能降温的具体办

11-06 xlueops.exe:无线网络安装向导

11-06 女斗合众国:win7系统cf与主机连接不稳定怎

11-06 0xc000022-[cf烟雾头]cf怎么调烟雾头

11-06 qizideyouhuo:应用程序无法正常启动0xc0000

11-06 ipz-185:win7系统vcf文件怎么打开

11-06 傻哥蹦迪:win10系统s4怎么打开usb调试

11-06 八神浩树gtaste:回收站清空了怎么恢复

11-06 妖尾之黑色守护:win10系统电脑没有1440x900

11-06 校园至尊魔王小说:win7系统浏览网页时字体

11-06 女斗合众国:win10系统访问共享文件夹提示请

11-06 tokyo hot n0654:恢复win7系统默认字体一招

11-06 雨酷仙境:设置win7系统转移临时文件夹腾出

11-06 阿穆纳伊之杖:win7系统开始菜单在右边还原

11-06 tunespotting:win10系统火狐flash插件总是

11-06 甘尔葛分析师：计谋网站seo关键词暴涨有什

11-06 蔡贵霖: 计谋网站seo关键词暴涨有什么秘密

11-06 博益网首页:ao3网页版进入不了解决方法

11-06 漏斗子专栏: 网站数据分析小白易懂精华篇

11-06 见证双虹怎么做:win7系统开启telnet命令的

11-06 颾狐蝶蜋:系统资源不足无法完成请求的服务

11-06 国光中学校歌:提交网站到alexa查询详细步骤

11-06 西安有情天:静态网页和动态网页的区别

11-06 红木雅尚斋:外部链接构造对网站的好处

11-06 前官礼遇：防止域名劫持–增强域安全性的10

11-06 密传二转答案: 中文分词算法有哪些

11-06 金泉家园邮编:百度快照劫持的表现及应对方

Python fftpack.fft函数代码示例发布时间：2022-05-27

Python fftpack.dct函数代码示例发布时间：2022-05-27

Python util.grid_equal函数代码示例

1 Python 入门教程

Python入门教程 Python 是一种解释型、面向对象、动态数据类型的高级程序设计语言。 P

阅读：13828|2022-01-22

2 Python wikiutil.getFrontPage函数代码示例

Python wikiutil.getFrontPage函数代码示例

阅读：10223|2022-05-24

3 Python 简介

Python 简介 Python 是一个高层次的结合了解释性、编译性、互动性和面向对象的脚本

阅读：4099|2022-01-22

4 Python tests.group函数代码示例

Python tests.group函数代码示例

阅读：4046|2022-05-27

5 Python util.check_if_user_has_permission

Python util.check_if_user_has_permission函数代码示例

阅读：3854|2022-05-27

6 Python 操练实例98

Python 练习实例98 Python 100例题目：从键盘输入一个字符串，将小写字母全部转换成大

阅读：3515|2022-01-22

7 Python 环境搭建

Python 环境搭建本章节我们将向大家介绍如何在本地搭建 Python 开发环境。 Py

阅读：3035|2022-01-22

8 Python output.darkgreen函数代码示例

Python output.darkgreen函数代码示例

阅读：2662|2022-05-25

9 Python 基础语法

Python 基础语法 Python 语言与 Perl，C 和 Java 等语言有许多相似之处。但是，也

阅读：2660|2022-01-22

10 Python 中文编码

Python 中文编码前面章节中我们已经学会了如何用 Python 输出 Hello, World!，英文没

阅读：2306|2022-01-22

客服电话

电子邮件

Python fftpack.diff函数代码示例

示例1: bench_random

示例2: test_period

示例3: test_int

示例4: test_random_odd

示例5: test_random_even

示例6: test_expr_large

示例7: test_expr

示例8: find_peaks

示例9: test_sin

示例10: test_zero_nyquist

示例11: time_diff

示例12: test_definition

示例13: print

请发表评论

全部评论

上一篇：

下一篇：

Python util.grid_equal函数代码示例

Python util.get_worker_name函数代码示例

Python util.get_webmention_target函数代

Python util.get_uuid函数代码示例

Python util.get_type_by_name函数代码示例

Python util.grid_equal函数代码示例

Python util.get_worker_name函数代码示例

Python util.get_webmention_target函数代

Python util.get_uuid函数代码示例

Python util.get_type_by_name函数代码示例

Python util.get_stdout函数代码示例

关于我们

产品与服务

解决方案

139-2527-9053