大地測(cè)量礦區(qū)控制網(wǎng)設(shè)計(jì)

1、學(xué)院： 班級(jí)： 姓名： 學(xué)號(hào)： 燈湖礦區(qū)控制網(wǎng)設(shè)計(jì)一、工程控制網(wǎng)設(shè)計(jì)任務(wù)書燈湖礦區(qū)控制網(wǎng)設(shè)計(jì)任務(wù)通知書礦大測(cè)繪隊(duì)：燈湖礦區(qū)位于XX省XX煤田東北域。本礦區(qū)已列入國(guó)家經(jīng)濟(jì)建設(shè)計(jì)劃，準(zhǔn)備進(jìn)行重點(diǎn)開發(fā)。根據(jù)建設(shè)計(jì)劃，首先需施測(cè)80平方公里的1:2000的地圖。所給圖為1:50000地形圖，虛線范圍內(nèi)為1:2000測(cè)圖區(qū)。為了加快建井速度，需在豎井間進(jìn)行對(duì)向掘進(jìn)，最大貫通距離為3km，礦井最大單翼長(zhǎng)度為5km，井深一般不超過(guò)800m。為了滿足礦山設(shè)計(jì)、建井、生產(chǎn)三階段測(cè)繪各種比例尺地形圖、井巷貫通以及工業(yè)場(chǎng)地施工測(cè)量的需要，應(yīng)在全礦區(qū)建立統(tǒng)一的具有足夠精度密度的平面控制網(wǎng)與高程控制網(wǎng)。要求你隊(duì)在任務(wù)通

2、知書下達(dá)后一周內(nèi)完成控制網(wǎng)技術(shù)設(shè)計(jì)，一月內(nèi)完成建網(wǎng)任務(wù)。 XX省煤炭管理局 年 月 日 二、測(cè)區(qū)自然地理?xiàng)l件及已有測(cè)繪成果資料（一）自然地理?xiàng)l件1.地理概況本測(cè)區(qū)中心位置為東經(jīng)11713，北緯3412。測(cè)區(qū)地面高程為+30+244m。位于XX省XX縣境內(nèi)，屬于華東丘陵地帶，南部及中部為山地。其上林木繁茂，不利通視，其余為平原，密布村莊和廠礦，建筑群遍地，影響觀測(cè)。燈湖位于測(cè)區(qū)東部，湖面開闊，中有大堤。2.交通情況測(cè)區(qū)內(nèi)有鐵路支線通過(guò)，公路四通八達(dá)，村間大道可行汽車，交通方便。3.氣候情況測(cè)區(qū)地處華東近海地區(qū)，氣候宜人。全年平均降雨量為1040mm，雨量集中于6、7、8三個(gè)月份。全年平均氣溫為+

3、15，夏季氣溫較高，一般為30左右，冬季有雪，但不寒冷，最低氣溫為-5，凍土線深度為0.1m。全年平均風(fēng)力為23級(jí)，夏季略受臺(tái)風(fēng)的影響。宜于野外作業(yè)時(shí)間為311月份，年平均作業(yè)時(shí)間利用率為21天/月。4.居民及居民點(diǎn)測(cè)區(qū)內(nèi)地少人多，人口稠密，多為漢族，少數(shù)為回族。測(cè)量作業(yè)所需人力、物力、材料及食宿均可就地解決。（二）已有測(cè)繪成果資料1.三角成果及其精度國(guó)家二等網(wǎng)點(diǎn)3個(gè)：點(diǎn)名 等級(jí) x（m） y（m） H（m） 玉 山 二 195.4廣具山 二 135.0太 山 二 244.0此點(diǎn)系1959年有XX測(cè)繪隊(duì)施測(cè)，作業(yè)依據(jù)的規(guī)范為一、二、三、四等三角測(cè)量細(xì)則(1958年)。三點(diǎn)標(biāo)石保存完好。坐標(biāo)系統(tǒng)

4、為1954年北京坐標(biāo)系，三度分帶，中央子午線為117。2.高程網(wǎng)成果及其精度國(guó)家等水準(zhǔn)路線由西向東橫穿測(cè)區(qū)北部。根據(jù)二等水準(zhǔn)路線圖，本測(cè)區(qū)內(nèi)及附近應(yīng)有三個(gè)二等水準(zhǔn)點(diǎn)，依點(diǎn)之記只找到其中的2個(gè)標(biāo)石保存完好：點(diǎn)號(hào) 等級(jí) x（m） y（m） H（m）-4 二35.000-2 二46.000該二等水準(zhǔn)路線系國(guó)家測(cè)繪隊(duì)于1979年施測(cè)。施測(cè)精度及埋石質(zhì)量均符合規(guī)范要求。高程為1956年黃海高程系統(tǒng)。3.地形圖資料1984年XX航測(cè)隊(duì)航測(cè)本地區(qū)1:50000地形圖，成圖效果良好，可供技術(shù)設(shè)計(jì)圖上選點(diǎn)使用。地形圖見(jiàn)成果圖。三、本次設(shè)計(jì)的要求1.首級(jí)測(cè)圖控制網(wǎng)可采用E級(jí)(四等)GPS網(wǎng)，點(diǎn)數(shù)為612個(gè)（含已知

5、點(diǎn)），平均邊長(zhǎng)為25km。采用GPS靜態(tài)相對(duì)定位測(cè)量，觀測(cè)時(shí)間為45分鐘，GPS接收機(jī)用3臺(tái)或4臺(tái)。要求設(shè)計(jì)出幾個(gè)同步時(shí)段構(gòu)成的異步環(huán)。 GPS網(wǎng)觀測(cè)與數(shù)據(jù)處理采用基線解算、無(wú)約束平差、約束平差，經(jīng)坐標(biāo)轉(zhuǎn)換最后得到GPS點(diǎn)的1954坐標(biāo)系坐標(biāo)。GPS點(diǎn)的水準(zhǔn)高程通過(guò)已知水準(zhǔn)高程公共點(diǎn)采用多項(xiàng)式擬合法求得。估算最弱中誤差。對(duì)于1:2000比例尺測(cè)圖，首級(jí)控制網(wǎng)精度估算后，相鄰點(diǎn)最弱點(diǎn)位誤差5cm。2.加密網(wǎng)為一級(jí)（或二級(jí)）導(dǎo)線網(wǎng)。每條導(dǎo)線長(zhǎng)度、平均邊長(zhǎng)導(dǎo)線點(diǎn)數(shù)按相應(yīng)要求等級(jí)布設(shè)（平均500m一點(diǎn)）觀測(cè)方法采用全站儀按規(guī)范要求測(cè)角、測(cè)邊。估算最弱點(diǎn)中誤差。對(duì)于1:2000比例尺測(cè)圖，加密控制網(wǎng)精度

6、估算后，相鄰點(diǎn)最弱點(diǎn)位誤差5cm。3.水準(zhǔn)網(wǎng)可布設(shè)為四等附合或閉合水準(zhǔn)網(wǎng)，平均12km一個(gè)水準(zhǔn)點(diǎn)。估算最弱點(diǎn)高程中誤差，其值應(yīng)2cm。四、本次設(shè)計(jì)成果 1.GPS網(wǎng)首先對(duì)測(cè)區(qū)的地形地貌有一個(gè)大概的了解，根據(jù)上述對(duì)GPS網(wǎng)的要求選擇合適的點(diǎn)，除三個(gè)已知點(diǎn)之外，我還重新選了九個(gè)點(diǎn)，一共是十二個(gè)觀測(cè)點(diǎn)，對(duì)已知點(diǎn)用-1，-2，-3,進(jìn)行標(biāo)記，其他電用1-9標(biāo)記。采用4臺(tái)GPS接收機(jī)的作業(yè)方法，分為6個(gè)時(shí)段，獨(dú)立基線為18條。然后對(duì)GPS網(wǎng)中的獨(dú)立基線進(jìn)行圖上測(cè)量，換算到實(shí)際距離，然后對(duì)每條獨(dú)立基線定權(quán)，設(shè)計(jì)基線的觀測(cè)精度為：10mm+10D(km)根據(jù)此式將每條獨(dú)立基線的權(quán)值確定。如下：基線圖上距離（

7、cm）實(shí)際距離(km)權(quán)值P-1 15.02 3.1375 0.8175 -1 -25.753.593750.6405 1 33.64 2.2751.4355 2 34.582.8650.9643 2 45.393.368750.7179 4 -32.051.281253.3559 3 64.452.781251.0148 3 -34.983.11250.82948 7 -34.492.806250.99889 7 95.403.3750.7155 7 62.801.752.18219 -24.512.818750.991 5 85.583.48750.6735 8 -23.822.38751

8、.3231 6 -26.323.950.53397 5 13.502.18751.5325 5 33.892.431251.2827 5 64.602.8750.9568然后應(yīng)用MATLAB軟件對(duì)其進(jìn)行精度估算，程序如下：B=zeros(36,18);JK=1,99;2,99;1,3;3,2;2,4;4,99;3,6;3,99;7,99;7,9;7,6;9,99;5,8;8,99;6,99;5,1;5,3;5,6;i=1;while 1; if i=19; break; endj=JK(i,1);k=JK(i,2); B(i-1)*2+1,(j-1)*2+1)=1; B(i-1)*2+2,(j

9、-1)*2+2)=1; if k=99; a=1; else B(i-1)*2+1,(k-1)*2+1)=-1; B(i-1)*2+2,(k-1)*2+2)=-1; end i=i+1;endQ=blkdiag(0.8175,0.8175,0.6405,0.6405,1.4355,1.4355,0.9643,0.9643,0.7179,0.7179,3.3559,3.3559,1.0148,1.0148,0.82948,0.82948,0.99889,0.99889,0.7155,0.7155,2.1821,2.1821,0.991,0.991,0.6735,0.6735,1.3231,1.3

10、231,0.53397,0.53397,1.5325,1.5325,1.2827,1.2827,0.9568,0.9568);P=inv(Q);BTPB=B*P*B;Qxx=inv(BTPB);得到的B,Qxx結(jié)果如下:經(jīng)觀察，7號(hào)和8號(hào)點(diǎn)的Qxx最大，為0.85697；其中誤差：7=8=29.7745 *=27.5.56mm50mm因此此次設(shè)計(jì)沒(méi)有超限，符合標(biāo)準(zhǔn)。2. 水準(zhǔn)網(wǎng)水準(zhǔn)網(wǎng)的設(shè)計(jì)和GPS原理大致相同，保持圖上距離在1.6-3.2cm以內(nèi)，設(shè)立的水準(zhǔn)點(diǎn)可以不完全通過(guò)GPS點(diǎn)，只需要經(jīng)過(guò)其中的幾個(gè)即可，以方便進(jìn)行高程擬合。我選擇了29個(gè)水準(zhǔn)點(diǎn)。定權(quán)結(jié)果如下：點(diǎn)號(hào)圖上距離權(quán)值P1,992.

11、15 1.1331,22.101.172,32.690.8313,243.020.69724,42.01.244,232.281.04923,182.101.1718,271621.5727,192.241.07419,992.560.894,52.201.15,61.921.36,71.641.557,252.720.8225,82.480.9348,222.620.8722,213.140.65521,201.781.4220,152.181.11315,992.311.0316,992.560.8916,172.700.8317,182.920.7348,91.801.49,102.18

12、1.11310,112.780.79111,122.221.08712,132.281.04913,143.00.70414,262.181.11326,153.00.7048,282.01.2428,291.821.38529,123.100.669然后方法如上，用MATLAB進(jìn)行解算，B=100000000000000000000000000001-100000000000000000000000000001-10000000000000000000000000000100000000000000000000-100000000-100000000000000000001000000001

13、000000000000000000-100000000000000000000000-10000100000000000000000000000100000000-100000000000000000000-10000000100000000000000000000100000000000001-100000000000000000000000000001-100000000000000000000000000001-10000000000000000000000000000100000000000000000-100000000000-100000000000000001000000000

14、0010000000000000-1000000000000000000000000000-1100000000000000000000000000-110000000000000000000000-10000100000000000000000000000100000000000000000000000000000100000000000000000000000000001-100000000000000000000000000001-10000000000000000001-100000000000000000000000000001-100000000000000000000000000

15、001-100000000000000000000000000001-100000000000000000000000000001-100000000000000000000000000001-10000000000000000000000000000100000000000-100000000000000000-100000000001000000000010000000000000000000-100000000000000000000000000001-100000000000-100000000000000001經(jīng)執(zhí)行后，其得到的Qxx為26*26的矩陣：經(jīng)觀察得：10號(hào)點(diǎn)的Qxx最大

16、，精度計(jì)算如下： =10* =18.04mm20mm所以水準(zhǔn)網(wǎng)設(shè)計(jì)合理，沒(méi)有超出精度要求。3. 導(dǎo)線的精度評(píng)定導(dǎo)線設(shè)計(jì)為一級(jí)加密導(dǎo)線網(wǎng)，在此只舉出其中的一條導(dǎo)線進(jìn)行精度評(píng)定。首先，在兩個(gè)GPS點(diǎn)之間設(shè)計(jì)一條導(dǎo)線，其圖上距離大概在0.8cm，然后用直尺量出兩點(diǎn)之間的距離以及X和Y，確定其權(quán)值，導(dǎo)線網(wǎng)的權(quán)值根據(jù)一下公式定權(quán)：則算出的導(dǎo)線點(diǎn)的近似坐標(biāo)和權(quán)值如下：點(diǎn)號(hào)X坐標(biāo)XY坐標(biāo)Y權(quán)P玉山#1-#2394-307.8.220#2-#3496.3.3-61.220#3-#4404.5.8-293.9.320#4-#54924.286.8.120#5#6449.4.6219.2.320#6-#7483.

17、6129.4.720#7-#8493.8.4-78.720#8-#9312.4.8-383.720#9-#10466.8.6-179.720#10-#11499.9.5-8.720#11-#12428.6.1-257.5.520#12-#13456.8.9-203.4.120#13#14498.8.7-34.9.220#14太山424.7-265.220將數(shù)據(jù)導(dǎo)入MATLAB的程序?yàn)椋哼\(yùn)行程序后所得的Qxx的矩陣為：B=zeros(33,28);JKH=-1,-2,1;1,-1,2;2,1,3;3,2,4;4,3,5;5,4,6;6,5,7;7,6,8;8,7,9;9,8,10;10,9,11

18、;11,10,12;12,11,13;13,12,14;14,13,-3;-3,14,-2;JK=-1,-2;-1,1;1,2;2,3;3,4;4,5;5,6;6,7;7,8;8,9;9,10;10,11;11,12;12,13;13,14;14,-3;-3,-2;XYJS=,;,;,;,;,;,;,;,;,;,;,;,;,;,;,;,;,;ro=2062.65; i=1; while 1 if i=17; break; endj=JKH(i,1);k=JKH(i,2);h=JKH(i,3); if j=0; j=14+abs(j); end if k=0; k=14+abs(k); end

19、if h=0; h=14+abs(h); end dxjk=XYJS(k,1)-XYJS(j,1);dyjk=XYJS(k,2)-XYJS(j,2); dxjh=XYJS(h,1)-XYJS(j,1);dyjh=XYJS(h,2)-XYJS(j,2); sjko=dxjk*dxjk+dyjk*dyjk; sjho=dxjh*dxjh+dyjh*dyjh;j=JKH(i,1);k=JKH(i,2);h=JKH(i,3); if j=0&k=0 A=0; end if j=0 B(i,k*2-1)=-ro*dyjk/sjko;B(i,k*2)=ro*dxjk/sjko; end if k=0 B(

20、i,j*2-1)=ro*dyjk/sjko;B(i,j*2)=-ro*dxjk/sjko; end if k=0&j=0 B(i,j*2-1)=ro*dyjk/sjko;B(i,j*2)=-ro*dxjk/sjko;B(i,k*2-1)=-ro*dyjk/sjko;B(i,k*2)=ro*dxjk/sjko; end if j=0&h=0 A=0; end if j=0 B(i,h*2-1)=ro*dyjh/sjho;B(i,h*2)=-ro*dxjh/sjho; end if h=0 B(i,j*2-1)=B(i,j*2-1)-ro*dyjh/sjho;B(i,j*2)=B(i,j*2)+r

21、o*dxjh/sjho; end if h=0&j=0 B(i,j*2-1)=B(i,j*2-1)-ro*dyjh/sjho;B(i,j*2)=B(i,j*2)+ro*dxjh/sjho;B(i,h*2-1)=ro*dyjh/sjho;B(i,h*2)=-ro*dxjh/sjho; end i=i+1; end while 2 if i=34; break; endj=JK(i-16,1);k=JK(i-16,2); if j=0; j=14+abs(j); end if k=0; k=14+abs(k); end dxjk=XYJS(k,1)-XYJS(j,1);dyjk=XYJS(k,2)

22、-XYJS(j,2); sjko=sqrt(dxjk*dxjk+dyjk*dyjk); j=JK(i-16,1);k=JK(i-16,2); if j=0&k=0 A=0; end if j=0 B(i,k*2-1)=dxjk/sjko;B(i,k*2)=dyjk/sjko; end if k=0 B(i,j*2-1)=-dxjk/sjko;B(i,j*2)=-dyjk/sjko; end if k=0&j=0 B(i,j*2-1)=-dxjk/sjko;B(i,j*2)=-dyjk/sjko;B(i,k*2-1)=dxjk/sjko;B(i,k*2)=dyjk/sjko; end i=i+1

23、; end P=blkdiag(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0.49,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,0.72); BTPB=B*P*B; Qxx=inv(BTPB);經(jīng)觀測(cè)得18號(hào)點(diǎn)位最弱點(diǎn)，誤差為5*0.9464=4.732小于5導(dǎo)線網(wǎng)設(shè)計(jì)合理，沒(méi)有超出精度范圍。圖見(jiàn)最后衛(wèi)星遙感影像圖設(shè)計(jì)一、控制網(wǎng)設(shè)計(jì)任務(wù)書對(duì)選定的一幅衛(wèi)星遙感圖片，面積約為100200平方千米，標(biāo)定一直控制點(diǎn)，設(shè)計(jì)測(cè)繪地形圖（1：5000）工程所需控制點(diǎn)，按課程這幾大綱要求進(jìn)行控制網(wǎng)設(shè)計(jì)。二、本次設(shè)計(jì)要求與提示 根據(jù)已知點(diǎn)確定比

24、例尺，估算面積，確定控制網(wǎng)登記和點(diǎn)數(shù)。如面積為150平方千米，平均邊長(zhǎng)為5km，則GPS點(diǎn)數(shù)約為8點(diǎn)（含已知點(diǎn)）?？稍O(shè)計(jì)首級(jí)為E級(jí)的GPS網(wǎng)，網(wǎng)形設(shè)計(jì)構(gòu)網(wǎng)后進(jìn)行精度估算，制定觀測(cè)綱要。GPS點(diǎn)的水準(zhǔn)高程設(shè)計(jì)方案應(yīng)有聯(lián)測(cè)的二等水準(zhǔn)路線，設(shè)計(jì)好后，進(jìn)行高程精度估算，制定觀測(cè)方案。加密控制點(diǎn)設(shè)計(jì)為一級(jí)導(dǎo)線，去其中2條E級(jí)邊進(jìn)行附和導(dǎo)線設(shè)計(jì)，進(jìn)行精度估算，制定觀測(cè)方案。三、本次設(shè)計(jì)成果 1.GPS網(wǎng)首先對(duì)測(cè)區(qū)有一個(gè)大概的了解，根據(jù)上述對(duì)GPS網(wǎng)的要求選擇合適的點(diǎn)，我選了總共12個(gè)點(diǎn)（含已知的3個(gè)點(diǎn)），分別用1-9對(duì)選定的未知點(diǎn)標(biāo)記，采用4臺(tái)GPS接收機(jī)額作業(yè)方法，分為6個(gè)時(shí)段，獨(dú)立基線為18條。然后對(duì)

25、GPS網(wǎng)中點(diǎn)在Google Earth 上面找到點(diǎn)的位置，把相應(yīng)的經(jīng)緯度轉(zhuǎn)換成國(guó)家80坐標(biāo)，利用X , Y值計(jì)算兩點(diǎn)之間的實(shí)際距離，最后對(duì)每條獨(dú)立基線定權(quán)。設(shè)計(jì)基線的觀測(cè)精度為：10mm+10D(km)根據(jù)此式將每條獨(dú)立基線的權(quán)值確定。然后應(yīng)用MATLAB軟件對(duì)其進(jìn)行精度估算，程序如下：B=zeros(36,18);JK=1,2;2,5;1,4;1,99;3,99;3,4;5,99;5,4;4,8;3,6;6,7;4,7;9,99;9,99;8,99;8,99;8,7;6,99;i=1;while 1; if i=19; break; endj=JK(i,1);k=JK(i,2); B(i-1

26、)*2+1,(j-1)*2+1)=1; B(i-1)*2+2,(j-1)*2+2)=1; if k=99; a=1; else B(i-1)*2+1,(k-1)*2+1)=-1; B(i-1)*2+2,(k-1)*2+2)=-1; end i=i+1;endQ=blkdiag(0.977,0.977,1.177,1.177,0.805,0.805,1.034,1.034,1.052,1.052,0.485,0.485,0.736,0.736,1.147,1.147,0.772,0.772,2.809,2.809,1.936,1.936,1.112,1.112,0.8615,0.8615,0.7

27、904,0.7904,1.493,1.493,0.7679,0.7679,1.8645,1.8645,0.6799,0.6799);P=inv(Q);BTPB=B*P*B;Qxx=inv(BTPB);得到的Qxx結(jié)果如下: 經(jīng)觀察，3號(hào)點(diǎn)的Qxx最大，為0.84098；其中誤差為：5*0.9170=4.585因此此次設(shè)計(jì)沒(méi)有超限，符合標(biāo)準(zhǔn)。2水準(zhǔn)網(wǎng)選取了31個(gè)未知點(diǎn)，利用已知點(diǎn)D027和四個(gè)GPS點(diǎn)來(lái)控制水準(zhǔn)網(wǎng)的精度。應(yīng)用MATLAB軟件對(duì)其進(jìn)行精度估算，程序如下：B=zeros(38,31);JK=1,99;1,2;2,3;3,4;4,5;5,99;7,99;7,8;8,9;9,99;11,

28、99;11,12;12,13;13,18;18,28;28,27;27,99;29,99;29,30;30,31;31,99;13,14;14,15;15,16;16,17;17,99;19,99;19,20;20,10;10,6;6,13;20,21;21,22;22,23;23,24;24,25;25,26;26,27;i=1;while 1; if i=39; break; endj=JK(i,1);k=JK(i,2); B(i,j)=1; if k=99; a=1; else B(i,k)=-1; end i=i+1;endQ=blkdiag(1,1,1,1,1,1,1,1,1,1,1

29、,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1);P=inv(Q);BTPB=B*P*B;Qxx=inv(BTPB);得到的Qxx結(jié)果如下:經(jīng)觀察得：23號(hào)點(diǎn)的Qxx最大為2.3393，精度計(jì)算如下： =10* =15.295mm20mm所以水準(zhǔn)網(wǎng)設(shè)計(jì)合理，沒(méi)有超出精度要求。3導(dǎo)線網(wǎng)選取了10個(gè)點(diǎn)作為加密控制點(diǎn)B=zeros(25,20);JKH=-2,-1,1;1,-2,2;2,1,3;3,2,4;4,3,5;5,4,6;6,5,7;7,6,8;8,7,9；9,8,10；10,9，-3；-3,10，-4;JK=-2,-1;1,-2

30、;2,1;3,2;4,3;5,4;-3,5;-4,-3;-3,-4;6,-3;7,6;8,7;9,8;10,9;11,10;-5,11;-6,-5;XYJS=92750,6438;92625,7188;91825,7188;91250,7937;90937,7188;89750,6438;89062,6125;88375,6125;88375,6937;88062,7500;ro=2062.65; i=1; while 1 if i=13; break; endj=JKH(i,1);k=JKH(i,2);h=JKH(i,3); if j=0; j=10+abs(j); end if k=0;

31、k=10+abs(k); end if h=0; h=10+abs(h); end dxjk=XYJS(k,1)-XYJS(j,1);dyjk=XYJS(k,2)-XYJS(j,2); dxjh=XYJS(h,1)-XYJS(j,1);dyjh=XYJS(h,2)-XYJS(j,2); sjko=dxjk*dxjk+dyjk*dyjk; sjho=dxjh*dxjh+dyjh*dyjh;j=JKH(i,1);k=JKH(i,2);h=JKH(i,3); if j=0&k=0 A=0; end if j=0 B(i,k*2-1)=-ro*dyjk/sjko;B(i,k*2)=ro*dxjk/sjko; end if k=0 B(i,j*2-1)=ro*dyjk/sjko;B(i,j*2)=-ro*dxjk/sjko; end if k=0&j=0 B(i,j*2-1)=ro*dyjk/sjko;B(i,j*2)=-r