|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
阶次分析:我的理解是主要处理变转速信号的。在频域内的变速信号,通过等角度采样之后变成“定速”信号。理论部分大致弄清楚了,只是代码部分不太明白
然后论坛上最热的关于阶次分析(有仿真例子)的帖子是这个http://forum.vibunion.com/forum. ... &fromuid=229660。代码如下:
clear;
L = 5000; % Length of signal
Fs = 1000; % Sampling frequency
T = 1/Fs; % Sample time
t = (0:L-1)*T;
x = 0.7.*sin(2.*pi.*t.*t);
y=resample(x,0.005); % 0.005 是重采样转数增量
L=length(y);
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(y,NFFT)/L;
Fs=1/0.005;
f = Fs/2*linspace(0,1,NFFT/2);
plot(f,2*abs(Y(1:NFFT/2)))
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')
resample略:
先说我的第一个问题吧。我不明白这个代码里面“0.005 是重采样转数增量”。。。这个0.005是怎么算出来的。然后之后算重采样之后的采样频率时“Fs=1/0.005;”这里的0.005是上面那个转速增量?这个是采样频率怎么来的?(后面还有一个比较长的,需要自己输入数据的,暂时没研究)
然后是impulse主任的这个帖子http://forum.vibunion.com/thread-123498-1-1.html。代码如下
clear
Ls=1024*64*4; %样本数
fs=4096*8; %采样频率
t0 =0:1/fs:(Ls-1)/fs;
t1=t0(length(t0));
f0=5; %起始频率
f1=30; %终止频率
x_1x = 1*chirp(t0,f0,t1,f1,[],280); %一倍频信号
x_5x = 0.65*chirp(t0,0.5*f0,t1,0.5*f1,[],115); %0.5倍频信号
x_2x = 0.75*chirp(t0,2*f0,t1,2*f1,[],20); %2倍频信号
x_noise=0.15*rand(1,length(x_1x)); %噪声
x=x_1x+x_2x+x_5x;%+x_noise; %合成信号
[maxtab, mintab] = peakdet(x_1x, 0.5); %通过基频信号(相当于转轴信号)检波,模拟后面的键相信号,函数附在附件中
t=maxtab(:,1)';
tacho=zeros(1,Ls);
tacho(t)=1; %模拟键相信号,假定1为键相脉冲位置
tacho(t+1)=1;
subplot(211)
plot(t0,x);
title('振动信号')
subplot(212)
plot(t0,tacho)
title('转速(键相)信号')
主任的例子只有“模拟产生变转速振动信号和键相信号”部分,里面用了这个函数“[maxtab, mintab] = peakdet(x_1x, 0.5);”。看了一下这个是算波峰和波谷的,然后这个里面是用它的波峰maxtab做键相信号。我在matlab里面运行了一下。。大概是140个点。想问一下这个之后该怎么处理才能得到帖子上面的图片。然后这个重采样之后采样频率一般是怎么计算的?
刚开始研究阶次分析,新人小白一个,如果问题有点幼稚,望各位朋友不吝赐教!多多留言!在此不胜感激!
|
|