我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图
' v F3 g: q% l要得出的图像跟故障应该有的图像相似" k2 u! }' o+ P( L3 s
并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征9 v2 b4 U3 L! a( T* z
我编的程序如下 你看看用得着不
! q& q( f% n8 B8 p' s& f1 r2 G4 C$ e
%***************1.bch0****************%
2 ]4 N. q7 u( ~) j6 O$ jclc+ ?5 F/ E* @$ Y2 e7 ]- C$ Y
clf( i4 u" s0 l+ n% Q
clear
1 d; [1 F6 p4 ?' u* P; r6 \bch0=load('F:\MATLAB7\lijiande\bch0.dat');
; T/ H+ o3 g4 g" ?fs=100;%设定采样频率5 q' v' @/ s1 V; H0 @# c8 k$ f% G; T
N=128;
5 _! M, e; S6 B* q. c0 In=0:N-1;' r1 U5 p. p. _, c) T1 [6 p2 y6 U
t=n/fs;
0 Y) N' X! n& m, a+ k, {0 Jbch0=bch0(n+1);" r' a$ z" |' d6 ?" i5 e# q- \
%采样 V$ I9 N, ]! j* b6 t/ X
figure(1);+ L4 V8 o. y: B; |8 y
subplot(211);; A& r2 _8 S" ~- ^' \( |
plot(t,bch0);%作正弦信号的时域波形
# P( o' t( z+ t0 C5 ~( H1 Qaxis([0,N/fs,-0.4,0.4]);9 C5 u, O. Z; f0 a ^0 Q3 f
xlabel('t');' x! ]) k4 H4 D! D" a0 e6 ~, f5 K
ylabel('y');( H' w$ U: |+ R& x
title('bch0采样t=[0,1.28]');$ n# O5 ~3 ~" f1 g) [4 w
grid;5 X9 x/ ~& f, }* p' ?! h
%进行FFT变换并做频谱图! ~7 S% m3 ?% A, ]
y=fft(bch0,N);%进行fft变换
+ X, E3 L5 K/ f9 Hmag=abs(y);%求幅值
* ]/ n: h, U) g9 b# mf=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换
# j% x* P8 X0 F; Y4 M5 wfigure(1);1 o+ t: D5 a* d& O u
subplot(212);
) f9 C% l3 H* P% t& @plot(f,mag);%做频谱图
1 E3 z& m& [9 C0 G; h/ `' b; Daxis([0,100,0,4]);
1 \1 L5 G P) C) |& I- axlabel('频率(Hz)');+ @, F, n2 r/ w" `4 i
ylabel('幅值');9 [+ t0 k2 A' J5 G; ^7 M9 _/ V
title('bch0幅频谱图N=128');$ L, n; f) C) G: _
grid;: P8 {, r a8 t ]( K- { W) A T
! l& z, G+ L5 ^9 W3 f V
不好意思我写在写字板粘贴过来的不知道为什么标点没了 现在改过来~ |