主题: [MATLAB基础] 为什么我这个方程没有解?
查看单个帖子
旧 2010-06-01, 21:50   #1
edongliang
初级会员
 
注册日期: 2010-06-01
年龄: 37
帖子: 6
声望力: 0
edongliang 正向着好的方向发展
默认 为什么我这个方程没有解?

for ii=0:0.05:1
x=vpa(double(solve(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))-ii)),8)
// x=vpa(x,8)表示对X小数点后保留8位
end
x
上面的红色部分是matlab解出来的

只有x=0的解下面就不动了,像死机了一样的!!
这是什么原因?
edongliang 当前离线   回复时引用此帖