我有Ab = c形式的问题,其中A是未知的,M×M矩阵,b是M×1矢量,c是M×1矢量。我有很多b和c的样本(超定问题),我想找到“最好的”A.任何可以轻松解决这种回归的库的建议?
非线性但多项式模型怎么样?也就是说,形式中的一些等式:
[c_0 c_1 ... c_M] = [1 b_0 b_1 ... b_M ... b_0 ^ k b_1 ^ k .. b_M ^ k] A
其中k是多项式阶数,A现在是M×k矩阵。
答案 0 :(得分:1)
找到了一个很好的解决方案:DLIB很好地解决了这个问题。这里有一个例子: http://dlib.net/least_squares_ex.cpp.html
注意:如果需要使用矩阵,DLIB仍然需要将数据存储为向量。也就是说,如果您正在求解Ab = c且A是3x3矩阵,则上面链接的示例中的parameter_vector必须为typedef matrix<double,9,1> parameter_vector;
,否则将无法编译。