爆破震动信号的时频分析
上传者:刘文丽|上传时间:2015-05-08|密次下载
爆破震动信号的时频分析
第27卷 第12期 岩 土 工 程 学 报 Vol.27 No.12 2005年 12月 Chinese Journal of Geotechnical Engineering Dec., 2005
爆破震动信号的时频分析
内容需要下载文档才能查看张义平
内容需要下载文档才能查看中南大学 资源与安全工程学院
左宇军
摘
内容需要下载文档才能查看要
内容需要下载文档才能查看进行时频分析比较研究
小波变换
细地反映了质点震动强度随时间的衰减起伏变化
关键词
张义平(1970– )男
非平稳随机信号
TD 235.1 文献标识码
STFT
内容需要下载文档才能查看STFT
内容需要下载文档才能查看使用的窗函数固定
DWT
其分析结果较详
HHT自适应性强
内容需要下载文档才能查看并且大部分能量都集中在有限的能量谱线上
主要从事岩土灾害控制
LI Xi-
内容需要下载文档才能查看bing
Abstract
WANG Wei-hua
(School of Resources and Safety Engineering, Central South University, Changsha 410083, China)
mechanics of explosion; non-stationary random signal; time-frequency analysis; Hilbert-Huang transform
0 引 言
爆破震动信号是一种非平稳随机信号
内容需要下载文档才能查看[1,2]
内容需要下载文档才能查看ÈçÏà¹Øº¯Êý
内容需要下载文档才能查看ÊÇʱ±äº¯Êý
内容需要下载文档才能查看[8,9]
Ôڱ任Öл¹»á²úÉú©йЧӦ
内容需要下载文档才能查看HHT
²¢ÇÒ
FFT
WT和HHT三种不同的时频分
析方法对爆破震动信号进行分析
对爆破震动信
号进行时频分析就是为了更全面而准确地得到信号能
内容需要下载文档才能查看量随时间
短时Fourier变换
内容需要下载文档才能查看WT
1 STFT方法
STFT是以Fourier变换为基础的最早和最常用的
───────
内容需要下载文档才能查看基金项目
50490274
内容需要下载文档才能查看10472134
内容需要下载文档才能查看2002CB412703
ÎÈ̬»ù´¡ÉϵĴ¿ÆµÓò·ÖÎö·½·¨
万方数据
第12期 张义平
内容需要下载文档才能查看对于非平稳随机信号
为使变换
内容需要下载文档才能查看具有时域局部性
内容需要下载文档才能查看就得到信号的短时谱
其表达式为
S(w,τ
内容需要下载文档才能查看)=∫f(t)g*(w?τ)e
?iwt
dt
R
³£½«STFT取模平方就可得到信号能量在时
频二维中的分布情况
即窗口面积
对信号的高频成
内容需要下载文档才能查看分
显微镜
的功能
内容需要下载文档才能查看[16]
内容需要下载文档才能查看 内容需要下载文档才能查看CWT
内容需要下载文档才能查看 内容需要下载文档才能查看MA
将ψ(t)∈
内容需要下载文档才能查看L2(R)
(2)
时称ψ
内容需要下载文档才能查看(t)为一个基本小波函数或母小波函数
对于连续的情
内容需要下载文档才能查看况
(3)
内容需要下载文档才能查看其中
b
内容需要下载文档才能查看为平移因子
小波
ψ
内容需要下载文档才能查看j/2j,k(t)=2?ψ(2?jt?k)
ÔÚʵ¼ÊÓ¦ÓÃ
ÖÐÐèÒ
内容需要下载文档才能查看ªÑ¡ÔñС²¨»ù
»¹ÐèÍâÐÎÓë±»·ÖÎöµÄ¶ÔÏóÓÐÒ»
¶¨µÄÏàËÆÐÔ
内容需要下载文档才能查看[16,17]
CWT
任给函数f
内容需要下载文档才能查看(t)∈L2
(R)
(5)
其中?f,ψa,b?为f(t)与ψ
内容需要下载文档才能查看a,b之内积
万方数据
f(t)的DWT系数
C∞
*j,k=∫
?∞
f(t)ψ
内容需要下载文档才能查看j,k(t)dt
³ß¶ÈaÊǰ´¶þ½øÖÆ
ÀëÉ¢µÄ
ÓëCWT相比
Empirical Mode
Decomposition
法和Hilbert变换两部分组成
下包络线的平均值求
内容需要下载文档才能查看再提取固有模态函数IMF必须满足2个条件
内容需要下载文档才能查看1
内容需要下载文档才能查看2
内容需要下载文档才能查看内容需要下载文档才能查看 内容需要下载文档才能查看
ͬÀíµÃµ½Ï°üÂçÏß
内容需要下载文档才能查看3
定义h
内容需要下载文档才能查看11
若经过k次筛选后的结果
h1k满足IMF
内容需要下载文档才能查看的定义
内容需要下载文档才能查看内容需要下载文档才能查看
ÖØ¸´ÒÔÉÏɸѡ¹ý³Ìc3
那么x(t)就可分解为n个IMF分量及rn的和 n
即 x(t)=
∑ci
(t)+rn
(t)
i=
内容需要下载文档才能查看1
¶Ôÿ¸öIMF
分量进行Hilbert
内容需要下载文档才能查看变换
综合所有IMF分量的瞬时频谱就可获得Hilbert
内容需要下载文档才能查看谱
(8)
幅值函数
a(t)=
内容需要下载文档才能查看c2(t)+H2[c(t)]
(10)
在此基础上求出瞬时频率
1474 岩 土 工 程 学 报 2005年
1dΦ(t)
2πdt
¾-Hilbert变换得到的震幅和频率都是f(t)=
内容需要下载文档才能查看ÒÀ´Î·Ö½â³öµÍƵ
解是信号本身所决定的一个自适应分解过程
ÔÙ
内容需要下载文档才能查看EMD分
时间的函数
表明
内容需要下载文档才能查看了分解的可靠性
(12)
式中除去了余差r
内容需要下载文档才能查看n
¾ÍµÃµ½
Hilbert边际谱
T
h(ω
内容需要下载文档才能查看)=∫H(ω,t)dt
内容需要下载文档才能查看
»òÄÜÁ¿
¿ÉÒÔ¶¨ÒåHilbert瞬时能量
如下:
IE
内容需要下载文档才能查看(t)=
∫H2(ω,t)dω
ω
Èç¹ûÕð
内容需要下载文档才能查看·ùµÄƽ·½¶Ôʱ¼ä»ý·Ö
(15)
Hilbert
内容需要下载文档才能查看能量谱提供了每个频率的能量计算式
利用四个余弦信号与
内容需要下载文档才能查看一个指数趋势项相加作为仿真信号
(16)
内容需要下载文档才能查看
图1 仿真信号与分解结果
Fig. 1 Simulative signal and EMD components
对该信号
其分解结果得到的4个分量c1
c4和一个趋势项c
内容需要下载文档才能查看5
内容需要下载文档才能查看b
万方数据
内容需要下载文档才能查看
从上看出HHT方法
内容需要下载文档才能查看是一种更具适应性的时频局部化分析方法
首次给出了IMF的定义
使信号分析更加灵活多变能够表达信号的局部特征
瞬时频率定义为相位函数的导数,
不需要整个波来定义局部频率
这比小波变换有了明显的进
步
内容需要下载文档才能查看[18]
(b)中的fft代表FFT所得的频谱)
内容需要下载文档才能查看内容需要下载文档才能查看
Ö÷ÒªÊܱ¬ÆÆ²ÎÊý
内容需要下载文档才能查看2
µ«Ã»ÓÐÈκÎʱ
内容需要下载文档才能查看ÓòÐÅÏ¢
图2 爆破震动信号及其频谱
Fig. 2 Blasting vibration signal and its spectrum
第12期 张义平
ʹÐźŵÄÄÜÁ¿¨CƵÂʨCʱ¼ä·Ö²¼ºÜÄÑ
(b)
内容需要下载文档才能查看为时频谱)
内容需要下载文档才能查看
Ò²²»ÄÜ
内容需要下载文档才能查看·´Ó³ÖʵãÕð¶¯ÄÜÁ¿ËæÊ±¼äµÄ±ä»¯×´¿ö
2
ƵÂÊÖ÷Òª¼¯ÖÐÔÚ0.01的
归一化频率线附近
时间分辨
率和频率分辨率需要折衷
分辨率单一
内容需要下载文档才能查看[15]
Òò´ËSTFT不适合类似爆破震
内容需要下载文档才能查看动信号的时频分析
相应的分解结果如
图4
内容需要下载文档才能查看
ÔÚʱÓòÉÏÐźÅÄÜÁ¿·Ö²¼ÔÚ100
内容需要下载文档才能查看ÓÖÖ÷Òª¼¯ÖÐÔÚ100
ÆäÖÐ100点和600
内容需要下载文档才能查看点处表现更为明显
CWT
内容需要下载文档才能查看存在一定的冗余
内容需要下载文档才能查看ÆäÖÐa8为低频分量
信号
的能量在时域上主要分布在1000采样点的范围内在频域上主要集中在
d4
内容需要下载文档才能查看内容需要下载文档才能查看
½Ï¾ßÌ嵨ÃèÊöÁË·ÇÆ½ÎÈÐźŵÄʱƵ·Ö²¼
内容需要下载文档才能查看ͬʱ
¶¨Á¿¸ø³ö
²»Í¬µÄС²¨»ù»áÓв»Í¬µÄ½á
内容需要下载文档才能查看¹û
ÔÚÕû
内容需要下载文档才能查看¸öÐźŷÖÎö¹ý³ÌÖÐÖ»ÄܲÉÓÃÕâͬһ¸öС²¨»ù
图4 信号的db3小波谱
内容需要下载文档才能查看Fig. 4 Wavelet spectrum of signal based on db3
图5 信号的DWT分量
Fig. 5 DWT components of signal
4.3 信号的HHT分析
利用EMD分解图2中的原信号后得到如图6所示的分量
内容需要下载文档才能查看15
和式
分别见图7
b所示
内容需要下载文档才能查看
·Ö½âËùµÃµÄÿһ¸öIMF分量都有不同
的震幅和频率
内容需要下载文档才能查看˲ʱÄÜÁ¿Æ×¸üÖ±¹ÛµØ±íÃ÷ÁËÄÜÁ¿Ö÷Òª
¼¯ÖÐÔÚ100
1200采样点区间内
1476 岩 土 工 程 学 报 2005年
际谱可看出能量集中在10 Hz
内容需要下载文档才能查看以内的范围
内容需要下载文档才能查看3
而边际
谱中某一频率仅代表有这样频率的信号存在的可能性
内容需要下载文档才能查看[10]
内容需要下载文档才能查看
WT和HHT都能表现
内容需要下载文档才能查看出良好的局部化特征
且主要集中在有
内容需要下载文档才能查看限的时频谱线上
具有较宽的频谱分布
是由
内容需要下载文档才能查看于小波基有限长并受到测不准原理的限制
从而发生能量泄漏
HHT法中也存在端点效
应等问题
内容需要下载文档才能查看[20]
内容需要下载文档才能查看Æä±ê×¼»ùÊÇÓɵȷù
ÔÚFourier变换中会产生能量泄漏和高频
内容需要下载文档才能查看放大
内容需要下载文档才能查看无法满足时频分析的要求
对爆破震动信号给出
了时频描述
对于爆破震动这类具
有复杂
要求信号在剧变时刻有较高的时间分辨率
STFT
内容需要下载文档才能查看的时间分辨率和频率分辨率面临折衷考虑
所以STFT无法两者
内容需要下载文档才能查看兼顾
(3) WT是一种以小波基为变换基础的多分辨率时频分析方法
而STFT分析得出的是较平缓的
时间衰减规律
是一种比较适合非平稳信号
内容需要下载文档才能查看分析的时频方法
小波谱的能量在频率范围内分布
较宽
EMD
内容需要下载文档才能查看依据信号本身的固有特性进行分解
具有自适应性强和高效
内容需要下载文档才能查看的优点
HHT较WT
内容需要下载文档才能查看更具适应性
有利于更好地指导爆破设计与研究
[1] Luy, Hao H, Ma G, et al. Simulation of structural response
under high frequency ground excitation[J]. Earthquake Engng Struct Dyn, 2001, 30(2): 307–325.
[2] Dowding C H. Monitoring and control of blasting effects[A].
Ming Engineering Handbook[M]. Prentice Hall. 1985:746–760. [3] 何 军
爆破震动信号的小波分析
内容需要下载文档才能查看[J]
下载文档
热门试卷
- 2016年四川省内江市中考化学试卷
- 广西钦州市高新区2017届高三11月月考政治试卷
- 浙江省湖州市2016-2017学年高一上学期期中考试政治试卷
- 浙江省湖州市2016-2017学年高二上学期期中考试政治试卷
- 辽宁省铁岭市协作体2017届高三上学期第三次联考政治试卷
- 广西钦州市钦州港区2016-2017学年高二11月月考政治试卷
- 广西钦州市钦州港区2017届高三11月月考政治试卷
- 广西钦州市钦州港区2016-2017学年高一11月月考政治试卷
- 广西钦州市高新区2016-2017学年高二11月月考政治试卷
- 广西钦州市高新区2016-2017学年高一11月月考政治试卷
- 山东省滨州市三校2017届第一学期阶段测试初三英语试题
- 四川省成都七中2017届高三一诊模拟考试文科综合试卷
- 2017届普通高等学校招生全国统一考试模拟试题(附答案)
- 重庆市永川中学高2017级上期12月月考语文试题
- 江西宜春三中2017届高三第一学期第二次月考文科综合试题
- 内蒙古赤峰二中2017届高三上学期第三次月考英语试题
- 2017年六年级(上)数学期末考试卷
- 2017人教版小学英语三年级上期末笔试题
- 江苏省常州西藏民族中学2016-2017学年九年级思想品德第一学期第二次阶段测试试卷
- 重庆市九龙坡区七校2016-2017学年上期八年级素质测查(二)语文学科试题卷
- 江苏省无锡市钱桥中学2016年12月八年级语文阶段性测试卷
- 江苏省无锡市钱桥中学2016-2017学年七年级英语12月阶段检测试卷
- 山东省邹城市第八中学2016-2017学年八年级12月物理第4章试题(无答案)
- 【人教版】河北省2015-2016学年度九年级上期末语文试题卷(附答案)
- 四川省简阳市阳安中学2016年12月高二月考英语试卷
- 四川省成都龙泉中学高三上学期2016年12月月考试题文科综合能力测试
- 安徽省滁州中学2016—2017学年度第一学期12月月考高三英语试卷
- 山东省武城县第二中学2016.12高一年级上学期第二次月考历史试题(必修一第四、五单元)
- 福建省四地六校联考2016-2017学年上学期第三次月考高三化学试卷
- 甘肃省武威第二十三中学2016—2017学年度八年级第一学期12月月考生物试卷
网友关注
- 民俗与旅游
- 餐饮服务与管理教学大纲
- 第一章艺术的起源`本质特征与种类
- 数字艺术概论4
- oracle物理设计原则
- 生物化学笔记
- 金属材料与热处理 金属材料基础知识介绍
- FPC的基础入门
- 生物化学试卷及答案
- 市场分析
- 第四讲 中世纪文明[资料]
- [宝典]数字电搅扰路抗设计绳尺
- 理论力学重修第二讲
- 健康管理的历史、现状和挑战
- [中医中药]中医养生学
- 1养生保健概论
- 12DNA的生物合成
- 薪酬管理
- 生物地理学课程讲义 (可编辑)
- 【精品】安全工程、工程经济类专业参考目录23
- 13复杂电路
- 薪酬管理
- 数字媒体艺术概论_李四达编著2
- 普通车床操纵实训引诱书
- 第四章 市场营销环境
- 徐州医学院 儿童营养和儿童保健原则课件
- 《社区护理(第二版)》第四章社区重点人群保健
- 大学生必备资料
- 基础写作
- 建立病原微生物实验室生物安体系
网友关注视频
- 七年级英语下册 上海牛津版 Unit3
- 沪教版牛津小学英语(深圳用) 六年级下册 Unit 7
- 【部编】人教版语文七年级下册《逢入京使》优质课教学视频+PPT课件+教案,安徽省
- 沪教版牛津小学英语(深圳用) 四年级下册 Unit 2
- 19 爱护鸟类_第一课时(二等奖)(桂美版二年级下册)_T502436
- 第19课 我喜欢的鸟_第一课时(二等奖)(人美杨永善版二年级下册)_T644386
- 青岛版教材五年级下册第四单元(走进军营——方向与位置)用数对确定位置(一等奖)
- 小学英语单词
- 冀教版小学英语五年级下册lesson2教学视频(2)
- 冀教版英语五年级下册第二课课程解读
- 【部编】人教版语文七年级下册《过松源晨炊漆公店(其五)》优质课教学视频+PPT课件+教案,江苏省
- 七年级英语下册 上海牛津版 Unit5
- 冀教版小学数学二年级下册第二周第2课时《我们的测量》宝丰街小学庞志荣
- 三年级英语单词记忆下册(沪教版)第一二单元复习
- 六年级英语下册上海牛津版教材讲解 U1单词
- 冀教版小学数学二年级下册第二单元《租船问题》
- 精品·同步课程 历史 八年级 上册 第15集 近代科学技术与思想文化
- 8.对剪花样_第一课时(二等奖)(冀美版二年级上册)_T515402
- 《空中课堂》二年级下册 数学第一单元第1课时
- 沪教版牛津小学英语(深圳用) 四年级下册 Unit 12
- 二年级下册数学第二课
- 冀教版英语四年级下册第二课
- 外研版英语三起5年级下册(14版)Module3 Unit1
- 北师大版小学数学四年级下册第15课小数乘小数一
- 【部编】人教版语文七年级下册《逢入京使》优质课教学视频+PPT课件+教案,辽宁省
- 冀教版小学英语四年级下册Lesson2授课视频
- 苏科版数学 八年级下册 第八章第二节 可能性的大小
- 冀教版小学数学二年级下册第二单元《有余数除法的简单应用》
- 二次函数求实际问题中的最值_第一课时(特等奖)(冀教版九年级下册)_T144339
- 苏科版八年级数学下册7.2《统计图的选用》
精品推荐
- 2016-2017学年高一语文人教版必修一+模块学业水平检测试题(含答案)
- 广西钦州市高新区2017届高三11月月考政治试卷
- 浙江省湖州市2016-2017学年高一上学期期中考试政治试卷
- 浙江省湖州市2016-2017学年高二上学期期中考试政治试卷
- 辽宁省铁岭市协作体2017届高三上学期第三次联考政治试卷
- 广西钦州市钦州港区2016-2017学年高二11月月考政治试卷
- 广西钦州市钦州港区2017届高三11月月考政治试卷
- 广西钦州市钦州港区2016-2017学年高一11月月考政治试卷
- 广西钦州市高新区2016-2017学年高二11月月考政治试卷
- 广西钦州市高新区2016-2017学年高一11月月考政治试卷
分类导航
- 互联网
- 电脑基础知识
- 计算机软件及应用
- 计算机硬件及网络
- 计算机应用/办公自动化
- .NET
- 数据结构与算法
- Java
- SEO
- C/C++资料
- linux/Unix相关
- 手机开发
- UML理论/建模
- 并行计算/云计算
- 嵌入式开发
- windows相关
- 软件工程
- 管理信息系统
- 开发文档
- 图形图像
- 网络与通信
- 网络信息安全
- 电子支付
- Labview
- matlab
- 网络资源
- Python
- Delphi/Perl
- 评测
- Flash/Flex
- CSS/Script
- 计算机原理
- PHP资料
- 数据挖掘与模式识别
- Web服务
- 数据库
- Visual Basic
- 电子商务
- 服务器
- 搜索引擎优化
- 存储
- 架构
- 行业软件
- 人工智能
- 计算机辅助设计
- 多媒体
- 软件测试
- 计算机硬件与维护
- 网站策划/UE
- 网页设计/UI
- 网吧管理