GEO_WU: lapack 线性方程组最小二乘求解函数dgels

1. DGELS函数简介

DGELS函数是利用最小二乘求解线性方程组（超定或欠定）AX=B或A'X=B的函数，处于driver level。利用LU分解或者QR分解进行求解。

2. DGELS函数调用格式

DGELS(TRANS, M, N, NRHS, A, LDA, B, LDB, WORK, LWORK, INFO)

输入参数：

TRANS（CHAR）：值为'N'，表示No transpose, 对应AX=B；值为'T',表示Transpose，对应A'X=B;

M,N（INT）: 系数矩阵A的维数，M>N,超定（下面只以超定方程组为例）,M<N，欠定；

NRHS（INT）：观测值向量B的列数（通常为1）；

A（Double）：系数阵；

LDA（INT）：同DGESV函数中的LDA，LDA>=max(1,M)，通常等于M；

B（Double）：观测值向量，当TRANS=‘N’时，B（M，NRHS），当TRANS=‘T’时，B（N，NRHS）；

LDB（INT）：同DGESV函数中的LDB，LDB>=max(1,M)，通常等于M；

LWORK（INT）：Length of Work, 数组WORK的长度，见http://www.netlib.org/lapack/lug/node118.html；

输出参数：

B（Double）：B(1:N)为计算的结果；B（N+1：M）误差信息，具体见函数说明；

A（Double）：M>N, QR分解矩阵信息；M<N, LU分解矩阵信息；

WORK（Double）：维数为LWORK，存储空间(WORKspace)信息，具体见说明http://www.netlib.org/lapack/lug/node120.html

INFO：=0，成功执行；=-i，第i个参数illegal；=i，A的第i个对角元素为0，A非满秩阵，无LS解。

3. DGELS函数举例：

4. 同类型函数

5. 参考资料

GEO_WU