我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图( ?0 P j K; J/ l3 N; B4 c2 u) ^
要得出的图像跟故障应该有的图像相似
' G: ?1 ^1 n/ V$ v5 S" a并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征
; A2 r5 m' W" {# P我编的程序如下 你看看用得着不4 p0 T# J/ \! r n
8 [; B8 u* p: } U
%***************1.bch0****************%
- y; a$ n( E7 bclc( B$ ` E+ F* N; q2 g7 K2 R; Q+ b
clf7 Q, ?. ~+ D) A6 }+ l- R: A
clear
( h& D5 ]9 A! Nbch0=load('F:\MATLAB7\lijiande\bch0.dat');5 [) q- K0 i& L( t
fs=100;%设定采样频率( m% `! i5 L3 t" ^. S% Y% u
N=128;5 \8 C" h, d) k; \& X5 Y* U; c
n=0:N-1;; E8 f! }1 }+ v# [# }
t=n/fs;
5 V4 G. r8 b3 A) }* Ybch0=bch0(n+1);
. w1 a& `! B! M; D9 r; Q# S%采样
- w3 e. Q4 h' ?8 p3 Afigure(1);
$ ]* e( P6 d# G8 [, _: X5 qsubplot(211);- O5 R7 e" `! Y* `3 Z
plot(t,bch0);%作正弦信号的时域波形
9 y- t: V" g. }& H* m" @axis([0,N/fs,-0.4,0.4]);
) V" @/ N" g# Q6 ]5 Sxlabel('t');
6 q$ f# |, e' \3 X* T4 a$ \6 kylabel('y');" K' j4 ~5 s2 y$ n" w
title('bch0采样t=[0,1.28]');
4 f9 ^' F5 v# W6 r4 f8 f4 Vgrid;/ m* `3 M* I$ J
%进行FFT变换并做频谱图
- i/ q1 k& t) U% R1 jy=fft(bch0,N);%进行fft变换3 n( s1 h* z I3 E+ ~
mag=abs(y);%求幅值
. v+ J1 k" Q% |! y% ` tf=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换
; E; \0 `+ @& ~figure(1);
# r0 _* [, r7 e' h7 x4 f- r, r/ ssubplot(212);
' k0 O# l2 ~" w# r9 V1 [plot(f,mag);%做频谱图
0 q, z- u: l* n, c5 N7 Taxis([0,100,0,4]);0 j# r: L: n" x5 s, u( |' l
xlabel('频率(Hz)');- w- n( B6 N0 U
ylabel('幅值');
. i( L9 W. H" O9 Ytitle('bch0幅频谱图N=128');
4 E2 K5 a! w$ n1 e8 |4 Dgrid;
+ W4 g# i: R0 n( `
6 [2 n8 e1 a' u4 E不好意思我写在写字板粘贴过来的不知道为什么标点没了 现在改过来~ |