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

超級小寶

我有實驗數(shù)據(jù)是dat格式的 要用matlab分析出頻譜圖跟時域圖要得出的圖像跟故障應(yīng)該有的圖像相似并且?guī)臀医忉屢幌?轉(zhuǎn)子不平衡的時域圖 頻譜圖各是什么特征我編的程序如下 你看看用得著不%***************1.bch0****************%clcclfclearbch0=load('F:\MATLAB7\lijiande\bch0.dat');fs=100;%設(shè)定采樣頻率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);%進行對應(yīng)的頻率轉(zhuǎn)換figure(1);subplot(212);plot(f,mag);%做頻譜圖axis([0,100,0,4]);xlabel('頻率(Hz)');ylabel('幅值');title('bch0幅頻譜圖N=128');grid;
謝謝大家了~�。。�

超級小寶

我有實驗數(shù)據(jù)是dat格式的 要用matlab分析出頻譜圖跟時域圖
要得出的圖像跟故障應(yīng)該有的圖像相似
并且?guī)臀医忉屢幌?轉(zhuǎn)子不平衡的時域圖 頻譜圖各是什么特征
我編的程序如下 你看看用得著不
" |/ ~( n, h" ~, i; @/ a, O' V
%***************1.bch0****************%
) M' S* n- e8 dclc
clf
clear
8 E5 Y( o% P- o+ S+ X8 D1 G1 O  Fbch0=load('F:\MATLAB7\lijiande\bch0.dat');
fs=100;%設(shè)定采樣頻率
N=128;
n=0:N-1;
t=n/fs;
4 C. w& k4 D! t; t  H6 b% y0 {bch0=bch0(n+1);
%采樣
' {+ O. R6 K* [1 ]2 G3 Xfigure(1);
. L9 c. d" Q5 ]" W4 I# }* G/ fsubplot(211);
plot(t,bch0);%作正弦信號的時域波形
( _4 r, k8 S2 G6 T; v  p' Jaxis([0,N/fs,-0.4,0.4]);
6 o" _$ e" L! H  M6 K' Q8 }xlabel('t');
/ R' }. y1 J! z2 Q$ Tylabel('y');
8 K0 o4 _- I: E2 }- {title('bch0采樣t=[0,1.28]');
grid;
%進行FFT變換并做頻譜圖
- z; d* g# G# J5 y4 C% t  t3 sy=fft(bch0,N);%進行fft變換
* s4 k9 F- S  Hmag=abs(y);%求幅值
f=(0:length(y)-1)'*fs/length(y);%進行對應(yīng)的頻率轉(zhuǎn)換
$ u# K, T4 Y' {% O- lfigure(1);
subplot(212);
plot(f,mag);%做頻譜圖
+ I  }4 u- Z$ s: P7 waxis([0,100,0,4]);
6 k3 N  d2 i' h. p% ^xlabel('頻率(Hz)');
ylabel('幅值');
title('bch0幅頻譜圖N=128');
grid;
# h1 e, a5 @+ V7 ?0 ~7 t
不好意思我寫在寫字板粘貼過來的不知道為什么標點沒了 現(xiàn)在改過來~

軸研小草

應(yīng)該和轉(zhuǎn)子模型以及試驗工裝模型有關(guān)，跨距等�？蓚饕辉囼灩ぱb模型圖，以便分析。

黑天鬼

你的轉(zhuǎn)子工作頻率多少

黑天鬼

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