我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图
$ i1 ^* S; s" j, C4 l要得出的图像跟故障应该有的图像相似, T' m/ i% Z: v$ F1 l. w& i y
并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征
' p( Z& t6 W9 X( i3 G, K( G+ i- Z H我编的程序如下 你看看用得着不0 ^' Y% D1 ]! L' T# D/ h+ r
3 S; v4 `: m$ G0 g: |; d9 w%***************1.bch0****************%
' ]5 l C8 J2 \7 Eclc2 p& N- z% Z0 I3 t$ m0 J6 u
clf4 U E: y; M% I2 q$ w$ H
clear+ m# @+ W6 o2 N
bch0=load('F:\MATLAB7\lijiande\bch0.dat');
$ n' w4 y/ m* j, C! v4 ^6 Ifs=100;%设定采样频率0 a3 H; [$ ] W) T
N=128;
1 [0 [5 e* D. B: u3 Nn=0:N-1;) i' A, `. W9 L% T# u. p% Q+ X" w
t=n/fs;8 T: q& a2 Z# ~4 S! R" S
bch0=bch0(n+1);4 i+ w6 C1 ^% ^" k' U
%采样( S- k3 ?5 y4 [ }2 p
figure(1);: g, _& B |: M9 y
subplot(211);
/ W2 t3 |- k3 {; {' ~0 ^+ Lplot(t,bch0);%作正弦信号的时域波形% h" k: y8 _6 _, {: o" g
axis([0,N/fs,-0.4,0.4]);. u/ o& T# o6 o9 J. N( h2 s+ ]6 @8 y
xlabel('t');- }; ~+ Y1 G$ F. [) N( \5 `( v2 L
ylabel('y');
7 Y$ z9 h* F6 i7 }0 P& r. _6 \+ o3 Ztitle('bch0采样t=[0,1.28]');9 t4 T( {8 U3 t
grid;( }1 u% {& i: d( d
%进行FFT变换并做频谱图1 x1 q* q( K" @
y=fft(bch0,N);%进行fft变换, ?5 g" E* Z2 B
mag=abs(y);%求幅值
5 A! X) _4 P- Xf=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换
$ i) L$ h+ K3 sfigure(1);6 p! b5 I+ p; h0 r
subplot(212);
- e b1 h* @2 j6 J( }plot(f,mag);%做频谱图
6 h/ C# o/ n! h8 H9 {3 Uaxis([0,100,0,4]);( y& C: ?% k/ ?- g/ ]9 f
xlabel('频率(Hz)');
, j9 P% S( {$ Q- ]( g/ rylabel('幅值');. C4 t& W+ n$ i% @# z
title('bch0幅频谱图N=128');7 f! {) z0 {& i
grid;- P0 G3 |2 n) J) }4 ^5 p: [5 b5 b
! ~$ Z. X: p8 c7 Y2 k5 M不好意思我写在写字板粘贴过来的不知道为什么标点没了 现在改过来~ |