现代信号处理第2章信号的时域分析.ppt

2020 2 18 机械工程及自动化研究所 现代信号处理技术及应用 第二章信号的时域分析 西安交通大学机械工程学院研究生学位课程 第二章信号的时域分析 2 1信号的预处理2 2信号的采样2 3时域统计分析2 4相关分析及应用 引言 以时间为自变量描述物理量的变化是信号最基本 最直观的表达形式 在时域内对信号进行滤波 放大 统计特征计算 相关性分析等处理 统称为信号的时域分析 通过时域分析方法 可以有效提高信噪比 求取信号波形在不同时刻的相似性和关联性 获得反映机械设备运行状态的特征参数 为机械系统动态分析和故障诊断提供有效的信息 第二章信号的时域分析 2 1信号的预处理2 2信号的采样2 3时域统计分析2 4相关分析及应用 2 1信号的预处理 传感器获取的信号往往比较微弱 并伴随着各种噪声 不同类型的传感器 其输出信号的形式也不尽相同 为了抑制信号中的噪声 提高检测信号的信噪比 便于信息提取 须对传感器检测到的信号进行预处理 所谓信号预处理 是指在对信号进行变换 提取 识别或评估之前 对检测信号进行的转换 滤波 放大等处理 2 1信号的预处理 常用的信号预处理方法信号类型转换应变测力传感器 热电阻传感器输出的信号均为电阻信号 为了便于后续处理常用电桥将电阻信号转变为电压信号信号放大常用的信号放大器包括 测量放大器 隔离放大器 可编程增益放大器等信号滤波 本节重点介绍 去除均值在计算信号的标准差等统计量时 需要去除信号均值去除趋势项常用的趋势项消除方法有滤波法 多项式拟合法 2 1 1信号的滤波处理 信号滤波处理是消除或减弱干扰噪声 保留有用信号的过程 把实现滤波功能的系统称之为滤波器 滤波器可分为两大类 即经典滤波器和现代滤波器 1 经典滤波器 定义 当噪声和有用信号处于不同的频带时 噪声通过滤波器将被衰减或消除 而有用信号得以保留分类根据幅频特性的不同 滤波器分为低通滤波器 高通滤波器 带通滤波器 带阻滤波器等类型 根据处理信号类型的不同 滤波器可分为模拟滤波器和数字滤波器 对于数字滤波器来说 根据滤波器的单位脉冲响应序列长度的无限和有限 数字滤波器可进一步分为无限冲击响应滤波器 IIR 和有限冲击响应滤波器 FIR 两类 1 经典滤波器原理 经典滤波概念和方法建立在频域分析基础上 滤波后的信号为 滤波器的传递函数或滤波器的频率响应函数 在噪声频带和有用成分频带分离的情况下 通过设计如下的滤波器函数 滤波器的单位脉冲响应函数 2 1 1 2 理想模拟滤波器 理想模拟滤波器是一个理想化的模型 对其讨论有助于进一步了解和改进实际滤波器的性能理想模拟滤波器的幅频特性曲线 图2 1 1 理想模拟滤波器 由于理想低通滤波器具有矩形幅频特性和线性相位特性 同时 理想高通 带通和带阻滤波器均可以由理想低通滤波器串联得到 因此 以后均以理想低通滤波器为例来说明理想低通滤波器的矩形幅频 相频特性 2 1 4 理想模拟滤波器 理想低通滤波器的单位脉冲响应函数 2 1 5 3 实际滤波器及其基本参数 实际的滤波器为了物理上可实现 通常在通带和阻带之间设置过渡带 通带边缘频率 阻带边缘频率 截止频率 滤波器的半功率点 实际滤波器及其基本参数 实际滤波器的幅频特性幅值在通带和阻带内一般不严格为1和0 它们分别允许的波动量分别为 波动的大小分别用通带和阻带内的衰减 实际滤波器的参数还有 波纹幅度 带宽 品质因数和倍频程选择性等 4 数字滤波器的设计 数字滤波器有无限冲击响应IIR型滤波器和有限冲击响应FIR型滤波器之分IIR型数字滤波器的传递函数是 FIR型数字滤波器的传递函数是 滤波器的单位脉冲响应函数 数字滤波器的设计方法简述 FIR型数字滤波器的设计方法主要是建立在对理想滤波器频率特性作某种近似的基础上的 这些近似方法有窗函数法 频率抽样法等 IIR型数字滤波器的设计属于间接设计法 IIR型数字滤波器目前最通用的设计方法是利用已经很成熟的模拟滤波器的设计方法来进行设计 而模拟滤波器的设计方法又有巴特沃斯 Butterworth 滤波器 切比雪夫 Chebyshev 和椭圆滤波器等不同的设计方法 数字滤波器的设计方法简述 不论是FIR型数字滤波器还是IIR型数字滤波器的设计都包括三个步骤 Matlab数字滤波器设计演示 FDATool fda SPTool fda 2 现代滤波器 当噪声频带和有用信号频带相互重叠时 经典滤波器就无法实现滤波功能现代滤波器也称统计滤波器 从统计的概念出发对信号在时域进行估计 在统计指标最优的意义下 用估计值去逼近有用信号 相应的噪声也在统计最优的意义下得以减弱或消除常用的统计滤波器有维纳滤波器和卡尔曼滤波器两类 1 维纳滤波器 20世纪40年代第二次世界大战期间 由于军事上的需要 Wiener提出并解决了平稳过程的最佳线性滤波问题采用线性最小均方误差估计准则 设计的最佳滤波器称为维纳 Wiener 滤波器 h t 脉冲响应函数 按最小均方误差准则确定 2 1 8 2 1 9 维纳滤波器 维纳滤波器可根据t时刻 及t以前时刻的观测值x t 实现以下三个方面的应用 2 卡尔曼滤波器 维纳滤波器由于计算量大 难以作实时处理 故不能广泛应用 同时它对非平稳信号的滤波也无能为力 60年代初由于航天事业发展的需要 卡尔曼 Kalman 和布西 Bucy 在解决非平稳 多输入输出随机序列的估计问题中引入了状态变量 在克服维纳滤波某些局限的基础上 提出了被后人称为卡尔曼滤波的新滤波方法 该方法在雷达 通信 控制 生物和勘探等领域得到了广泛的应用 卡尔曼滤波器 卡尔曼滤波是线性最小均方误差滤波器的另一种处理方法 卡尔曼滤波建立在已知随机信号模型的基础上 包括模型的阶次 模型的参数和激励白噪声的统计特性 其原理是用信号前一时刻的估计值与测量值的误差项的加权平均作为当前时刻的估计值 2 1 10 信号前一时刻估计值 测量值的误差项 常数a和c分别表示参数模型和测量模型的参数 用Gn表示n时刻的加权系数 维纳滤波器与卡尔曼滤波器对比 卡尔曼滤波的特点是把信号的先验知识用信号的模型表达出来 在时域上引入状态变量法进行处理 采用递推型的线性最小均方误差算法 卡尔曼滤波和维纳滤波都是在应用随机信号和观测噪声的前二阶矩的统计特性 以线性最小均方估计解决随机信号的滤波问题 维纳滤波需要给出随机信号和噪声的有理谱形式 卡尔曼滤波则要求把随机信号规定为白噪声驱动的线性系统的输出 维纳滤波理论适应于平稳随机过程 卡尔曼滤波适用于有限初始时间的非平稳随机过程 第二章信号的时域分析 2 1信号的预处理2 2信号的采样2 3时域统计分析2 4相关分析及应用 2 2信号的采样 将连续信号转换成离散的数字序列过程就是信号的采样 它包含了离散和量化两个主要步骤 本节主要介绍采样过程中采样与混频 量化与误差 采样长度与分辨率及窗函数与泄露四方面的内容 2 2 1采样与混频 设模拟信号为x t 间距为 t的采样脉冲函数为p t 采样过程是x t 和p t 相乘 得到离散信号x n t 设x t 的傅里叶变换为X 采样脉冲函数p t 的傅里叶变换为 2 2 1 2 2 2 模拟信号 采样脉冲函数及其频谱 离散信号及其频谱 根据频域卷积定理可知 则式 2 2 1 所示离散信号x n t 的傅里叶变换为 离散信号及其频谱 频谱混叠 不能由频谱准确地恢复原信号 2 2 3 采样定理 采样定理 为避免混叠 采样频率 s必须不小于信号中最高频率 max的两倍 即或 实际中采样频率的选取往往留有余地 一般选取采样频率 s为处理信号中最高频率的2 5 4倍 另外 由于测量信号中的高频部分往往是由干扰引起的噪声或我们不感兴趣的频谱 因此采样前须先对信号进行低通滤波 又称抗混滤波 然后再根据滤波后信号的最高频率 max设定采样间隔 p28 2 2 2量化与误差 量化是对信号采样点取值进行数字化转换的过程 量化结果以一定位数的数字近似表示信号在采样点的取值 由于模 数转换器的位数有限 模 数转换器的输出只能表达一系列具有一定间隔的电平 当模拟信号在采样点上的取值落在两个相邻电平之间时 就要舍入到相近的一个电平上 我们把这一过程称之为量化 若设模 数转换器的位数为N 采用二进制编码 转换器转换的电压范围为 V 则相邻电平之间的增量为 2 2 4 2 2 3窗函数与泄漏 理论上信号的长度是无限的 但任何观测信号都是在有限时间段内进行观测的 因此 信号采样过程须使用窗函数 将无限长信号截断成为有限长度的信号 从理论上看 截断过程就是在时域将无限长信号乘以有限时间宽度的窗函数 窗函数与泄漏 最简单的窗函数是矩形窗 无限带宽泄漏 窗函数与泄漏 泄漏与截断长度 所使用的窗函数等有关 不进行信号截断就没有泄漏误差 另外 使用不同的窗函数泄漏大小也不同 泄漏取决于窗函数频谱的旁瓣 如果窗函数的旁瓣小 相应的泄漏也小 其它窗函数 三角窗 其它窗函数 汉宁窗 2 2 4采样长度与分辩率 数字信号的分辨率包括时间分辨率和频率分辨率 数字信号的时间分辨率即采样间隔 t 它反映了数字信号在时域中取值点之间的细密程度 数字信号的频率分辨率为 2 T 其中T N t为数字信号的时间跨度 N为数字信号的长度 频率分辨率表示了数字信号的频谱在频域中取值点之间的细密程度 因此 当采样频率或采样间隔确定后 增大采样点数就可增加信号的时间长度和频率分辨率 第二章信号的时域分析 2 1信号的预处理2 2信号的采样2 3时域统计分析2 4相关分析及应用 2 3时域统计分析 信号的时域统计分析是指对信号的各种时域参数 指标的估计或计算 常用的时域参数和指标包括 1 均值 2 均方值 3 均方根值 4 方差 5 标准差 6 概率密度函数 7 概率分布函数 8 联合概率密度函数等 本节先介绍常见参数的概念 然后给出它们的应用 2 3 1时域指标参数 1 均值 2 均方值 方差 P32 时域指标参数 3 概率密度函数 随机信号的取值落在区间内的概率可用下式表示 概率密度函数定义为 概率分布函数的定义为 时域指标参数 4 有量纲参数指标 有量纲参数指标包括方根幅值 平均幅值 均方幅值和峰值四种 若随机信号符合平稳 各态历经条件 且均值为零 概率密度函数为p x 则有量纲参数指标的定义如下 方根幅值 平均幅值 均方幅值 峰值 时域指标参数 上述的有量纲参数指标也可在时域定义 p33 时域指标参数 5 无量纲参数指标有量纲参数指标不但与机器的状态有关 且与机器的运动参数如转速 载荷等有关 而无量纲参数指标具有对信号幅值和频率变化均不敏感的特点 这就意味着理论上它们与机器的运动条件无关 只依赖于概率密率函数的形状 所以无量纲参数指标是一种较好的机器状态监测诊断参数 无量纲参数指标包括了波形指标 峰值指标 脉冲指标和裕度指标 定义为 2 3 10 无量纲参数指标 当上式中的l和m取不同值时 就得到如下指标 无量纲参数指标 偏斜度指标S表示信号概率密度函数的中心偏离正态分布的程度 反映信号幅值分布相对其均值的不对称性 峭度指标K则表示信号概率密度函数峰顶的陡峭程度 反映信号波形中的冲击分量的大小 2 3 2参数指标的应用 1 利用概率密度函数和概率分布函数进行产品质量控制 研究材料的强度和控制设备的工作稳定性当图 a 代表一批零件的加工尺寸 则根据图 b 图 c 可判断加工过程的质量高低 进而可评价或判断机床工具是否应该调整 操作工人的技术熟练程度等 参数指标的应用 概率密度函数用于机器状态判断新变速箱噪声的概率密度曲线如图2 3 3 a 所示 旧变速箱噪声的概率密度曲线如图2 3 3 b 所示 参数指标的应用 机器状态的时域参数指标判断方法 诊断参数指标一般应满足如下要求 易于测量和计算 所需计算机存储量小 能敏锐地反映和预报机器的早期故障 不受机器运行状态 如负载 转速等变化的影响 能够指示故障的存在 以便及时排查故障 参数指标的应用 图2 3 4表示了28只汽车后桥齿轮在不同运行状态下 由振动加速度信号计算得到的无量纲参数指标 波形指标K的变化很小 没有足够的诊断能力 而峰值指标C和脉冲指标I可以作为齿轮运行状态的优良诊断指标 第二章信号的时域分析 2 1信号的预处理2 2信号的采样2 3时域统计分析2 4相关分析及应用 2 4 1相关的概念 所谓相关 就是指变量之间的线性联系或相互依赖关系 根据前面的讨论 变量之间的联系可通过反映变量的信号之间的内积或投影大小来刻画 设有实信号x t 和y t 它们的内积可写成 如果信号和随自变量时间的取值相似 内积结果就大 反之亦然 因此 通过式 2 4 1 可定义信号和的相关性度量指标 2 4 1 2 4 2 2 4 2自相关函数及其应用 信号x t 的自相关函数和自相关系数定义为 自相关函数及其应用 几种常见信号的自相关函数 自相关函数及其应用 几种常见信号的自相关函数 自相关函数及其应用 可以看出信号中的周期性分量在相应的自相关函数中不会衰减 且保持了原来的周期 因此 自相关函数可从被噪声干扰的信号中找出周期成分 在用噪声诊断机器运行状态时 正常机器噪声是由大量 无序 大小近似相等的随机成分叠加的结果 因此正常机器噪声具有较宽而均匀的频谱 当机器状态异常时 随机噪声中将出现有规则 周期性的信号 其幅度要比正常噪声