我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图: ^" ~6 x/ D$ h; x 要得出的图像跟故障应该有的图像相似4 y% G1 U I' F g9 I# @ 并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征4 ]1 S3 C8 i! ?7 P+ I! h 我编的程序如下 你看看用得着不 + V' O7 t9 s- ^' v) g6 M: h* ?7 V. ]) X% _; T# e %***************1.bch0****************%0 N' W8 a8 f& i, |0 u/ O clc ; \/ Z; v) i, e0 e% F, i8 Eclf! n) }% V) d/ L6 W clear / S; i# v9 B4 v7 v Cbch0=load('F:\MATLAB7\lijiande\bch0.dat');, J m* C) y, u5 p8 v' J fs=100;%设定采样频率 ]3 Z7 Q1 n. B5 G! y3 a6 k0 `6 M' D- GN=128; " ^. \+ E2 U* }7 u% En=0:N-1; 2 I! O, S4 o3 p$ w+ h' }t=n/fs; ! \" V) w+ T, y/ N% Wbch0=bch0(n+1);. j# u* v2 W! d: P# ` %采样3 C# }; E) L& r figure(1); % z* D F5 d: S# G" h6 M- ]* [* rsubplot(211); 3 n: H5 x) b' U- hplot(t,bch0);%作正弦信号的时域波形+ H+ t9 w+ g, v) z5 V r axis([0,N/fs,-0.4,0.4]); ' U2 f# K( R9 l2 X% a5 i' ~( m0 Wxlabel('t'); 0 m. t2 T* u. J$ U- Wylabel('y');/ M' `0 `* B+ n9 ~' k0 q& c A5 ] title('bch0采样t=[0,1.28]'); 8 ]' F7 w, T% Qgrid;- N8 _. w9 M& ^% h6 b %进行FFT变换并做频谱图 9 i6 a; W# E) O% f. {y=fft(bch0,N);%进行fft变换 ; L6 e$ |/ O$ U. {! P) k' y, lmag=abs(y);%求幅值& h p: |. a. c9 y f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换4 Y0 C0 X% t% r figure(1); 9 |1 K! X$ i) @$ u2 G3 {subplot(212);4 D! X6 \; N e" y2 b4 F) Z6 |! k plot(f,mag);%做频谱图 : `3 ~& q- Z+ {" Oaxis([0,100,0,4]);3 }5 W) b; @ ~6 ^. Q xlabel('频率(Hz)'); 1 ~3 Z) m6 j+ M6 zylabel('幅值'); ! x" m& A6 [8 i: mtitle('bch0幅频谱图N=128');, e# D6 U6 p4 S6 J grid; : Z! N) ?6 P; Y g ; a# \+ K# A$ z) F不好意思我写在写字板粘贴过来的不知道为什么标点没了 现在改过来~ |