1

我得到了案例P(t) = a * A(t) - b* B(t),每个 17281x1 加倍。

现在我想使用曲线拟合来获取变量ab.

配件 1

配件箱,这里 a 和 b 为 1 和 1

我知道fitand fittype,但在这种情况下它们似乎不起作用。

任何想法,如何解决这个问题?

4

1 回答 1

1

使用最小二乘法怎么样?如果我理解正确,您的问题可以表示为 P(t) = [A(t), B(t)] * [a; -b]

让和[a; -b] = x_[A(t), B(t)] = YP(t) = P

现在最小二乘解决方案将是:

x = ((Y'*Y)^-1)*Y'*P;

在 Matlab 中,您还可以在这种情况下使用“反斜杠运算符”:

x = Y\P;

为此,您可以在此处找到文档:mldivide

作为参考:

维基百科

数学作品

我希望这有帮助。

编辑:

这是我的测试代码:

A = [1;2;3]
B = [4;5;6]
P = [7;8;9]

Y = [A, -B]

disp('------- regular least squares formula -------')
x = ((Y'*Y)^-1)*Y'*P

a = x(1)
b = x(2)

disp('------- mldivide -------')
x = Y\P

a = x(1)
b = x(2)
于 2017-07-11T07:09:19.430 回答