我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图 & I8 R9 w9 z+ I; {5 t" j( H/ {' W要得出的图像跟故障应该有的图像相似, \+ g( T7 {/ X 并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征5 h: N: _' R% \ 我编的程序如下 你看看用得着不0 `1 K, o9 J4 `. a
& f9 `! t1 i# [%***************1.bch0****************% ' Y; a9 d6 [% p7 ?: Fclc 6 p8 j) W: C+ g; \+ Z% N/ u( Cclf 5 L( v* K+ ]# x c: Y8 b+ uclear / a" m8 c" m1 j( V! Ubch0=load('F:\MATLAB7\lijiande\bch0.dat');% `& a( f! U6 Z7 a fs=100;%设定采样频率 % X" b4 R5 i/ _ lN=128;1 ]$ j e, c9 Z3 c8 l: k1 f- U n=0:N-1;5 Z# m# X7 d$ ]) b- F* L t=n/fs;3 V* s+ }& L) N+ X' `6 ]' u9 J. y bch0=bch0(n+1); ' N5 Z* A7 M$ R. b. T/ Q%采样 , W9 f: Q+ e* wfigure(1);7 s0 I4 U4 Y$ H# x. N subplot(211);& g8 y0 Z) v2 ^' i- M1 N1 p plot(t,bch0);%作正弦信号的时域波形7 G% ?7 ^9 f7 g( w. T9 O axis([0,N/fs,-0.4,0.4]);' Z* r1 ] A+ h+ g xlabel('t'); , m. m) E k5 ^' W2 |( Hylabel('y');- Y% ^9 _% |: \+ f8 v title('bch0采样t=[0,1.28]');; y& d9 U( i5 F Y S grid; 2 c; _2 ^" l4 O1 G& }+ l& m%进行FFT变换并做频谱图 : O! [& Z8 k* uy=fft(bch0,N);%进行fft变换4 K# h/ x- ?4 {9 s2 \* A Q mag=abs(y);%求幅值 $ e! h, k% Q* y+ z; ^; jf=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换* h0 q+ ]& { b+ i figure(1); 4 T8 }2 I' \+ F" y, F) t6 m/ Z$ Dsubplot(212); W }% L m$ w* K( f2 Wplot(f,mag);%做频谱图 # Z- B; [9 l Faxis([0,100,0,4]);7 Y5 b8 G+ O* N# c+ R1 n xlabel('频率(Hz)'); 6 Z* y* F) I& t* a, w. E7 ^4 S$ tylabel('幅值');T6 s9 m q7 A0 ^, N0 i8 `, Y } title('bch0幅频谱图N=128');2 d; W" }4 ~ `* f grid; F9 {5 E' Q6 L" w3 y9 V0 N) r# l- n% c5 f+ I% w 不好意思我写在写字板粘贴过来的不知道为什么标点没了 现在改过来~ |