我有實驗數據是dat格式的 要用matlab分析出頻譜圖跟時域圖! e3 X' }/ B# b: i; [
要得出的圖像跟故障應該有的圖像相似
2 W; Z, ?$ ?: X" e% H& c并且幫我解釋一下 轉子不平衡的時域圖 頻譜圖各是什么特征7 c; ^/ h. J8 G- p7 K
我編的程序如下 你看看用得著不
# g% G8 e; [ Z. c& U# e4 T2 U
: C& X7 h P( f( D0 _4 A9 c6 d%***************1.bch0****************%* D: z/ ^4 Z q! R2 D7 E- ^! o
clc
: T! ~" @) C! h. I& H* u) P/ r# _8 a7 ^clf# D. E( t$ e0 e5 Y( V* i- g, ?
clear
. a+ x- Y9 w: e$ hbch0=load('F:\MATLAB7\lijiande\bch0.dat');$ @9 g2 A% t' R" ]( X4 E6 F9 q
fs=100;%設定采樣頻率
2 J3 S! L/ Y4 \6 Q4 W. l- y* lN=128;$ d/ \( _5 f( j
n=0:N-1;
8 I+ ?( e e* M' N4 pt=n/fs;
5 S) o# z6 {: `bch0=bch0(n+1);
# _: U! R$ ^' Z5 ^6 |9 ]%采樣
3 U; ]3 k# X: r/ Y R: Xfigure(1);; o( b2 e+ f, |# c e! E3 L) D
subplot(211);
$ r! F! U& J5 \8 l" |6 Zplot(t,bch0);%作正弦信號的時域波形
5 u& _% q% H8 T2 daxis([0,N/fs,-0.4,0.4]);- s* A% Z( V; w7 O2 Y
xlabel('t');$ `$ w+ `# {0 Z% u
ylabel('y');( B9 T y- U1 m. m w; D
title('bch0采樣t=[0,1.28]');, M D E" i1 g9 p" E
grid;5 B% U+ {. G& Y
%進行FFT變換并做頻譜圖
; S: y9 {+ t* T+ Q \* H: Py=fft(bch0,N);%進行fft變換 C% h3 {; c3 v/ |$ r t
mag=abs(y);%求幅值& z/ C( _) ^. e R: i J
f=(0:length(y)-1)'*fs/length(y);%進行對應的頻率轉換
% N& ?7 L/ p: N5 K5 |figure(1);* V5 e- v8 g4 E3 H. A3 P0 B
subplot(212);
0 V! }* V5 z1 wplot(f,mag);%做頻譜圖
! f5 u! E% P' U/ Aaxis([0,100,0,4]);
8 y0 _! k; k6 p3 F' L: bxlabel('頻率(Hz)');
' H5 |& y0 c! K, R$ ?) a2 ]ylabel('幅值');
5 X: m& g: q7 ]5 k% i* D# htitle('bch0幅頻譜圖N=128');
3 s9 s9 I, b1 u. v$ Lgrid;$ d/ V- z, W4 m2 r, X+ ~
/ z9 G5 w. W/ X6 J* C5 ~不好意思我寫在寫字板粘貼過來的不知道為什么標點沒了 現在改過來~ |