matlab處理實驗數據 實驗：利用轉子試驗臺測試轉子不平衡以及渦動狀態 故障診斷

超級小寶

我有實驗數據是dat格式的 要用matlab分析出頻譜圖跟時域圖要得出的圖像跟故障應該有的圖像相似并且幫我解釋一下 轉子不平衡的時域圖 頻譜圖各是什么特征我編的程序如下 你看看用得著不%***************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;
謝謝大家了~！！！

超級小寶

我有實驗數據是dat格式的 要用matlab分析出頻譜圖跟時域圖
! Q8 o2 M8 k1 \要得出的圖像跟故障應該有的圖像相似
. l* c$ x; n! U+ W并且幫我解釋一下 轉子不平衡的時域圖 頻譜圖各是什么特征
我編的程序如下 你看看用得著不
% w  s# ~9 [2 S" w% K
! l/ W# u- B% d$ v: J, h6 g6 {%***************1.bch0****************%
% W. p. q* J; |) |; Gclc
clf
clear
6 s  r' C& a8 w. C9 }bch0=load('F:\MATLAB7\lijiande\bch0.dat');
9 v2 i/ s) J( t5 ~fs=100;%設定采樣頻率
# q: w, p4 H* MN=128;
n=0:N-1;
8 a! z, O) l7 it=n/fs;
8 k- }0 d7 B5 w5 i% X3 H: t' J# Wbch0=bch0(n+1);
%采樣
/ `1 G6 O% i" Z+ l* h2 T# Ofigure(1);
& i* r; x! O1 h$ [, ~, gsubplot(211);
plot(t,bch0);%作正弦信號的時域波形
: P/ Q/ W9 z; o' z& c. d( saxis([0,N/fs,-0.4,0.4]);
( d% I- g6 ~; @xlabel('t');
ylabel('y');
title('bch0采樣t=[0,1.28]');
& `8 W& x2 a2 O; `1 R. M4 J4 q+ s( Zgrid;
%進行FFT變換并做頻譜圖
7 P( w0 U1 C7 ], u6 |y=fft(bch0,N);%進行fft變換
) }7 u5 `7 R) b. C2 b1 Y1 W, Kmag=abs(y);%求幅值
% ?$ m" b0 ^2 ~/ w$ ~$ l% Af=(0:length(y)-1)'*fs/length(y);%進行對應的頻率轉換
figure(1);
subplot(212);
/ |0 Z8 n" v0 Bplot(f,mag);%做頻譜圖
axis([0,100,0,4]);
xlabel('頻率(Hz)');
ylabel('幅值');
title('bch0幅頻譜圖N=128');
0 f9 C+ h1 z/ l: c. Ogrid;
  B. I. B: V& W( |7 U( X$ m
不好意思我寫在寫字板粘貼過來的不知道為什么標點沒了 現在改過來~

軸研小草

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

黑天鬼

你的轉子工作頻率多少

黑天鬼

把數據下下來看看再說吧