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;
) f2 H) R, L, q% e謝謝大家了~！！！

超級小寶

我有實驗數(shù)據(jù)是dat格式的 要用matlab分析出頻譜圖跟時域圖
1 c$ [; ?, |* v4 M要得出的圖像跟故障應該有的圖像相似
并且?guī)臀医忉屢幌?轉子不平衡的時域圖 頻譜圖各是什么特征
我編的程序如下 你看看用得著不

1 q, G  R" }7 U%***************1.bch0****************%
clc
clf
clear
bch0=load('F:\MATLAB7\lijiande\bch0.dat');
fs=100;%設定采樣頻率
) S9 j! n/ Q1 p, c) S, @/ Y: Z. ?N=128;
. A5 V* i% p$ m, j: `* j. d( F7 fn=0:N-1;
t=n/fs;
1 u7 l2 V4 I7 y8 n' T" r6 j3 q' Abch0=bch0(n+1);
  o& c. r! q& h, B# g% C6 S* i%采樣
figure(1);
subplot(211);
plot(t,bch0);%作正弦信號的時域波形
axis([0,N/fs,-0.4,0.4]);
xlabel('t');
( c6 f- |# j, j& X; X, p. T4 Q: zylabel('y');
% a4 ~$ H" }% r# ytitle('bch0采樣t=[0,1.28]');
' ]7 y. |7 [( }3 G2 r# \0 P) Kgrid;
%進行FFT變換并做頻譜圖
+ d' ]% R/ m# H7 O. l& K# y7 ky=fft(bch0,N);%進行fft變換
mag=abs(y);%求幅值
3 u! m- e: q/ s7 |% J- O, {( }( K2 Hf=(0:length(y)-1)'*fs/length(y);%進行對應的頻率轉換
+ c, z9 S: U' h, t/ w7 a& cfigure(1);
subplot(212);
plot(f,mag);%做頻譜圖
' |2 m# `* @8 F, A4 ?7 }: L  P. c5 D" @axis([0,100,0,4]);
xlabel('頻率(Hz)');
ylabel('幅值');
title('bch0幅頻譜圖N=128');
grid;

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

軸研小草

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

黑天鬼

你的轉子工作頻率多少

黑天鬼

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