function stabledesk2
clc
clear
x0=rand(2,1);  %x,y
ub1=[1-eps;1-eps];
lb1=[eps;eps];
a=fmincon(@fun1,x0,[],[],[],[],lb1,ub1)
%a=x0;
fun1(a)
end

function ra=fun1(x)

[A1,C1]=solve(fun2([1;x(1)-1],[x(2)-1;1]),'A','C');
ro1=dis(A1,C1,0,0);
rp1=dis(A1,C1,1,1);

if x(1)<=1e-6
    ro2=1;
    rp2=0;
else
[A2,C2]=solve(fun2([1;x(1)-1],[-1;-1]),'A','C');
ro2=dis(A2,C2,0,0);
rp2=dis(A2,C2,1,-1);
end

if x(2)<=1e-6
    ro3=1;
    rp3=0;
else
    [A3,C3]=solve(fun2([-1;-1],[x(2)-1;1]),'A','C');
    ro3=dis(A3,C3,0,0);
    rp3=dis(A2,C2,-1,1);
end

ra=double(max(max(double(rp1/ro1),double(rp2/ro2)),double(rp3/ro3)));
end


function lineq=fun2(x,y)
syms A C
[x,y];
lineq=[A*x(1)+x(2)+C;A*y(1)+y(2)+C];

end

function r=dis(A,C,x0,y0)
r=abs(A*x0+y0+C)/sqrt(A^2+1);
end

2楼: Originally posted by cheersend at 2012-06-28 12:43:28
难道不是应该先给出“稳定”的定义的吗？