MATLAB爱好者论坛-LabFans.com

MATLAB爱好者论坛-LabFans.com (https://www.labfans.com/bbs/index.php)
-   MATLAB论坛 (https://www.labfans.com/bbs/forumdisplay.php?f=6)
-   -   [MATLAB基础] 为什么我这个方程没有解? (https://www.labfans.com/bbs/showthread.php?t=11677)

edongliang 2010-06-01 21:50

为什么我这个方程没有解?
 
for ii=0:0.05:1
x=vpa(double(solve([COLOR="Red"]9*u^2/((1/3*(9*(3^(1/2)-2/3*3^(1/2)*u)^2+6*u^2)^(1/2)+1/3*(9+36*(1/2*3^(1/2)-2/3*3^(1/2)*u)^2+24*u^2)^(1/2))*(9*(3^(1/2)-2/3*3^(1/2)*u)^2+6*u^2)*(1/3*(9+36*(1/2*3^(1/2)-2/3*3^(1/2)*u)^2+24*u^2)^(1/2)-1/3*(9*(3^(1/2)-2/3*3^(1/2)*u)^2+6*u^2)^(1/2)))^(1/2)/(1/4*3^(1/2)+2*((1/2+1/6*(9*(3^(1/2)-2/3*3^(1/2)*u)^2+6*u^2)^(1/2)+1/12*(9+36*(1/2*3^(1/2)-2/3*3^(1/2)*u)^2+24*u^2)^(1/2))*(-1/2+1/6*(9*(3^(1/2)-2/3*3^(1/2)*u)^2+6*u^2)^(1/2)+1/12*(9+36*(1/2*3^(1/2)-2/3*3^(1/2)*u)^2+24*u^2)^(1/2))*(1/2-1/6*(9*(3^(1/2)-2/3*3^(1/2)*u)^2+6*u^2)^(1/2)+1/12*(9+36*(1/2*3^(1/2)-2/3*3^(1/2)*u)^2+24*u^2)^(1/2))*(1/2+1/6*(9*(3^(1/2)-2/3*3^(1/2)*u)^2+6*u^2)^(1/2)-1/12*(9+36*(1/2*3^(1/2)-2/3*3^(1/2)*u)^2+24*u^2)^(1/2)))^(1/2)+1/2*((1/2+1/6*(9+36*(1/2*3^(1/2)-2/3*3^(1/2)*u)^2+24*u^2)^(1/2))*(-1/2+1/6*(9+36*(1/2*3^(1/2)-2/3*3^(1/2)*u)^2+24*u^2)^(1/2)))^(1/2))-[/COLOR]ii)),8)
// x=vpa(x,8)表示对X小数点后保留8位
end
x
上面的红色部分是matlab解出来的

只有x=0的解下面就不动了,像死机了一样的!!
这是什么原因?


所有时间均为北京时间。现在的时间是 15:06

Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.