我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图, F7 B2 Y' F$ y9 U
要得出的图像跟故障应该有的图像相似
0 @0 e1 U$ R8 [- L# w% G并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征
: Y" x1 d8 U7 F: d g6 B: t0 k我编的程序如下 你看看用得着不( v% ] E1 a& x
$ ~4 A! x3 P6 B$ y6 V* k- S%***************1.bch0****************%. \1 Z: Y ]7 {. ]0 q5 }4 S6 u
clc) r+ J3 T! [ n
clf
. ~ u, x- q& }clear; P' f* q$ y* K6 F: @' @# M$ y$ `
bch0=load('F:\MATLAB7\lijiande\bch0.dat');6 b- h' u6 D- J- T9 n$ C
fs=100;%设定采样频率
; F" l( f4 y+ X6 o. i6 hN=128;& h \! l1 o8 {" ]5 S. K; H2 l" w4 _
n=0:N-1;
3 P+ g) q/ s% \, J, Ct=n/fs;
- s8 P- i& P9 @; h' V( Hbch0=bch0(n+1);' K* f, T5 y$ O+ V) M
%采样' k! X* e ], R" y; V2 e% P- ?! k
figure(1);8 e- [1 G, P! V2 V8 b a
subplot(211);
! c/ u4 o! r& g% I2 n1 ^plot(t,bch0);%作正弦信号的时域波形
: Z4 P. R# g! A9 T4 c; \6 waxis([0,N/fs,-0.4,0.4]);; N t; C0 K7 J" O: f3 N1 R3 l8 B" X- \( |
xlabel('t');
! o1 b& d+ E) n' b8 \ylabel('y');0 m2 T; \% s) _
title('bch0采样t=[0,1.28]'); }0 Z c6 p9 x6 y$ Q
grid;4 A& p# d0 Z& Z
%进行FFT变换并做频谱图 `4 S9 G5 r4 ^ d+ V3 ]
y=fft(bch0,N);%进行fft变换
4 P' X) o0 X( B* [mag=abs(y);%求幅值
$ z h. ?% R7 l0 ef=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换
/ c6 ~. `( C! ^8 L$ E; nfigure(1);
2 i0 U/ @4 [) i5 |# Esubplot(212);
X/ B8 Z. o$ n. k5 jplot(f,mag);%做频谱图
0 k; C! K/ A7 Kaxis([0,100,0,4]);; E7 }# q# p: a3 m
xlabel('频率(Hz)');
# `9 p5 {# x* B7 lylabel('幅值');/ [# \% l% J' J) t5 b/ V6 L
title('bch0幅频谱图N=128');3 p1 [0 c# A2 U$ ~+ q% U
grid;
* m4 j0 H d6 Y# m( \, u5 k, H2 b
5 g" h8 X( H. P9 N& x! @; m$ j不好意思我写在写字板粘贴过来的不知道为什么标点没了 现在改过来~ |