功能:有理式连乘拟合

格式:
[a, er, f] = RationalMulFit(x, y, n)
[a, er, f] = RationalMulFit(x, y)

x : 矩阵变量存储的拟合自变量,每一行对应一个样本
y : 矩阵变量存储的拟合因变量,每一行对应一个样本, 这里y的行数必须与x一致, 且y的列数必须为1
n : 有理式个数, 此个数越多, 拟合参数越多, 一般也会更加精确, 此参数默认为3

a : 返回的拟合参数
er: 返回的误差平方和
f : 返回的拟合表达式

说明:
1、本函数主要拟合如下表达式
$$
\begin{cases}
   g_i &= \dfrac{a_j + a_{j+1}x_1 + a_{j+2}x_2 + ... a_{j+m}x_m}{1 + a_{j+m+1}x_1 + a_{j+m+2}x_2 + ...}
      &,j = (i - 1)(m + m + 1) + 1
\\
   y &= g_1 \times g_2 \times ... \times g_n
\end{cases}
$$
   
2、本函数内部使用梯度优化窗口功能进行优化.

例子:
//下面演示使用本函数去拟合1到8区间的正弦函数
x = 1:0.1:8;
x = x';
y = Sin(x);
[a,er,f] = RationalMulFit(x,y,4)//回车得到如下结果, f当中的x(,1)表示x变量的第1列数据
a =
[-1.70656254966939
  0.55038620890590
 -0.01070430402376
 -2.38281449335200
  0.27849860983223
 -0.00741210426056
  0.31126826581556
 -0.56487767410679
 -0.00738471295414
 -1.39193402124384
  0.22010031991974
 -0.01424793665337 ]
er =
[ 0.08621304576281 ]
f = "(a(1) + a(2)*x(,1))/(1 + a(3)*x(,1)) * (a(4) + a(5)*x(,1))/(1 + a(6)*x(,1)) * (a(7) + a(8)*x(,1))/(1 + a(9)*x(,1)) * (a(10) + a(11)*x(,1))/(1 + a(12)*x(,1))"
