[求助]如何快速的在求解析解后代入数值计算？

[felonwan]

用matlab求得一个高次方程的解析解（sym数组），通常很复杂，有时没完全必要管它到底什么结果，只是想把原来的sym变量转化为double变量（数组），代入求数值解。

例如：由德瓦耳斯方程在给定已知a，b，μ，T，P的情况下求ρ
syms a b r u p R T;
rho=simple(solve('a*r^3-a*b/u*r^2-u*p*r-b*p-R*T','r'));

除是对应的n×1 double型外，其它未知参数都是1×1 double。
方程有3个解，所以rho是3×1 sym变量。

之后该怎么快速地代入数据，求解T，P对应的ρ？……

[watcher]

在获得解析表达式之后，对符号变量赋“数值”，然后用eval函数求解析表达式的值。