機(jī)械振動(dòng)信號(hào)分析及故障報(bào)警_課程設(shè)計(jì)(存儲(chǔ)版)
【正文】 師簽字 基層教學(xué)單位主任簽字 說(shuō)明： 此表一式四份，系、指導(dǎo)教師、學(xué)生各一份，報(bào)送院教務(wù)科一份。提取時(shí)域波形指標(biāo)如均值、峰峰值、峭度、偏度、脈沖因 數(shù)等。 信號(hào)分析的主要目的是尋找一種簡(jiǎn)單有效的信號(hào)變換方法，使信號(hào)所包含的重要信息能顯現(xiàn)出來(lái)。 傅里葉變換 有限長(zhǎng)序列可以通過(guò)離散傅里葉變換 (DFT)將其頻域也離散化成有限長(zhǎng)序列 .但其計(jì)算量太大 ,很難實(shí)時(shí)地處理問(wèn)題 ,因此引出了快速傅里葉變換 (FFT). 1965 年， Cooley 和 Tukey 提出了計(jì)算離散傅里葉變換（ DFT）的快速算法，將 DFT 的運(yùn)算量減少了幾個(gè)數(shù)量級(jí)。這樣變換以 后， 總的運(yùn)算次數(shù)就變成 N+2（ N/2） 2=N+N2/2。 第四章 MATLAB界面 設(shè)計(jì) 兩個(gè)畫(huà)圖按鈕簡(jiǎn)化成了一個(gè)下拉菜單，坐標(biāo)圖一和二分別顯示時(shí)域和頻域波形，計(jì)算按鈕進(jìn)行計(jì)算并把數(shù)值放到下面四個(gè)文本框中。樣本序號(hào) n/個(gè) 39。 %第一層細(xì)節(jié)信號(hào)的包絡(luò)譜 y=hilbert(d1)。,)。 fid=fopen(39。 N=1000。 xdata=(xdatamean(xdata))/std(xdata,1)。 % 峰峰值 set(,39。%偏度 。string39。gui_Name39。gui_LayoutF39。 else gui_mainf(gui_State, varargin{:})。 % 最大值 C=min(xdata)。,num2str(E))。 isequal(get(hObject,39。)。BackgroundColor39。), get(0,39。amp。white39。)) set(hObject,39。 =2400。 =39。 clear ss function pushbutton3_Callback(hObject, eventdata, handles) global F global J global FF if (FF=) J=0。), get(0,39。 load fid=fopen(39。 xdata=(xdatamean(xdata))/std(xdata,1)。 ylabel(39。,c,l,39。 d2=wrcoef(39。db1039。CurrentAxes39。)。)。樣本序號(hào) n/個(gè) 39。 d4=wrcoef(39。db1039。d39。 p=abs(fft(ydata,nfft))。 ylabel(39。頻率 f/Hz39。 y=ymean(y)。,2)。d39。db1039。 plot(xdata)。 xdata=fread(fid,N,39。 ylabel(39。 p=abs(fft(ydata,nfft))。d39。db1039。 d4=wrcoef(39。樣本序號(hào) n/個(gè) 39。)。)。 isequal(get(hObject,39。 fwrite(ss,J)。 =39。139。), get(0,39。BackgroundColor39。)。 isequal(get(hObject,39。defaultUicontrolBackgroundColor39。,39。 function edit1_Callback(hObject, eventdata, handles) function edit1_CreateF(hObject, eventdata, handles) if ispc amp。%峭度 set(,39。,num2str(A))。 ischar(varargin{1}) = str2func(varargin{1})。gui_OutputF39。 參考文獻(xiàn) 1）《微型計(jì)算機(jī)控制系統(tǒng)》 賴壽宏，機(jī)械工業(yè)出版社（教材） 2）《單片機(jī)及應(yīng)用》李大友，高等教育出版社（教材） 3）《信號(hào)處理原理及應(yīng)用》 謝平 等 機(jī)械工業(yè)出版社（教材） 4)《 Matlab 程序設(shè)計(jì)及其在信號(hào)處理中的應(yīng)用》 聶祥飛等 西南交通大學(xué)出版社 附錄一 matlab 總程序 function varargout = keshe0628(varargin) gui_Singleton = 1。 else J=1。,num2str(E))。 % 最大值 C=min(xdata)。)。39。功率譜 p/w39。 set(gcf,39。db1039。 plot(xdata)。 根據(jù)采樣定律，一個(gè)頻帶有限的信號(hào)，可以對(duì)它進(jìn)行時(shí)域采樣而不丟失任何信息， FFT 變換則說(shuō)明對(duì)于時(shí)間有限的信號(hào) (有限長(zhǎng)序列 )，也可以對(duì)其進(jìn)行頻域采樣，而不丟失任何信息。 設(shè) x(n)為 N項(xiàng)的復(fù)數(shù)序列，由 DFT 變換，任一 X（ m）的計(jì)算都需要 N次復(fù)數(shù)乘法和 N1次復(fù)數(shù)加法，而一次復(fù)數(shù)乘法 等于四次實(shí)數(shù)乘法和兩次實(shí)數(shù)加法，一次復(fù)數(shù)加法等于兩次實(shí)數(shù)加法，即使把一次復(fù)數(shù)乘法和一次復(fù)數(shù)加法定義成一次“運(yùn)算”（四次實(shí)數(shù)乘法和四次實(shí)數(shù)加法），那么求出 N 項(xiàng)復(fù)數(shù)序列的 X（ m） ,即 N點(diǎn) DFT變換大約就需要 N^2次運(yùn)算。它的主要特
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1