_{程序Untitled52.m的DGP1，GLS去势p=2，a=-13时。}

_{程序Untitled61.m的DGP2，GLS去势p=1，a=-7时。}

_{程序Untitled62.m的DGP2，GLS去势p=2，a=-13时。}

_二、数据生成系统存在线性时间趋势项时，单位根检验犯第I类错误的概率

3、DGP3：y_t = α+ δt + y_t-1+ ε_t

4、DGP4： y_t = α +δt + u_t，u_t = u_t-1+ ε_t

α=0.5，δ=0.5

_{这时也共4个程序， 程序Untitled71.m的DGP3，GLS去势p=1，a=-7时，模拟的样本可以自己选定，根据不同的样本计算不同的概率，这时的概率计算的临界值由我给的上一个下载里的https://bbs.pinggu.org/thread-198766-1-1.html&page=3。}

matlab的DF-GLS的分位数为5%时的单侧临界值模拟程序