clear;clc;
%% 求不规则脉冲之前的时间间隔,进而求出电机转速
fid =textread('E:\毕业论文\实验计划\velocity test\velocity test\v-3.txt');
x=(fid(1:36000))';%导入采集到的脉冲数据,电机转一圈一个脉冲
fs=12000;
N=36000;
t=1/fs:1/fs:N/fs;%采样时间段
[maxv,tk]=findpeaks(x,'minpeakheight',64); %maxv峰峰值点 1000最小高度
%% 找出峰值maxv对应的 时间点序列tkt
tkt=[];
for i=1:length(tk)
tkt1=t(tk(i));
tkt=[tkt;tkt1];
end
%% 求出相邻峰值之间的时间间隔,然后倒数即使瞬时速度
vt=[];%时间间隔
for i=1:length(tk)-1
vt1=tkt(i+1)-tkt(i);
vt=[vt;vt1];
end
v=1./vt;%瞬时速度
figure
subplot(2,1,1);
plot(t,x);xlabel('Time(sec)');ylabel('Amplitude');title('Original signal')
% plot(t,envelope);xlabel('Time(sec)');ylabel('Amplitude'