求一个matlab代码， 公式不懂也没关系。其实蛮简单，就是一组时间序列数据，可看成向量n行1列，套下面DFC的公式，注意我要算的是taot=1,2,3.....n， 也就是第1行到n行的求和是一个数，下一次是第2行到到n行一个数， 总共有n个结果，求其图像，求其中最大的结果。

求DFC方程1的Matlab代码， Yt可看成一组向量数据，t可以看成是一组向量的顺序. taot 可以看成 1，2，3......

也求方程2(方差)的matlab代码， t>taoT 就是1， 反之0.  方程2中减号后面的变量是ols的方程(21)估计量，这个我可以求出。
若有大神告诉我如何蒙特卡罗求critical value 在这个布朗运动的分布组合， 我也十分感激。