我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图& q( V7 \3 g+ R: e3 x8 F
要得出的图像跟故障应该有的图像相似. ^5 T+ Z2 o" a( V; g
并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征5 J& l6 s. `6 l/ A. [1 G
我编的程序如下 你看看用得着不. U7 E: V: i3 ~" `
8 e3 i, l v2 j6 `! v/ y0 K
%***************1.bch0****************%9 Y+ j/ A* O" W4 k; j7 Q& J9 p
clc0 D% \3 E6 K, L% u% `
clf4 I; s8 \, i9 G! {9 {+ j( D4 Y& p7 \
clear
% }) h0 F; B. y+ P$ tbch0=load('F:\MATLAB7\lijiande\bch0.dat');# W! c, t3 J3 ~8 h; G- M+ W
fs=100;%设定采样频率5 r0 b [# Z& }2 X/ v' L
N=128; ~8 y* l! S. P8 h' ?
n=0:N-1;
, V' Z+ V4 T7 e! O4 S" Vt=n/fs;
( p( H: B$ N, l% f" s+ sbch0=bch0(n+1);
% b, S+ n' ~- v) u7 u%采样' E; u8 x. V- H
figure(1);$ Z* x" d4 E- u1 \" y7 @8 i8 ?" N
subplot(211);
. r! }/ |1 w* n! Vplot(t,bch0);%作正弦信号的时域波形9 j5 W1 q" L1 l; E
axis([0,N/fs,-0.4,0.4]);
& @4 |- B- ?0 Pxlabel('t');9 Q5 \+ j, F' R! J8 c" o! x( V
ylabel('y');. `' g H6 A0 T2 R
title('bch0采样t=[0,1.28]');1 _( q* @0 e0 V" [. ?1 g6 I9 s
grid;1 Z8 ]$ Q% f# i1 [. d# O
%进行FFT变换并做频谱图
- s8 N% D# b* My=fft(bch0,N);%进行fft变换
9 c7 U1 z( x% X0 U, N2 S& Smag=abs(y);%求幅值
, F6 p1 G" D$ f7 k3 @' ]f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换, f0 s$ M& T- A3 U- H% \
figure(1);
6 l* V1 i' t+ `, z: }: f0 U# Bsubplot(212);
5 N% g& [$ K0 S7 m! p9 u& |plot(f,mag);%做频谱图
7 I& V) B' E" _, J! h3 o4 Yaxis([0,100,0,4]); ?* l$ ]7 S, }3 q v; v k, ~
xlabel('频率(Hz)');
! b8 y+ w n5 G; T* [ylabel('幅值');
3 }: v) y9 N0 a0 _2 \# ttitle('bch0幅频谱图N=128');
( ~7 l' S. t w9 N4 jgrid;
6 Q8 a8 _& K6 r+ p+ A2 y6 ]( o$ N& a! d
& ^! |: \! |6 h q1 t不好意思我写在写字板粘贴过来的不知道为什么标点没了 现在改过来~ |