24小时热门版块排行榜

返回列表

ldoop

铁杆木虫 (著名写手)

应助: 7 (幼儿园)
金币: 3485.8
散金: 3151
红花: 16
沙发: 1
帖子: 2225
在线: 285.5小时
虫号: 533053
注册: 2008-03-25
专业: 群落生态学

[求助] matlab 程序问题

各位大虾们帮忙看段程序，每次运行总是出现这样的错误：
Improper assignment with rectangular empty matrix.

Error in improved2 (line 56)
                     c(z)=find(divergentsp==b(z));
此时的值 b=4, c=[]
请问怎么调试，另外如何简化这段程序？万分感谢

源程序如下：
if abundance(species)==1
         if species~=S
               b=find(divergentsp==species);
               if isempty(b)~= 1
                  divergentsp(b)=b(1)-1;
                  divergentsp(b(1))=divergentsp(species);
                  for z=1:length(b)
                     c(z)=find(divergentsp==b(z));
                     divergentsp(divergentsp==b(z))=divergentsp(divergentsp==b(z))-1;
                     if isempty(c(z))~=1
                           for x=1:length(c(z))
                              d(x)=find(divergentsp==c(x));
                              divergentsp(divergentsp==c(x))=divergentsp(divergentsp==c(x))-1;
                              if isempty(d(x))~=1
                                 for y=1:length(d(x))
                                       e(y)=find(divergentsp==d(y));
                                       divergentsp(divergentsp==d(y))=divergentsp(divergentsp==d(y))-1;
                                       if isempty(e(y))~=1
                                          for w=1:length(e(y))
                                             f(w)=find(divergentsp==e(w));
                                             divergentsp(divergentsp==e(w))=divergentsp(divergentsp==e(w))-1;
                                             if isempty(f(w))~=1
                                                   g(u)=find(divergentsp(divergentsp==f(u)));
                                                   divergentsp(divergentsp==f(u))=divergentsp(divergentsp==f(u))-1;
                                                   if isempty(g(u))~=1
                                                      for q=1:length(g(u))
                                                         h(q)=find(divergentsp(divergentsp==g(q)))
                                                         divergentsp(divergentsp==g(q))=divergentsp(divergentsp==g(q))-1;
                                                         if isempty(h(q))~=1
                                                               for l=1:length(h(q))
                                                                  divergentsp(divergentsp==h(l))=divergentsp(divergentsp==h(l))-1;
                                                               end
                                                         end
                                                      end
                                                   end
                                             end
                                          end
                                       end
                                 end
                              end
                           end
                     end
                  end
               end
         end

回复此楼