sedumi([1; 1],1,[1 1])
close allclear allclc
Pi=[3.4309 1.4534; 1.4534 3.3298];VAe=[-0.1946 -0.0045 -0.1897 -0.0186 -0.1793 -0.0483 -0.1622 -0.0964 -0.1408 -0.1542 -0.1218 -0.2039 -0.1100 -0.2340 -0.1046 -0.2477 -0.0112 -0.2018 -0.0000 -0.1961 0.0000 -0.1961 0.0003 -0.1958 0.0031 -0.1932 0.0112 -0.1853 0.0278 -0.1693 0.0546 -0.1431 0.0896 -0.1084 0.1265 -0.0714 0.1574 -0.0399 0.1785 -0.0183 0.1901 -0.0063 0.1950 -0.0012 0.1961 -0.0000 0.1961 0.0000 0.1946 0.0045 0.1897 0.0186 0.1793 0.0483 0.1622 0.0964 0.1408 0.1542 0.1218 0.2039 0.1100 0.2340 0.1046 0.2477 0.0112 0.2018 0.0000 0.1961 -0.0000 0.1961 -0.0003 0.1958 -0.0031 0.1932 -0.0112 0.1853 -0.0278 0.1693 -0.0546 0.1431 -0.0896 0.1084 -0.1265 0.0714 -0.1574 0.0399 -0.1785 0.0183 -0.1901 0.0063 -0.1950 0.0012 -0.1961 0.0000 -0.1961 -0.0000];Vbe = [0.9809 0.9817 0.9826 0.9820 0.9780 0.9714 0.9660 0.9632 0.9794 0.9806 0.9806 0.9806 0.9812 0.9826 0.9852 0.9882 0.9901 0.9894 0.9867 0.9838 0.9818 0.9808 0.9806 0.9806 0.9809 0.9817 0.9826 0.9820 0.9780 0.9714 0.9660 0.9632 0.9794 0.9806 0.9806 0.9806 0.9812 0.9826 0.9852 0.9882 0.9901 0.9894 0.9867 0.9838 0.9818 0.9808 0.9806 0.9806];
sdpvar sqrtalphaBB=sqrtalpha*Pi^(-1/2);cost=-logdet(BB);%log(det(Binv)), we know det(Binv)>=0cond=[];for i=1:length(Vbe) cond=[cond norm(BB*VAe(i,:)',2)<=Vbe(i)];endoptimize(cond,cost)%,sdpsettings('debug',1));
alpha=value(sqrtalpha)^2;
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++| Test| Solution| Solver message|+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++| Core functionalities| N/A| Successfully solved (YALMIP)|| LP| Correct| Successfully solved (GLPK-GLPKMEX-CC)|| LP| Correct| Successfully solved (GLPK-GLPKMEX-CC)|| QP| Correct| Successfully solved (QUADPROG)|| QP| Correct| Successfully solved (QUADPROG)|| SOCP| N/A| Unknown problem in solver (Turn on 'debug' in sdpsettings) (Undefined function or variable "y_s".)|| SOCP| N/A| Unknown problem in solver (Turn on 'debug' in sdpsettings) (Undefined function or variable "y_s".)|| SOCP| N/A| Unknown problem in solver (Turn on 'debug' in sdpsettings) (Undefined function or variable "y_s".)|| SDP| N/A| Unknown problem in solver (Turn on 'debug' in sdpsettings) (Undefined function or variable "y_s".)|| SDP| N/A| Unknown problem in solver (Turn on 'debug' in sdpsettings) (Undefined function or variable "y_s".)|| SDP| N/A| Unknown problem in solver (Turn on 'debug' in sdpsettings) (Undefined function or variable "y_s".)|| SDP| N/A| Unknown problem in solver (Turn on 'debug' in sdpsettings) (Undefined function or variable "y_s".)|| MAXDET| N/A| Unknown problem in solver (Turn on 'debug' in sdpsettings) (Undefined function or variable "y_s".)|| MAXDET| N/A| Unknown problem in solver (Turn on 'debug' in sdpsettings) (Undefined function or variable "y_s".)|| Infeasible LP| N/A| Infeasible problem (GLPK-GLPKMEX-CC)|| Infeasible QP| N/A| Infeasible problem (QUADPROG)|| Infeasible SDP| N/A| Unknown problem in solver (Turn on 'debug' in sdpsettings) (Undefined function or variable "y_s".)|| Moment relaxation| N/A| Unknown problem in solver (Turn on 'debug' in sdpsettings) (Undefined function or variable "y_s".)|| Sum-of-squares| N/A| Unknown problem in solver (Turn on 'debug' in sdpsettings) (Undefined function or variable "y_s".)|| Bilinear SDP| N/A| No suitable solver|+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
which maxstep -all
>> which maxstep -allC:\work\solvers\sedumi\maxstep.mC:\Program Files\MATLAB\R2015a\toolbox\dsp\filterdesign\@adaptfilt\maxstep.m % adaptfilt methodC:\Program Files\MATLAB\R2015a\toolbox\dsp\dsp\maxstep.m % Shadowed