各位stata or matlab的高手，求助阿！
有一个模型其实是很简单的，但是我不知道怎么编程，麻烦各位帮帮我啊！

P(t)=f(t)+a*P(t-1)-a*f(t-1)

f(t) 很简单，多项式：f(t)=b0+b1*t+c*sin(2pi/12)t +d*cos(2pi/12)t

所以其实就是：

P(t)=b0+b1*t+c*sin(2pi/12)t +d*cos(2pi/12)t +a*P(t-1)-a*[b0+b1*(t-1)+c*sin(2pi/12)(t-1) +d*cos(2pi/12)(t-1)]

P(t) 和t 是数据，做参数估计：a，b0, b1, c, d 和常数

麻烦各位帮帮忙，stata或者matlab ！谢谢谢谢