Matlab振動程序-代碼作業(yè)

Matlab業(yè)-CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN一、課題任務(wù)要求 隨著機(jī)械工業(yè)不斷向自動化、高精度、智能化等方向的發(fā)展，在機(jī)械設(shè)備運(yùn)行及生產(chǎn)過程中進(jìn)行參量測試、分析與診斷等處理過程已成為必要環(huán)節(jié)，許多信號處理方法如時域統(tǒng)計分析、相關(guān)分析、相干分析、頻譜分析等已經(jīng)被廣泛被應(yīng)用與機(jī)械工程測試領(lǐng)域。本文為機(jī)械測試信號的時域和頻域分析，其中時域分析包括對信號最大FFTMATLAB域分析和頻域分析，并對繪出的頻譜進(jìn)行。、要FFT隨機(jī)數(shù)rand之和來擬合振動信號，用MATLAB中的函數(shù)mea、mi、max、st對序中的值、最大值、最小值、差等時域性能進(jìn)行分析，過用函數(shù)ffypsyrcey對振動信號進(jìn)行頻域的性能分析。在設(shè)計過程中的理DF、譜的和以及倒譜的和。DFT的定和DFTX(k)=DFT[x(n)]=

N1(nWkn

0≤k≤N-1n0 NX(n)=

1 N1X(kWkn) N Nk0DFT的DFT的在信號中的頻的分量的振動值的過DFT法FFT很出來。譜的和:X(Ω)=DFT[x(m)x(n)為N則X*Ω=DFT[N(-m)]1DFT[R(M)]=N1

DFT[x(m)]DFT[xN(-m)]1^Nx(Ω)=N

X(Ω)X*(Ω)=

|X(Ω)|^2綜上所述，先用FFTDFT，再計算幅頻特性的平方，再除以N，得到該隨機(jī)信號的功率譜估計。功率譜的意義在于它可以對信號中的周期成分進(jìn)行分析。③倒譜的概念與意義：倒譜可以分析復(fù)雜頻譜圖上的周期成分，分離和提取在密集泛頻信號中的成分。倒頻譜對于整個譜的形狀不敏感。三、實現(xiàn)程序A1=3;A2=4;f1=100;f2=250;fs=1000;t=0:1/fs:2;R=rand(1,N);Y=X1+X2+R;title('振動信號的波形');xlabel('時間/秒');ylabel('幅度');grid;holdon;%時域分析m=mean(Y);disp(m);mi=min(Y);disp(mi);mx=max(Y);disp(mx);st=std(Y);disp(st);%l=length(Y);figure(2);plot(f,abs(r));grid; holdfigure(3);psd(Y,2048,1000,kaiser(512,5),0,;figure(4);yc=rceps(Y);plot(yc);、運(yùn)行結(jié)果(1)時序列的平均值為序列的最小值為序列的最大值為序列的標(biāo)準(zhǔn)差為(2)。上FFT，從該頻譜中可以看100Hz，250Hz100Hz2