我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图
# w T8 }9 K- [# W: y0 V% a要得出的图像跟故障应该有的图像相似
! _7 `. t1 r& K8 U) z9 C并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征8 W, f7 O0 D; _6 R( L4 W
我编的程序如下 你看看用得着不
1 h' J1 S# h: a P I) k. i* d- j, Z/ `3 n+ o S
%***************1.bch0****************%0 f# {/ ?' m) o$ S; Z
clc
7 M6 g- D' K. rclf: D5 m7 L9 Y# f; h: Q
clear
' M ~: k4 r6 i4 \( g5 e4 s$ R- i8 |bch0=load('F:\MATLAB7\lijiande\bch0.dat');9 D; B% Y7 v: D2 [! s& D
fs=100;%设定采样频率
9 l0 Z1 r/ s, r: TN=128;& V7 x* p( _! w
n=0:N-1;1 H; v9 {# q" i
t=n/fs;- t5 g2 b3 m- l; Y' r2 v5 E; ~5 k
bch0=bch0(n+1);
: E4 K% b& k+ E& t0 |' h* L' q%采样
4 \0 e( z1 B8 z! ^figure(1);2 s h( y( ?6 \* L3 y! t- U3 ]. I
subplot(211);
3 f# Q& @; `( \' ^2 J2 [plot(t,bch0);%作正弦信号的时域波形
- _& d$ ^1 o% E' Laxis([0,N/fs,-0.4,0.4]);
; o1 ?" h0 ]* T5 O3 axlabel('t');
5 H, d0 H4 K5 yylabel('y');7 U, C# O- {) F {
title('bch0采样t=[0,1.28]');; r% `# W$ \* [: E/ B2 Y5 M
grid;/ Q, v, Y! L3 G9 O5 x
%进行FFT变换并做频谱图
* b7 `- ]5 c) q% s- i7 O, B" ny=fft(bch0,N);%进行fft变换7 ? @( x2 K9 Z
mag=abs(y);%求幅值
X5 y* w5 V# r7 w% ^: cf=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换& I: G" T+ G3 J/ J
figure(1);
4 ^; q% [" K8 W" F# P6 K1 @% v: msubplot(212);9 p. A- ?4 [6 Q- h+ E
plot(f,mag);%做频谱图
3 e+ k T! u6 r) l! \ I6 a* E* n8 Paxis([0,100,0,4]); ]0 b1 s( J( O6 O% F" {1 X
xlabel('频率(Hz)');
) E+ i0 k; g0 \6 Vylabel('幅值');% ]6 e7 T% c$ i
title('bch0幅频谱图N=128');
8 E* o5 U) i+ M2 H* S wgrid;$ h \- ?- n2 C `6 m; S1 T6 @. Y# ]
" @4 V! |" ~" b4 h5 {" I1 Y
不好意思我写在写字板粘贴过来的不知道为什么标点没了 现在改过来~ |