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

tuoni

运行报错信息如下：

, V  G1 D, u. d无法从 sym 转换为 logical。

出错 M (line 65)
    while differ>eps

+ V( X9 [0 C/ f+ R5 M+ x/ {
代码如下：
clc
% U6 r4 L9 v/ z* I7 x* \+ [+ Kclear
%Newton iteration-Multi Points


syms x1 x2 x3 x4 x5 x6 sum_n n_NaF n_AlF3;
syms RR;
7 w( Q1 g+ x! t4 ~& L- D: lsyms N
RR=ones(1,7)/10e11;


7 L0 G4 _7 m6 ^5 Z9 i8 f

/ x7 ?9 A( ]9 A8 A+ m, H& i%****************
3 Q( @/ D) f. r( J: sT=['n_AlF3  ','x1  ', 'x2  ', 'x3  ', 'x4  ', 'x5  ', 'x6  ','n_sum  '];        %
disp(T);
: X6 o6 E% @4 I4 Q5 {# X
6 `# T; g2 b; m5 L# W* m
N=0;
* x- _, N3 H+ Q" D2 I% m9 j6 k& A' Efor n_AlF3=0.01:0.01:0.5;                %%
" [8 ?6 q4 }% h* A1 p3 k) j1 M   
0 K; N; n0 q' \4 r    disp(n_AlF3);
* Y9 s, ~4 _& n1 a& [2 U    disp('--------');
* N, d5 V" v6 q* Q! g   
! Q7 [7 i: h- K+ s! \    if N==0
# R# Q+ D8 P" Q        s1=1-3*n_AlF3;
        s2=0;
$ Z6 M: X5 I$ n, H, l& h4 ~6 K; u        s3=0;
        s4=n_AlF3/3;
' M- l5 V8 E" u% U+ B! ^% m        s5=0;
        s6=1-n_AlF3;
        s7=0;
' F+ \* n1 w0 ^, |* s' l/ u) v1 ?) b- f$ c    elseif N>0
        s1=RR(1);
        s2=RR(2);
        s3=RR(3);
% D& p9 O2 }4 p5 Q) Q2 e        s4=RR(4);
% |) o5 @% m: c3 {) Z        s5=RR(5);
        s6=RR(6);
( E8 h# Z. s% w0 M  k: b        s7=RR(7);
    end
6 a9 f4 |4 p1 ~% k+ J. T# d7 a3 y   
    eps=1.0e-6;
( e3 e/ |( D( h. n( X    differ=1;
% |/ F" R" b5 [: d' e, }# x    n=1;
" N, j! q+ T1 w* r$ H' [   
    f1=x1*x3-0.2471*sum_n*x2; %0.247146 改为0.2471
6 A: m- E3 k3 U$ U    f2=x1*x4-0.1352*sum_n*x3;%0.135198 改为0.1352
3 A6 w; ~- ]) D& @4 G+ J    f3=x1*x5-0.0003*x4^2; %0.000299 改为0.0003
3 S) Q( S/ H% n: ]    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;
0 n0 ~4 K' I3 }. @; F    f6=x1+x2+x3+x4+x5+x6-sum_n;
    f7=x4^2*x2-0.0000008*x6*sum_n^2;
8 k0 K1 O. }3 ^& W! O. n    F=[f1 f2 f3 f4 f5 f6 f7];
- i+ _8 R- u" `2 `/ S3 z   
" p0 ?+ H8 c' i! ]9 T    dF=transpose([diff(F,'x1');diff(F,'x2') ;diff(F,'x3');diff(F,'x4');diff(F,'x5'); diff(F,'sum_n'); diff(F,'x6')]);
   
* Y1 O( e2 n, \" v+ T    x0 = [s1,s2,s3,s4,s5,s6,s7];
   
2 s2 t/ T0 g8 f1 u. T    Fx = subs(F,[x1 x2 x3 x4 x5 sum_n x6],x0);
- T$ S9 S+ T9 s9 Z: i$ {/ T) J3 m7 ^    dFx = subs(dF, [x1 x2 x3 x4 x5 sum_n x6],x0);
3 Y5 d2 C7 \- H" N8 V    %disp('*******************');
& c1 a$ j" c$ p5 G! y+ p. X    %disp(vpa(Fx));
    %disp(vpa(dFx));
   
+ O$ l1 r5 i- r  a% J1 P    r= transpose(x0)- inv(dFx)*transpose(Fx);
6 ?0 s; a* t3 x  q3 B; P- B6 T    %disp(vpa(r));
: m5 |6 U" J- n$ i   
    while differ>eps
6 J0 @- J4 o& }5 o* v        %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));
2 [- a- x' C) T3 H) Y' x  [' X' ~) k        %disp(vpa(Fx));
        %disp(vpa(dFx));
        
        differ=norm( inv(dFx)*transpose(Fx));
        
        r=r-inv(dFx)*transpose(Fx);
7 T7 u8 h0 r. I% `( s        
        %disp(vpa(r));
        
        n =n+1;
        if(n>100);
            disp('iteration exceed! Not converged!');
6 e+ G: P' k0 G- f1 w8 Z            return;
# a8 a) b& t% v9 Z- ~        end;
/ J) `: A- q) `6 N7 |! K% D    end;
   
3 V/ t- t2 X8 l+ g) x) H% ^) c" M    %r=[n_AlF3,r];
2 A( ^( Q1 Y! `4 u    RR=vpa(transpose(r));
    disp(RR);
  o- L6 S/ r7 |2 y    N=N+1;
end;


% v5 {& x5 Z0 B; x3 E2 Y

shelby

这个没有见过。

tuoni

已解决。

aplerll

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