xsj 基于小波变换的碰磨故障信号的特征提取

基于小波变换的碰磨故障信号的特征提取，可以画出信号原图，轴心轨迹，频谱图以及多层小波变换的重构信号

clear;
close all;%关闭所有
clc;%清除命令窗口 
figure（1）;
set（gcf‘Color‘‘W‘）;
fidu=fopen（‘c:\MATLAB7\work\xsj\Omega8=5\u1.dat‘）;
fidv=fopen（‘c:\MATLAB7\work\xsj\Omega8=5\u3.dat‘）;
fidt=fopen（‘c:\MATLAB7\work\xsj\Omega8=5\tol.dat‘）;

U=fscanf（fidu‘%g‘）; 
V=fscanf（fidv‘%g‘）;
C=fscanf（fidt‘%g‘）;
status=fclose（fidu）;
status=fclose（fidv）;
status=fclose（fidt）;

A=U*1000;
B=V*1000; 
E=C/2/5/pi/18;

subplot（221）;% x 方向波形
plot（EA‘b‘）;
title（‘X轴信号图‘）;
xlabel（‘时间t/s‘）;
ylabel（‘位移x/mm‘）; 
axis（[20 20.4 -Inf Inf]）;
%set（gca‘ytickLabel‘[-2 -1 0 1 2 ]）;
%set（gca‘xtickLabel‘[0 0.001 0.002 0.003 0.004]）;

subplot（222）;
plot（A（1:50000）B（1:50000）‘b‘）;
title（‘轴心轨迹图‘）;
xlabel（‘位移x/mm‘）; 
ylabel（‘位移y/mm‘）; 
%axis（[-4 4 -4 4]）;

Y=fft（A131072）;
N=length（Y）;
M1=abs（Y（1:N/2

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2623  2010-06-09 10:00  xsj\five.m

     文件       2637  2010-06-13 16:14  xsj\fivepfive.m

     文件        692  2010-06-03 13:50  xsj\MaxMode.m

     文件    9026193  2004-07-22 14:44  xsj\Omega8=5.5\tol.dat

     文件   13622356  2004-07-22 14:43  xsj\Omega8=5.5\u1.dat

     文件   13197988  2004-07-22 14:43  xsj\Omega8=5.5\u3.dat

     文件        190  2004-07-22 13:50  xsj\Omega8=5.5\新建 文本文档.txt

     文件    9026193  2004-07-22 14:46  xsj\Omega8=7\tol.dat

     文件   13639093  2004-07-22 14:46  xsj\Omega8=7\u1.dat

     文件   13197378  2004-07-22 14:46  xsj\Omega8=7\u3.dat

     文件        189  2004-07-22 13:53  xsj\Omega8=7\新建 文本文档.txt

     文件       2644  2010-06-13 16:03  xsj\onepeight.m

     文件       2635  2010-06-09 09:38  xsj\onepfive.m

     文件       2518  2010-06-10 19:02  xsj\oneptwo.asv

     文件       2947  2010-06-13 15:58  xsj\oneptwo.m

     文件        418  2010-06-10 17:58  xsj\r1.asv

     文件        192  2010-06-10 17:58  xsj\r1.m

     文件       2637  2010-06-13 16:11  xsj\three.m

     文件       2637  2010-06-13 16:12  xsj\threepfive.m

     文件       2635  2010-06-13 16:05  xsj\two.m

     文件       2635  2010-06-13 16:09  xsj\twopeight.m

     文件       2596  2010-06-09 09:16  xsj\twopfive.m

     文件       2635  2010-06-13 16:07  xsj\twopone.m

     文件       2635  2010-06-09 09:53  xsj\twopseven.m

     文件       2619  2010-06-09 09:51  xsj\twopthree.m

     文件       2707  2010-06-09 10:08  xsj\twoptwo.m

     文件        187  2010-06-07 22:47  xsj\新建 文本文档.txt

     目录          0  2010-06-01 08:41  xsj\Omega8=5.5

     目录          0  2010-06-01 08:41  xsj\Omega8=7

     目录          0  2010-06-10 19:02  xsj

............此处省略3个文件信息