?? computeconstraintresiduals.m
字號:
function res = computeconstraintresiduals(p,x)
res= [];
if ~isempty(p.F_struc)
vecres = p.F_struc*[1;x];
if p.K.f>0
res = -abs(vecres(1:p.K.f));
end
if p.K.l>0
res = [res;vecres(p.K.f+1:p.K.f+p.K.l)];
end
if p.K.q(1)>0
top = 1+p.K.f+p.K.l;
for i = 1:length(p.K.q)
n = p.K.q(i);
X = vecres(top:top+n-1);top = top+n;
res = [res;X(1)-norm(full(X(2:end)))];
end
end
if p.K.s(1)>0
top = 1+p.K.f+p.K.l+p.K.q;
for i = 1:length(p.K.s)
n = p.K.s(i);
X = reshape(vecres(top:top+n^2-1),n,n);top = top+n^2;
res = [res;min(eig(X))];
end
end
end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -