大型實驗周MATLAB處理與分析

1、大型實驗周MATLAB辦理與解析大型實驗周MATLAB辦理與解析大型實驗周MATLAB辦理與解析實驗一柴油機(jī)狀態(tài)監(jiān)測與故障診斷綜合實驗本實驗以低速第一通道為例進(jìn)行解析。繪制時域波形及頻域波形的程序以下：clcFs=24000m=load（C:UsersAdministrator.LBDZ-20130910XVDesktopdaxingshiyanzh柴油機(jī)故障診斷實驗故障狀態(tài)低速1119189Waveform.TXT）fori=2:11figure(i):plot(m(1:1024,i);title(通道,int2str(i-1)A=mean(m(:,i);m(:,i)=m(:,i)-A;en

2、dfori=2:11y(:,i)=psd(m(:,i);figure(i+10);plot(y(:,i)title(通道,int2str(i-1)end圖像以下：故障頻域:故障頻域:正常時域：正常頻域:900求時頻域特色值程序：Cif;fs=24000;N=1024;lm=1024;%采樣頻率和數(shù)據(jù)點數(shù)n=0:N-1;t=n/fs;%時間序列z=fft(y);%對信號進(jìn)行快速Fourier變換mag=abs(z);%求得Fourier變換后的振幅f=n*fs/N;%頻序次列A=mean(x)%均值B=sum(x-A).A2)/lm%方差y=x-A;C=sqrt(sum(y.A2)/lm)%均方

3、根D=sum(y.A4)/lm%峭度E=sum(abs(y)/lm;%均對均值F=C/E%波形指標(biāo)G=max(y);%最大值H=min(y);%最小值匸G-H%峰值G=D/(CA4)%峭度指標(biāo)K=sum(f*mag)/sum(mag)%頻譜重心運行結(jié)果：正常時頻域特色值以下：A=-0.8026B=140.5535C=11.8555D=5.4590e+05F=2.0430I=200.4649G=27.6329K=1.1969e+04故障時頻域特色值以下:=-0.5562=32.4501=5.6965=2.5621e+04F=1.6648I=94.6457G=24.3311K=1.1954e+04

4、實驗二傳動系統(tǒng)狀態(tài)監(jiān)測與故障診斷綜合實驗此次實驗以電機(jī)頻率18HZ第五通道為例。繪制時域波形及頻域波形的程序：clcFs=24000m=load(C:UsersAdministrator.LBDZ-20130910XVDesktopdaxingshiyanzh柴油機(jī)故障診斷實驗故障狀態(tài)低速1119189_Waveform.TXT)fori=2:11figure(i):plot(m(1:1024,i);title(通道,int2str(i-1)A=mean(m(:,i);m(:,i)=m(:,i)-A;endfori=2:11y(:,i)=psd(m(:,i);figure(i+10);plot

5、(y(:,i)title(通道,int2str(i-1)end圖像以下：正常時域：正常頻域：故障頻域:故障頻域:求時頻域特色值程序以下:Cif;fs=24000;N=1024;lm=1024;%采樣頻率和數(shù)據(jù)點數(shù)n=0:N-1;t=n/fs;%時間序列z=fft(y);%對信號進(jìn)行快速Fourier變換mag=abs(z);%求得Fourier變換后的振幅f=n*fs/N;%頻序次列A=mean(x)%均值B=sum(x-A).A2)/lm%方差y=x-A;C=sqrt(sum(y.A2)/lm)%均方根D=sum(y.A4)/lm%峭度E=sum(abs(y)/lm;%均對均值F=C/E%波

6、形指標(biāo)G=max(y);%最大值H=min(y);%最小值I=G-H%峰值G=D/(CA4)%峭度指標(biāo)K=sum(f*mag)/sum(mag)%頻譜重心運行結(jié)果：正常：A=0.0198B=0.4638=0.6810=0.7675F=1.2770I=5.2093G=3.5681K=1.1986e+04故障：A=-0.0372B=1.6258C=1.2751D=9.2125F=1.2909I=8.7846G=3.4853K=1.1986e+04實驗三金屬箔式應(yīng)變片電橋性能及應(yīng)用實驗單臂電橋：MATLAB程序以下：U=3.36.59.913.116.419.723.026.229.532.8p=p

7、olyfit(W,U,1)%進(jìn)行曲線擬合U1=polyval(p,W);plot(W,U,*,W,U1,r)%作岀數(shù)據(jù)點和擬合曲線的圖形title(單臂電橋電壓(mv)與質(zhì)量(g)的關(guān)系)e=U1-U;e1=abs(e);e2=max(e1)f=e2/32.8運行結(jié)果：p=0.16390.0067e2=0.0642圖像以下：同理半橋電路：運行結(jié)果：p=0.40840.02670.0012圖像以下:全橋電路運行結(jié)果：p=f=0.2137則矯捷度為0.4135，非線性誤差為21.37%。圖像以下：實驗四差動變壓器性能實驗及振動測試實驗表2-1MATLAB程序：(這里為方便運算用y代替V，用x代替X

8、)1.121.381.601.852.032.24x=5.05.566.577.588.599.510.010.511.011.512.012.513.013.514.014.5p1=polyfit(x(1:10),y(1:10),1);%分段擬合p2=polyfit(x(11:end),y(11:end),1);Y1=vpa(poly2sym(p1,x),6)Y2=vpa(poly2sym(p2,x),6)x1=min(x(1:11):0.1:max(x(1:11);x2=min(x(11:end):0.1:max(x(11:end);y1=polyval(p1,x1);y2=polyval

9、(p2,x2);plot(x,y,*,x1,y1,x2,y2)運行結(jié)果：Y1=4.91824-0.463758*xY2=0.464727*x-4.46091圖像以下:吐X曲線下面計算矯捷度和非線性誤差：量程為1mm時：輸入MATLAB序：y=0.20.40.62x=10.010.511.0p=polyfit(x,y,1)y1=polyval(p,x);e=y1-y;e1=abs(e);e2=max(e1);f=e2/0.62運行結(jié)果：p=f=0.0108則矯捷度為0.4200，非線性誤差為同理1.08%。量程為-1mm時程序：y=0.720.460.2x=99.510.0p=polyfit(x

10、,y,1)y1=polyval(p,x);e=y1-y;e1=abs(e);e2=max(e1);f=e2/0.72運行結(jié)果：p=-0.52005.4000f=1.4649e-15量程為3mm時程序：y=0.20.40.620.881.121.381.60 x=10.010.511.011.512.012.513.0p=polyfit(x,y,1)y1=polyval(p,x);e=y1-y;e1=abs(e);e2=max(e1);f=e2/1.38運行結(jié)果：p=f=0.0202量程為-3mm時程序：y=1.711.481.230.980.720.460.2x=77.588.599.510.

11、0p=polyfit(x,y,1)y1=polyval(p,x);e=y1-y;e1=abs(e);e2=max(e1);f=e2/1.71運行結(jié)果：p=-0.50575.2671f=0.0100實驗五電渦流傳感器的位移特色及應(yīng)用實驗表3-1MATLAB程序：X=2.53.54.55.56.57.58.59.510.511.5U=0.060.620.981.181.291.351.381.401.421.43title(U-X曲線)plot(X,U,*)點擊工具基本擬合保形插值，得U-X曲線以下:計算量程為1mm、3mm及5mm時矯捷度和線性度同實驗四對線性度矯捷度的求法。量程為1mm時運行結(jié)

12、果:量程為3mm時運行結(jié)果:P=P=f=f=1.1639e-150.0814量程為5mm時運行結(jié)果：P=f=0.1739表3-2量程3mm時運行結(jié)果：量程1mm時運行結(jié)果:P=P=f=f=0.09121.5209e-16表3-3量程1mm時運行結(jié)果：P=f=1.748e-16量程3mm時運行結(jié)果：P=f=0.1070表3-6用MATLAB作出振動頻率和輸出峰值曲線以下表3-7用MATLAB作出V-n曲線以下:實驗六溫度控制與測溫傳感器溫度特性測試實驗表4-1MATLAB程序：T=50556065707580859095100105110;40.545.049.554.1;p=polyfit(T

13、,U,1)U1=polyval(p,T);plot(T,U,*,T,U1,R)title(Uo2-T曲線)e=U1-U;e1=abs(e);e2=max(e1);f=e2/54.1運行結(jié)果：f=0.0177圖像如右圖4-1所示。表4-2MATLAB程序：T=5055606570758085120;28.649.451.4;p=polyfit(T,U,1)U1=polyval(p,T);plot(T,U,*,T,U1,R)title(Uo2-T曲線)e=U1-U;e1=abs(e);e2=max(e1);f=e2/51.4運行結(jié)果：f=0.0089圖像如右圖4-2所示。Ua2-T曲線圖6-2表4

14、-2曲線實驗原理同上：表5-1運行結(jié)果：k=1.0e+03*2.50700.0160f=0.0330圖像如圖7-1所示。表5-2運行結(jié)果:k=1.0e+03*-2.81730.24830.2227圖像如圖7-2實驗七擴(kuò)散硅壓阻式壓力傳感器的壓力測試實驗P1（増大卜U成塔大）應(yīng)表5-3運行結(jié)果：k=1.0e+03*1.8937-0.2230f=0.0585圖像如右圖：表5-4運行結(jié)果：圖7-3表5-3曲線k=1.0e+03*-1.9780-0.0333f=0.0352圖像如右圖：P1倔?。└pP2（鬲小圖7-4表5-4曲線實驗八霍爾傳感器特色及應(yīng)用實驗表6-3MATLAB程序：V=46810121

15、62024;U1=polyval(k,V);plot(V,RPM,*,V,U1,r)title(V-RPM)運行顯示圖像：表6-1量程為2mm時運行結(jié)果：p=f=0.0860量程為4mm時運行結(jié)果：p=f=0.0733量程為6mm時運行結(jié)果：p=0.5269圖像以下:實驗九扭矩傳感器靜態(tài)性能實驗逐漸加砝碼時MATLAB程序：300U=21.042.564.685.5103.6123.0p=polyfit(W,U,1)y1=polyval(p,W);plot(W,U,*,W,y1,r)title(U-W曲線)e=y1-U;e1=abs(e);e2=max(e1);f=e2/123運行結(jié)果：0.40811.9467f=0.0157圖像以下:曲絨從100mm移到內(nèi)端MATLAB程序：L=10090807060U=123.0111.296.380.465.6p=polyfit(L,U,1)y1=polyval(p,L);plot(L,U,*,L,y1,r)title(U-L曲線)e=y1-U;e1=abs(e);e2=max(e1);f=