无法从 sym 转换为 logical怎么办？

tuoni · 发表于 2020-5-7 11:19

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

运行报错信息如下：

无法从 sym 转换为 logical。

出错 M (line 65)
while differ>eps

代码如下：
clc
clear
%Newton iteration-Multi Points

syms x1 x2 x3 x4 x5 x6 sum_n n_NaF n_AlF3;
syms RR;
syms N
RR=ones(1,7)/10e11;

%****************
T=['n_AlF3  ','x1  ', 'x2  ', 'x3  ', 'x4  ', 'x5  ', 'x6  ','n_sum  '];       %
disp(T);

N=0;
for n_AlF3=0.01:0.01:0.5;             %%

disp(n_AlF3);
disp('--------');

if N==0
      s1=1-3*n_AlF3;
      s2=0;
      s3=0;
      s4=n_AlF3/3;
      s5=0;
      s6=1-n_AlF3;
      s7=0;
elseif N>0
      s1=RR(1);
      s2=RR(2);
      s3=RR(3);
      s4=RR(4);
      s5=RR(5);
      s6=RR(6);
      s7=RR(7);
end

eps=1.0e-6;
differ=1;
n=1;

f1=x1*x3-0.2471*sum_n*x2; %0.247146 改为0.2471
f2=x1*x4-0.1352*sum_n*x3;%0.135198 改为0.1352
f3=x1*x5-0.0003*x4^2; %0.000299 改为0.0003
f4=x1+6*x2+5*x3+4*x4+7*x5+14*x6-( n_NaF+3*n_AlF3);
f5=x2+x3+x4+2*x5+3*x6-n_AlF3;
f6=x1+x2+x3+x4+x5+x6-sum_n;
f7=x4^2*x2-0.0000008*x6*sum_n^2;
F=[f1 f2 f3 f4 f5 f6 f7];

dF=transpose([diff(F,'x1');diff(F,'x2') ;diff(F,'x3');diff(F,'x4');diff(F,'x5'); diff(F,'sum_n'); diff(F,'x6')]);

x0 = [s1,s2,s3,s4,s5,s6,s7];

Fx = subs(F,[x1 x2 x3 x4 x5 sum_n x6],x0);
dFx = subs(dF, [x1 x2 x3 x4 x5 sum_n x6],x0);
%disp('*******************');
%disp(vpa(Fx));
%disp(vpa(dFx));

r= transpose(x0)- inv(dFx)*transpose(Fx);
%disp(vpa(r));

while differ>eps
      %disp('*******************');
      Fx = subs(F, [x1 x2 x3 x4 x5 sum_n x6], transpose(r));
      dFx = subs(dF, [x1 x2 x3 x4 x5 sum_n x6], transpose(r));
      %disp(vpa(Fx));
      %disp(vpa(dFx));

      differ=norm( inv(dFx)*transpose(Fx));

      r=r-inv(dFx)*transpose(Fx);

      %disp(vpa(r));

      n =n+1;
      if(n>100);
         disp('iteration exceed! Not converged!');
         return;
      end;
end;

%r=[n_AlF3,r];
RR=vpa(transpose(r));
disp(RR);
N=N+1;
end;

shelby · 发表于 2020-5-7 13:50

这个没有见过。

tuoni · 发表于 2020-5-7 15:27

已解决。

aplerll · 发表于 2020-5-8 17:37

楼主怎么解决的能不能分享一下啊。

帐号		自动登录	找回密码
密码			注册

无法从 sym 转换为 logical怎么办？

EDA365欢迎您登录！

推荐内容 /1