|
2#
樓主 |
發(fā)表于 2011-5-26 11:47:56
|
只看該作者
我有實(shí)驗(yàn)數(shù)據(jù)是dat格式的 要用matlab分析出頻譜圖跟時(shí)域圖
5 S; h$ l+ I! m要得出的圖像跟故障應(yīng)該有的圖像相似
1 W* J" A/ l3 P/ q. h) n并且?guī)臀医忉屢幌?轉(zhuǎn)子不平衡的時(shí)域圖 頻譜圖各是什么特征# L4 X% k; p7 y2 O) m
我編的程序如下 你看看用得著不
4 _, i) h& M& \* G' S m+ U7 \3 x# T$ u) n- ]+ s
%***************1.bch0****************%
0 n5 _. P* r" f' @5 A0 D/ V6 w+ jclc& Q Z# s3 k$ V
clf
4 A1 ]9 R0 A- G% e3 U' W. O' Iclear. g0 x& a. @8 Q$ n! M6 d
bch0=load('F:\MATLAB7\lijiande\bch0.dat');3 ^' O! e. h' M+ ? f( w7 E
fs=100;%設(shè)定采樣頻率
# s8 y$ U8 p& h$ l/ U: Q5 gN=128;
2 m4 [0 R9 z& w. P+ on=0:N-1;* ?2 {7 ^% p/ x6 j- B' ] d
t=n/fs;
' q8 u; J$ M6 e" v7 Q' Xbch0=bch0(n+1);
& ^4 U4 q$ u, Y/ e7 h%采樣( e. o' u4 G3 }# J
figure(1);$ Z( F1 f7 r5 z, `! ?
subplot(211);
+ |- [4 x x* M+ Q! vplot(t,bch0);%作正弦信號(hào)的時(shí)域波形1 e1 {% w$ a0 X- L/ ~& ]
axis([0,N/fs,-0.4,0.4]);
8 e2 H+ k, Y5 O8 r: P$ Xxlabel('t');
* x) A' M4 @* S* P0 Pylabel('y');
) b! V! r+ j0 E( q0 u' Atitle('bch0采樣t=[0,1.28]');* s5 w5 e9 n' e. ?6 P- t* b+ ~
grid;7 R9 s& `+ V- ]& a( P
%進(jìn)行FFT變換并做頻譜圖% t$ L7 z2 o7 Z1 K% |9 ^/ Q
y=fft(bch0,N);%進(jìn)行fft變換
5 q$ \/ f3 I& J) W/ jmag=abs(y);%求幅值! S* Y1 a! c5 R2 {/ g
f=(0:length(y)-1)'*fs/length(y);%進(jìn)行對(duì)應(yīng)的頻率轉(zhuǎn)換* m8 ^3 f4 k, \
figure(1);
! E2 \% P7 h/ ?; g9 i3 g) csubplot(212);" X' h+ q3 l( x' A3 A: {7 u# J
plot(f,mag);%做頻譜圖& ]; a+ H) _ w3 e) b( u+ z
axis([0,100,0,4]);
" n0 s( _4 c0 W: zxlabel('頻率(Hz)');
F |% t f. l8 s4 ^' ? @6 zylabel('幅值');2 |( r1 B" b- O; m$ \, Q" _
title('bch0幅頻譜圖N=128');8 f; M) ~. e' p5 N0 o+ J; H* H
grid;) K/ T& S5 {/ ~ e
% s; Z9 _# m- S9 M& w
不好意思我寫在寫字板粘貼過(guò)來(lái)的不知道為什么標(biāo)點(diǎn)沒(méi)了 現(xiàn)在改過(guò)來(lái)~ |
|