给定联合自回归具体的模型以后做数值模拟产生出具体的数值，利用所给的数字计算格兰杰因果值，并且在频域上分析画出频域格兰杰因果值的图。能不能给出一个模拟的例子和相关的matlab代码