matlab求解复杂方程为零的符号解,d为要求的变量,代码如下,
clear;
clc;
syms a b c o k p v O B1 B2 B r eo f t d
g2 = (7*(a-b*c)^2*k)/(32*b*(p+v)) - (eo*f*(a-b*c)*k)/(4*(p+v));
g1 = (g2*o+f)/(p+O);
m1 = ((a-b*c)^2*o*k)/(8*b*(p+v)*(p+O));
m2 = ((a-b*c)^2*k)/(8*b*(p+v));
r1 = ((a-b*c)^2*o*k)/(16*b*(p+v)*(p+O));
r2 = ((a-b*c)^2*k)/(16*b*(p+v));
tao2 = ((r^2 * m1*exp(v*d))*(1 - exp(-O * (t - d))))/(O * B1);
G2 =(B^2*r2*(1-exp(-v*t)))/(B2*v)+(o*r^2*m1*exp(O*d)*(1-exp(-v*(t-d))))/(O*B1*v)+(o*r^2*m1*exp(O*d)*((exp(-v*(t-d))-exp(-O*(t-d)))))/(O*B1*(v-O));
m32 = (m1^2*r^2*exp(O*d)*(1-(exp(O*d))/2))/(p*B1) + (B^2*m2*r2)/(p*B2);
r32 = (r1*m1*r^2*exp(O*d))/(p*B1) + (B^2*r2^2)/(2*p*B2);
[d]=solve([m1*tao2 + m2*G2 + m32-((m2*r2*B^2*(p+v))/(B2*p*v))],[d]);
av = 100;bv = 0.5;cv = 2;ov = 1;kv = 0.4;pv = 0.2;vv = 1;Ov = 0.5;B1v = 15;B2v = 13;Bv = 2;rv = 2;eov = 10;fv = 0.7;tv=10;
dv = subs(d,{a, b, c, o, k, p, v, O, B1, B2, B, r, eo, f, t},[av bv cv ov kv pv vv Ov B1v B2v Bv rv eov fv tv]);
为什么求不出来呢??
显示
d =
Empty sym: 0-by-1
【相关推荐】
因为最近 MATLAB 被卡脖子了,我在上大学的时候也使用 MATLAB ,确实是因为好用。
但是发现脱离不开这个环境,而且也不方便部署吧。按照的软件也是超级的大,很麻烦呢。
现在好了,看来天下真的没有免费的午餐呢。让人家卡脖子了。
说明solve函数找不到解