matlab解非线性方程组解不出来:代码如下:
来源:学生作业帮 编辑:大师作文网作业帮 分类:综合作业 时间:2024/11/11 07:43:35
matlab解非线性方程组解不出来:代码如下:
global b T temp;
b=1.2;
T=15;
temp=1;
syms a R;
[a,R]=solve('a-(R-b)/2+0.5*((R-b)^2-0.281*R*T)^0.5=0','13*temp-13*(a+b)*a-0.912*T*R=0')
运行结果如下:
a =
-(0.25*(0.281*R*T + (R - 2.0*a)^2 - 1.0*R^2))/a
R =
0.070153846153846153846153846153846*R*T + a*(a - (0.25*(0.281*R*T + (R - 2.0*a)^2 - 1.0*R^2))/a)
global b T temp;
b=1.2;
T=15;
temp=1;
syms a R;
[a,R]=solve('a-(R-b)/2+0.5*((R-b)^2-0.281*R*T)^0.5=0','13*temp-13*(a+b)*a-0.912*T*R=0')
运行结果如下:
a =
-(0.25*(0.281*R*T + (R - 2.0*a)^2 - 1.0*R^2))/a
R =
0.070153846153846153846153846153846*R*T + a*(a - (0.25*(0.281*R*T + (R - 2.0*a)^2 - 1.0*R^2))/a)
=1.2;T=15;temp=1;
syms a R
s=solve(a-(R-b)/2+0.5*((R-b)^2-0.281*R*T)^0.5,13*temp-13*(a+b)*a-0.912*T*R);
>> vpa(s.R)
ans =
-0.50580706462237194771798754220257
>> vpa(s.a)
ans =
-1.9755961125729077963914624034358
syms a R
s=solve(a-(R-b)/2+0.5*((R-b)^2-0.281*R*T)^0.5,13*temp-13*(a+b)*a-0.912*T*R);
>> vpa(s.R)
ans =
-0.50580706462237194771798754220257
>> vpa(s.a)
ans =
-1.9755961125729077963914624034358