matlab處理實驗數(shù)據(jù) 實驗：利用轉子試驗臺測試轉子不平衡以及渦動狀態(tài) 故障診斷

超級小寶

我有實驗數(shù)據(jù)是dat格式的 要用matlab分析出頻譜圖跟時域圖要得出的圖像跟故障應該有的圖像相似并且?guī)臀医忉屢幌?轉子不平衡的時域圖 頻譜圖各是什么特征我編的程序如下 你看看用得著不%***************1.bch0****************%clcclfclearbch0=load('F:\MATLAB7\lijiande\bch0.dat');fs=100;%設定采樣頻率N=128;n=0:N-1;t=n/fs;bch0=bch0(n+1);%采樣figure(1);subplot(211);plot(t,bch0);%作正弦信號的時域波形axis([0,N/fs,-0.4,0.4]);xlabel('t');ylabel('y');title('bch0采樣t=[0,1.28]');grid;%進行FFT變換并做頻譜圖y=fft(bch0,N);%進行fft變換mag=abs(y);%求幅值f=(0:length(y)-1)'*fs/length(y);%進行對應的頻率轉換figure(1);subplot(212);plot(f,mag);%做頻譜圖axis([0,100,0,4]);xlabel('頻率(Hz)');ylabel('幅值');title('bch0幅頻譜圖N=128');grid;
3 E8 m- @6 I, u6 ^- D9 `& _$ b謝謝大家了~！！！

超級小寶

我有實驗數(shù)據(jù)是dat格式的 要用matlab分析出頻譜圖跟時域圖
4 s8 g0 F8 C; W/ Q6 k9 Y# X! V要得出的圖像跟故障應該有的圖像相似
并且?guī)臀医忉屢幌?轉子不平衡的時域圖 頻譜圖各是什么特征
3 Y" A6 `: n' @9 a6 w! @  g我編的程序如下 你看看用得著不

' P6 r) i) _. [! ~0 f! V%***************1.bch0****************%
clc
clf
clear
bch0=load('F:\MATLAB7\lijiande\bch0.dat');
9 ~: f) n6 w% t. i& P& f* Afs=100;%設定采樣頻率
* {  h" s$ X, e  G3 yN=128;
n=0:N-1;
t=n/fs;
3 t; p. ~% I; }3 jbch0=bch0(n+1);
%采樣
, ^% Q8 H7 X& e3 d' f9 f$ L9 _figure(1);
, q( L! P" L4 nsubplot(211);
- c3 k$ z3 R  u( L0 Mplot(t,bch0);%作正弦信號的時域波形
axis([0,N/fs,-0.4,0.4]);
' C7 q* F* d2 h* [' Q+ `* M- w1 e' R& Uxlabel('t');
ylabel('y');
title('bch0采樣t=[0,1.28]');
grid;
%進行FFT變換并做頻譜圖
y=fft(bch0,N);%進行fft變換
- ?, h3 g1 h( d& ]mag=abs(y);%求幅值
f=(0:length(y)-1)'*fs/length(y);%進行對應的頻率轉換
% d, N  H' h+ m* l8 z+ M" Sfigure(1);
% f9 [3 [: \* J' ksubplot(212);
plot(f,mag);%做頻譜圖
axis([0,100,0,4]);
1 X  [0 B  E) M( C! C7 ]* b& ]xlabel('頻率(Hz)');
ylabel('幅值');
title('bch0幅頻譜圖N=128');
& E! E! ?: M/ L$ u& Dgrid;

不好意思我寫在寫字板粘貼過來的不知道為什么標點沒了 現(xiàn)在改過來~

軸研小草

應該和轉子模型以及試驗工裝模型有關，跨距等。可傳一試驗工裝模型圖，以便分析。

黑天鬼

你的轉子工作頻率多少

黑天鬼

把數(shù)據(jù)下下來看看再說吧