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');title('Original signal')

% hold on;                                                  %绘制原波形

% for i=1:length(tk)

% plot(tkt(i),maxv(i),'r.'); %绘制最大值点

% text(tkt(i),maxv(i),sprintf(' (%g)',tkt(i)));

% end

subplot(2,1,2);

plot(v);xlabel('Time(sec)');ylabel('Amplitude');title('Original signal')

%% 用以上方法求出来的时间间隔不对，峰值点数也不对

%% 想尝试用上升沿，或者下降沿做，但是失败

% hold on;                                                  %绘制原波形

% plot(maxv,'*','color','R');                %绘制最大值点

% for i=1:length(tk)

% plot(tk(i),maxv(i),'r.'); %绘制最大值点

% text(tk(i),maxv(i),sprintf(' (%g)',tk(i)));

% end

% %% 充分利用上升沿，下降沿

% threshold=1.5;

% if X1(1)

% flag=0;

% else

% flag=1;

% end

% old=X1(1);

% new=X1(2);

% count=0;

% for t=2:1:236%236是随便取的，看你要统计多少时间而定

% new=X1(t);

% %上升沿

% if(new>threshold)&&(old<=threshold)&&(flag==0)

% flag=1;

% count=count+1;

% end

% %下降沿

% if(new=threshold)&&(flag==1)

% flag=0;

% end

% %

% old=new;

% end

% count

1544795795(1).jpg

(36.4 KB, 下载次数: 0)

2018-12-14 22:03 上传

2018-12-14 21:52 上传

点击文件名下载附件

375.89 KB, 下载次数: 7