如下:
a1 = [1 2 3 4 6 8];
b1 = [10 9 8 7 6 6];
a2 = [1 2 4 6];
b2 = [5 4 3 2];
a1是包含a2的,而且b1是由a1得到的,b2是由a2得到的。我的需求是:求出a1和a2相等的部分,然后得到各自相等的那部分在两个向量中的下标,这样我就可以得到相等那部分在b1中是哪些数,而这些数重新组成的向量跟b2是等长的。
这个需求只需要用intersect这个函数就可以了。
intersect是求出两个向量之间相等的部分,还可以得到在对应向量中相等部分的下标。
其用法如下:
[C,IA,IB] = intersect(a,b);intersect其实顾名思义就是交集的意思,C就是两个向量之间的交集,而且会被自动升序排序,而IA是相等的部分在a中的下标,IB是相等的部分在b中的下标。
注意:对于a中重复的值,intersect只会取第一个值,如果a和b里都有重复的值,则都取重复值的第一个。
我的应用:
d1是已知的15分钟下的日期,d2是已知的30分钟下的日期。M1是已知的15分钟下的收盘价计算得到的MACD,M2是已知的30分钟下的收盘价计算得到的MACD。为了在30分钟下得到15分钟的MACD。因为必须保证在30分钟下15分钟的MACD是和30分钟的MACD等长的。之前想到的办法就是用for循环判断15分钟的日期值是否等于30分钟的日期值,等于的时候则取出来放到一个新的数组里。但知道了intersect之后,几行代码搞定。如下:
[C,IA,IB] = intersect(d1,d2);
newM = M1(IA); newM就是和30分钟等长的而且同步的15分钟的MACD。