基于优化的直线拟合
7视图(30天)
显示旧的评论
你好!
我用以下优化脚本我拟合后的曲线
两点之间
和
我的问题是,我怎么能解决同样的通过添加约束优化问题
吗?
clc;
清晰;
抽搐
%的数据:时间和响应数据
10 t = [0]”;
7 y = [1];
%看数据
次要情节(2,1,1)
情节(t y‘*’,“MarkerSize”,10)
网格在
包含(“时间”)
ylabel (“响应”)
持有在
%曲线符合
E =[(大小(t)) exp (- t))
%解决约束线性最小二乘问题
% cNew = lsqlin (E, y, [], [], [1], y(1),[],[],[],选择)% Solver-based方法
p = optimproblem;
c = optimvar (“c”2);
p。ObjectiveSense =“最小化”;
p。目标=总和((E * bxcy)。^ 2);
%约束的例子:p.Constraints.intercept c = c (1) + (2) = = 0.82
索尔=解决(p);
cNew = sol.c;
tf = (0:0.1:10) ';
Ef =[(大小(tf) exp (- tf)];
yhatc = Ef * cNew;
%绘制曲线\
次要情节(2,1,2)
情节(t y‘*’,“MarkerSize”,10)
网格在
包含(“时间”)
ylabel (“响应”)
持有在
yhatc情节(tf)
标题(y (t) = c₁+ c_2e ^ {- t}’)
toc