数据采集及处理.ppt

x(t)y(t)时延器乘法器y(t-τ)X(t)y(t-τ)积分器Rxy(τ)Fori=0ToNr(i)=0Forj=0ToNr(i)=r(i)+x(j)*y(j+i)NextNextFori=0ToNr(i)=r(i)/NNext计算公式：第95页,共102页，星期六，2024年，5月相关分析的工程应用案例：机械加工表面粗糙度自相关分析被测工件相关分析?提取出回转误差等周期性的故障源。第96页,共102页，星期六，2024年，5月案例：地下输油管道漏损位置的探测案例：AGV小车定位，声位笔定位第97页,共102页，星期六，2024年，5月实验：自相关分析第98页,共102页，星期六，2024年，5月3.信号频率成分直接估计算法f=？Fs=?N=1024dt=1.0/Fspi=3.1415926XR=0XI=0Forn=0ToN-1XR=XR+x(n)*cos(2*pi*f*n*dt)*dtXI=XI+x(n)*sin(2*pi*f*n*dt)*dtNextA=sqr(XR*XR+XI*XI)Q=atn(XI/XR)离散傅里叶变换的计算公式为：第99页,共102页，星期六，2024年，5月第100页,共102页，星期六，2024年，5月3、数字滤波利用离散时间系统特性对输入信号进行加工处理，把输入序列x(n)变换成一定的输出序列y(n)，从而达到改变信号频率构成的目的。x(n)h(n)y(n)y(n)=x(n)*h(n)设数字滤波器的脉冲响应序列为{h0,h1,h2,…,hm},则,展开：y(k)=h0x(k)+h1x(k+1)+h2x(k+2)+…+hmx(k+m)k=0,1,........第101页,共102页，星期六，2024年，5月将数字滤波器系数带入公式中就可以对信号进行数字滤波。Fori=0ToN－7Y(i)=-0.064X(i-0)+0.041X(i-1)+0.301X(i-2)+0.454X(i-3)+0.301X(i-4)+0.041X(i-5)-0.064X(i-6)Next第102页,共102页，星期六，2024年，5月4．离散时间与离散频率的傅里叶变换—离散傅里叶级数…………第63页,共102页，星期六，2024年，5月例：从连续周期信号的抽样得到第64页,共102页，星期六，2024年，5月第一个域离散函数?第二个域周期函数连续函数非周期函数?且易证：一个域中的周期函数的周期第65页,共102页，星期六，2024年，5月5.4.3从离散傅立叶级数(DFS)到离散傅立叶变换(DFT)效仿连续周期信号有傅立叶级数，记作：离散周期序列也有傅立叶级数，记作：周期性以N为周期第66页,共102页，星期六，2024年，5月离散周期序列的傅立叶级数(DFS)的正负运算对周期序列的基频是是K次谐波分量，谐波系数是谐波成分中只有N个是独立的，是周期的第67页,共102页，星期六，2024年，5月第68页,共102页，星期六，2024年，5月有限长序列是周期序列的一个周期有限长序列x(n)只有的N个值x(n)可看成是周期序列的主值序列，记作周期序列当叫做的主值周期，记作有限长序列的以N为周期的周期延拓第69页,共102页，星期六，2024年，5月的主值序列也是周期性的，相当于有限长序列周期延拓当时，其主值序列 相当于一个有限长序列第70页,共102页，星期六，2024年，5月和都取主值周期，得到离散傅立叶变换(DFT)对第71页,共102页，星期六，2024年，5月周期为N和周期为2N的不同当主值周期为0—N-1时，N点的DFT为当主值周期为0—2N-1时，2NDFT（接下页）第72页,共102页，星期六，2024年，5月第73页,共102页，星期六，2024年，5月小