|
2#

樓主 |
發(fā)表于 2011-5-26 11:47:56
|
只看該作者
我有實(shí)驗(yàn)數(shù)據(jù)是dat格式的 要用matlab分析出頻譜圖跟時(shí)域圖* B6 V* I, O. I+ L* C E
要得出的圖像跟故障應(yīng)該有的圖像相似
5 \8 q% X- U+ q+ l. ~6 i并且?guī)臀医忉屢幌?轉(zhuǎn)子不平衡的時(shí)域圖 頻譜圖各是什么特征
/ Z" {% S2 W$ l5 i2 ^. N我編的程序如下 你看看用得著不
9 M, w! J, i6 b0 `1 V7 U/ j. h/ o0 X) H1 Q! M* S; K
%***************1.bch0****************%( P) P( h0 ]% m/ N9 J4 a
clc
. d# c* W f9 I7 gclf8 K# e u5 h' N( E- a, p
clear
) d" v; T4 X& |* a+ z' F0 ]bch0=load('F:\MATLAB7\lijiande\bch0.dat');
+ o9 d* F9 @( E2 b) F* Mfs=100;%設(shè)定采樣頻率
K, R/ u) `5 A$ ^: AN=128;
0 V; K; d7 H0 _- |1 i Pn=0:N-1;
2 f+ f; q9 \. R2 }8 ut=n/fs;# B E6 z* w! O4 D) C( J. ~
bch0=bch0(n+1);+ e+ L. M" f: o6 ?2 @( G9 x
%采樣; Z% _" h/ V: _/ {
figure(1);
% K) f/ @" k; L1 w3 bsubplot(211);
4 H) ^) H, U! Z8 b1 _plot(t,bch0);%作正弦信號(hào)的時(shí)域波形" _, Q0 k$ b7 ^
axis([0,N/fs,-0.4,0.4]);/ K' g. p9 a- Z" C' O) ~
xlabel('t');& t) ^, E7 p$ g9 B
ylabel('y');
3 u3 k8 ~$ @) Z8 P& ]title('bch0采樣t=[0,1.28]');
1 S# R" X$ L" {7 H5 jgrid;
. P6 Y# V2 `( |1 R7 y5 E& i%進(jìn)行FFT變換并做頻譜圖: k/ y9 u$ N+ q- z3 E9 W8 g
y=fft(bch0,N);%進(jìn)行fft變換
/ |+ P, C* ^" P! E+ W! a2 D. smag=abs(y);%求幅值, K. a3 @4 b3 ]# L$ u3 O, C* B% ?( P8 M" F
f=(0:length(y)-1)'*fs/length(y);%進(jìn)行對(duì)應(yīng)的頻率轉(zhuǎn)換) |- `# L# J8 i: ]% m- {! q
figure(1);) `; W- T4 S7 l- v6 r
subplot(212);5 `3 K& x4 S/ b. Q0 W
plot(f,mag);%做頻譜圖0 Z6 o6 O- c& v2 h0 M" Z( V" G% q
axis([0,100,0,4]);
9 |/ B4 L$ @0 `- Pxlabel('頻率(Hz)');* y3 u" f$ s. z$ K
ylabel('幅值');& N% S6 A2 i- a
title('bch0幅頻譜圖N=128');
/ s0 e- l9 e0 f5 t8 b" Pgrid;
3 b; b* J2 R* q# l5 K
. ?* s6 t- g$ y1 T不好意思我寫在寫字板粘貼過來的不知道為什么標(biāo)點(diǎn)沒了 現(xiàn)在改過來~ |
|