我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图
L) s$ N* O+ n2 ^8 g8 S: V# w要得出的图像跟故障应该有的图像相似; z) q- u: Y1 ]
并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征7 L2 f) [7 i8 D" t1 Y6 \" f
我编的程序如下 你看看用得着不
3 W& @$ k* g1 v8 }- g$ H& h6 E. e0 o# n+ ?' r" g: M! y+ S6 S
%***************1.bch0****************%
% [9 E3 j0 E: Z Z- b" y lclc% H( ~' A# S: s5 _% {
clf
' R; Q6 G, D+ t ~( t' Kclear
# _0 m9 u. x5 W6 N8 w9 d/ [bch0=load('F:\MATLAB7\lijiande\bch0.dat');
) w/ u1 u5 X7 @fs=100;%设定采样频率
. c/ _; q( }8 ~N=128;
: u H$ t' o5 p7 e( U, Un=0:N-1;- u* C) h { q9 a/ ~ Y" c
t=n/fs;
3 y2 x& D3 @1 U/ u. Z4 c5 fbch0=bch0(n+1);
J! L! }! F' J0 g%采样/ w: C: t, D/ o9 V' W* L
figure(1);2 e* E+ P ], X' K
subplot(211);! L" W+ ^2 S) ?: Q; C: C H% O
plot(t,bch0);%作正弦信号的时域波形
3 _5 \! b. N, O" x) Oaxis([0,N/fs,-0.4,0.4]);
2 B5 | U \7 R* u, M% Wxlabel('t');
. O. I2 @' f* d' h8 I6 c; \" K" L, Wylabel('y');
4 \4 Y8 |, G8 z% l, gtitle('bch0采样t=[0,1.28]');) J' T1 `$ ?( l+ u; u, m
grid;% H7 S4 k1 A0 u# ?8 s+ Q
%进行FFT变换并做频谱图9 u' |+ {, i9 R0 L$ W7 }2 O6 P
y=fft(bch0,N);%进行fft变换; f9 B5 A9 v9 A8 J& _3 v, }, n
mag=abs(y);%求幅值
$ v, c5 k+ S C4 X; df=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换6 Z" R+ a) `! S4 @/ x6 p+ i( [' P
figure(1);
9 ]* p; P, `/ B2 k$ X' F- A3 ~subplot(212);
0 ?, D# `1 z, O9 e2 ]plot(f,mag);%做频谱图7 i% a# h- Q1 \
axis([0,100,0,4]);6 G+ ]5 B5 L, b, O
xlabel('频率(Hz)');
! U6 O7 R$ I; R2 Q: H1 Tylabel('幅值');
; P7 G3 C* P4 w2 L/ ?title('bch0幅频谱图N=128'); t# V( G, c" Q# u8 p
grid;: b) T/ n3 h" b* M
; I. `% f" g8 G3 h9 S
不好意思我写在写字板粘贴过来的不知道为什么标点没了 现在改过来~ |