最小二乘约束在绝对值
23日视图(30天)
显示旧的评论
答案(2)
Torsten
2023年6月16日14:32
编辑:Torsten
2023年6月16日20:04
rng (“默认”)
n = 32;
y =兰德(n - 1) + 1我*兰德(n, 1);
我=兰德(n) + 1 *兰德(n);
x0 =兰德(n, 1) + 1我*兰德(n, 1);
x0 =[真实(x0);图像放大(x0)];
x0 = x0. / (√x0 1: n)。^ 2 + x0 (n + 1:2 * n) ^ 2);√6 (x0 (1: n)。^ 2 + x0 (n + 1:2 * n) ^ 2)];
有趣= @ (x) (B * (x (1: n) + 1 * x (n + 1:2 * n)) - y)”* (B * (x (1: n) + 1 * x (n + 1:2 * n)) - y);
乐趣(x0)
nonlcon = @ (x)协议([],x (1: n)。x ^ 2 + (n + 1:2 * n) ^ 2 (n - 1));
索尔= fmincon(有趣,x0, []、[] [], [], [], [], nonlcon, optimset (“MaxFunEvals”,10000,“TolFun”1 e-12“TolX”1 e-12))
乐趣(sol)
索尔(1:n) ^ 2 +索尔(n + 1:2 * n) ^ 2 (n, 1)