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

tuoni

运行报错信息如下：
5 J9 b1 h0 R# e* X+ ?/ F
0 y# [) M! r4 X$ c: U+ K! `无法从 sym 转换为 logical。

. I- Z/ z! C8 m/ e6 i出错 M (line 65)
    while differ>eps

/ _( K( b6 S) l: |4 a) ?
& Z/ n8 w% R% m: ?7 H代码如下：
# g- q: ^. o4 H7 w  qclc
clear
) e2 Z( j# x! i%Newton iteration-Multi Points

6 f7 @& B" _5 `& n) h
syms x1 x2 x3 x4 x5 x6 sum_n n_NaF n_AlF3;
syms RR;
syms N
RR=ones(1,7)/10e11;



8 u6 w( q* L' V/ f- O3 t; I4 M
% l+ h0 |) R2 D" E%****************
0 g% k2 k7 S7 Y8 ~T=['n_AlF3  ','x1  ', 'x2  ', 'x3  ', 'x4  ', 'x5  ', 'x6  ','n_sum  '];        %
$ q! g5 N/ H& j/ Ddisp(T);
% Z8 U7 {2 u/ s) H5 l. p

N=0;
8 P( x8 D6 u- g9 tfor n_AlF3=0.01:0.01:0.5;                %%
4 e! U+ ]2 U& ~8 y   
* n6 s9 W* ]' i3 n3 p6 B3 R0 Z    disp(n_AlF3);
$ v: r. t9 V7 X4 \1 e    disp('--------');
4 Y4 O  k! q9 ?+ o  l  b   
! m2 i2 r" d% ^, Y" G. w    if N==0
2 S& p/ L* l, H& |3 |        s1=1-3*n_AlF3;
- Q2 _, A4 V1 B, @) p( A2 {% ]. q0 f        s2=0;
        s3=0;
- W* S- X2 I# \+ ~4 r7 c        s4=n_AlF3/3;
+ q+ X$ {: g1 Z8 A" q        s5=0;
        s6=1-n_AlF3;
# ~1 S  H( S  `) p/ y        s7=0;
1 z4 e; z3 ?( X    elseif N>0
9 G9 i- v; T- H- m! s        s1=RR(1);
        s2=RR(2);
2 U8 ~0 H; e" d4 j$ Q        s3=RR(3);
2 K4 ~' @: E! m/ c4 p3 l        s4=RR(4);
) P8 x: M& M" M# z8 N( v        s5=RR(5);
5 d, M) S: i- {8 b. U1 e9 I        s6=RR(6);
& B9 e3 L( T5 s3 h        s7=RR(7);
    end
; f% Q9 i& H3 ]" k3 `8 N# V7 ~   
    eps=1.0e-6;
    differ=1;
    n=1;
   
    f1=x1*x3-0.2471*sum_n*x2; %0.247146 改为0.2471
3 G" a8 X# Q$ p- a+ Q, I. B    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);
$ Y4 g% b8 _! T& [    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')]);
8 y1 s7 Q+ \9 j6 V/ v: Z7 h   
    x0 = [s1,s2,s3,s4,s5,s6,s7];
) W! Q: k: ^/ y& [1 z   
    Fx = subs(F,[x1 x2 x3 x4 x5 sum_n x6],x0);
! e$ O2 x+ J. L8 n0 A" |2 S. ?9 ?    dFx = subs(dF, [x1 x2 x3 x4 x5 sum_n x6],x0);
/ e8 }) n5 B* Y3 Y/ A/ U    %disp('*******************');
6 I! U2 M6 L5 f9 a% n    %disp(vpa(Fx));
    %disp(vpa(dFx));
! D! C- v% Q7 M! ?4 L   
    r= transpose(x0)- inv(dFx)*transpose(Fx);
    %disp(vpa(r));
' w7 M' {: b' n7 J% W; s   
    while differ>eps
        %disp('*******************');
        Fx = subs(F, [x1 x2 x3 x4 x5 sum_n x6], transpose(r));
" L8 @3 O1 ^3 M- C        dFx = subs(dF, [x1 x2 x3 x4 x5 sum_n x6], transpose(r));
        %disp(vpa(Fx));
* y2 F* S0 e5 u# ^7 [/ f' v9 f        %disp(vpa(dFx));
        
2 M% p! @- o; c, t7 Q6 r# s        differ=norm( inv(dFx)*transpose(Fx));
; M' X; L' N: r6 z. A1 D7 Q        
        r=r-inv(dFx)*transpose(Fx);
" v0 ^" {! ?0 }+ |! n; O; B/ ~+ w- y/ m        
5 _* I( p8 U. p        %disp(vpa(r));
! Y' L* D$ A# R' k        
* A! `6 U$ M( g0 `        n =n+1;
( b5 H' ^% h( \; e1 M        if(n>100);
; R2 u- S: G0 f$ M            disp('iteration exceed! Not converged!');
            return;
        end;
. n! e+ z. x- I3 M7 u    end;
   
5 H' I9 [3 u# b% A5 h2 t) [' h    %r=[n_AlF3,r];
- ?" d+ C( d; q    RR=vpa(transpose(r));
7 r) B. R; H* N: _* E    disp(RR);
    N=N+1;
end;

" B) [. g0 t; ~2 J- h4 }7 ~# \

shelby

这个没有见过。

tuoni

已解决。

aplerll

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