52phm社区,工业互联网人的技术交流社区

2022-01-13 17:40:13    互联网    130    当前专栏:数字信号处理    分类:算法开发    本站官网:www.52phm.cn   

扫码进群 ...

【数字信号处理】--数字信号分析与处理基础知识

文章目录

  • 一、信号的时域分析
    • [1.1信号波形参数识别]
    • [1.2检测方法 过零检测法]
    • [1.3数字信号微分与数字信号积分]
  • 二、信号的频域分析
    • [2.1周期信号的频谱分析]
    • [2.2 信号的频谱分析]
    • [2.3数字信号的频谱计算方法]
  • 三、信号的时差域相关分析
    • [3.1信号的相关系数]
    • [3.2 相关应用]
    • [3.3 数字滤波器和模拟滤波器的区别]
  • 四、信号的幅值域分析
    • [4.1概率密度曲线与概率分布曲线的应用]
  • 五、信号的数字滤波
    • [5.1滤波器概念]
    • [5.2频率域滤波:]
    • [5.3时域滤波]
    • [5.4 时域FIR滤波器]
    • [5.5 IIR滤波器]
    • [5.6 其他滤波器]
    • [5.7 数字音乐均衡器的设计]
  • 六、信号的时频分析简介
    • [6.1 信号小波分析]
    • [6.2四种分析方法时频分析精度对比]
    • [6.3连续小波]
    • [6.5离散小波变换]
    • [6.6几种小波变换的关系]

一、信号的时域分析

采样定理:F>2Fs; F>5Fs

频率混迭:采样频率不满足采样定理时,信号中的高频成分会被不正确的采样成低频成分,这一现象称为频率混迭。

1.1信号波形参数识别

  • 周期、频率
  • 峰值、双峰值:max(singal)- min(singal)
  • 初相位
  • 均值:反映信号变化的中心趋势,称为直流分量:mean(singal)
  • 均方值:表示信号强度,其平方根为有效值(RMS)
  • 方差:反映信号绕均值的波动程度:std(singal)

1.2检测方法 过零检测法

过零点有最大的斜率

上升沿过零点条件:x(k)<=0<x(k1)

相邻两个过零点时间差就是周期,第一个过零点位置与周期比例是初相位。

峰值检测原理:目前峰值检测方法的核心是在某一个或某几个尺度内搜索小波变换模极大极小值之间的过零点,在此基础上对过零点位置相对于峰值位置的误差进行一定程序的修正。

过零检测MATLAB代码:

2020-04-09_120206

1.3数字信号微分与数字信号积分

数字信号先微分再积分,可还原信号,但与原信号相比,丢失初始值信息。

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-gFbyBQLL-1590892967028)(https://s1.ax1x.com/2020/04/09/G4gpGQ.png)]

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-eW2Fx5dp-1590892967031)(https://s1.ax1x.com/2020/04/09/G4g92j.png)]

利用信号的峰值进行索道钢丝断裂检测

二、信号的频域分析

2.1周期信号的频谱分析

x(t)=x(tnT);基频f0=1/T。该信号可以用傅里叶级数分解。

信号的正交分解与合成:



利用三角正交函数集对信号进行分解。
三角函数集{cos(2×pi×i×f0×t),sin(2×pi×i×f0×t),i=1,2,3,4,5} 中任意两个函数相乘,积分值为零。



信号的傅里叶级数分解:

周期信号在三角函数集下的投影:x(t)=x(tnT) => {cos(2×pi×i×f0×t),sin(2×pi×i×f0×t)}





2.2 信号的频谱分析

频谱图



周期信号频谱图的特点:

  • 离散性
  • 谐波性
  • 收敛性

连续域上的积分求解映射到离散域上的求和操作。

MATLAB中使用yn=fft(x,n)函数,数据长度n必须是2的幂次方,结果yn是投影系数实部an和虚部bn组成。
Yn=anj*bn,n=1,2,3…Y0表示直流分量,Y1表示一次谐波f=f0时的值,Y2表示二次谐波f=2f0时的值.

fft注意事项

  • 不显示负频率部分 plot(f,A(1:N/2));
  • X坐标换为频率 f=linspace(0,fs/2,N/2);
  • 幅值量纲还原 A=abs(y)/(N/2); #还原幅频谱图的量纲
  • 幅频谱图:A1=abs(y)/(N/2);
  • 功率谱图:A2=A1.^2;
  • 对数功率谱图:P2=20*log10(A2);

2.3数字信号的频谱计算方法

1.信号截断:对信号进行截断,看做以W长度为大小的周期信号进行处理。由于信号截断的原因,产生了能量泄漏,即使信号频率与频谱离散取样点不相等,也能得到该频率分量的一个近似值。
信号整周期截断。

2.能量泄漏



  • 信号->信号截断、周期延拓–>离散周期信号
  • 信号截断–>能力泄漏- FFT–>栅栏效应;
    加窗–>控制能量泄漏–>减小栅栏效应误差;

信号截断函数:能量泄漏分主瓣泄漏和旁瓣泄漏,信号加窗截断是用窗函数与信号相乘,然后截断,减小周期延拓边界的跳变;抑制旁瓣能量泄漏。信号加窗减小栅栏效用。

平顶窗用的多

3.离散傅里叶变换





按照三角函数的展开式求出a0 an bn这三个傅里叶级数的系数

数字信号->信号截断->周期延拓->截断的周期信号

4.FFT

三、信号的时差域相关分析

3.1信号的相关系数

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-2RjyHiIk-1590892967051)(https://s1.ax1x.com/2020/04/14/GzAKCq.png)]
工程上通常只关心分子部分。使用FFT计算自相关函数

相关函数的性质:

  • 自相关函数是偶函数;
  • 当t=0,自相关函数具有最大值;
  • 周期信号的自相关函数仍然是同频率的周期信号,但不保留相位信息;
  • 两个非同频率的周期信号互不相关
% 自相关滤波
N=1024;
T=0.2;
x=linspace(0,T,N);
y=sin(2*pi*50*x);
figure;
plot(x,y);
s=xcorr(y,'unbiased'); %加修正的相关系数
x1=linspace(-T,T,2*N-1);
figure;
plot(x1,s);

3.2 相关应用

  • 声波传播速度测量
  • 超声波探伤
  • 地下输油管道漏损位置的探测
  • 自相关滤波,对信号进行自相关去噪

3.3 数字滤波器和模拟滤波器的区别

- 影响采样频率;

- 影响有效A/D位数;

  • 数字滤波器由数字乘法器、加法器和延时单元组成的一种算法或装置。
  • 模拟滤波器有有源和无源的,有源滤波器主要是有运放,或者跨到运放,及电阻,电容构成。
  • 从实现手段上看,模拟滤波器一般用电容,电感这些模拟器件搭建的,数字滤波器可以通过软件或者数字芯片来实现。
  • 模滤用H(S)表示,而数滤用H(Z)表示.模滤是以幅频特性的逼近为主要,而数字滤波器则可以实现相位的匹配;
  • 数字滤波器有两类:IIR,FIR,IIR主要有脉冲响应不变法与双线性变换法两种,而FIR可以实现线性相位,有窗口法,频率采样法,最优化方法.可以用FFT来实现.

四、信号的幅值域分析

幅值域分析称为信号的概率密度曲线和概率分布曲线分析。





MATLAB相关函数:

  • hist() %计算直方图
  • pdf=hist(y)/length(y) %
  • cdfplot() %计算概率分布函数

4.1概率密度曲线与概率分布曲线的应用

  • 齿轮故障诊断
  • 水泵故障诊断

五、信号的数字滤波

滤波的条件:有用信号与噪声信号在频率域上(频带)不重叠。

当采样信号中有较高频率成分信号时,首先使用模拟滤波器进行低通滤波,再进行信号的预处理。这样可以降低信号的采样频率,减小运算量 。

在信号去噪处理中,判断去噪效果的好坏主 要有信噪比 (SNR) 和信号重构均方误差(MSE) 2 个标准。一般来说, SNR 越大, MSE 越小,则去 噪效果越好。 均方根误差的定义为:

5.1滤波器概念





滤波器评价指标:

  • 过渡区陡峭
  • 脉冲区应紧支集:滤波器做傅里叶变换后,宽度应尽可能窄

5.2频率域滤波:







长数据信号的FFT频域滤波:对长信号进行帧划分,比如划成1024点,再进行FFT滤波;再采样帧数据重叠技术,对重迭部分进行掐头去尾处理,比如舍去64点。

5.3时域滤波




时域数字滤波就是将信号与滤波器的单位脉冲响应(Z域滤波器系数)进行卷积;寻找所需的滤波器单位脉冲响应的过程称为滤波器设计,信号卷积过程称为滤波。

求取系数的过程就是数字滤波器的设计过程。

  • Z变换应用于离散系统分析;频域转换到Z域
  • 拉普拉斯变换应用于连续系统分析; 时域滤波的关键是设计出Z变换形式的滤波器,就可以快速写出滤波公式。

  • 数字差分=简单的高通滤波器

  • 数据平滑=简单的低通滤波器
  • 高斯平滑滤波器



5.4 时域FIR滤波器

MATLAB函数:fir1函数()、fir2()函数、firls()函数
使用数组定义滤波器

滤波后的数据出现头部暂态现象,滤波器的系数(阶数越高)越长,头部暂态现象越长。

设计方法:

  • 脉冲响应函数法:
  • 窗函数设计法: 减小滤波器的纹波

5.5 IIR滤波器

设计方法:先借助模拟滤波器,再将模拟滤波器转为数字滤波器;

  • 脉冲响应不变法,可能会导致频率混叠
  • 模拟->数字双线性变换法,bilinear(b,a,Fs)1. 定义模拟滤波器频率同频带1. 映射到采样频率范围内1. 离散采样



5.6 其他滤波器

高通滤波器可以消除信号中的缓慢波动的低频现象。传感器晃动可能产生低频干扰。



低通滤波器可以消除信号中的细小快速(细小毛刺)波动的高频现象。



陷波滤波器(notch filter):属于带阻滤波器,其阻带很窄。可以快速衰减某一频率点信号,已达到阻止该频率信号的滤波效果。消除工频50Hz干扰

MATLAB中的使用:iirnotch(wo,bw)
[H,W]=freq(b,a):画出滤波器的频率特性

尖峰滤波器(peak filter):属于带通滤波器,其通带很窄。可以只让某一频率点信号通过,阻止其他频率成分的过滤效果。

MATLAB中的使用:iirpeak(wo,bw)

梳妆滤波器(comb filter):



半带滤波器:通带和阻带对称

希尔伯特全通滤波器90°相移滤波器:用于包络信号的检波处理。希尔伯特变换可以提取波形的包络

MATLAB中的使用:hilbert(x);



幅度调制与解调amplitude modulatiom

调幅、调频、调相



高频信号的传输距离更远。

MATLAB中的使用:

  • 调幅与解调:ammod(x,Fc,Fs);amdemod(y,Fc,Fs);
  • 调频与解调:fmmod(x,Fc,Fs,200);fmdemod(y,Fc,Fs,200);
  • 调相与解调(通过波形很难观察):pmmod(x,Fc,Fs,3.14);pmdemod(y,Fc,Fs,3.14);

5.7 数字音乐均衡器的设计

方法 1. FFT滤波:将信号进行FFT变换,根据均衡器设计不同频域的滤波器,两者相乘,再做IFFT,得到原信号1. 滤波器组:设计不同的带通滤波器即可,FIR、IIR都可以1. FIR2/FIRLS

六、信号的时频分析简介

6.1 信号小波分析




小波分析:可变长度窗口 dt*df=C

针对信号的不同频率特征,采用不同的窗口大小,如下图所示。低频信号使用宽窗口,高频信号使用窄窗口。



每个小波变换都有一个母小波,同时还有一个尺度函数,也被成为父小波。任何小波变换的小波基函数,其实就是对母小波和父小波缩放和平移后得到的集合。小波变换允许更加精确的局部描述以及信号特征的分离。一个傅立叶系数通常表示某个贯穿整个时间域的信号分量,因此,即使是临时的信号,其特征也被强扯到了整个时间周期去描述。而小波展开的系数则代表了对应分量它当下的自己,因此非常容易诠释。

窗口越宽,时间分辨率越低,频率分辨率越高;
窗口越窄,时间分辨率越高,频率分辨率越低; - 低频信号,时间精度调低,频率精度调高,用宽窗口- 高频信号,频率精度调低,时间精度调高,用窄窗口

6.2四种分析方法时频分析精度对比

  • 波形分析:没有频率分辨率
  • 频谱分析:没有时间分辨率
  • STFT分析(短时傅里叶分析):固定的时间–频率分辨率,等精度分析
  • 小波分析:可变时间–频率分辨率

6.3连续小波

连续小波变换的数学定义:信号与展缩小波函数的卷积分得到小波系数CWT。



通过展缩因子调节小波变换窗口的范围和频域轴的位置。 - 尺度因子s的概念:尺度意味着对一个信号的压缩或者展开。尺度s越大,信号展开的越大。- 平移因子τ的概念:平移将一个信号在时间轴上移动。 小波变换的计算



图示:信号与小波函数卷积



几种常见的小波函数 1. Haar小波



2. Mexican Hat小波



3. Morlet小波



MATLAB中的小波函数:wavelet

  • dwt() 单层小波分解
  • CWT()函数
[C,L]=wavedec(x,3,'db3');   %多层一维小波分解
wpt=wpdec(x,3,'db1','shannon'); %小波包分解:
[cd1,cd2,cd3]=detcoef(C,L,[1,2,3]); %获取细节系数
ca3=appcoef(C,L,'db3',3);   %获取近似系数 
A=wrcoef('a',C,L,'db3',3);  %对信号的低频部分重构,A表示信号的最小频率
D1=wrcoef('d',C,L,'db3',3); %对信号的高部分重构,D1表示信号的最大频率
D2=wrcoef('d',C,L,'db3',3); %对信号的高部分重构,D1表示信号的第二大频率
D3=wrcoef('d',C,L,'db3',3); %对信号的高部分重构,D1表示信号的第三大频率

小波变换结果与选择的小波函数有关

STFT谱阵和小波尺度谱阵的关系

谱阵STFT,等精度,代表信号代表在时–频面上的能量分别;Y轴为时间轴

小波尺度谱阵,变精度,尺度越小,精度越高;尺度越大,频率精度越低。代表在时–频面上的能量分别;X轴为时间轴
横纵坐标轴的方向不一样

尺度与频率的关系



6.5离散小波变换

时移因子直接按信号的采样间隔进行离散化,尺度因子则按照指数规律离散化

离散小波分解–>信号的多尺度滤波模型

离散二进小波变换:s=s^j-1,s={1,2,4,8…},
增加数据量

正交小波分解dwt:每一级正交小波分解后数据量不增加。



要求:滤波器对{hj,gj}是正交的。

MATLAB函数:wavedec()分解 waverec()合成
滤波:将该频域下小波系数置为零。

小波包变换:wpdec()分解 wprec()合成

6.6几种小波变换的关系



  • CWT信号:自相似性分析、信号跳变检测
  • DWT信号:信号多尺度分解
  • PWT信号:图像压缩

    【参考】:B站:【公开课】数字信号分析理论与实践(基于matlab) - 华中科技大学

    原文链接:https://blog.csdn.net/wxq_1993/article/details/106452297

免责声明:信息来源于互联网聚合整理或者作者原创,非原创文章标示"互联网",原创文章标示"作者昵称",本文仅代表文章作者的个人观点,且文中陈述文字和内容的时效性、真实性与本站无关,请读者仅作参考,如发现有害或侵权内容,请联系我们将在第一时间进行核实处理。


[推荐] 数字信号处理期末总复习

2021-12-14 10:44:43    互联网    219    分类:算法开发    专栏:数字信号处理   


[推荐] 信号调制基础理论知识

2021-12-14 22:00:52    互联网    277    分类:算法开发    专栏:数字信号处理   


[推荐] 数字信号处理的主要内容

2022-01-13 17:36:49    互联网    156    分类:算法开发    专栏:数字信号处理   



转发此文章到社区


进群

让志同道合读者学习交流



旋转机械故障诊断学基础知识

1. 机械故障诊断涉及哪些学科?做哪方面的科学研究?包含了哪些技术?形成了哪些方法?解决什么工程问题?(1) 涉及的学科有:机械、力学、电子、计算机、信号处理、人工智能等。(2) 机械故障诊断是研究机器或机组运行状态的变化在诊断信息中的反映,因此包括信号获取与传感技术、故障机理与征兆联系、信号处理与特征提取、识别分类与智能决策等方面的研究,根据基础和关键科学问题又可细分为机械系统运行状态下故障动态演化机理、机械系统动态信号处理的内积匹配原理与微弱信号特征增强机制、故障定量识别和剩余寿命预测原理、

2021-12-13 14:30:26    互联网    282    分类:算法开发    专栏:故障诊断   


基于Python的频谱泄露分析

1、频谱泄露  对于频率为fs的正弦序列,它的频谱应该只是在fs处有离散谱。但是,在利用DFT求它的频谱时,对时域做了截断,结果使信号的频谱不只是在fs处有离散谱,而是在以fs为中心的频带范围内都有谱线出现,它们可以理解为是从fs频率上“泄漏”出去的,这种现象称 为频谱“泄漏”。2、代码分析如果我们波形不能在fft_size个取样中形成整数个周期的话会怎样呢?将上篇博客中的采样对象...

2021-12-14 14:06:09    互联网    241    分类:算法开发    专栏:数字信号处理   


随机振动信号分析方法总结(时域、频域、时频域分析)

随机振动信号分析方法总结信号处理(信号滤波、时频域分析、神经网络、寿命预测)一、时域分析时域分析特征包括均值、方差、峭度、峰峰值等;振动信号降噪结果分析:对于去噪效果好坏的评价,常用信号的信噪比(SNR)、估计信号同原信号的均方根误差(RMSE)来判断。SNR 越高则说明混在信号里的噪声越小,否则相反。RMSE的计算值越小则表示去噪效果越好。信噪比定义:均方根误差定义:二、频域分析三、 时频联合域分析(Joint Time-Frequency Analysis,JTFA)即时频分析,

2021-12-14 18:49:46    互联网    376    分类:算法开发    专栏:振动信号预处理   


数据挖掘--风电机组异常数据识别与清洗

一、赛题背景(一)背景风能是一种环境友好且经济实用的可再生能源。中国是世界排名第一的风力发电国家、新装风力发电设备装机容量最大的国家,并且保持快速增长。由于风力发电正处于飞速发展阶段,风电场数量和规模不断扩大,然而受地理条件和环境因素限制,风电场多位于偏僻遥远的平原、山区或海上,因此为风电公司引入SCADA系统(数据采集与监视控制系统)对风电场群的日常运行进行集中监控、调度和管理,但风电机组受设备、环境、运行状态等因素影响,SCADA系统实时采集的风机运行数据会存在有大量异常值和缺失值,这些“脏数据”的

2021-12-19 01:00:15    互联网    254    分类:算法开发    专栏:风电机组   


提升预测精度:时间序列探索性数据分析-EDA

提升预测精度:探索性数据分析-EDAsales_train_validation.csv参考:Python版本EDA+传统时间序列方法:<li>整体聚合销量时间序列。 <ul>- 销售额普遍上升,这对沃尔玛来说是个好消息。我们可以找出一些每年的季节性,在圣诞节下降,这是一年中唯一的商店关闭的日子

2021-12-21 14:15:31    互联网    229    分类:算法开发    专栏:时间序列预测   


倒谱分析与倒谱求法

倒谱分析可检测频谱中的重复模式,使其对区分多个故障非常有用,该故障在不同的主要频谱(即FFT、阶次、包络和增强频谱)中很难看到。    最重要的行业应用与机械诊断相关,如齿轮箱分析,以及其他应用,如:1 回声检测和去除2 以及语音分析    在以下网站中:http://www.bksv.cn/Products/analysis-software/signal-a

2021-12-26 13:31:12    互联网    234    分类:算法开发    专栏:振动信号预处理   


语音信号的预处理

概述语音信号是一种非平稳的时变信号,它携带着大量信息。在语音编码、语音合成、语音识别和语音增强等语音处理中,都需要提取语音中包含的各种信息语音处理的目的对语音信号进行分析,提取特征参数,用于后续处理加工语音信息,如语音增强和语音合成中的应用根据所分析的参数类型,语音信号可以分成:时域分析最简单、最直观直接对语音信号的时域波形进行分析特征参数:语音的短时能力、平均幅度、短时...

2022-01-13 16:53:53    互联网    208    分类:算法开发    专栏:语音信号预处理   


语音信号的数据分析

数据集和代码均已上传到Github中,欢迎大家下载使用。Github地址:https://github.com/JasonZhang156/Sound-Recognition-Tutorial如果这个教程对您有所帮助,请不吝贡献您的小星星Q^Q.数据分析本节针对ESC-10数据集进行基本的数据分析,包括数据样本数,数据类别数,每类声音样本数等信息。并且对每类样本的声音波形,功率谱进...

2022-03-23 22:57:45    互联网    83    分类:算法开发    专栏:语音信号预处理   


计算机视觉-OpenCV视频处理

opencv视频处理主要设计了实时摄像头读取视频,本地视频读取并处理两个大的方面, 主要介绍了视频读取及处理的基本操作,内容概述如下: * 摄像头视频读入、本地视频读入 * 图像及视频写入保存 * 视频流单帧单帧处理 * 每一帧视频流的基本属性,如帧率、视频流宽度高度等

2022-03-31 21:39:22    博客笔记    73    分类:算法开发    专栏:图像预处理   


基于CNN实现细胞识别分类

本数据集用于检测疟疾的细胞图像,细胞图像分为两类: 一类是感染的细胞,另一类是未感染的细胞目标需求1:要基于卷积神经网络CNN来识别哪些细胞已经感染、哪些细胞还未感染 目标需求2:可视化模型随着迭代次数的训练集与测试集损失值的变化情况 目标需求3:可视化模型随着迭代次数的训练集与测试集准确率的变化情况

2022-03-31 21:45:25    博客笔记    59    分类:算法开发    专栏:未分组