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

超級(jí)小寶

我有實(shí)驗(yàn)數(shù)據(jù)是dat格式的 要用matlab分析出頻譜圖跟時(shí)域圖要得出的圖像跟故障應(yīng)該有的圖像相似并且?guī)臀医忉屢幌?轉(zhuǎn)子不平衡的時(shí)域圖 頻譜圖各是什么特征我編的程序如下 你看看用得著不%***************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);%作正弦信號(hào)的時(shí)域波形axis([0,N/fs,-0.4,0.4]);xlabel('t');ylabel('y');title('bch0采樣t=[0,1.28]');grid;%進(jìn)行FFT變換并做頻譜圖y=fft(bch0,N);%進(jìn)行fft變換mag=abs(y);%求幅值f=(0:length(y)-1)'*fs/length(y);%進(jìn)行對(duì)應(yīng)的頻率轉(zhuǎn)換figure(1);subplot(212);plot(f,mag);%做頻譜圖axis([0,100,0,4]);xlabel('頻率(Hz)');ylabel('幅值');title('bch0幅頻譜圖N=128');grid;
謝謝大家了~�,。,�,！

超級(jí)小寶

我有實(shí)驗(yàn)數(shù)據(jù)是dat格式的 要用matlab分析出頻譜圖跟時(shí)域圖
要得出的圖像跟故障應(yīng)該有的圖像相似
; J$ F1 i; G/ c; ~! N3 g并且?guī)臀医忉屢幌?轉(zhuǎn)子不平衡的時(shí)域圖 頻譜圖各是什么特征
我編的程序如下 你看看用得著不
$ ^/ e; u4 a7 h/ t
%***************1.bch0****************%
6 Z4 Q. K) h$ g$ a* Y, [8 mclc
- T7 [0 M- ~' @& q* T; Xclf
clear
' w, k8 u( K) Cbch0=load('F:\MATLAB7\lijiande\bch0.dat');
fs=100;%設(shè)定采樣頻率
$ y# w: D/ t1 DN=128;
n=0:N-1;
t=n/fs;
bch0=bch0(n+1);
; N0 \6 C( A( D" R%采樣
3 |& ~/ U' U8 s5 E! D# N$ S0 A/ E- xfigure(1);
subplot(211);
plot(t,bch0);%作正弦信號(hào)的時(shí)域波形
axis([0,N/fs,-0.4,0.4]);
, X5 x( L* [) Q# y* ]* Xxlabel('t');
  p: f9 l- Q9 o: v3 K4 q. mylabel('y');
title('bch0采樣t=[0,1.28]');
' A8 J/ I3 w8 a/ ygrid;
%進(jìn)行FFT變換并做頻譜圖
y=fft(bch0,N);%進(jìn)行fft變換
) n1 B' l; J5 ]+ w! U& bmag=abs(y);%求幅值
f=(0:length(y)-1)'*fs/length(y);%進(jìn)行對(duì)應(yīng)的頻率轉(zhuǎn)換
figure(1);
subplot(212);
plot(f,mag);%做頻譜圖
4 l$ v% ~& B. C: d$ ?# Q: d2 eaxis([0,100,0,4]);
xlabel('頻率(Hz)');
ylabel('幅值');
title('bch0幅頻譜圖N=128');
grid;
, Z, X' {* |) t8 Z# _, ]
不好意思我寫在寫字板粘貼過來的不知道為什么標(biāo)點(diǎn)沒了 現(xiàn)在改過來~

軸研小草

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

黑天鬼

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

黑天鬼

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