【文章內(nèi)容簡(jiǎn)介】 由各邊的坐標(biāo)方位角和邊長(zhǎng)，按公式、計(jì)算坐標(biāo)增量。 （3） 導(dǎo)線點(diǎn)坐標(biāo)計(jì)算按公式、計(jì)算各導(dǎo)線點(diǎn)坐標(biāo)[3]。 東華理工大學(xué)長(zhǎng)江學(xué)院畢業(yè)設(shè)計(jì) 基于Matlab的導(dǎo)線計(jì)算程序設(shè)計(jì)3. 基于Matlab的導(dǎo)線計(jì)算程序設(shè)計(jì) 計(jì)算程序概述 計(jì)算程序特點(diǎn)相對(duì)于手工計(jì)算，坐標(biāo)計(jì)算程序計(jì)算的主要特點(diǎn)是計(jì)算速度快、精度高、數(shù)據(jù)處理自動(dòng)化，從而把人從繁重的計(jì)算工作中解放出來(lái)。從程序設(shè)計(jì)的角度看，程序設(shè)計(jì)與平差計(jì)算相對(duì)獨(dú)立。在平差手工計(jì)算時(shí)，我們總是面對(duì)需要計(jì)算的具體問(wèn)題，所以其數(shù)據(jù)是特定的，計(jì)算過(guò)程由人實(shí)時(shí)控制；在計(jì)算機(jī)程序計(jì)算中，在程序設(shè)計(jì)時(shí)數(shù)據(jù)是抽象的，必須考慮到實(shí)際計(jì)算中問(wèn)題的多樣性，以及數(shù)據(jù)計(jì)算過(guò)程的自動(dòng)化，所以在程序設(shè)計(jì)時(shí)必須考慮需要處理的所有問(wèn)題的普遍性和規(guī)律性。另外，相對(duì)于手工計(jì)算，在程序計(jì)算時(shí)，選擇平差方法的依據(jù)不同。在手工計(jì)算時(shí)，我們通常希望盡量降低計(jì)算工作量。當(dāng)必要觀測(cè)數(shù)t大于多余觀測(cè)數(shù)r時(shí)，我們可以選擇條件平差；當(dāng)必要觀測(cè)數(shù)t小于多余觀測(cè)數(shù)r時(shí)，我們可以選擇間接平差，這樣，可以降低平差計(jì)算量。然而，在計(jì)算機(jī)程序計(jì)算時(shí)，由于計(jì)算機(jī)計(jì)算的快速高效性，我們不是很關(guān)心計(jì)算量的問(wèn)題，而把主要精力集中于方法實(shí)現(xiàn)的現(xiàn)實(shí)性方面，也就是要求該方法具有較強(qiáng)的規(guī)律性，便于程序設(shè)計(jì)的技術(shù)實(shí)現(xiàn)。在坐標(biāo)計(jì)算程序設(shè)計(jì)中，使用間接平差，對(duì)于一般控制網(wǎng)，誤差方程形式統(tǒng)一、規(guī)律性強(qiáng)、便于程序設(shè)計(jì)；而使用條件平差，誤差方程形式多樣、規(guī)律性差，不利于程序設(shè)計(jì)。所以，在本課程中，我們主要使用間接平差方法進(jìn)行程序設(shè)計(jì)。總之，我們?cè)谶x擇數(shù)學(xué)模型的時(shí)候，一定要考慮算法同計(jì)算機(jī)程序設(shè)計(jì)的特點(diǎn)相統(tǒng)一。 計(jì)算程序的基本要求坐標(biāo)計(jì)算程序設(shè)計(jì)與其它程序設(shè)計(jì)相同，應(yīng)當(dāng)滿足一定的要求：（1）程序邏輯結(jié)構(gòu)簡(jiǎn)單，清晰易讀，符合結(jié)構(gòu)化程序設(shè)計(jì)要求，便于擴(kuò)展；（2）運(yùn)算速度快，占用內(nèi)存小，內(nèi)外存之間的交換不宜過(guò)于頻繁；（3）數(shù)學(xué)模型及計(jì)算方法正確、先進(jìn)，計(jì)算結(jié)果精度高；（4）適應(yīng)性強(qiáng)，便于移植，充分考慮各種可能形式，盡量滿足不同要求與需要。 計(jì)算程序的設(shè)計(jì)步驟(1) 結(jié)構(gòu)總體設(shè)計(jì)；(2) 數(shù)據(jù)結(jié)構(gòu)設(shè)計(jì)；(3) 確定軟件各組成部分的算法及數(shù)據(jù)組織；(4) 選定某種表達(dá)式來(lái)描述各種算法；(5) 程序編寫；(6) 程序調(diào)試。 坐標(biāo)計(jì)算處理過(guò)程（1）數(shù)據(jù)輸入；（2）角度閉合差計(jì)算及分配；（3）方位角推算；（4）坐標(biāo)增量計(jì)算；（5）坐標(biāo)增量閉合差分配；（6）坐標(biāo)計(jì)算。 附合導(dǎo)線計(jì)算程序附合導(dǎo)線包括一個(gè)連接角的附合導(dǎo)線、兩個(gè)連接角的附合導(dǎo)線和無(wú)連接角的附合導(dǎo)線等，依據(jù)本文介紹的附合導(dǎo)線計(jì)算方法和基本公式，依據(jù)matlab語(yǔ)言對(duì)附合導(dǎo)線進(jìn)行編程,主要的程序代碼如下（具體程序代碼見附錄1）： %觀測(cè)數(shù)據(jù)輸入 n=input(39。請(qǐng)輸入導(dǎo)線邊數(shù)39。)。 string1=(char((39。請(qǐng)輸入第39。)39。,a,(39。條導(dǎo)線邊長(zhǎng)39。)39。))39。 string2=(char((39。請(qǐng)輸入第39。)39。,s,(39。個(gè)轉(zhuǎn)折角39。)39。))39。 %起算數(shù)據(jù)輸入 %起算方位角的計(jì)算 %角度閉合差計(jì)算 fb=fwjABfwjCD+sum(jiao1)n*180。 %各觀測(cè)角改正后的角度 jiao2=jiao1+(fb/6)*(1)。 %真方位角的計(jì)算 fwjz(k)=fwjAB+jiao2(k)。 fwjz(k)=fwjz(k1)+jiao2(k)180。 %坐標(biāo)閉合差計(jì)算 xz(i)=bian(i)*cosd(fwjz(i))。 yz(i)=bian(i)*sind(fwjz(i))。 fx=xA+sum(xz)xC。 fy=yA+sum(yz)yC。 %坐標(biāo)閉合差分配 vx(i)=(1)*fx/sum(bian)*bian(i)。 vy(i)=(1)*fy/sum(bian)*bian(i)。 %坐標(biāo)計(jì)算 x(1)=xA+xz(1)+vx(1) y(1)=yA+yz(1)+vy(1)。 x(i)=x(i1)+xz(i)+vx(i)。 y(i)=y(i1)+yz(i)+vy(i)。 閉合導(dǎo)線計(jì)算程序閉合導(dǎo)線就是已知一條邊，測(cè)量若干個(gè)邊長(zhǎng)和夾角后又閉合到已知邊的導(dǎo)線測(cè)量方法。通過(guò)計(jì)算平差后，可計(jì)算得到經(jīng)過(guò)的未知點(diǎn)的平面坐標(biāo)。依據(jù)本章介紹的閉合導(dǎo)線計(jì)算方法和基本公式，依據(jù)matlab語(yǔ)言對(duì)閉合導(dǎo)線進(jìn)行編程，主要的程序代碼如下（具體程序代碼見附錄2）： %觀測(cè)數(shù)據(jù)的輸入 n=input(39。請(qǐng)輸入閉合導(dǎo)線的點(diǎn)數(shù)39。)。 string1=(char((39。請(qǐng)輸入第39。)39。,a,(39。條導(dǎo)線的邊長(zhǎng):39。)39。))39。 string2=(char((39。請(qǐng)輸入第39。)39。,b,(39。個(gè)轉(zhuǎn)折角:39。)39。))39。 %起算數(shù)據(jù)的輸入 string3=(char(39。請(qǐng)輸入起始已知方向的起算方位角39。))。 %角度閉合差的計(jì)算 fb=sum(jiao2)(n2)*180。 %真方位角計(jì)算 fwjz(1)=fwj2+jiao2(1)+v(1)180。 %坐標(biāo)閉合差計(jì)算 xzq=cosd(fwj2)*bian(1)。 yzq=sind(fwj2)*bian(1)。 xz(i)=cosd(fwjz1(i))*bian(i+1)。 yz(i)=sind(fwjz1(i))*bian(i+1)。 fx=xzq+sum(xz)。 fy=yzq+sum(yz)。 fs=sqrt(fx^2+fy^2)。 k=fs/sum(bian)。 %坐標(biāo)計(jì)算 X(1)=x+xzq+vx(1)。 Y(1)=y+yzq+vy(1)。 X(i)=X(i1)+xz(i1)+vx(i)。 Y(i)=Y(i1)+yz(i1)+vy(i)。 支導(dǎo)線計(jì)算程序由已知控制點(diǎn)出發(fā)，不附合、不閉合于任何已知點(diǎn)的導(dǎo)線。依據(jù)本章介紹的支導(dǎo)線計(jì)算方法和基本公式，依據(jù)matlab語(yǔ)言對(duì)支導(dǎo)線進(jìn)行編程，主要的程序代碼如下（具體程序代碼見附錄3）： %觀測(cè)數(shù)據(jù)輸入 n=input(39。please the 待定點(diǎn)個(gè)數(shù):39。)。 s1=char(39。please input the 邊長(zhǎng) of a:39。)。 s2=char(39。please input the 轉(zhuǎn)折角 of a:39。)。 jiao2(i)=degree1+minute1/60+second1/3600。 %已知數(shù)據(jù)輸入 %方位角推算 fwj(j)=f0+jiao2(j)180。 fwj(j)=fwj(j1)+jiao2(j)180。 %坐標(biāo)增量計(jì)算 xzl(k)=S(k)*cosd(fwj(k))。 yzl(k)=S(k)*sind(fwj(k))。 %坐標(biāo)推算 x(l)=x0+xzl(l)。 y(l)=y0+yzl(l)。 x(l)=x(l1)+xzl(l)。 y(l)=y(l1)+yzl(