我有實驗數據是dat格式的 要用matlab分析出頻譜圖跟時域圖
% c. H8 Q$ s" p( E4 [0 t要得出的圖像跟故障應該有的圖像相似
J( d* n! c8 H: P) P并且幫我解釋一下 轉子不平衡的時域圖 頻譜圖各是什么特征: ?, F9 O4 Y4 L: o* A% n" A
我編的程序如下 你看看用得著不
) G9 \) g6 v' o0 M( @
) o8 \/ b! \0 y4 G%***************1.bch0****************%" v9 X0 D) T8 h- Q* {4 D: m
clc
8 H5 \9 h1 j; w* g A8 C$ T- Zclf' e+ O3 o0 H- {
clear5 F+ i. ^( n% O& K9 F
bch0=load('F:\MATLAB7\lijiande\bch0.dat');
) r/ Y6 ^& z: b. R. sfs=100;%設定采樣頻率- @9 L. @6 S( x# Y( q: z
N=128;: Y0 b' @8 ^) k& u! |/ }' z& o6 d6 T
n=0:N-1;0 m# n1 t3 S9 H% j; U
t=n/fs;
# Q/ L% h! A- m( a- T1 U; n7 fbch0=bch0(n+1);2 ?& d- K- t2 V) X6 y
%采樣
) x7 Q% P+ M# Afigure(1);
9 I4 U6 _8 l6 v0 Hsubplot(211);) p$ H6 q8 s; `& p5 R
plot(t,bch0);%作正弦信號的時域波形3 o% S7 z) ?0 u: R+ l" L3 |7 m& _
axis([0,N/fs,-0.4,0.4]);. ]$ C0 M- d% |( h/ a
xlabel('t');+ B5 A O# {: p" V" l
ylabel('y');) n5 H& B! s9 M5 _# Z( I+ k
title('bch0采樣t=[0,1.28]');
4 r8 Y+ Z2 }# U; ngrid;) |5 U5 _8 e& _: x
%進行FFT變換并做頻譜圖 R$ S( o n0 s9 _1 f, s8 n
y=fft(bch0,N);%進行fft變換
( B& m; S/ Z: P( u! w0 G4 t" _mag=abs(y);%求幅值0 u/ X- y% L+ S; [! `+ n9 j
f=(0:length(y)-1)'*fs/length(y);%進行對應的頻率轉換
" V' x: R4 G L' l: T" yfigure(1);
* s; X. U6 s4 W# f2 U/ Ksubplot(212);9 W5 `6 w" `; K# i' ~. M
plot(f,mag);%做頻譜圖
# e4 n* r; w3 P8 U1 U( S7 eaxis([0,100,0,4]);+ B% m3 c* U6 \3 i* O
xlabel('頻率(Hz)');6 \ L f6 d- F3 L: j
ylabel('幅值');8 B0 M; f7 C* `' i6 v# V
title('bch0幅頻譜圖N=128');- P- Q( |4 W+ B. h
grid;0 U6 O2 o' f" h8 ~# O' J: j1 A
4 j1 Z2 Q: z3 Q2 m9 s. q. q( y, D
不好意思我寫在寫字板粘貼過來的不知道為什么標點沒了 現在改過來~ |