航空器檢測(cè)與診斷技術(shù)試驗(yàn)指導(dǎo)書

1、鄭州航空工業(yè)管理學(xué)院航空器檢測(cè)與診斷技術(shù)課程實(shí)驗(yàn)指導(dǎo)書航空工程系2015-11 TOC o 1-5 h z HYPERLINK l bookmark0 o Current Document 安全注意事項(xiàng) 1 HYPERLINK l bookmark2 o Current Document 實(shí)驗(yàn)基本要求2 HYPERLINK l bookmark4 o Current Document 實(shí)驗(yàn)一智能故障診斷系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn) 3飛行力學(xué)課程實(shí)驗(yàn)指導(dǎo)書安全注意事項(xiàng)為保證實(shí)驗(yàn)安全進(jìn)行，請(qǐng)同學(xué)們一方面要嚴(yán)格遵守實(shí)驗(yàn)室制度，另一方面在 使用實(shí)驗(yàn)設(shè)備時(shí)注意如下事項(xiàng)：一、檢查各儀器電源線是否插緊插好，各儀器是否

2、可靠接地，以防觸電。二、所有儀器設(shè)備工作過程中發(fā)現(xiàn)異常應(yīng)立即斷電，并上報(bào)老師處理。飛行力學(xué)課程實(shí)驗(yàn)指導(dǎo)書實(shí)驗(yàn)基本要求1、安全第一，嚴(yán)格遵守實(shí)驗(yàn)室要求2、請(qǐng)?zhí)崆邦A(yù)習(xí)實(shí)驗(yàn)內(nèi)容3、按照實(shí)驗(yàn)指導(dǎo)書要求做實(shí)驗(yàn)4、注意觀察實(shí)驗(yàn)現(xiàn)象，并認(rèn)真記錄實(shí)驗(yàn)數(shù)據(jù)5、對(duì)實(shí)驗(yàn)現(xiàn)象和數(shù)據(jù)進(jìn)行科學(xué)、詳細(xì)分析6、將實(shí)驗(yàn)裝置收拾好7、撰寫實(shí)驗(yàn)報(bào)告飛行力學(xué)課程實(shí)驗(yàn)指導(dǎo)書實(shí)驗(yàn)一智能故障診斷系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)一、實(shí)驗(yàn)?zāi)康?、了解常用智能故障診斷系統(tǒng)的基本原理。2、學(xué)會(huì)使用matlab進(jìn)行智能故障診斷系統(tǒng)的設(shè)計(jì)實(shí)現(xiàn) 二、實(shí)驗(yàn)內(nèi)容.基于matlab的灰度關(guān)聯(lián)分析法利用版色美或分析進(jìn)行除合if價(jià)的步景是；。根據(jù)i i價(jià)目的確定評(píng)價(jià)淅歷體系.收

3、集評(píng)價(jià)教據(jù) 設(shè)數(shù)據(jù)序的形成如卜辯牌：修仁，與f力其中 m為指標(biāo)的個(gè)數(shù).掙為叁評(píng)人演出中X = fx,ah若* x,（H） /,?）確定賽號(hào)數(shù)據(jù)，考數(shù)據(jù)到應(yīng)該星一個(gè)理個(gè)的比較標(biāo)準(zhǔn).可以以存指標(biāo)的優(yōu)1B破蛆劣值）構(gòu)成卷考數(shù)據(jù)刷.也 可根據(jù)評(píng)價(jià)目的選擇其它參照值紀(jì)作茶gn，,看便匕冷瓜”分送個(gè),w自個(gè)被泮階對(duì)索布標(biāo)序列a七較序列）與多學(xué)序列對(duì)象足案的能時(shí)有值.即I冷- .* = t. bbb. , / = 1. ； n口句定 mtim si ica t kt -1*1 rl, fWmiiRRM霹 * t k . III #15）計(jì)算兌取系數(shù)III卜式.分別計(jì)算毋個(gè)比較序列與參考序列財(cái)應(yīng)疣素的關(guān)聯(lián)系毆

4、P maxniax 即4J x： fki i ri式中為分辯第數(shù).在8J）內(nèi)取值.若P越小.關(guān)注系敷間差異越大.區(qū)分能力趙遢通常o瓶。&如 肇fXom/為呆優(yōu)值數(shù)席列越大及好方 & g 例最ftfflttffi列 越大型不好6）計(jì)算關(guān)會(huì)序q若評(píng)價(jià)時(shí)象世:蛟序。1）分別H-其一指標(biāo)與器序列時(shí)應(yīng)尤第的“兼系粒的均值.以反映樸群 價(jià)均象與參考序列的關(guān)紙關(guān)系.井稱其為美麻序.妃為：芯=-L y（it*e7）儂據(jù)各察附室的美聯(lián)序.1匕統(tǒng)合講一結(jié)果仿真程序clc;%清屏clear; % 清參數(shù)clear allx(1,:)=1 1.0684 1.2388 1.4878 1.6827;x(2,:)=1 1.

5、0951 1.2307 1.4220 1.6135;飛行力學(xué)課程實(shí)驗(yàn)指導(dǎo)書x(3,:)=1 1.0164 0.9743 0.9635 0.9992;%列出各數(shù)值，可修改x(4,:)=1 1.1063 1.3318 1.6009 1.7871; x(5,:)=1 0.9470 0.8996 0.8107 0.7402; x(6,:)=1 0.9643 0.8977 0.7561 0.7352; x(7,:)=1 0.9914 0.9600 0.9514 0.9314; x(8,:)=1 0.9938 0.9631 0.9477 0.9323; x(9,:)=1 1.0539 1.2146 1.3

6、184 1.3455; x(10,:)=1 0.9986 1.0068 1.0069 0.9986; x(11,:)=1 0.9663 0.9401 0.8677 0.8621; x(12,:)=1 0.8304 0.8232 0.8129 0.7995; x(13,:)=1 1.0217 1.1382 1.2040 1.2433; x(14,:)=1 1.0075 0.9747 0.9646 0.9503; x(15,:)=1 0.7468 1.1340 0.8214 0.9461; x(16,:)=1 0.9741 0.9121 0.9617 0.9032; m=16; n=5; x0=1

7、 1.0732 1.2123 1.2295 1.3364; for j=1:m for i=1:n A(j,i)=abs(x(j,i)-x0(i);%求出序列差矩陣 Aend end a=max(max(A); %求 A 的最大值 ab=min(min(A);%求 A 的最小值 bfor j=1:m for i=1:n f(j,i)=(b+0.5*a)/(A(j,i)+0.5*a); %求關(guān)聯(lián)系數(shù)矩陣 fend end for j=1:m c(j)= mean(f(j,:); %求關(guān)聯(lián)度end figure(1); plot(c); figure(2) bar(c); cc,ss=sort(b

8、);.灰色預(yù)測(cè)法飛行力學(xué)課程實(shí)驗(yàn)指導(dǎo)書function GM1_1(X0)%format long ;m,n=size(X0);X1=cumsum(X0); % 累加X2=;for i=1:n-1X2(i,:)=X1(i)+X1(i+1);endB=-0.5.*X2 ;t=ones(n-1,1);B=B,t ;% 求 B矩陣YN=X0(2:end);P_t=YN./X1(1:(length(X0)-1) %對(duì)原始數(shù)據(jù)序列X0進(jìn)行準(zhǔn)光滑性檢驗(yàn),%序列 X0 的光滑比 P(t)=X0(t)/X1(t-1)A=inv(B.*B)*B.*YN.;a=A(1)u=A(2)c=u/a ;b=X0(1)-c

9、 ;X=num2str(b),exp,(,num2str(-a),k,),num2str(c);strcat(X(k+1)=,X)%syms k;for t=1:length(X0)k(1,t)=t-1;endk飛行力學(xué)課程實(shí)驗(yàn)指導(dǎo)書Y_k_1=b*exp(-a*k)+c;for j=1:length(k)-1Y(1,j)=Y_k_1(j+1)-Y_k_1(j);endXY=Y_k_1(1),Y % 預(yù)測(cè)值CA=abs(XY-X0) ; % 殘差數(shù)列Theta=CA % 殘差檢驗(yàn)絕對(duì)誤差序列XD_Theta= CA ./ X0 % 殘差檢驗(yàn) 相對(duì)誤差序列AV=mean(CA); %殘差數(shù)列平均

10、值R_k=(min(Theta)+0.5*max(Theta)./(Theta+0.5*max(Theta) ;% P=0.5 R=sum(R_k)/length(R_k) % 關(guān)聯(lián)度Temp0=(CA-AV)A2 ;Temp1=sum(Temp0)/length(CA);S2=sqrt(Temp1) ; % 絕對(duì)誤差序列的標(biāo)準(zhǔn)差 %AV_0=mean(X0); % 原始序列平均值Temp_0=(X0-AV_0).A2 ;Temp_1=sum(Temp_0)/length(CA);S1=sqrt(Temp_1) ;% 原始序列的標(biāo)準(zhǔn)差TempC=S2/S1*100;% 方差比C=strcat(

11、num2str(TempC),%) % 后驗(yàn)差檢驗(yàn) 方差比 %SS=0.675*S1 ;Delta=abs(CA-AV);飛行力學(xué)課程實(shí)驗(yàn)指導(dǎo)書TempN=find(Delta=SS);N1=length(TempN);N2=length(CA);TempP=N1/N2*100;P=strcat(num2str(TempP),%) %后驗(yàn)差檢驗(yàn)計(jì)算小誤差概率調(diào)用例子：X0=2,874,3.278,3.337,3.39,3.679;GM1_1(X0).神經(jīng)網(wǎng)絡(luò)法P=2.50 0.11 10.40 12.60 45.5;0.24 0.36 10.30 0.30 48.2;1.81 0.27 11.

12、72 11.50 46.3;1.77 0.22 10.33 14.43 50.0;1.75 0.23 10.40 9.80 49.9;1.82 0.25 10.33 12.50 72.2;1.84 0.26 10.10 11.20 18.5;P=P;T=0.9 0.1 0.1 0.1 0.1 0.1 0.1;0.1 0.9 0.1 0.1 0.1 0.1 0.1;0.1 0.1 0.9 0.1 0.1 0.1 0.1;0.1 0.1 0.1 0.9 0.1 0.1 0.1;0.1 0.1 0.1 0.1 0.9 0.1 0.1;0.1 0.1 0.1 0.1 0.1 0.9 0.1;這里的矩陣

13、0.1 0.1 0.1 0.1 0.1 0.1 0.9;net=newff(minmax(P),11 7,tansig,logsig,trainlm);%飛行力學(xué)課程實(shí)驗(yàn)指導(dǎo)書大小范圍可用 minmaxnet.trainParam.epochs=15;net.trainParam.goal=0.01;LP.lr=0.1;net=train(net,P,T);4.knnfunction target=KNN(in,out,test,k)% in:training samples data,n*dmatrix% out:training samples class label,n*1% test:

14、testing data% target: class label given by knn% k: the number of neighborsClassLabel=unique(out);c=length(ClassLabel);n=size(in,1);% target=zeros(size(test,1),1);dist=zeros(size(in,1),1);for j=1:size(test,1)cnt=zeros(c,1);for i=1:ndist(i)=norm(in(i,:)-test(j,:);endd,index=sort(dist);for i=1:kind=fin

15、d(ClassLabel=out(index(i);cnt(ind)=cnt(ind)+1;飛行力學(xué)課程實(shí)驗(yàn)指導(dǎo)書endm,ind=max(cnt);target(j)=ClassLabel(ind);End4.距離函數(shù)法算法步驟：step.1-初始化距離為最大值step.2-計(jì)算未知樣本和每個(gè)訓(xùn)練樣本的距離diststep.3-得到目前K個(gè)最臨近樣本中的最大距離maxdiststep.4-如果dist小于maxdist ,則將該訓(xùn)練樣本作為K-最近鄰樣本step.5-重復(fù)步驟2、3、4,直到未知樣本和所有訓(xùn)練樣本的距離都算完step.6-統(tǒng)1t K-最近鄰樣本中每個(gè)類標(biāo)號(hào)出現(xiàn)的次數(shù)step.

16、7-選擇出現(xiàn)頻率最大的類標(biāo)號(hào)作為未知樣本的類標(biāo)號(hào)KNN的matlab簡(jiǎn)單實(shí)現(xiàn)代碼function target=KNN(in,out,test,k)% in:training samples data,n*d matrix% out: training samples class label,n*1% test:testing data% target:class label given by knn% k:the number of neighborsClassLabel=unique(out);飛行力學(xué)課程實(shí)驗(yàn)指導(dǎo)書c=length(ClassLabel);n=size(in,1);% target=zeros(size(test,1),1);dist=zeros(size(in,1),1);for j=1:size(test,1)cnt=zeros(c,1);for i=1:ndist(i)=norm(in(i,:)-test(j,:);endd,index=sort(dist);for i=1:kind=find(ClassLabel=out(index(i);cnt(ind)=cnt(ind