資源描述：
《信號時域與頻域分析》由會員上傳分享，免費在線閱讀，更多相關內容在工程資料-天天文庫。
1、信號時域與頻域分析實驗報告姓名:楊班級：機械學號：213實驗數據中，電機轉速為1200r/min,釆樣頻率為1280Hz。Hz3為X位移振幅數據，Hz4為Y位移振幅數據，Hz5為速度振幅數據。Matlab中信號特征對應函數編程ma=max(Hz)mi=min(Hz)me=mean(Hz)pk=ma-miva=var(Hz);st=std(Hz);ku=kurtosis(Hz);rm=rms(Hz);%最人值%最小值%平均值%峰■峰值%方差%標準差%峭度%均方根、X軸位移測量分析plot(Fs3,Hz3)時域圖:ma=52.0261mi=56.7010me=1.8200pk=108.7271v
2、a=1.3870e+03st=37.2431ku=1.5462rm=37.2693頻域圖：fs=1280;x=Hz3;N=length(Hz3);df=fs/N;f=O:df:N*df-df;y二fft(x);y=abs(y)*2/N;figure(l);Plot(ty);xlabelC頻率/Hz1)ylabef幅健)頻譜幅值取得最大值51.9847um,頻率為20Hz,與電機轉速對應頻率一致,應為電機軸未動平衡所致；二倍頻處有較大振幅，可能為軸承間隙過大所致。二、Y軸位移測量分析plot(Fs4,Hz4)時域圖:ma=61.3987mi=-74.6488me=-1.1948pk=136.0
3、475av=42.6109va=2.2428e+03st=47.3582ku=1.5135rm=47.3501頻域圖:fs=1280;x=Hz4;N=length(Hz4);df=fs/N;f=0:df:N*df-df;y二fft(x);y=abs(y)*2/N;figure(l);Plot(ty);xlabelC頻率/Hz')ylabef幅值J頻譜幅值取得最大值66.6319um,頻率為20Hz,與電機轉速對應頻率一致,應為電機軸未動平衡所致；二倍頻處有較大振幅，可能為軸承間隙過大所致。三、速度測量分析plot(Fs5,Hz5)時域圖：ma=1.0907mi=0.1652me=0.6323
4、pk=0.9255av=0.6323va=0.0397st=0.1993ku=1.9936rm=0.6630頻域圖：fs=1280;x=Hz5;N=length(Hz5);df=fs/N;f=0:df:N*df-df;y=fft(x)；y=abs(y)*2/N;figure(l);plot(ty);xlabelC頻率/Hz1)ylabel(?幅值')速度頻譜幅值取得最大值1.2647mm/s,頻率為20Hz,與電機轉速對應頻率一致。振幅主要集中在20Hz附近，高頻較為穩(wěn)定。速度自相關函數圖:fs=1280;x=Hz;N=length(Hz);df=fs/N;f=O:df:N*df-df;[a
5、,b]=xcorr(Hz,'unbiased');plot(b*df,a);titled速度自相關函數匕xlabel('t');、軸心軌跡測量分析figure(l);plot(Hz3,Hz4)gridon;五、思考題分析1.(1)A,B的波形圖:fs=1024;t=O:l/fs:l;x=10*sin(60*pi*t+pi/4)+5*sin(120*pi*t+pi/3);y=12*sin(60*pi*t+pi/2)+4*sin(120*pi*t+pi/6);subplot(2,l,l);plot(t’x);title('A的波形圖')xlabef時間/sjylabef幅直)subplot(2,
6、l,2);Plot(t’y);title('B的波形圖')xlabef時間/sjylabef幅直)figure(2)；plot(x,y);title(1軸心軌跡1)；(2)B的幅值譜，相位譜，功率譜：fs=1024;df=l/fs;t=O:df:l;N=length(t);y=12*sin(60*pi*t+pi/2)+4*sin(120*pi*t+pi/6);f=fft(y,N);mag=abs(f)*2/N;phase=angle(f);power=mag.A2;F=O:df*N:(N-l)/N/df;figure(l);plot(F(l:(N+l)/2)zmag(l:(N+l)/2));
7、titleCB的幅值譜');figure(2);plot(F(l:(N+l)/2),phase(l:(N+l)/2));title('B的相位譜')；figure(3);plot(F(l:(N+l)/2),power(l:(N+l)/2));B的幅值譜:B的titleCB的功率譜')；譜:B的功率譜:(3)A和B的互相關函數，A的自相關函數:fs=1024;df=l/fs;t=O:df:l;N=length(
