function?[resultpoint1point2]=bestpath（A1ij）

if?nargin==0
????A?=[nan?nan?5?2;3?nan?nan?1;nan?6?nan?3]????????%用來測試閉合回路的算法
????x=3;y=1;
elseif?nargin==3
????A=A1;x=i;y=j;
end
move=[1?0;-1?0;0?1;0?-1];????????????????????????%指針的移動方向，下、上、右、左
[mn]=size（A）;???????????????????????????????????%m、n分別為A的行數和列數

R=[];
%------------行的可行性標記
[marksijrowcolB]=mark（Axy）;??????????????%初始化可行性路徑
prev=[0?0];?????????????????????????????????????%用于保存上一次移動時的移動方向
[~R~~~]=seekPath（B0marksmoveRprev0ijrowcol）;??%查詢比和路徑

[resultresult1point1point2]=res（ArowcolRmn）;????????????????????????%查詢結果，用二進制表示路徑
function?[isfindRtimesrowcolprev]=seekPath（AisfindmarksmoveRprevtimesijrowcol）
??if?（i==row?&&j==col）???????????????????%判斷是否又回到查找起點
????????times=times+1;
????????if?times==2??????????????????????%此時回到起點，查找結束
????????????isfind=1;
????????????return
????????end
??end
??x=i;y=j;
??if?isnan（A（ij））??????????????%該數據點為空所執行的算法
??????next=prev;????????????????%如果數據格為空則只能直走?????
??????u=i-prev（1）;v=j-prev（2）;??%按照前一次的方向移動
??????if?prev（1）
??????????mrow=2;
??????else
??????????mrow=1;
??????end
??????if?marks（u+1v+1mrow）==1
?????????[isfindRtimesrowcol]=seekPath（AisfindmarksmoveRnexttimesuvrowcol）;
?????????if?isfind==1
??????????????R=[R;i?j];
??????????????return
?????????end
??????end
??else
????????for?k=1:4
????????????if?k==1||k==2
????????????????mrow=2;
????????????else
????????????????mrow=1;
????????????end
????????????u=i+move（k1）;
????????????v=j+move（k2）;
????????????
????????????if?prev（1）~=move（k1）||prev（2）~=move（k2）???????????????%判斷是否原路返回
????????????????if?marks（u+1v+1mrow）==1
????????????????????next=[-move（k1）-move（k2）];
????????????????????[isfindRtimesrowcol]=seekPath（AisfindmarksmoveRnexttimesuvrowcol）;
???????????????????????if?isfind==1
???????????????????????????R=[R;i?j];
???????????????????????????return
???????????????????????end
????????????????end
???????????end
????????end
??end
????function?[resultresult1point1point2]=res（ArowcolRmn）
????????[m1n1]=size（R）;
????????point=[];
????????num=0;
????????
????????A（rowcol）=0;???????????????????%為了以后記錄方便
????????result=zeros（mn）;??????????????%只記錄轉折點的值
????????result1=zeros（mn）;??
????????result2=zeros（mn）;?
?
????????result2（rowcol）=1;
????????
????????for?i=1:m1
????????????if?~isnan（A（R（i1）R（i2）））?
????????????????result2（R（i1）R（i2））=1;
????????????end
????????end
????????
????????%只保留拐角處的結點，
????????for?i=1:m
????????????T1=find（result2（i:）==1）;
????????????if?length（T1）>1
????????????????result（iT1（1））=1;
????????????????result（iT1（length（T1）））=1;
????????????????
????????????????result1（iT1（1））=A（iT1（1））;
????????????????result1（iT1（length（T1）））=A（iT1（length（T1）））;
????????????end
????????end
????????for?j=1:n
????????????T2=find（result2（:j）==1）;
????????????if?length（T2）>1
?????????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????4668??2013-04-23?17:52??matlab?運輸問題表上作業\bestpath.m
?????文件????????8252??2013-04-24?12:45??matlab?運輸問題表上作業\vogel.m
?????目錄???????????0??2013-04-25?21:11??matlab?運輸問題表上作業\