matlab编程方面的一个疑问

wanggongming · 发表于 2005-10-25 22:45:44

function [src_x,src_y,des_x,des_y,dis,no]=fly_control(px,py,n) tx=ones(1,n); %%¶¨Òå¶þÎ¬¿Õ¼äµÄÁ½¸ö»®·Öp1ºÍp2 ty=ones(1,n); tx1=ones(1,n); ty1=ones(1,n); no=-1; for i=1:n %%p1ºÍp2×ø±ê³õÊ¼»¯ tx(i)=Inf; ty(i)=Inf; tx1(i)=Inf; ty1(i)=Inf; end if (n==1) dis=inf; %Ö»ÓÐÒ»¸öµãÊ±,¾àÀëÎªÎÞÇî´ó src_x=px(1); %¼ÇÂ¼¸ÃµãµÄ×ø±ê src_y=py(1); no=1; return; %±êÖ¾,±íÃ÷¸Ã¿Õ¼äÖÐÖ»ÓÐÒ»¸öµã elseif (n==2) %ÓÐÁ½¸öµãÊ±,¾àÀëÎª¶þÕßÖ®¼äµÄ¾àÀë dis=sqrt((px(1)-px(2))^2+(py(1)-py(2))^2); src_x=px(1); %¼ÇÂ¼ÕâÁ½¸öµãµÄ×ø±ê src_y=py(1); des_x=px(2); des_y=py(2); no=2; %±êÖ¾,±íÃ÷¸Ã¿Õ¼äÖÐÖ»ÓÐÁ½¸öµã return; else %¶à¸öµãÊ±µÄÌØÊâÇé¿ö m=mean(px); %ÇóÖÐÎ»µã n1=0; %»®·ÖÁ½¸ö¿Õ¼äp1ºÍp2 n2=0; for i=1:n if (px(i)<=m) n1=n1+1; tx1(n1)=px(i); ty1(n1)=py(i); else n2=n2+1; tx2(n2)=px(i); ty2(n2)=py(i); end; end; end; %end if src_x1=0; src_y1=0; des_x1=0; des_y1=0; dis1=0; no1=0; %µÝ¹éÔÚ·ÖÀí³öµÄÁ½¸ö¿Õ¼äÖÐÀûÓÃ·ÖÖÎ·¨Çó½â [src_x1,src_y1,des_x1,des_y1,dis1,no1]=fly_control(tx1,ty1,n1); [src_x2,src_y2,des_x2,des_y2,dis2,no2]=fly_control(tx2,ty2,n2); if (dis1<dis2) %ÕÒµ½Á½¸ö¿Õ¼ä½âÖÐµÄ×îÐ¡Õß distemp=dis1; %¾àÀë src_xtemp=src_x1; %×ø±ê src_ytemp=src_y1; des_xtemp=des_x1; des_ytemp=des_y1; else distemp=dis2; %¾àÀë src_xtemp=src_x2; %×ø±ê src_ytemp=src_y2; des_xtemp=des_x2; des_ytemp=des_y2; end %ºÏ²¢Ç°½á¹û³õÊ¼»¯ discombine=Inf; %¾àÀë³õÊ¼»¯ src_xcombine=Inf; %×ø±ê³õÊ¼»¯ src_ycombine=Inf; des_xcombine=Inf; des_ycombine=Inf; %ºÏ²¢ for i=1:n1 if (i<n2) %×öÒ»ÏÂÅÐ¶Ï,·ÀÖ¹p2ÖÐµÄ×ø±êÔ½½ç k=i; else k=n2; end while_flag=1; while(abs(ty2(k)-ty1(i))>distemp) if ((ty2(k)-ty1(i))>0) %%È·¶¨É¨Ãè¿ªÊ¼µã k=k-1; else k=k+1; end if ((k<1)||(k>n2)) %%ÏÂ±êÔ½½çÔòÍË³ö while_flag=0; break; end if ((ty2(k)-ty1(i))>distemp)&&((ty1(i)-ty2(k-1))>distemp) %%·ÀÖ¹É¨ÃèÖ¸Õë×÷ÎÞÇîÕñµ´ while_flag=0; break; end end if (while_flag==0) %%¶Ô¸ø¶¨µÄp1(i)¶øÑÔ,p2ÖÐÃ»ÓÐÈÎºÎµãºÍËü¿ÉÄÜ¹¹³É×î½Ó½üµã¶Ô,½áÊø±¾´ÎÑ»·,¿¼²ìµãp1(i+1) continue; end %¿ªÊ¼ºÏ²¢²éÕÒ while(abs(ty2(k)-ty1(i))<=distemp) %É¨ÃèÖ¸ÕëÔ¼Êø distance=sqrt((tx1(i)-tx2(k))^2+(ty1(i)-ty2(k))^2); %Çó¾àÀë if (distance<discombine) %ÕÒµ½Ò»¸ö¸üÐ¡µÄ¾àÀë discombine=distance; %¾àÀë src_xcombine=tx1(i); %×ø±ê src_ycombine=ty1(i); des_xcombine=tx2(k); des_ycombine=ty2(k); end end end %%end for if (discombine<distemp) dis=discombine; src_x=src_xcombine; src_y=src_ycombine; des_x=des_xcombine; des_y=des_ycombine; else dis=distemp; src_x=src_xtemp; src_y=src_ytemp; des_x=des_xtemp; des_y=des_ytemp; end

wanggongming · 发表于 2005-10-25 22:52:38

<

>上述程序是我用分治法求解空中飞行管理问题的源程序，从原理上说是没有什么问题的，但是好像调试通不过！！我发现是matlab6好像不支持双重递归（一个程序中有两个递归函数），我这个函数的函数体是 function [src_x,src_y,des_x,des_y,dis,no]=fly_control(px,py,n)，里面用了两个递归[src_x1,src_y1,des_x1,des_y1,dis1,no1]=fly_control(tx1,ty1,n1); 和 [src_x2,src_y2,des_x2,des_y2,dis2,no2]=fly_control(tx2,ty2,n2); 如果只有一个递归体，函数可以执行（我用只含一个递归体的排序和阶乘函数试过，没有什么问题），但是为什么有多个递归体就不可以了呢？？？ 希望高手给予指导！！！！！！

wanggongming · 发表于 2005-10-25 22:53:48

不知道什么原因，我的代码中的汉字注释不能显示出来！！！！

		自动登录	找回密码
密码			注-册-帐-号