我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图
8 Y, b7 {2 r& V要得出的图像跟故障应该有的图像相似
/ U# | V" w+ e \- c" J/ Y; R; D并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征
$ _- G5 B& Q; X: n" h我编的程序如下 你看看用得着不6 J, v( i% _& c# I5 Y' P# z
& r; {0 y& W& S7 Q6 J
%***************1.bch0****************%( v) J% W5 }0 V# ~1 t; E( o0 N
clc- F3 n0 N+ N. n8 I2 J! R3 `
clf
2 x' B. ]# k. X( nclear' E. _( H2 Q( m! X) z3 k
bch0=load('F:\MATLAB7\lijiande\bch0.dat');5 i* f6 D' @3 H. @! X' g
fs=100;%设定采样频率
6 F- V% q- w. l* `N=128;
# o. ~! o' j; [4 j4 Pn=0:N-1;2 `1 n6 L1 ^$ i) K
t=n/fs;
6 e% P6 Y6 X' Bbch0=bch0(n+1);7 q# _1 @: G4 |0 V
%采样0 {# z) r) F/ D( J! i' {1 K6 h
figure(1);9 [/ ^2 W2 c0 V+ `. K) d3 \
subplot(211);( f! d! g: l" w
plot(t,bch0);%作正弦信号的时域波形
: D# z# Q; [3 N3 C8 a. }# F- \axis([0,N/fs,-0.4,0.4]);
- o, H. J$ \( Sxlabel('t');4 Y2 \" z0 i5 V) e9 o" O% I
ylabel('y');" W$ ~9 @% |5 `( W5 `& {
title('bch0采样t=[0,1.28]');
1 l+ ?3 m, z7 Y/ Q2 [% H" hgrid;
2 O2 o# u& B. w. p4 R* P- G; l%进行FFT变换并做频谱图
! K: ^, h4 R' n" f! Jy=fft(bch0,N);%进行fft变换
% ^0 |1 l5 H E' a! hmag=abs(y);%求幅值1 R, a. \3 t Y3 D- I4 r5 d
f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换
, V I+ p! I* l) _' Y+ `figure(1);* O7 A' ~) ]% ]/ t$ d, A
subplot(212);" X p( f$ G4 ^3 k" N1 \# A1 _# ^; r
plot(f,mag);%做频谱图
/ E; c6 j( v' T9 o2 P: W* Iaxis([0,100,0,4]);, Q; G3 U/ [+ B1 i3 ^
xlabel('频率(Hz)');9 {+ _9 N! C, Q4 E1 ^
ylabel('幅值');% I9 t2 n" R0 C X: D) e
title('bch0幅频谱图N=128');) K% ?% a: f. k( H
grid;
1 M, @' i$ c! a! _. ^; N- d3 X: G
不好意思我写在写字板粘贴过来的不知道为什么标点没了 现在改过来~ |