目前,我正在尝试找出如何使用模型来解决非线性系统。模型是:
m2(x) = a*x + 2*exp(b*x)
我必须使用数据解决a和b的问题:
x = [0.00,0.25,0.50,0.75,1.00];
y = [2.10,3.70,6.26,10.03,16.31];
使用它,我解决了下面的成本函数:
(a,b) = sum(y(i) - m2(x(i))^2
我已经尝试了一段时间,但我想不出如何从m2(x(i))中提取a和b。fsolve对我不起作用,我也尝试过常规的solve功能。我不是在寻找答案,而是在寻找一个方向,因为我遇到了一个完全的障碍。
编辑:该模型捕获了y对自变量x的依赖关系。我运行此函数来获得5个方程。
function [f] = func(x,y)
syms a b
f = a*x + 2 * exp(b*x) == y;
end
使用驱动程序:
for i = 1:5
f(i) = func(x(i),y(i));
end
disp(f');
发布于 2021-03-22 20:07:08
您的未知数是a,b
,所以我在模型中将它们参数化为w
,而不是您所做的。我使用了lsqnonlin
,解决方法如下:
x = [0.00,0.25,0.50,0.75,1.00];
y = [2.10,3.70,6.26,10.03,16.31];
m2 = @(w) w(1)*x + 2*exp(w(2)*x);
cost_func = @(w) m2(w)-y;
v = lsqnonlin(cost_func,[1 1]); % normal least squares
plot(x,y,'x',x,m2(v));
https://stackoverflow.com/questions/66752278
复制