【fj102】在https://bbs.pinggu.org/thread-125006-1-1.html给出的答案：
写的方法就是把(1-L)^d展开就可以了
可以用loop or do while
你可以truncate 500 or 1000, or user define
至于simulate
你完全可以用gauss or matlab先simulate一组normal的data
然后进行用fix的d arfima算就可以了
还有这个帖子：https://bbs.pinggu.org/thread-1135195-1-1.html
用R或SPLUS做