我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图 N2 y6 v5 c1 R0 n6 }# u, e要得出的图像跟故障应该有的图像相似 ; `1 \$ a5 E: S7 F, N/ A并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征3 c2 I6 [ ^6 O 我编的程序如下 你看看用得着不) ~9 S) b6 N) h" O" P
6 y+ h3 s+ k% F' j- F. E5 w%***************1.bch0****************% ' ?+ u5 @, Y' f8 O/ H& Q* gclc5 v: k |) Q% L( w& N% O clf' h( d9 A5 }; P6 [. X' [, k clear3 G/ z$ [+ S" X& J U D bch0=load('F:\MATLAB7\lijiande\bch0.dat'); - P+ f' X& l% @, l+ a7 ^, ~: Ufs=100;%设定采样频率 % n, p+ j2 V9 `5 P8 L' xN=128;1 S$ _/ l$ x _) S n=0:N-1; " T1 Q& V( k' k& T; Y9 _3 a5 [t=n/fs; , q# I3 }# J1 ?; obch0=bch0(n+1); 9 Y. ]- q. d6 {2 ~* }& K$ e, u%采样 ; w, m( J" v6 h3 nfigure(1); 6 e' E9 g8 E) @8 o, D5 O) S& qsubplot(211); # P/ N" @6 |/ r& Vplot(t,bch0);%作正弦信号的时域波形# `; j4 E8 K, P8 Q# n axis([0,N/fs,-0.4,0.4]);9 J" Y( I. L( \" c xlabel('t');p+ m7 f0 K/ a2 d ylabel('y'); 8 @* ]' I$ J G& I! ptitle('bch0采样t=[0,1.28]'); # c) \7 u, H. zgrid;, O! v. _7 T _ %进行FFT变换并做频谱图! E" C5 L- A: I0 f3 H# K y=fft(bch0,N);%进行fft变换 ' O E, U8 D* ^8 R; D/ Dmag=abs(y);%求幅值 # U- J, N7 j: T7 ~) z% Ff=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换 # s7 H8 L9 d$ a* Wfigure(1);% k0 P% b" f9 T: G( O/ g subplot(212); $ q. r0 H( _' y/ U4 L; g S# Q3 @plot(f,mag);%做频谱图1 C2 |9 X& Z+ U9 ]+ w/ z5 n3 u axis([0,100,0,4]);; Y2 k( w# M6 n# U; B8 T% U xlabel('频率(Hz)'); X+ D# h- d5 Z& k( ^ylabel('幅值'); : Q$ W! H; P: J2 r2 utitle('bch0幅频谱图N=128');5 }3 c# |6 ^- K! K grid; + P; n- h( E/ N* R; A4 M. h6 w- t/ N& h# e- Y& s5 f: G 不好意思我写在写字板粘贴过来的不知道为什么标点没了 现在改过来~ |