matlab处理实验数据 实验：利用转子试验台测试转子不平衡以及涡动状态 故障诊断

超级小宝

我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图要得出的图像跟故障应该有的图像相似并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征我编的程序如下 你看看用得着不%***************1.bch0****************%clcclfclearbch0=load('F:\MATLAB7\lijiande\bch0.dat');fs=100;%设定采样频率N=128;n=0:N-1;t=n/fs;bch0=bch0(n+1);%采样figure(1);subplot(211);plot(t,bch0);%作正弦信号的时域波形axis([0,N/fs,-0.4,0.4]);xlabel('t');ylabel('y');title('bch0采样t=[0,1.28]');grid;%进行FFT变换并做频谱图y=fft(bch0,N);%进行fft变换mag=abs(y);%求幅值f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换figure(1);subplot(212);plot(f,mag);%做频谱图axis([0,100,0,4]);xlabel('频率(Hz)');ylabel('幅值');title('bch0幅频谱图N=128');grid;
; k+ c( w& C+ C谢谢大家了~！！！

超级小宝

我有实验数据是dat格式的 要用matlab分析出频谱图跟时域图
$ _: g% s( R% x# _; U' B( W要得出的图像跟故障应该有的图像相似
3 ?( B: N) O7 m/ t' O并且帮我解释一下 转子不平衡的时域图 频谱图各是什么特征
7 \) g: E, m2 T- I0 @1 |+ Q我编的程序如下 你看看用得着不
: K' M  X) I. C; q
%***************1.bch0****************%
clc
9 w5 u5 ]1 U# `7 q, y) P# V( _. gclf
% _5 w. n- m# T' L  F/ z: bclear
bch0=load('F:\MATLAB7\lijiande\bch0.dat');
fs=100;%设定采样频率
! T2 j& c  p" J8 S, iN=128;
) T/ h+ _6 u6 a0 Cn=0:N-1;
& f: `9 b# {* p  f1 m7 f4 ot=n/fs;
bch0=bch0(n+1);
%采样
figure(1);
! a+ \7 I% h$ c, P, asubplot(211);
plot(t,bch0);%作正弦信号的时域波形
# n; |  B' L8 F5 m, j0 Raxis([0,N/fs,-0.4,0.4]);
% c4 e$ c' N/ Y6 l) n% N0 ?. Qxlabel('t');
: \1 S$ J  ^# |0 l* Qylabel('y');
2 ~" X  ]  Y$ d* B4 G) Ltitle('bch0采样t=[0,1.28]');
grid;
( H9 r  s6 U2 V+ R# K; p2 o+ w, T%进行FFT变换并做频谱图
y=fft(bch0,N);%进行fft变换
  p. l2 p$ K" f7 V) N8 `' l( Gmag=abs(y);%求幅值
f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换
figure(1);
subplot(212);
plot(f,mag);%做频谱图
axis([0,100,0,4]);
3 v. r, K  n# {# H; E$ b/ yxlabel('频率(Hz)');
ylabel('幅值');
, U, q, f. ?+ A- Q1 p0 wtitle('bch0幅频谱图N=128');
1 w: J6 k" F& E/ Fgrid;

8 c2 |. G3 x# A% E- j: @不好意思我写在写字板粘贴过来的不知道为什么标点没了 现在改过来~

轴研小草

应该和转子模型以及试验工装模型有关，跨距等。可传一试验工装模型图，以便分析。

黑天鬼

你的转子工作频率多少

黑天鬼

把数据下下来看看再说吧