r=[1 2 3];
r=subs( solve( 'v^3-8.314*T0*v^2/P+aV*v/P-bV*8.314*T0*v/P-bV^2*v-aV*bV/P==0' ,'v') ); %解一元三次方程
d=max(r); %求三个根中的最大值
Undefined function 'max' for input arguments of type 'sym'.
望高手指点一二、·~
谢谢你的回答,参数的值是已知的,关于v的一元三次方程~解到的三个根是0.0012800607434096664775264967576924,
- 0.000031029871704833238763248378846181 -0.0076485053529296446744317748441282*i,和
- 0.000031029871704833238763248378846181 + 0.0076485053529296446744317748441282*i
需要找到最大的正实根,用什么命令执行呢?~~谢谢
先用isreal判断,看有几个实根,如果多于1个,才有必要继续寻找最大的
你这个只有3个根,且只有1个实根,实在是没必要继续找了
clear all;clc;
x1=0.0012800607434096664775264967576924;
x2=-0.000031029871704833238763248378846181...
-0.0076485053529296446744317748441282*i;
x3=-0.000031029871704833238763248378846181...
+0.0076485053529296446744317748441282*i;
xx1=isreal(x1);
xx2=isreal(x2);
xx3=isreal(x3);
我举的是个特例哈,改变参数T0的值时v的根要发生变化,所以想用一般方法找最大正实根~~谢~~