是显示上的细化.
另外经过FFT后1024离散点变成的是一个什么样的矩阵呢?
下面是程序
%==========================================================%
%额定工作转速1486转/分,。齿轮箱为两级斜齿轮减速机构,第一级主动齿轮直接安装在电机输出轴上,齿数为24,被动齿轮齿数68
%D:\mylab\柳亦兵实验数据\gear\LOCAL\TIME\LF:齿轮故障实验,
%LF,低频采样FS=1600hz;
%中频MF,中频采样FS=6400hz;
%高频hF,高频采样FS=25600hz;
%HL高负荷;30.7Nm
%中负荷ML;17.9Nm
%低负荷LL;5.1Nm
%===========================================================
clear;
%%Z1,点蚀面积较小,.000无故障;.024,共24齿故障
[fid]=fopen('e:\work\TIME\MF\KHLMF3.002');%
t11=fread(fid,125952,'float32');fclose(fid);
clf;%擦除当前图形
fs=1600;%低频采样FS=1600hz;
subplot(2,1,1);
plot(t11);
N=length(t11);%用于计算矩阵长度N=1024
nf=0:N/2-1;%nf=511
xf=nf*fs/N;%xf=511*1600/1024=798 (频率轴数据)
y=fft(t11);
yf=abs(y)/N*2;%abs:对绝对值或复数求模
subplot(2,1,2);
plot(xf,yf(1:512));%画出频谱图
xlabel('Hz'); |