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

tuoni

. `5 [) Z  {! J& }运行报错信息如下：

无法从 sym 转换为 logical。

出错 M (line 65)
  ~* S: ^  _- k( V! h6 m6 F- _    while differ>eps

+ [  s4 J% E# b3 W
- ?9 E4 @9 v/ k  E! I6 `1 D代码如下：
clc
7 b; _( @2 P8 p5 ?! S, Y: s1 Y* pclear
3 A0 w% Z- G$ p: B9 |$ Z' K6 g%Newton iteration-Multi Points
3 k% G+ N5 k6 E: _( A- W; ?1 g

. h( Z) @$ q7 ^% ]( ssyms x1 x2 x3 x4 x5 x6 sum_n n_NaF n_AlF3;
0 s+ W5 A- G& p1 g  Jsyms RR;
syms N
& L% X7 q# K! K8 k* V5 g; F% P. tRR=ones(1,7)/10e11;

7 @' ^0 f' a+ k* }* l. ^
( u# O$ V4 o. I1 q8 S

%****************
- ]  ]1 |6 S% n" ET=['n_AlF3  ','x1  ', 'x2  ', 'x3  ', 'x4  ', 'x5  ', 'x6  ','n_sum  '];        %
disp(T);


N=0;
for n_AlF3=0.01:0.01:0.5;                %%
   
4 L2 a  N# Y6 ?& L    disp(n_AlF3);
    disp('--------');
   
3 z$ S: s. f2 ^; j    if N==0
        s1=1-3*n_AlF3;
        s2=0;
0 u& _; ]1 Q& P" n- a8 Y/ `0 j9 A        s3=0;
        s4=n_AlF3/3;
        s5=0;
; J& j' B. c8 B; E  l. K        s6=1-n_AlF3;
        s7=0;
0 i- `6 N* F# \9 t  q0 M2 E    elseif N>0
        s1=RR(1);
        s2=RR(2);
        s3=RR(3);
        s4=RR(4);
        s5=RR(5);
        s6=RR(6);
* j' k  q7 |% H% b7 v        s7=RR(7);
    end
   
9 j4 B  ~  V3 r! X. |    eps=1.0e-6;
5 w5 A2 ]2 v9 \    differ=1;
1 C9 a. U  z  r" L2 k9 i    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
8 t! X+ P6 Y* t* d6 i    f3=x1*x5-0.0003*x4^2; %0.000299 改为0.0003
1 E3 ]2 y# u% X8 F    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;
" ]  q4 E6 U: ^; i1 J6 d9 c    f6=x1+x2+x3+x4+x5+x6-sum_n;
2 E$ i2 d$ j6 E4 P1 M, v9 \    f7=x4^2*x2-0.0000008*x6*sum_n^2;
    F=[f1 f2 f3 f4 f5 f6 f7];
4 S( s* W4 n" {# |+ i& ~   
4 |2 H4 E/ m+ @8 Y9 e$ B    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 {- k4 `# H! l, 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);
    %disp('*******************');
' b- }6 G& J% C1 x( s; r% T6 o    %disp(vpa(Fx));
8 p) F1 ?% V4 ]; [# A* ?  I    %disp(vpa(dFx));
  f' {) n- U3 D8 l/ b$ f" w, M   
# H  R# L( f# @3 o2 X    r= transpose(x0)- inv(dFx)*transpose(Fx);
    %disp(vpa(r));
   
    while differ>eps
        %disp('*******************');
2 `+ L* B& \6 S& [. _" w) p( j        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));
% V. G0 `0 p! n3 J5 }' a9 V) L        
$ F9 Z# F& c& ?0 f: x        differ=norm( inv(dFx)*transpose(Fx));
  i* \2 ?7 n* o# V        
        r=r-inv(dFx)*transpose(Fx);
        
9 \& `$ k4 e; {8 X% h+ H        %disp(vpa(r));
7 f7 ]/ W+ N' }% N        
        n =n+1;
! \0 B$ U5 |- y9 Z        if(n>100);
            disp('iteration exceed! Not converged!');
' }5 K1 G( N5 e            return;
        end;
    end;
   
    %r=[n_AlF3,r];
    RR=vpa(transpose(r));
    disp(RR);
    N=N+1;
end;
0 @$ N: N( z% T

% i, s9 s- d- [

aplerll

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

tuoni

已解决。

shelby

这个没有见过。