機(jī)械故障大作業(yè)-資料下載頁(yè)
2025-03-26 01:41本頁(yè)面
【正文】 39。)。ylabel(39。振動(dòng)加速度A/v39。)。nfft=16384。S=psd(sig,nfft)。 %對(duì)sig做功率譜subplot(212)。 %繪制功率譜plot((0:nfft/21)/nfft*fs,S(1:nfft/2))。xlabel(39。頻率f/Hz39。)。ylabel(39。功率譜P/W39。)。[c,l]=wavedec(sig,4,39。db239。)。 %利用db2對(duì)sig進(jìn)行3級(jí)小波分解c3=wrcoef(39。a39。,c,l,39。db239。,3)。%重構(gòu)第13層細(xì)節(jié)d1d3和第3層蓋帽c3d3=wrcoef(39。d39。,c,l,39。db239。,3)。d2=wrcoef(39。d39。,c,l,39。db239。,2)。d1=wrcoef(39。d39。,c,l,39。db239。,1)。figure。subplot(414)。plot(t,c3)。ylabel(39。c339。)。%繪制c3subplot(413)。plot(t,d3)。ylabel(39。d339。)。%繪制d3subplot(412)。plot(t,d2)。ylabel(39。d239。)。%繪制d2subplot(411)。plot(t,d1)。ylabel(39。d139。)。%繪制d1y=hilbert(d1)。%對(duì)d1進(jìn)行Hilbert變換，得yydata=abs(y)。%ydata=|y|ydata=ydatamean(ydata)。對(duì)ydata去均值（目的是去除幅度較大的直流分量）P=psd(ydata,nfft)。%ydata的功率譜為Pfigure。plot((0:nfft/21)/nfft*fs,P(1:nfft/2))。%繪出d1的Hilbert包絡(luò)譜xlabel(39。頻率f/Hz39。)。P=P(1:nfft/2)。%故障頻率f1為包絡(luò)譜中幅度最大處的頻率[M,f1]=max(P)。f1=f1*fs/nfft1。 完美DOC格式整理
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1