已知变量x，y满足约束条件1≤x+y≤4，-2≤x-y≤2，若目标函数z=ax+y（a＞0）仅在点(3，1)处取得最大值，则实数a的取值范围为（    ）。