我有一个优化问题，我想通过符号化的方式来优化厚度t1, t2, t3。解决方案应大致为：

t1 = 1.72053; 
t2 = 1.91395;
t3 = 0.781725;

%Given
M_y             = 2000000;  %[Nmm]
b               = 100;        %[mm]
E               = 70000;      %[N/mm^2]
k_Druck         = 3.62;
k_Biegung       = 21.7;
sigma_zul       = 100;        %[N/mm^2]

syms t1 t2 t3 positive

% center of gravity
zGes =(b^2*t1+b^2*t3)/(t1*b+t2*b+2*t3*b);

zSeite = zGes - b/2;
zUnten = b - zGes;
zOben= -zGes;

%area moment of inertia
I_y = 2*((t3*b^3/12)+zSeite^2*t3*b)+ ((b*t1^3/12)+zUnten^2*t1*b)  +  ((b*t2^3/12)+zOben^2*t2*b);

sigma_Oben       = -M_y / I_y * zOben;           %Stabilitätsversagen (druckbelastete Oberseite)
sigma_Unten      = M_y / I_y * zUnten;          %Festigkeitsversagen (zugbelasteten Unterseite)

sigma_krit_Oben      =k_Druck * E * (t2/b)^2;       %kritische Spannung Oberseite
sigma_krit_Seiten    =k_Biegung * E * (t3/b)^2;   %kritische Spannung Seitenwand

%% condition
Beulen_Oben = sigma_Oben == sigma_krit_Oben;
Beulen_Seiten= sigma_Oben == sigma_krit_Seiten;
Festig_Unten = sigma_Unten==sigma_zul;

%Solver
eqns = [Beulen_Oben,Beulen_Seiten,Festig_Unten];
vars = [t1, t2, t3];
[t1, t2, t3]=vpasolve(eqns,vars)

我试图使用Matlab中的vpasolve函数来解决这个问题，但是在positive条件下，函数总是输出Empty sym: 0-by-1。当我输入给定的厚度解时，给定的条件几乎满足。我假设他们是一个有点小康由于小数位数，我输入。