我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图) V( s9 N6 w3 K: c1 m: s8 i! K0 L
要得出的图像跟故障应该有的图像相似8 R0 @* p. o. i" L/ x) G% X/ r
并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征! M9 P" ^6 @; D3 C" ~2 T2 o( `0 B# E
我编的程序如下 你看看用得着不
% r- Z" Q0 l4 x, W7 N$ i- J; V8 n# y' D5 q8 k$ F$ n
%***************1.bch0****************%$ y+ @; M7 M& R( r b
clc
+ ]7 K: m; Z" Q* _1 m gclf
/ P' `5 V8 P8 x. J" L- Oclear! H) X$ H+ R8 g
bch0=load('F:\MATLAB7\lijiande\bch0.dat');
' `" J& x9 \) h @4 K5 Ufs=100;%设定采样频率# \# @% ?: t5 S6 h8 b5 z
N=128;) x- Q i' }8 [% u% A
n=0:N-1;! t& ^; d m& v' }& q% V& x: }
t=n/fs;' g) [- @1 S6 y0 C. {. ^( u4 V$ c
bch0=bch0(n+1);
1 h: F/ W$ d/ Y5 \%采样6 l- w9 @9 T% r3 g6 G
figure(1);. U# a5 K# z, z
subplot(211);
9 d$ y" Y, f, k, d2 ?plot(t,bch0);%作正弦信号的时域波形7 g, q1 x* I: T& C
axis([0,N/fs,-0.4,0.4]);
1 O# _. r8 m |xlabel('t');
% D5 C' G) G" c& u7 R- Gylabel('y');0 s9 H b9 b! ^9 \8 j3 m
title('bch0采样t=[0,1.28]');
* t+ S* f2 {6 T% x- X% Dgrid;# d. y# i/ `4 F. z1 N
%进行FFT变换并做频谱图
% a8 i+ f, h7 s3 {& I/ @y=fft(bch0,N);%进行fft变换4 |+ c1 a: u! \) x2 F, a
mag=abs(y);%求幅值
9 i0 e" y1 z: P4 p. Q0 k7 ff=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换
) @; @: h9 P! mfigure(1);) l3 K0 Y5 l# o# }$ h8 r
subplot(212);( i! e( J# U% }
plot(f,mag);%做频谱图0 d& p% o6 W- q- P; H3 f0 C
axis([0,100,0,4]);( V* f6 E' o* p& m
xlabel('频率(Hz)');
9 ~4 ~/ ^! n! b* uylabel('幅值');
5 A% B, j% @/ l$ o' jtitle('bch0幅频谱图N=128');
+ @% y8 D+ U9 q$ U* R0 X. }! Kgrid;5 V, w, L# @. I- l2 s! i0 |3 h
( V D, ^. Z. P不好意思我写在写字板粘贴过来的不知道为什么标点没了 现在改过来~ |