若奇函数f(x)满足f(x+2)=-f(x),x∈[0,2]时,f(x)=3^x+1,求f(4001)的值