请教：这种情况下极大似然法如何实现

翻身仗

R为n维观测值向量（比如：R = 31, 25, 52, 36, .....）

X Y 均为已知的 n 维向量

T = a*t^2 +b*t+c （t=1,2，.....，n； a、b、c为待求的系数，均为标量）

R = Rg + E (E也是n维向量，是残差项,服从均值为0，标准差为 d 的正态分布；Rg是R的估计值)

Rg= f(X,Y,T) （ f 是函数，形如: (X/Y)^T ，但比这个稍复杂一点点）

已知 R 服从均值为 f(X,Y,T) ，标准差为 d 的正态分布。

问题是：要用极大似然法求 a b c  d 的值。请问，matlab里如何实现？