Question

matlab解方程

程序如下，请大家帮忙看看怎么解s，一直运行出错

clc
clear
syms s;
a=0.16; b=1.55;
d1=0.3; d2=0.3;
u=(-(a+b-1)+sqrt((a+b-1).^2+4.*a))./2; v=u;
a11=1-2.*u-(a.*b.*v)/((u+a).^2); a12=-(b.*u)./(u+a);

d12=2;
A=u.*(1-2.*u-(a.*b.*v)./((u+a).^2)+b.*v./(u+a))
B1=s.*(d1+d12.*v+d12.*u)-(1-2.*u-(a.*b.*v)./((u+a).^2)).*d2
C1=a11.*(-s)-a12.*s;
A1=4.*d1.*u.*C1;
A2=4.*C1.*(d1.*d2+d2.*d12.*v+d1.*u.*B1./A);
e=(A1+sqrt(A1.^2+4.*A.^2.*A2))./(2.*A.^2)
s=solve('B1./A+e=6')

结果是这样
??? Error using ==> solve>getEqns at 202
' B1./A+e=6 ' is not a valid expression or equation.

Error in ==> solve at 67
[eqns,vars] = getEqns(varargin{:});

Error in ==> Untitled at 16
s=solve('B1./A+e=6')

各位大神帮帮忙

Answer 1

你的错误，就是这一句出了问题。即  s=solve('B1./A+e=6') 

解决方法很简单，把solve函数里的表达式的单引号去掉就可以了。即为

s=solve(1./A+e==6) 

说明，对于低版本的应这样改，s=solve(1./A+e-6) 

更改后运行，可以得到

s=0.47628

建议，在solve()函数命令前，再套一个命令，即

s=double(solve(B1./A+e==6))   %数值为双精度