先介绍配置,matlab2010a,maple15,装了Sedumi和sostools两个工具包
碰到问题,上结果
??? Error using ==> sprintf
Function is not defined for sparse inputs.
Error in ==> int2str at 16
s = sprintf('%.1f',x); % .1 to avoid precision loss on hpux
Error in ==> mysympower at 39
charexpr = [charexpr,char(vartable(i)),'^',int2str(Z(j,i)),'*'];
Error in ==> monomials>oldconstructZ at 96
Z = mysympower(vartable,Z);
Error in ==> monomials at 69
Z = [Z; oldconstructZ(vartable,i)];
这是代码,谁给分析下结果,先谢谢啦(为什么别人可以出结果)
clear; clc
syms x y;
qu1=[0.9,1];tqu=[0 0.4];
Program1 = sosprogram([x;y]);
VEC = monomials([x; y],[0 1 2 3 4]);
[Program6,v] = sospolyvar(Program1,VEC,'wscoeff');
Program6 = sosineq(Program6,subs(v,y,0.4)+x-1);
Program6 = sosineq(Program6,-diff(v,y)-0.5*0.09*x^2*diff(diff(v,x),x));
Program6= sossetobj(Program6,subs(v,{x,y},{1,0}));
Program6= sossolve(Program6);
SOLcoeff_1 = sosgetsol(Program6,coeff_1);
SOLcoeff_2 = sosgetsol(Program6,coeff_2);
SOLcoeff_3= sosgetsol(Program6,coeff_3);
SOLcoeff_4= sosgetsol(Program6,coeff_4);
SOLcoeff_5 = sosgetsol(Program6,coeff_5);
SOLcoeff_6 = sosgetsol(Program6,coeff_6);
SOLcoeff_7 = sosgetsol(Program6,coeff_7);
SOLcoeff_8= sosgetsol(Program6,coeff_8);
SOLcoeff_9 = sosgetsol(Program6,coeff_9);
SOLcoeff_10 = sosgetsol(Program6,coeff_10);
SOLcoeff_11 = sosgetsol(Program6,coeff_11);
SOLcoeff_12 = sosgetsol(Program6,coeff_12);
SOLcoeff_13 = sosgetsol(Program6,coeff_13);
SOLcoeff_14 = sosgetsol(Program6,coeff_14);
SOLcoeff_15 = sosgetsol(Program6,coeff_15);