我得到了案例P(t) = a * A(t) - b* B(t)
,每个 17281x1 加倍。
现在我想使用曲线拟合来获取变量a
和b
.
配件箱,这里 a 和 b 为 1 和 1
我知道fit
and fittype
,但在这种情况下它们似乎不起作用。
任何想法,如何解决这个问题?
我得到了案例P(t) = a * A(t) - b* B(t)
,每个 17281x1 加倍。
现在我想使用曲线拟合来获取变量a
和b
.
配件箱,这里 a 和 b 为 1 和 1
我知道fit
and fittype
,但在这种情况下它们似乎不起作用。
任何想法,如何解决这个问题?
使用最小二乘法怎么样?如果我理解正确,您的问题可以表示为
P(t) = [A(t), B(t)] * [a; -b]
。
让和[a; -b] = x
_[A(t), B(t)] = Y
P(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)