我有實驗數(shù)據(jù)是dat格式的 要用matlab分析出頻譜圖跟時域圖6 U( D: e" g) z: L* [
要得出的圖像跟故障應該有的圖像相似
1 {/ l- T9 o+ H并且?guī)臀医忉屢幌?轉子不平衡的時域圖 頻譜圖各是什么特征, ?2 B% Q- n, @# s5 i
我編的程序如下 你看看用得著不/ i# U- @8 k' q7 W
+ ~: F, z: a! \, V: n! m6 @1 g
%***************1.bch0****************% ?+ z9 a8 k3 Z& L8 }
clc* s4 T# x' B" ~4 n7 T, h% A. r. |
clf* w# F' e y- g+ l' _) B
clear$ b: o9 D) E2 X( d& C) B
bch0=load('F:\MATLAB7\lijiande\bch0.dat');! N( m6 {! |- i+ D( p2 z
fs=100;%設定采樣頻率
0 @; K* a, b1 o. n$ ^N=128;- h! _9 y: l( k" ]
n=0:N-1;
" `5 K1 k6 J4 B2 p5 Tt=n/fs;
. {' O. Z- ]3 ]- lbch0=bch0(n+1);" d4 G5 `* K( x/ @8 o& e9 P4 f
%采樣; {$ _ M, I, n7 t5 T d6 m( R
figure(1);0 `( ^ d, P8 v6 F
subplot(211);! S; p* r I! C& B
plot(t,bch0);%作正弦信號的時域波形
; m* y: T" V" R2 k& X/ Z9 G2 o" C/ @axis([0,N/fs,-0.4,0.4]);
* `1 i( S* X) T N$ exlabel('t');- i; p' a7 Z7 h: @4 g
ylabel('y');* d' K4 c3 ]9 s
title('bch0采樣t=[0,1.28]');6 o2 F7 {# l/ V2 n5 `
grid;
0 F5 S; D" j( A3 ^. e) l%進行FFT變換并做頻譜圖$ z3 u3 C8 M, k. U, E& d
y=fft(bch0,N);%進行fft變換& z- g, e" X" A! |% T4 y d
mag=abs(y);%求幅值4 N4 J. S# E1 v- ?3 |$ z1 l
f=(0:length(y)-1)'*fs/length(y);%進行對應的頻率轉換9 z- B* ^0 e, x! e
figure(1);4 H5 b9 B7 }; s2 b0 g
subplot(212);. F! b: b4 ?1 N! `
plot(f,mag);%做頻譜圖# M8 c' z6 r5 j7 j, O3 l: ` V
axis([0,100,0,4]);
0 d. |7 [+ O( J" H" Xxlabel('頻率(Hz)');* ?9 J" n ~) t) t/ F
ylabel('幅值');# C0 ~2 f5 o% {# ]
title('bch0幅頻譜圖N=128');' G7 t1 ~7 b) g% L
grid;
# d, ~ c) w0 R2 J
; z4 ~' R. t2 u不好意思我寫在寫字板粘貼過來的不知道為什么標點沒了 現(xiàn)在改過來~ |