本人有一组数据,会用程序求出其中一行的时间差。程序如下:
v=rate
a=v(500,:)%第500行
k1 = find(a>=0.001,1) %找出第一个>=0.001的数
k2=find(a>=0.001,1,'last') %找出最后一个>=0.001的数
t1=(k1-1)*0.05 %把点转化成时间用(i-1)* dt
t2=(k2-1)*0.05
tend=t2-t1 %结束的时间-开始的时间为延时
但是这组数据有200800行,我想算出每一行的tend,于是写了个for循环,按理说应该得到一个200800x1的数据,现在得到的是200800x0的空矩阵,想问下如何修改?
或者有更好地方法来算每一行的时间差(tend)吗?困惑了一下午,希望有人可以指导一下,非常感谢!!
for循环语句如下:
for i=1:200800 %有200800行
a=v(i,:)
k1 = find(a>=0.001,1) %找出第一个>=0.001的数
k2=find(a>=0.001,1,'last') %找出最后一个>=0.001的数
t1=(k1-1)*0.05 %把点转化成时间用(i-1)* dt
t2=(k2-1)*0.05
y(i,:)=t2-t1;
end