【正文】 更好的算法。 對PTV算法的展望PTV匹配算法的研究一直是PTV技術的關鍵研究課題之一，如何優(yōu)化其匹配算法，更快更有效率地進行計算，也是國內(nèi)外學者們共同追求的目標。在本文中，匹配幾率法在研究二維規(guī)則流場時得到的結(jié)果近乎完美。但是除了本文介紹的幾種常見的算法，還存在很多非常新穎的算法，比如Spring Model算法，這個算法的出現(xiàn)為三維PTV技術的發(fā)展打下了基礎，人工神經(jīng)網(wǎng)絡算法由于其低錯誤率和完全自動進行的優(yōu)點，具有良好的發(fā)展前景，但是還處于初步發(fā)展階段，應用于實際生產(chǎn)生活的時機還不成熟。隨著計算機水平的飛速發(fā)展，各種算法的不斷完善，相信在不遠的將來，PTV技術將會日趨完善，更加廣泛地應用于各個領域的研究。參考文獻參考文獻[1] 陳釗，郭永彩，高潮．三維PIV原理及其實現(xiàn)方法 [J] ．實驗流體力學，2006,20（4）：7783．[2] 由長福,祁海鷹，等. 采用PTV技術研究循環(huán)流化床內(nèi)氣固兩相流 [J] . 應用力學學報，2004，21（4）：16．[3] 張洋，王元， [J] . 空氣動力學報，2010，28（3）：250254．[4] 阮曉東,吳鋒,楊化勇， [J] .機械工程學報，1990，78（4）：678689．[5] 周云龍，宋連壯，周紅娟. 基于圖像處理的氣固流化床中的氣泡行為的分析[J] .化工自動化及儀表，2011，38（1）：6064．[6] 張洋,王元， [J] . 西安交通大學學報，2009，43（5）：119123．[7] 靳斌,楊冠玲，何振江，郭永采，謝利利. 一種利用示蹤粒子群體運動特征的PTV方法 [J] . 光學技術，2000，26（1）：1618．[8] 周云龍，李洪偉，范振儒. 基于PTV法對油氣水三相流流場的測定[J]. 化工學報，2008,59（10）：25052510．[9] , on ascending and descending motion statistics of particles in turbulent air flow by PTV Int Symp on Applications of Laser Techniques to Fluid Mechanics, Lisbon, Portugal, July 0912,2012.[10] , . A new twoframe particle tracking algorithm using match in Fluids, 1996, 22:2332.[11] ,. Feature point tracking and trajectory analysis for video imaging in cell biology. Journal of Structural Biology, 2005,151:182195． 附 錄1匹配幾率法代碼：function zx(f,Tm,Tn,Tp)n=300。v=2*f。B1=zeros(n,v)。for i=0001:f[M2,MT]=mpretrack_init(39。C:\MATLAB7\work\39。,2,1,3,5,1,i)。MT(:,1)。MT(:,2)。dd = length(MT(:,1))。for j=1:dd B1(j,2*i1)=MT(j,1)。 B1(j,2*i)=MT(j,2)。endendTx=5。c=1。cmax=f1。while c=cmaxf1=zeros(n,2)。f2=zeros(n,2)。f=zeros(n,4)。 for j=1:n f1(j,1)=B1(j,2*c1)。 f1(j,2)=B1(j,2*c)。 f2(j,1)=B1(j,2*c+1)。 f2(j,2)=B1(j,2*c+2)。 endn1max=5。n2max=10。displx=5*Tm*ones(n,1)。disply=5*Tm*ones(n,1)。p=zeros(n+1,n+1)。q=zeros(n+1,n+1)。T=zeros(1,n)。R=zeros(1,n)。pmax=zeros(n,1)。 for k=1:n if f1(k,1)~=0 for j=1:n if f2(j,1)~=0 d=sqrt((f2(j,1)f1(k,1)).^2+(f2(j,2)f1(k,2)).^2)。 if d=Tm T(1,j)=d。 end end end can=find(T)。 if length(can)n1max [g,is]=sort(T(can))。 can=can(is(1:n1max))。 end ind{k}=can。 ns(k)=length(can)。 if ns(k)=1 p(k,can(:))=1/(ns(k)+1)。 p(k,n+1)=1/(ns(k)+1)。 else p(k,n+1)=1。 end T=zeros(1,n)。 clear can。 endendci=0。cmax=10。 A=。 B=4。 C=2。 nn=zeros(n,1)。 for k=1:n if f1(k,1)~=0 for j=1:n if f1(j,1)~=0amp。amp。f1(j,1)~=f1(k,1) d=sqrt((f1(1:n,1)f1(k,1)).^2+(f1(1:n,2)f1(k,2)).^2)。 if d=Tn R(1,j)=d。 end end end aneb=find(R)。 if length(aneb)n2max [disti,is]=sort(R(aneb))。 aneb=aneb(is(1:n2max))。 end nn(k)=length(aneb)。 while cicmax if ns(k)0 for j=1:ns(k) dispx=f2(ind{k}(j),1)f1(k,1)。 dispy=f2(ind{k}(j),2)f1(k,2)。 if (dispxdisplx(k))^2+(dispydisply(k))^2Tx^2 p(k,ind{k}(j))=p(k,ind{k}(j))*C。 end q(k,ind{k}(j))=0。 if nn(k)~=0 for i=1:nn(k) l=find((f2(:,1)f1(aneb(i),1)dispx).^2+(f2(:,2)f1(aneb(i),2)dispy).^2=Tp^2)。 if length(l)0 x=find(p(aneb(i),l)~=0)。 if length(x)0 x=min(x)。 q(k,ind{k}(j))=q(k,ind{k}(j))+p(aneb(i),l(x))。 end end end end if q(k,ind{k}(j))~=0 p(k,ind{k}(j))=A*p(k,ind{k}(j))+B*q(k,ind{k}(j))。 end endp(k,ind{k}(1:ns(k)))=p(k,ind{k}(1:ns(k)))/(sum(p(k,ind{k}(1:ns(k))))+p(k,n+1))。 p(k,n+1)=p(k,n+1)/(sum(p(k,ind{k}(1:ns(k))))+p(k,n+1))。 end clear aneb dispx dispy 。 R=zeros(1,n)。 ci=ci+1。 end end [z,y]=max(p(k,:))。 if y~=301 f(k,1)=f1(k,1)。 f(k,2)=f1(k,2)。 f(k,3)=f2(y,1)。 f(k,4)=f2(y,2)。 end end figure plot([0,1000],[0,0]) hold on plot(1000,0,39。39。) plot([0,0],[1000,0]) xlabel(39。m39。) ylabel(39。m39。) for i=1:n if f(i,1)~=0amp。amp。f(i,3)~=0 plot(f(i,1),f(i,2),39。r.39。) plot(f(i,3),f(i,4),39。b.39。) line([f(i,1),f(i,3)],[f(i,2),f(i,4)]) end end c=c+1。endend附 錄2英文文獻：