我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图
: N3 w5 G# h: C$ k4 E要得出的图像跟故障应该有的图像相似
- }1 e4 ?$ s5 M* U并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征# t9 W$ N& D0 @& W: L& p$ Z( u% X
我编的程序如下 你看看用得着不/ I2 d2 f/ t2 G4 f7 F
! o9 S+ [/ E; q0 h%***************1.bch0****************%. @8 c2 s# \: [% \
clc/ m4 ]. R; w- G1 s( F4 F
clf1 x1 E" H; u9 _' [! f5 ?
clear" `+ R3 P* M4 L: y) i: z- V$ z) K
bch0=load('F:\MATLAB7\lijiande\bch0.dat');
+ v+ N; _6 \7 o7 Z6 `0 hfs=100;%设定采样频率2 `' U: S* ]5 S' ]$ i1 T
N=128;
/ v5 p1 _$ a) N1 cn=0:N-1;
( W; k! u8 `# _- c( K2 @ O) Bt=n/fs;
2 Y! r* p5 J7 k: l; o* E( ybch0=bch0(n+1);
! U+ S& \" i- i%采样
# q/ a5 m! k* l6 E; {7 P4 F9 l- nfigure(1);
) O$ m/ ]4 y; u0 ~; `/ W/ ]$ @subplot(211);
' U& @% X: [2 i7 W* Y. `4 x1 I# Hplot(t,bch0);%作正弦信号的时域波形7 r% i! U. A" P! e7 |$ G
axis([0,N/fs,-0.4,0.4]);, i! _4 X: S) ?# `" r
xlabel('t');" P* _2 L) W7 l' n% E- R% S
ylabel('y');. J, G5 _/ H6 ?$ e3 C% Y
title('bch0采样t=[0,1.28]');3 \4 O/ A# q, f8 P0 h
grid;% }6 l9 r! p/ X1 I+ j# O
%进行FFT变换并做频谱图: y1 m& k& b. a; M$ e
y=fft(bch0,N);%进行fft变换8 G6 x- b$ ?! e$ d( u
mag=abs(y);%求幅值
' {0 ~- N" v2 ?1 q5 q- Wf=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换- u s# u3 N5 @9 O* q3 s8 @6 D
figure(1);$ \: r! d& @+ z7 h2 Q5 g' i
subplot(212);$ r1 R3 j# Z% O6 F: V& a4 J
plot(f,mag);%做频谱图1 K$ i3 ?7 H$ D2 M9 i- A
axis([0,100,0,4]);
1 N5 e, @6 v2 B! O- K8 k% hxlabel('频率(Hz)');% f" Q% W; H2 V: ~- F9 j' r
ylabel('幅值');6 e7 a1 y# v9 [" ~
title('bch0幅频谱图N=128');
# N/ C& u. o* c. `* B' ggrid;+ q0 q0 b& }. K' T/ M) u, C: R7 H
8 }/ y3 o9 k: X0 b4 U不好意思我写在写字板粘贴过来的不知道为什么标点没了 现在改过来~ |