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

tuoni

& P$ m% C3 O! C运行报错信息如下：

无法从 sym 转换为 logical。
5 |) v5 I; x$ A; I" h0 Z
- O5 I' k" t; e4 V: ]7 F" G出错 M (line 65)
& U  T% m5 t" X& f6 i  q    while differ>eps
) y' b, G6 G- B* x# Y

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


1 n) c4 O+ g2 X# asyms x1 x2 x3 x4 x5 x6 sum_n n_NaF n_AlF3;
0 _! N2 W! p2 Y' f9 @9 nsyms RR;
9 _! r. E) c) A- W6 m* |! J. `syms N
RR=ones(1,7)/10e11;
4 v( X  N5 ?4 J0 {$ L

; G+ T) [; f  M8 \" }

) w: ]  E; I5 D( |3 O, `5 ?3 s/ j%****************
. R9 W2 \  [! ET=['n_AlF3  ','x1  ', 'x2  ', 'x3  ', 'x4  ', 'x5  ', 'x6  ','n_sum  '];        %
) o4 H- [" i& j! t% ]2 a6 Zdisp(T);
: R% P& t! Y4 q( {% f( Y! B
" M; ?; C2 X- k1 S8 a' K0 m  j. i
N=0;
, n* Y. t: O2 }: ^& rfor n_AlF3=0.01:0.01:0.5;                %%
   
    disp(n_AlF3);
: k/ r0 h4 l4 H  r3 G    disp('--------');
   
6 [2 a2 {8 `* L  y  M/ C    if N==0
        s1=1-3*n_AlF3;
        s2=0;
  H. \! z8 ]5 k% D# p& O1 o! F5 W( N        s3=0;
5 g1 ~+ Y0 e/ k; S5 m) f) b        s4=n_AlF3/3;
        s5=0;
        s6=1-n_AlF3;
        s7=0;
    elseif N>0
; u3 _, ]1 n* L( ~% e; F* t        s1=RR(1);
        s2=RR(2);
        s3=RR(3);
        s4=RR(4);
        s5=RR(5);
$ O" r/ k  q; Z        s6=RR(6);
        s7=RR(7);
    end
3 b2 A0 e5 U1 ]: `$ X   
    eps=1.0e-6;
    differ=1;
    n=1;
   
5 n  {1 o% g6 V" l1 k8 X. l    f1=x1*x3-0.2471*sum_n*x2; %0.247146 改为0.2471
1 n+ h8 ~. z5 o    f2=x1*x4-0.1352*sum_n*x3;%0.135198 改为0.1352
    f3=x1*x5-0.0003*x4^2; %0.000299 改为0.0003
. j# z# \% [' ?' o+ ^, u    f4=x1+6*x2+5*x3+4*x4+7*x5+14*x6-( n_NaF+3*n_AlF3);
+ Y' {3 c# }5 p! ]8 y6 x# k    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];
; a0 b3 Y6 d2 ~8 L   
    dF=transpose([diff(F,'x1');diff(F,'x2') ;diff(F,'x3');diff(F,'x4');diff(F,'x5'); diff(F,'sum_n'); diff(F,'x6')]);
   
0 o4 H. I+ u7 z4 O3 p: p    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);
# D" l1 ~7 }( }( P3 d8 p6 X    %disp('*******************');
2 p* W/ |0 k9 F: m    %disp(vpa(Fx));
" }# m' q. p+ ]) g2 ^    %disp(vpa(dFx));
   
' ?; H9 t5 N: I' {1 ]9 K  j0 a: W    r= transpose(x0)- inv(dFx)*transpose(Fx);
+ r) n' ^$ x) @- S, e& D# K    %disp(vpa(r));
   
0 s7 G# H9 @$ r. t8 v# R' X7 G    while differ>eps
3 r2 o4 G2 C3 A, O$ T& x5 @        %disp('*******************');
        Fx = subs(F, [x1 x2 x3 x4 x5 sum_n x6], transpose(r));
# F1 _& S7 M$ ]        dFx = subs(dF, [x1 x2 x3 x4 x5 sum_n x6], transpose(r));
        %disp(vpa(Fx));
! C, \/ J% f. k; c        %disp(vpa(dFx));
) @; K7 g# t1 s5 Y9 e/ R        
        differ=norm( inv(dFx)*transpose(Fx));
9 f$ f5 h* h) d8 Z3 N" {% L+ G. |        
9 @2 u( x$ \; ?( l4 `        r=r-inv(dFx)*transpose(Fx);
        
5 W) p& w' m' w        %disp(vpa(r));
        
1 X: z3 D% N/ N        n =n+1;
3 w) p, B, K7 j- r( @0 N1 f0 V) x        if(n>100);
            disp('iteration exceed! Not converged!');
$ C8 q6 L: P% ~1 [# Q6 X9 Q9 Z            return;
3 j2 ]1 a0 V* J8 l1 R        end;
    end;
   
% b- e& j; l: ?3 R1 l, x9 X0 e1 v, J3 A    %r=[n_AlF3,r];
  [; a1 {# U/ B$ f    RR=vpa(transpose(r));
7 `9 |/ s8 A0 L, V4 g  B- ?# K9 V    disp(RR);
    N=N+1;
end;
! A; A7 D7 b" G
- O1 Z: u4 e8 u/ I8 R) K/ R

shelby

这个没有见过。

tuoni

已解决。

aplerll

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