鉆井布局設(shè)計(jì)(數(shù)學(xué)建模論文)

1、課 程 設(shè)計(jì)（論文）任務(wù)書(shū)專(zhuān)業(yè)班級(jí)：  應(yīng)用數(shù)學(xué)0801 學(xué)生姓名：    指導(dǎo)教師（簽名）：      一、課程設(shè)計(jì)（論文）題目鉆井布局模型二、本次課程設(shè)計(jì)（論文）應(yīng)達(dá)到的目的課程設(shè)計(jì)（論文）是實(shí)踐教學(xué)環(huán)節(jié)的重要組成部分，其目的是通過(guò)課程設(shè)計(jì)（論文）加深學(xué)生對(duì)本課程基本知識(shí)的理解，提高綜合運(yùn)用知識(shí)的能力；掌握本課程的主要內(nèi)容、工程設(shè)計(jì)或撰寫(xiě)論文的步驟和方法；提高制圖能力，學(xué)會(huì)應(yīng)用有關(guān)設(shè)計(jì)資料進(jìn)行設(shè)計(jì)計(jì)算和理論分析的方法，以提高學(xué)生獨(dú)立分析問(wèn)題、解決問(wèn)題的能力，逐步增強(qiáng)實(shí)際工程訓(xùn)練。三、本次課程設(shè)計(jì)（論文）任務(wù)的主要內(nèi)容和要

2、求（包括原始數(shù)據(jù)、技術(shù)參數(shù)、設(shè)計(jì)要求等）              勘探部門(mén)在某地區(qū)找礦。初步勘探時(shí)期已零散地在若干位置上鉆井，取得了地質(zhì)資料。進(jìn)入系統(tǒng)勘探時(shí)期后，要在一個(gè)區(qū)域內(nèi)按縱橫等距的網(wǎng)格點(diǎn)來(lái)布置井位，進(jìn)行“撒網(wǎng)式”全面鉆探。由于鉆一口井的費(fèi)用很高，如果新設(shè)計(jì)的井位與原有井位重合（或相當(dāng)接近），便可利用舊井的地質(zhì)資料，不必打這口新井。因此，應(yīng)該盡量利用舊井，少打新井，以節(jié)約鉆探費(fèi)用。 四、應(yīng)收集的資料及主要參考文獻(xiàn)：    

3、60;       收集資料：網(wǎng)格遍歷求解、窮舉算法、MATLAB編程 等參考文獻(xiàn)：1周承高, 廖園.優(yōu)化方法及應(yīng)用程序設(shè)計(jì)M. 北京: 中國(guó)鐵道出版社, 1989.2 張圣勤編著.MATLAB7.0實(shí)用教程M.北京：機(jī)械工業(yè)出版社，20093劉來(lái)福, 曾文藝.數(shù)學(xué)模型與數(shù)學(xué)建模M. 北京:  北京師范大學(xué)出版社, 1997, 8.4 謝兆鴻編著.數(shù)學(xué)建模技術(shù)M.北京：中國(guó)水利水電出版社，2003五、審核批準(zhǔn)意見(jiàn)教研室主任（簽字）        

4、    目錄摘要 31 問(wèn)題重述 42 模型建設(shè) 53 相關(guān)變量及其說(shuō)明 54 模型的建立 54.1 問(wèn)題1） 54.2 問(wèn)題2） 64.3 鉆井布局算法實(shí)現(xiàn) 74.4 對(duì)問(wèn)題3）的討論 105 參考文獻(xiàn) 126 致謝 137 附錄 14鉆井布局模型摘要勘探部門(mén)在某地區(qū)找礦。要在一個(gè)區(qū)域內(nèi)按縱橫等距的網(wǎng)格點(diǎn)來(lái)布置井位，進(jìn)行“撒網(wǎng)式”全面鉆探。由于鉆一口井的費(fèi)用很高，如果新設(shè)計(jì)的井位與原有井位重合（或相當(dāng)接近），便可利用舊井，不必打新井。因此應(yīng)該盡量利用舊井，少打新井，以節(jié)約鉆探費(fèi)用。本文根據(jù)鉆井布局的實(shí)際情況，對(duì)鉆進(jìn)布局問(wèn)題做了深入的研究。力求找到變化中的不變量。

5、討論了在一個(gè)區(qū)域內(nèi)按縱橫等距的點(diǎn)來(lái)布置井位的問(wèn)題，通過(guò)移動(dòng)網(wǎng)格使舊井利用最多，采用網(wǎng)格遍歷求解的方法得到很好的計(jì)算精度。最后得到問(wèn)題1)的解為：I=2,4,5,10時(shí)的舊井可以利用；經(jīng)分析可知問(wèn)題1）是問(wèn)題2）的特殊形式，故問(wèn)題2)可采用類(lèi)似的方法，求得的解為：I=1, 6, 7, 8, 9, 11的舊井可以利用。關(guān)鍵詞：網(wǎng)格遍歷法，結(jié)點(diǎn)，移動(dòng)，覆蓋1 問(wèn)題重述勘探部門(mén)在某地區(qū)找礦。要在一個(gè)區(qū)域內(nèi)按縱橫等距的網(wǎng)格點(diǎn)來(lái)布置井位，進(jìn)行“撒網(wǎng)式”全面鉆探。由于鉆一口井的費(fèi)用很高，如果新設(shè)計(jì)的井位與原有井位重合（或相當(dāng)接近），便可利用舊井，不必打新井。因此應(yīng)該盡量利用舊井，少打新井，以節(jié)約鉆探費(fèi)用。設(shè)

6、平面上有n個(gè)點(diǎn)pi，其坐標(biāo)為(ai,bi)，I=1，2，n，表示以有的n個(gè)井位。新布置的井位是一個(gè)正方形網(wǎng)格N的所有結(jié)點(diǎn)（所謂“正方形網(wǎng)格”是指每個(gè)格子都是正方形的網(wǎng)格；結(jié)點(diǎn)是指縱線(xiàn)和橫線(xiàn)的交叉點(diǎn)）。假定每個(gè)格子的邊長(zhǎng)都是1單位。整個(gè)網(wǎng)格是可以在平面上任意移動(dòng)的。若一個(gè)已知點(diǎn)pi與某個(gè)網(wǎng)格結(jié)點(diǎn)Xi的距離不超過(guò)給定誤差=(0.05單位)，則認(rèn)為pi處的舊井資料可用，不必在結(jié)點(diǎn)Xi處打新井。為進(jìn)行輔助決策，勘探部門(mén)要求我們研究以下問(wèn)題：1） 假定網(wǎng)格的橫向和縱向是固定的（比如東西南北），并規(guī)定兩點(diǎn)間的距離為其橫向距離（橫坐標(biāo)之差絕對(duì)值）及縱向距離（縱坐標(biāo)之差絕對(duì)值）的最大值。在平面上平行移動(dòng)網(wǎng)格N

7、，使可利用的舊井?dāng)?shù)盡可能大。試提供數(shù)值計(jì)算方法，并對(duì)下面的數(shù)值例子用計(jì)算機(jī)進(jìn)行計(jì)算。2） 在歐氏距離的誤差意義下，考慮網(wǎng)格的橫向和縱向不固定（可以旋轉(zhuǎn)）的情形，給出算法及計(jì)算結(jié)果。3） 如果有n口舊井，給出判定這些井均可利用的條件和算法（你可以任意選定一種距離）。數(shù)值例子 n=12個(gè)點(diǎn)的坐標(biāo)如下表所示：i123456789101112ai0.501.413.003.373.404.724.725.437.578.388.989.50bi2.003.501.503.515.502.006.244.102.014.503.410.802 模型假設(shè)1） 正方形網(wǎng)格N可以覆蓋n個(gè)井位；2) 每個(gè)格子的

8、邊長(zhǎng)為1；3) 舊井位點(diǎn)pi(ai,bi)初始坐標(biāo)位于第一象限，網(wǎng)格N可以任意移動(dòng)；4) 把舊井位點(diǎn)看作質(zhì)點(diǎn)；初始位置網(wǎng)格的某條橫線(xiàn)和縱線(xiàn)與坐標(biāo)軸重合，并記交點(diǎn)Q與原點(diǎn)O重合3 相關(guān)變量及其說(shuō)明1) pi(ai,bi) I=1,2,n：舊井位點(diǎn)的坐標(biāo)；2) Dpixi：pi到結(jié)點(diǎn)xi的距離；3) ：誤差；4) (x1i,y1i),(x2i,y2i)：pi所在方格對(duì)角線(xiàn)頂點(diǎn)坐標(biāo)；（x1i<x2i,y1i<y2i）;5) B：網(wǎng)格每次移動(dòng)的步長(zhǎng) C：每次旋轉(zhuǎn)的角度；6) m,n：橫向與縱向移動(dòng)的次數(shù)（橫向向右，縱向向上）；7) s：旋轉(zhuǎn)的次數(shù)（逆時(shí)針）；8) Q(p,q)：變化后Q點(diǎn)的

9、坐標(biāo)；9) 如果第pi個(gè)舊井可利用；則Ti=1否則Ti=0；10) Nj=：第j個(gè)位置可利用舊井的總數(shù)；如果平移，則j=1,2,mxn;如果旋轉(zhuǎn)，則j=1,2,m x n x s;4 模型的建立4.1 問(wèn)題1）1 計(jì)算移動(dòng)后pi所在方格對(duì)角線(xiàn)頂點(diǎn)坐標(biāo),移動(dòng)后，原對(duì)角線(xiàn)頂點(diǎn) 坐標(biāo)變?yōu)椋▁1i + m.B , y1i + n.B），（x2i + m.B , y2i + n.B）（其中:m.B,n.B<=1）； 如果（x1i+mB>=ai&&y1i+nB>=bi），則(x1i +mB-1,y1i+nB)è(x2i,y2i).即將新坐標(biāo)作為pi所在方格的(x

10、2i,y2i)同理：如果（x1i+mB>ai&&y1i+nB<bi）則(x1i +mB-1,y1i+nB）è(x1i,y1i).如果（x1i+mB<ai&&y1i+nB>bi）則(x1i +mB,y1i+nB-1）è(x1i,y1i).如果（x1i+mB <ai&&y1i+nB<bi）則（x1i +mB,y1i+nB）è(x1i,y1i).2. 求Dpixi.令d1=|x1i-ai|,d2=|ai-x2i|,d3=|y1i-bi|,d4=|bi-y2i|.則Dpixi=mc(d1,

11、d2,d3,d4) mc(di)指di中次小的.3. 判斷誤差并確定pi是否可用 如果Dpixi<=,Ti=1如果Dpixi>,Ti=04. 計(jì)算NjNj= j=1,2,.,m×n5. N=max(Nj),同時(shí)記下m,n值及Q(mB,nB).根據(jù)上述步驟，利用計(jì)算機(jī)計(jì)算（見(jiàn)附錄）可得最后結(jié)果為：可利用點(diǎn)為第2, 4, 5,10 點(diǎn)4.2 問(wèn)題2）同一坐標(biāo)系下，將網(wǎng)格N旋轉(zhuǎn)相當(dāng)于將pi反方向旋轉(zhuǎn)每次旋轉(zhuǎn)的角度c=/r,r=max()所以 c 并且 s*c<=/21. 執(zhí)行問(wèn)題1）中1后，以逆時(shí)針旋轉(zhuǎn)，旋轉(zhuǎn)后點(diǎn)pi的坐標(biāo)變?yōu)椋╝i',bi'）其中ai&#

12、39;=*cos(arctan(bi/ai)+s*c),bi'=*sin(arctan(bi/ai)+s*c)2. 求旋轉(zhuǎn)后pi所在方格的對(duì)角線(xiàn)頂點(diǎn)坐標(biāo)如果（ai'>=x1i&&bi'<=y2i）,則（x1i,y1i）è(x1i,y1i)如果（ai'<x1i&&bi'<y2i）,則（x1i-1,y1i）è(x1i,y1i)如果（ai'<x1i&&bi'>y2i）,則（x1i-1,y1i+1）è(x1i,y1i)3. 求 Dpi

13、xid1'=,d2'=d3'=,d4'=則Dpixi=min(d1',d2',d3',d4')4. 執(zhí)行問(wèn)題1）的3，4，5根據(jù)以上步驟，過(guò)程中取B=c=0.05，用計(jì)算機(jī)編程求得可利用點(diǎn)為6個(gè)，即第1, 6, 7, 8, 9, 11 點(diǎn).4.3 鉆井布局算法實(shí)現(xiàn)：我們研究的問(wèn)題是：已知舊井的位置，確定可以重用舊井的最多的個(gè)數(shù)。假設(shè)有一網(wǎng)格可以在井分布的區(qū)域移動(dòng)，這樣只要舊井與網(wǎng)格的某個(gè)節(jié)點(diǎn)距離<=0.05那么認(rèn)為這口舊井是可以重用的，具體計(jì)算機(jī)實(shí)現(xiàn)時(shí)我們考慮到網(wǎng)格的每個(gè)小方格是全等的，這樣我們?cè)谝苿?dòng)網(wǎng)格時(shí)只需整體上下移動(dòng)1

14、單位即可，另外我們是以一個(gè)參考點(diǎn)Q(開(kāi)始在原點(diǎn))在一單位的小方格中移動(dòng)來(lái)處理的。移動(dòng)步長(zhǎng)我們?cè)O(shè)定為0.01，即共要移動(dòng)100*100次。在第二問(wèn)，考慮旋轉(zhuǎn)的情況下，便于計(jì)算，我們移動(dòng)的是原來(lái)的舊井而不是方格，這樣便于處理，同時(shí)在第一問(wèn)的基礎(chǔ)上每移動(dòng)一次，我們都要對(duì)原來(lái)的舊井點(diǎn)坐標(biāo)逆時(shí)針旋轉(zhuǎn)90度（我們步長(zhǎng)暫設(shè)定為1度，根據(jù)計(jì)算機(jī)實(shí)現(xiàn)情況，當(dāng)步長(zhǎng)為1時(shí)運(yùn)行時(shí)間已達(dá)到15分鐘左右），計(jì)算可以重用的舊井的個(gè)數(shù)。到了第三問(wèn)，我們還是考慮到網(wǎng)格的全等對(duì)稱(chēng)性，考慮將所有的舊井移動(dòng)到近量小的范圍，我們用逆向思維考慮：將網(wǎng)格放在適當(dāng)?shù)奈恢?，如果舊井都可以重用的話(huà)，一定可以將這些井的橫縱坐標(biāo)加上m,n單位（n,

15、m可以小于零）將所有的舊井移動(dòng)到某個(gè)節(jié)點(diǎn)的周?chē)?，并且新位置井位與這個(gè)節(jié)點(diǎn)的距離<=0.05。這樣只要將原來(lái)的節(jié)點(diǎn)移動(dòng)到點(diǎn)（1，1）周?chē)?，做一個(gè)盡量小且能包圍所有新位置井位的圓，如果圓的直徑<=2*0.05那么這些舊井都可以重用。計(jì)算機(jī)實(shí)現(xiàn)流程圖如下：井點(diǎn)坐標(biāo)數(shù)據(jù)值顯示井點(diǎn)坐標(biāo)舊井點(diǎn)數(shù)數(shù)據(jù)值輸入程序入口處理結(jié)果數(shù)據(jù)處理選擇是否需要考慮旋轉(zhuǎn)可進(jìn)行相應(yīng)處理紀(jì)錄未移動(dòng)前pi所在方格的定點(diǎn)坐標(biāo)紀(jì)錄可重用井點(diǎn)變量加1,JNk1:=JNk1+1用MD與0.05比較，若MD<0.05認(rèn)為此舊井可以重用函數(shù)MinD具體實(shí)現(xiàn)過(guò)程count具體實(shí)現(xiàn)求出pi到四個(gè)頂點(diǎn)的最短距離MD求出網(wǎng)格移動(dòng)后pi

16、I所在方格的定點(diǎn)坐標(biāo)x1,y1,x2,y2紀(jì)錄s,t的值s,t循環(huán)10000 返回循環(huán)s,t繼續(xù)移動(dòng)網(wǎng)格的位置結(jié)果輸出：沒(méi)移動(dòng)網(wǎng)格一次時(shí)，可重用舊井的個(gè)數(shù)及相應(yīng)網(wǎng)格的位置（Q點(diǎn)的坐標(biāo)(s*0.01,t*0.01)）以上是在沒(méi)有考慮旋轉(zhuǎn)時(shí)的流程圖，對(duì)于考慮旋轉(zhuǎn)時(shí)的情況，只需改變一下count過(guò)程和MinD函數(shù)即可，還需要增加一層循用來(lái)考慮旋轉(zhuǎn)90度的情況，解決解決方案與不考慮旋轉(zhuǎn)時(shí)大同小異。 用程序解決問(wèn)題一，我們將運(yùn)行結(jié)果存到文件jieguo.xls中經(jīng)分析知在我們的精確度范圍內(nèi)最多可以重用井的個(gè)數(shù)是4個(gè)，相應(yīng)網(wǎng)格的位置是：祥見(jiàn)jieguo.xls。在算法優(yōu)化的基礎(chǔ)上我們還可以對(duì)這兩個(gè)點(diǎn)的周?chē)?/p>

17、區(qū)域在進(jìn)行循環(huán)查找看看是否有更多的舊井可以重用。 最后我們分析一下我們的算法復(fù)雜度，這對(duì)于優(yōu)化算法，提高程序的執(zhí)行效率是非常重要的。解決問(wèn)題一時(shí)總的運(yùn)行次數(shù)大約是n*10000次，若考慮旋轉(zhuǎn)時(shí)運(yùn)行次數(shù)更是龐大的驚人：n*900000，因此程序執(zhí)行時(shí)間是比較長(zhǎng)的，大約需要8s左右，分析輸出的數(shù)據(jù)知，很多結(jié)果是沒(méi)用的或者不需要程序循環(huán)，因此，我們可以這樣優(yōu)化程序：開(kāi)始的時(shí)候我們的循環(huán)步值適當(dāng)大些，這樣可以粗略的找到合適的網(wǎng)格位置，我們?cè)僭谶@個(gè)位置上進(jìn)行細(xì)化最終得到最終最優(yōu)結(jié)果，這樣處理我們的程序循環(huán)次數(shù)會(huì)減少到至少0.1倍，同樣在考慮旋轉(zhuǎn)的情況下 我們也作這樣的處理以便優(yōu)化程序、提高運(yùn)行效率、減少

18、運(yùn)行時(shí)間。4.4 對(duì)問(wèn)題 3）的討論1 根據(jù)題意，n口舊井均可利用的條件是由n口舊井的位置決定的，也就是說(shuō)，n口舊井在滿(mǎn)足一定位置關(guān)系的條件下都是可以重用的。我們要計(jì)算的就是這種關(guān)系。2. 如果舊井pi(ai,bi)可利用，則pi'(ai+m,bi+n)也可利用。由此，可把所有的舊井位點(diǎn)移到某個(gè)點(diǎn)的周?chē)?。令d=max(pi pj)（問(wèn)題1）中的距離）。 pi pj:表示移動(dòng)后舊井位點(diǎn)任意兩點(diǎn)的距離，如果d<=2 則滿(mǎn)足條件，即n口舊井均可利用。3. 對(duì)具體坐標(biāo)取值得討論。我們考慮將所有的舊井點(diǎn)移動(dòng)到點(diǎn)（1，1）周?chē)，F(xiàn)假設(shè)舊井點(diǎn)坐標(biāo)的小數(shù)部分是(xi,yi)，現(xiàn)在分4種情況對(duì)舊井

19、點(diǎn)坐標(biāo)移動(dòng)到點(diǎn)（1，1）的周?chē)笮伦鴺?biāo)的值進(jìn)行討論。(1)、xi>=0.5并且yi>=0.5 那么新坐標(biāo)是(xi,yi)è(xi,yi)；(2)、xi<0.5并且yi>=0.5 那么新坐標(biāo)是(xi+1,yi)è(xi,yi)；(3)、xi>=0.5并且yi<0.5 那么新坐標(biāo)是(xi,yi+1)è(xi,yi)；(4)、xi<0.5并且yi<0.5 那么新坐標(biāo)是(xi+1,yi+1)è(xi,yi)；4. 接下來(lái)要求pipj的最大距離d（問(wèn)題1）中的距離）。令U1=min(xi),U2=max(xi);V1

20、=min(yi),V2=max(yi);（I=1,2,n）則d=max(U2-U1,V2-V1)由此可得所有舊井均可利用的條件：d<=2。5 參考文獻(xiàn)1周承高, 廖園.優(yōu)化方法及應(yīng)用程序設(shè)計(jì)M. 北京: 中國(guó)鐵道出版社, 1989.2 3劉來(lái)福, 曾文藝.數(shù)學(xué)模型與數(shù)學(xué)建模M. 北京:  北京師范大學(xué)出版社, 1997, 8.4 謝兆鴻編著.數(shù)學(xué)建模技術(shù)M.北京：中國(guó)水利水電出版社，20036 致謝本篇論文的創(chuàng)作，離不開(kāi)老師耐心而正確的指導(dǎo)，離不開(kāi)各類(lèi)參考資料給予的啟發(fā)和幫助。在此，我表示鄭重的感謝。同時(shí)對(duì)在本篇論文創(chuàng)作過(guò)程中給予過(guò)我無(wú)私幫助的學(xué)長(zhǎng)和各位同學(xué)表示誠(chéng)摯的感謝。另外

21、，由于時(shí)間緊迫和水平有限，本篇論文難免有所錯(cuò)誤和疏漏，請(qǐng)您批評(píng)指正，不吝賜教。7 附錄：7.1 問(wèn)題1)的源程序：（1）、自定義過(guò)程及函數(shù) function MinD(k:integer): real;var secondD:array1.4 of real; i,j:integer; L:real;/做交換中間變量begin/求的Pi點(diǎn)與所在方框四個(gè)點(diǎn)即距離 SecondD1:=abs(x1k-Pk1);SecondD2:=abs(y1k-Pk2); SecondD3:=abs(x2k-pk1);SecondD4:=abs(y2k-Pk2); for i:=1 to 3 do for j:=

22、i+1 to 4 do begin if SecondDi<SecondDj then begin L:=SecondDi; SecondDi:=SecondDj; SecondDj:=L; end; end; result:=SecondD3;end;function MinDK(k:integer): real;var secondD:array1.4 of real;/存儲(chǔ)四個(gè)點(diǎn)的即距離 i,j:integer; L:real;begin /求的Pi點(diǎn)與所在方框四個(gè)點(diǎn)即距離SecondD1:=sqrt(Pk1-x1k)*(Pk1-x1k)+(Pk2-y1k)*(Pk2-y1k);S

23、econdD2:=sqrt(Pk1-x1k)*(Pk1-x1k)+(Pk2-y2k)*(Pk2-y2k);SecondD3:=sqrt(Pk1-x2k)*(Pk1-x2k)+(Pk2-y1k)*(Pk2-y1k);SecondD4:=sqrt(Pk1-x2k)*(Pk1-x2k)+(Pk2-y2k)*(Pk2-y2k);/對(duì)距離進(jìn)行排序，同時(shí)返回最小值 for i:=1 to 3 do for j:= 1 to 4 do begin if SecondDi<SecondDj then begin l:=SecondDi; SecondDi:=SecondDj; SecondDj:=L;

24、end; end; result:=SecondD3;end;procedure TForm1.count(a:real;b:real;x:real;y:real;k:integer);varboo1,boo2:boolean;begin if x>a then boo1:=true else boo1:=false; if y>b then Boo2:=true else boo2:=false; / 求出坐標(biāo) if boo1 and boo2 then begin x2k:=x1k; y2k:=y1k; x1k:=x1k-1;y1k:=y1k-1; end; if boo1 a

25、nd not boo2 then begin x2k:=x1k-1; y2k:=y1k+1; end; if not boo1 and boo2 then begin x2k:=x1k+1; y2k:=y1k-1; end;end;procedure TForm1.countK(a:real;b:real;x:real;y:real;k:integer);varboo1,boo2:boolean;begin boo1:=a<x;boo2:=b>y; / 求出坐標(biāo) if boo1 and boo2 then begin end; if not boo1 and boo2 then b

26、egin x2k:=x1k+1; y2k:=y1k+1; end; if not boo1 and not boo2 then begin x2k:=x1k+1; y2k:=y1k-1; end; if boo1 and not boo2 then begin x2k:=x1k-1; y2k:=y1k-1; end;end;（2）、變量定義部分 var Form1: TForm1; P,P1:array of ar;/存儲(chǔ)舊井的坐標(biāo) JN:array1.10000 of ar1; /存儲(chǔ)結(jié)果 x1,y1,x2,y2:array of real; /存儲(chǔ)某個(gè)P點(diǎn)所在的方格對(duì)角線(xiàn)的倆頂點(diǎn)的坐標(biāo) J

27、W:array1.10000 of string; /不考慮旋轉(zhuǎn)的情況結(jié)果 i:integer; n:integer; /舊井的個(gè)數(shù) s,t:integer;/s,t網(wǎng)格移動(dòng)部數(shù)（3）、程序?qū)崿F(xiàn)部分 var k:integer; MD:real;/存儲(chǔ)最短距離 mess:string;begin for s:=0 to 99 do for t:=0 to 99 do begin for k:=1 to n do begin x1k:=int(Pk1); y1k:=int(Pk2); x2k:=x1k+1+0.01*s; y2k:=y1k+1+0.01*t; x1k:=x1k+0.01*s; y

28、1k:=y1k+0.01*t; count(Pk1,Pk2,x1k,y1k,k);/求得i所在的方格的對(duì)郊縣的兩頂點(diǎn)的坐標(biāo) MD:=MinD(k); if MD<0.05 then begin JNs*100+t+11:=JNs*100+t+11+1;/如果最小距離小與0.05 那么可用井的個(gè)數(shù)JNs*t+11加1； JWs*100+t+1:=JWs*100+t+1+' '+inttostr(k); end; end; JNs*100+t+12:=s; JNs*100+t+13:=t; mess:=floattostr(JNs*100+t+11)+'個(gè)'+char(9)+floattostr(JNs*100+t+12*0.01)+char(9)+floattostr(JNs*100+t+13*0.01)+char(9)+JWs*100+t+1; end;7.2 問(wèn)題2)的源程序： var k,o:integer;/k,循環(huán)n次計(jì)算n各界點(diǎn)；o 循環(huán)90次 MDk:real;/存儲(chǔ)最短距離 mess:string; f_x,f_y,bi:real;beginfor k:=1 to n do for o:=1 to 2 do P1ko:=Pko; /初始化P1k 用于初始化 for s:=0 to 99 do