极大似然函数估计Matlab子函数

模块化编程便于程序的阅读、理解和二次开发，本文将《系统辨识与自适应控制Matlab仿真》书中递推极大似然参数估计算法改写为Matlab子函数的形式，方便以后调用。

function [theta] = LG_RML(u,y,oa,ob,od)
% Writer: GeoLee
% Date: 2015.9.25
% 参数说明:
% CARMA模型：A[z^(-1)]*y(k) = B[z^(-1)]*u(k) + D[z^(-1)]*e(k)
% A[z^(-1)] = 1 + a1*z^(-1) + ... + a_na*z^(-na);
% B[z^(-1)] = b0 + b1*z^(-1) + ... + b_nb*z^(-nb);
% D[