双层规划问题求解小试

文章中的优化模型：

分别单独求解

clc,clear

f1 = [-1,-2];
A1 = [-2 3
    1 1];
b1 = [12; 14];
[x1 y1] = linprog(f1,A1,b1);

f2 = [0 -1];
A2 = [-3 1;3 1];
b2 = [-3; 30];
[x2 y2] = linprog(f2,A2,b2);

找到最优解：

clc,clear
x0 = 3;
y0 = 3;
%使用迭代解双层规划问题
x_old = x0;
y_old = y0;
x = upper_func(y0);
y = lower_func(x0);
flag=0;
while abs(x-x_old)>1e-3&&abs(y-y_old)>1e-3
    x = upper_func(y);
    y = lower_func(x);
    flag = flag+1;
    x_old = x;
    y_old = y;
end
disp(strcat('x1:',num2str(x)))
disp(strcat('y1:',num2str(y)))
disp(strcat('-x1-2*y1:',num2str(-x-2*y)))

其中两个函数upper_func.m和lower_func.m

function x1 = upper_func(y)
% y in single number
f = -1;
A = [-2;1];
b = [12-3*y;14-y];
x1 = linprog(f,A,b);

function y1 = lower_func(x)

f =-1;
A = [1;1];
b = [-3+3*x;30-3*x];
y1 = linprog(f,A,b);

输出

x1:8
y1:6
-x1-2*y1:-20

实测x0, y0取别的初值的时候不可以。但当

x0 = 3.666666;
y0 = x0;

的时候，迭代收敛，输出

x1:6
y1:12
-x1-2*y1:-30

现在面临两个问题：

1. 怎么在算数的层面上证明可以取到的这个值$-30$为全局的最小值？
   带入前述的第一个数，可知这个数已经不在可行域内。
2. 怎么判断迭代初值收敛的条件，并求解初值收敛的范围？