大型实验周MATLAB处理与分析.doc

实验一 柴油机状态监测与故障诊断综合实验 本实验以低速第一通道为例进行分析。 绘制时域波形及频域波形的程序如下： clc Fs=24000 m=load(C:\Users\Administrator.LBDZV\Desktop\daxingshiyanzhou\柴油机故障诊断实验\故障状态\低速\1119189_Waveform.TXT) for i=2:11 figure(i):plot(m(1:1024,i)); title([通道,int2str(i-1)]) A=mean(m(:,i)); m(:,i)=m(:,i)-A; end for i=2:11 y(:,i)=psd(m(:,i)); figure(i+10);plot(y(:,i)) title([通道,int2str(i-1)]) end 正常时域： 故障频域： 求时频域特征值程序： clf; fs=24000; N=1024; lm=1024;%采样频率和数据点数 n=0:N-1; t=n/fs;%时间序列 z=fft(y);%对信号进行快速Fourier变换 mag=abs(z);%求得Fourier变换后的振幅 f=n*fs/N;%频率序列 A=mean(x)%均值 B=sum((x-A).^2)/lm%方差 y=x-A; C=sqrt(sum(y.^2)/lm)%均方根 D=sum(y.^4)/lm%峭度 E=sum(abs(y))/lm;%均对均值 F=C/E%波形指标 G=max(y);%最大值 H=min(y);%最小值 I=G-H%峰值 G=D/(C^4)%峭度指标 K=sum(f*mag)/sum(mag)%频谱重心 A = -0.8026 B = 140.5535 C = 11.8555 D = 5.4590e+05 F = 2.0430 I = 200.4649 G = 27.6329 K = 1.1969e+04 实验二 传动系统状态监测与故障诊断综合实验 此次实验以电机频率18HZ第五通道为例。 绘制时域波形及频域波形的程序： clc Fs=24000 m=load(C:\Users\Administrator.LBDZV\Desktop\daxingshiyanzhou\柴油机故障诊断实验\故障状态\低速\1119189_Waveform.TXT) for i=2:11 figure(i):plot(m(1:1024,i)); title([通道,int2str(i-1)]) A=mean(m(:,i)); m(:,i)=m(:,i)-A; end for i=2:11 y(:,i)=psd(m(:,i)); figure(i+10);plot(y(:,i)) title([通道,int2str(i-1)]) end 正常时域： 故障频域： 求时频域特征值程序如下： clf; fs=24000; N=1024; lm=1024;%采样频率和数据点数 n=0:N-1; t=n/fs;%时间序列 z=fft(y);%对信号进行快速Fourier变换 mag=abs(z);%求得Fourier变换后的振幅 f=n*fs/N;%频率序列 A=mean(x)%均值 B=sum((x-A).^2)/lm%方差 y=x-A; C=sqrt(sum(y.^2)/lm)%均方根 D=sum(y.^4)/lm%峭度 E=sum(abs(y))/lm;%均对均值 F=C/E%波形指标 G=max(y);%最大值 H=min(y);%最小值 I=G-H%峰值 G=D/(C^4)%峭度指标 K=sum(f*mag)/sum(mag)%频谱重心 A = 0.0198 B = 0.4638 C = 0.6810 D = 0.7675 F = 1.2770 I = 5.2093 G = 3.5681 K = 1.1986e+04 实验三 金属箔式应变片电桥性能及应用实验 单臂电桥： MATLAB程序如下： W=[20 40 60 80 100 120 140 160 180 200] U=[3.3 6.5 9.9 13.1 16.4 19.7 23.0 26.2 29.5 32.8] p=polyfit(W,U,1)%%进行曲线拟合 U1=polyval(p,W); plot(W,U,*,W,U1,r)%作出数据点和拟合曲线的