我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图
3 Y8 L7 [! R, S# f) ~要得出的图像跟故障应该有的图像相似
+ |% j q4 u/ L0 P并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征
# }) A4 L M" I, [ i我编的程序如下 你看看用得着不
; ?, W2 M& e7 P0 g! S' U; @: M6 ]) {$ `$ [; V; ?8 |$ Z# h" t
%***************1.bch0****************%
, e7 R4 B3 F2 Q+ K% e7 Nclc$ w, x d- @3 x) r7 E8 F( |
clf
1 q! b# w3 F) T0 ?$ lclear" d% Z, v/ K/ x& g
bch0=load('F:\MATLAB7\lijiande\bch0.dat');
- G3 L! g! Y% Afs=100;%设定采样频率1 r N9 @/ S$ a0 h& _
N=128;/ Z: _0 @4 M% F% C7 `
n=0:N-1;& S( `5 e2 X$ z9 b- I
t=n/fs;
$ ]1 a; J, v, m& U d9 Fbch0=bch0(n+1);
! }+ G, o* K% J& W' q9 K- Y%采样+ g5 k( d/ `3 f* p; z: V7 Y# d
figure(1);
+ s% H5 ?4 U& [( ^subplot(211);
; ?# f% B7 ?6 `3 X; Nplot(t,bch0);%作正弦信号的时域波形
$ ~4 G8 _0 v/ Q0 o2 Eaxis([0,N/fs,-0.4,0.4]);
2 `3 Z5 o5 }9 c7 ^ L; [xlabel('t');' A) a. T3 b6 ~) G. N/ ?$ r
ylabel('y');
$ [$ f+ \ T/ G! z# k) wtitle('bch0采样t=[0,1.28]');
0 P0 g/ d9 d/ D' B9 x+ `, ggrid;: S ~) n7 N2 L: W6 P: V/ B
%进行FFT变换并做频谱图
1 E7 }3 v8 o6 R5 B% ?' ^$ }y=fft(bch0,N);%进行fft变换$ t/ ] z7 ^1 C: A
mag=abs(y);%求幅值
% Y. G7 A7 o' B+ Cf=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换" }% k) h$ z, m3 M9 q
figure(1);- ~! ^0 I Y* U9 Z; n; G& Y8 _
subplot(212);- v; E, F0 b3 G: h* T" E I
plot(f,mag);%做频谱图
1 x' B+ o7 B/ v' X+ Yaxis([0,100,0,4]);1 N% B2 B: \9 A8 l' d& G
xlabel('频率(Hz)');4 c0 Z: J' ^; i0 ?( A
ylabel('幅值');- ^# k' v8 P5 D7 H: y1 r
title('bch0幅频谱图N=128');: V5 A, t1 u& S* X
grid;* r, B: @0 Z! \3 W* {, l% b2 b
2 ^ H+ W9 a4 P( j6 z) [
不好意思我写在写字板粘贴过来的不知道为什么标点没了 现在改过来~ |