请教高手

yanyun2929

for N1=1:10;
    s(N1)=0;
end     
a1=0.01;a2=0.02;a3=0.04;
r1=0.43;r2=0.5;r3=0.5;
N2=4;N3=1;
for  N1=1:10;
      
        for T=1:5;
            n1=binornd(N1,r1);
            n2=binornd(N2,r2);
            n3=binornd(N3,r3);
       if n1==0,n2==0,n3==0;s=0;
      elseif n1>0,n2==0,n3==0;
            if a1*r1*N1<a1*n1;b1=a1*r1*N1;
            else b1=0;end
            if b1~=0;
          d1=(a1*r1*N1)/b1;wlr1=max(r1*a1*N1,a1*n1)-a1*n1;
          wlr=wlr1+1-0.01-(a1*r1*N1);t1=(wlr*d1+a1*r1*N1)/a1;
          nt1=binornd(n1,p1);
          else t1=r1*N1;p1=t1/n1;nt1=binornd(n1,p1);end
         if  nt1*a1<1-0.01;
           s=nt1;end
           
      elseif n1>0,n2>0,n3==0;
         if r1*a1*N1<a1*n1;
             b1=r1*a1*N1;
            else b1=0;end
         if  r2*a2*N2<a2*n2;
               b2=r2*a2*N2;
           else b2=0;end
         if  (b1+b2)~=0;
             d1=(a1*r1*N1)/(b1+b2);d2=(a2*r2*N2)/(b1+b2);
             wlr1=max(r1*a1*N1,a1*n1)-a1*n1;wlr2=max(r2*a2*N2,a2*n2)-a2*n2;
             wlr=wlr1+wlr2+1-0.01-(a1*r1*N1+a2*r2*N2);
             t1=(wlr*d1+a1*r1*N1)/a1;t2=(wlr*d2+a2*r2*N2)/a2;
             p1=t1/n1;p2=t2/n2;
             nt1=binornd(n1,p1);nt2=binornd(n2,p2);
         else t1=r1*N1;t2=r2*N2;
             p1=t1/n1;p2=t2/n2;
             nt1=binornd(n1,p1);nt2=binornd(n2,p2);end
          if nt1*a1+nt2*nt2<1-0.01;
              s=nt1+nt2;end
              
          elseif n1==0,n2>0,n3>0;
                 if r2*a2*N2<a2*n2;
                    b2=r2*a2*N2;
                 else b2=0;end
                 if r3*a3*N3<a3*n3;
                     b3=r3*a3*N3
                   else b3=0;end
                  if b2+b3~=0;
                      d2=(a2*r2*N2)/(b2+b3);d3=(a3*r3*N3)/(b2+b3);
                      wlr2=max(r2*a2*N2,a2*n2)-a2*n3;wlr3=max(r3*a3*N3,a3*n3)-a3*n3;
                      wlr=wlr2+wlr3+1-0.01-(a2*r2*N2+a3*r3*N3);      
                      t2=(wlr*d2+a2*r2*N2)/a2;t3=(wlr*d3+a3*r3*N3)/a3;
                      p2=t2/n2;p3=t3/n3;
                      nt2=binornd(n2,p2);nt3=binornd(n3,p3);
                  else    t2=r2*N2;t3=r3*N3;p2=t2/n2;p3=t3/n3;nt2=binornd      (n2,p2);nt3=binornd(n3,p3);
                  if    nt2*a2+nt3*a3<1-0.01;
                        s=nt2+nt3;end
  matlab说我这里的elseif违法使用?? elseif n1==0,n2==0,n3>0;
                   if r3*a3*N3<a3*n3;
                     b3=r3*a3*N3;
                    else b3=0;end
                   if b3~=0;
                      d3=(a3*r3*N3)/b3;
                      wlr3=max(r3*a3*N3,a3*n3)-a3*n3;wlr=wlr3+1-0.01-(a3*r3*N3);
                      t3=(wlr*d3+a3*r3*N3)/a3;p3=t3/n3; nt3=binornd(n3,p3);
                  else   t3=r3*N3; p3=t3/n3;  nt3=binornd(n3,p3);
                  if   nt3*a3<1-0.01;
                      s=nt3;end
                     
               elseif n1>0,n2==0,n3>0;
                      if r1*a1*N1<a1*n1;
                      b1=r1*a1*N1;
                      else b1=0;end
                      if r3*a3*N3<a3*n3;
                     b3=r3*a3*N3;
                    else b3=0;end
                    if (b1+b3)~=0;
                        d1=(a1*r1*N1)/(b1+b3);d3=(a3*r3*N3)/(b1+b3);wlr1=max(r1*a1*N1,a1*n1)-a1*n1;wlr3=max(r3*a3*N3,a3*n3)-a3*n3; wlr=wlr1+wlr3+1-0.01-(a1*r1*N1+a3*r3*N3);t1=(wlr*d1+a1*r1*N1)/a1;t3=(wlr*d3+a3*r3*N3)/a3;p1=t1/n1;p3=t3/n3;nt1=binornd(n1,p1);nt3=binornd(n3,p3);
                    else  t1=r1*N1;t3=r3*N3;p1=t1/n1;p3=t3/n3;nt1=binornd(n1,p1); nt3=binornd(n3,p3);
                   if   nt1*a1+nt3*a3<1-0.01;
                       s=nt1+nt3;end
                       
                elseif n1==0,n2>0,n3==0;
                       if  r2*a2*N2<a2*n2;
                         b2=r2*a2*N2;
                       else b2=0;end
                       if b2~=0;
                           d2=(a2*r2*N2)/b2;wlr2=max(r2*a2*N2,a2*n2)-a2*n2; wlr=wlr2+1-0.01-(a2*r2*N2);t2=(wlr*d2+a2*r2*N2)/a2;p2=t2/n2;nt2=binornd(n2,p2);
                       else   t2=r2*N2;p2=t2/n2; nt2=binornd(n2,p2);
                           if nt2*a2<1-0.01;
                               s=nt2;end
                              
                       else  if r1*a1*N1<a1*n1;
                                   b1=r1*a1*N1;
                              else b1=0;end
                              if  r2*a2*N2<a2*n2;
                                  b2=r2*a2*N2;
                             else b2=0;end
                           if r3*a3*N3<a3*n3;
                               b3=r3*a3*N3;
                            else b3=0;end
                            if (b1+b2+b3)~=0;
                                d1=(a1*r1*N1)/(b1+b2+b3);
                               d2=(a2*r2*N2)/(b1+b2+b3);
                               d3=(a3*r3*N3)/(b1+b2+b3);
                            wlr1=max(r1*a1*N1,a1*n1)-a1*n1;
                           wlr2=max(r2*a2*N2,a2*n2)-a2*n2;
                           wlr3=max(r3*a3*N3,a3*n3)-a3*n3;
                           wlr=wlr1+wlr2+wlr3+1-0.01-(a1*r1*N1+a2*r2*N2+a3*r3*N3);
                          t1=(wlr*d1+a1*r1*N1)/a1;
                          t2=(wlr*d2+a2*r2*N2)/a2;
                         t3=(wlr*d3+a3*r3*N3)/a3;
                        p1=t1/n1;p2=t2/n2;p3=t3/n3;
                           nt1=binornd(n1,p1);
                           nt2=binornd(n2,p2);
                         nt3=binornd(n3,p3);
                        else t1=r1*N1;t2=r2*N2;t3=r3*N3;
                      t3=r3*N3;p2=t2/n2;p3=t3/n3;
                          nt1=binornd(n1,p1);
                          nt2=binornd(n2,p2);
                            nt3=binornd(n3,p3);
                                end
                     if nt1*a1+nt2*a2+nt3*a3<1-0.01;
                         s=nt1+nt2+nt3;
                          end
                      end
                      s(N1)=s(N1)+s;
                  end
                  s(N1)= s(N1)/5;
              end
                 plot(s);
             end
            
                     
                        

redlobe

当然是啦
正确的用法是：if
              elseif 表达式
              elseif 表达式
                .
                .
              else 表达式
              end