我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图
0 z o' |9 B, J' d3 `要得出的图像跟故障应该有的图像相似
* [. [. I5 T+ b5 K& D5 e并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征
/ y" G6 ?: M$ O+ R# R% y3 G7 y4 i我编的程序如下 你看看用得着不# ~$ ~( g( S4 O$ }( r
3 U$ y5 Q* l8 j2 M/ d% X%***************1.bch0****************%7 v2 N5 B$ [- n8 w2 Z
clc
. m) ^. B$ ?# K' H- l# x! q6 b0 ~+ |clf
, [: P/ t" k L) e. ]8 qclear0 \- M& `8 v% a) `& H
bch0=load('F:\MATLAB7\lijiande\bch0.dat');
1 K: M2 \* @ |fs=100;%设定采样频率
, L% K$ L' W+ q# E0 Y8 D: ]N=128;2 w& p& D* N2 J1 ^ V2 O5 u
n=0:N-1;
7 S1 D; c M* }1 _2 N: W) Nt=n/fs;+ c# p" L/ S; g- S
bch0=bch0(n+1);9 w$ V: } N; ]4 T
%采样; c/ P [% \) K& e" C4 n; i- E
figure(1);1 y# U" D4 d( P: o" V: z* s* R: ?( S
subplot(211);
1 O R/ V! q4 [3 |plot(t,bch0);%作正弦信号的时域波形6 G+ j: ?3 Z% ~3 |1 N
axis([0,N/fs,-0.4,0.4]);
{; H8 C" ^9 P! k4 r3 h5 Jxlabel('t');2 \: P- ]6 ]% d4 L. H
ylabel('y');
% u9 z% J; B0 z/ }+ ]- |) K; b( u& P$ ?# Ititle('bch0采样t=[0,1.28]');
4 l H4 o9 {0 i$ S% f8 s% J7 bgrid;! R' S* L0 `) S# h! U; ~$ e
%进行FFT变换并做频谱图
0 e. J% I+ _* o& m! m0 @: }y=fft(bch0,N);%进行fft变换
/ [& v+ ~( o" h( ^mag=abs(y);%求幅值
* f6 ~3 e# @; I! v: N: cf=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换
7 h) u8 j _% S. j8 R6 z6 Dfigure(1);" N0 k2 e- q& Z4 F, O
subplot(212);
: S7 A# ]9 f1 t" ^7 Xplot(f,mag);%做频谱图
$ T. s. @5 Y4 h" raxis([0,100,0,4]);! O7 T. h# F( ?. {& p
xlabel('频率(Hz)');
$ X6 ?- Q3 n3 U2 y3 fylabel('幅值');
2 d0 m- { Z- W- m7 [/ btitle('bch0幅频谱图N=128');) H( i r. \8 O& r
grid;
2 s$ {3 p6 l$ O0 {0 K
/ B1 @9 i; s5 i不好意思我写在写字板粘贴过来的不知道为什么标点没了 现在改过来~ |