我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图, f- n8 H8 a2 V
要得出的图像跟故障应该有的图像相似
2 P: P9 u" V. Z并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征
! N$ L+ ?, X# x3 Z我编的程序如下 你看看用得着不
# ~; ]/ {( P5 x
( B, C5 U- j8 S" c3 u%***************1.bch0****************%
9 E9 {+ e( Y# s* k0 eclc
3 k$ ]) e$ {2 U8 Q" `clf
% G& K. D5 M7 [: Y- Oclear" [4 Y! y6 v6 |/ }# I+ u# |
bch0=load('F:\MATLAB7\lijiande\bch0.dat');
: P5 O# @! Q5 B3 o( Yfs=100;%设定采样频率, {. _- n% P$ S7 J* Z9 r4 j3 v
N=128;; Z' j# {% x; h0 ^+ H
n=0:N-1;
2 A" K! e8 U1 I( w7 lt=n/fs;. b0 R5 P; P- P8 M; y
bch0=bch0(n+1);# ^( P1 I# j9 [9 g" b! N
%采样
' R5 w1 g" n' k vfigure(1);$ V- ?* q7 _" X, b( ^$ g# ?" ]
subplot(211);7 O9 A/ e+ Z! l# j
plot(t,bch0);%作正弦信号的时域波形- {; ]) S+ ^" r
axis([0,N/fs,-0.4,0.4]);
/ x( A2 p" q( Z0 x( W. ^xlabel('t');( c+ @/ d1 @; s
ylabel('y');
, s; m+ E2 n% O& g2 Htitle('bch0采样t=[0,1.28]');3 K7 W/ {* u. f
grid;
i2 v0 x" o7 B) m" [1 z7 K%进行FFT变换并做频谱图* w. n0 }) M+ k7 Q2 b# g
y=fft(bch0,N);%进行fft变换
) i7 Q3 T% M6 e% qmag=abs(y);%求幅值4 M4 m- a+ |1 k! { @4 T; F
f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换
4 ` R! A0 f, i g8 H: ]5 Wfigure(1);: {, n% L. U: }0 e
subplot(212);
! b' R# N7 F$ fplot(f,mag);%做频谱图& P9 C1 O6 T# [; ^& o4 h
axis([0,100,0,4]);' _$ z7 P1 h* f+ o! U. L- d
xlabel('频率(Hz)');
0 x& ^2 d) v* F7 Xylabel('幅值');# `/ |' z; S5 x0 O
title('bch0幅频谱图N=128');" T/ ^6 \4 R# t+ U7 R
grid;
. S. K$ d9 b" X8 y7 K! {% Q9 c) d6 D( I. l
不好意思我写在写字板粘贴过来的不知道为什么标点没了 现在改过来~ |