【正文】 謝在這篇論文完成之際，我衷心地感謝我的指導(dǎo)老師，李秀海老師在我大學(xué)的最后學(xué)習(xí)階段——畢業(yè)設(shè)計(jì)階段給自己的指導(dǎo)，李老師對(duì)我們認(rèn)真負(fù)責(zé)，及時(shí)傳遞相關(guān)信息，并在我們做課題期間給予了極大的幫助和支持。從最初的選題，到收集資料，到課題研究、編程序、撰寫、修改，到論文最終定稿，他給了我細(xì)心的指導(dǎo)和幫助。李老師用自己的工作和休息時(shí)間，指導(dǎo)我們的畢業(yè)論文，他的這種無私奉獻(xiàn)的敬業(yè)精神令人欽佩，在此我表示我誠(chéng)摯的謝意。同時(shí)，感謝所有任課老師，是他們耐心而又專業(yè)的講授，使我獲得豐富的基礎(chǔ)理論知識(shí)和一定的項(xiàng)目實(shí)踐經(jīng)驗(yàn)，為順利完成畢業(yè)設(shè)計(jì)打下了堅(jiān)實(shí)的基礎(chǔ)。我還要感謝我的各位同學(xué)，很榮幸和他們成為同學(xué)和朋友，大家在一起互相幫助，互相鼓勵(lì)，共同進(jìn)步，這段美好的學(xué)習(xí)和生活經(jīng)歷將成為我一生最珍貴的回憶。在此向他們表示我由衷的謝意，并祝所有的老師培養(yǎng)出越來越多的優(yōu)秀人才，桃李滿天下，祝愿各位同學(xué)前程似錦！感謝在寫論文階段給予我莫大幫助的同學(xué)。謝謝他們?cè)谖覍懻撐钠陂g陪伴我走過的那段充實(shí)而快樂的時(shí)光。感謝黑龍江工程學(xué)院的領(lǐng)導(dǎo)，老師對(duì)我的支持與幫助!感謝我的父母，他們的鼓勵(lì)與支持是我前進(jìn)的動(dòng)力!感謝所有關(guān)心，幫助我的老師同學(xué)和朋友！附 錄附錄A當(dāng)進(jìn)行一次觀測(cè)時(shí)的程序：X=[ ]Y=[ ]disp(39。當(dāng)觀測(cè)一次時(shí)的權(quán)為39。)P=[ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ]l=[ ]for j=1 for k=2:6 DX(j,k)=X(k)X(j)。 DY(j,k)=Y(k)Y(j)。 S(j,k)=sqrt((X(j)X(k))^2+(Y(j)Y(k))^2) endendfor j=2 for k=3:6 DX(j,k)=X(k)X(j)。 DY(j,k)=Y(k)Y(j)。 S(j,k)=sqrt((X(j)X(k))^2+(Y(j)Y(k))^2) endendfor i=1 j=input(39。j=39。)。k=input(39。k=39。) a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k) B1=[a b a b]endfor i=2 j=input(39。j=39。)。k=input(39。k=39。) a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。 B2=[a b 0 0]endfor i=3 j=input(39。j=39。)。k=input(39。k=39。) a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。 B3=[a b 0 0]endfor i=4 j=input(39。j=39。)。k=input(39。k=39。) a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。 B4=[a b 0 0]endfor i=5 j=input(39。j=39。)。k=input(39。k=39。) a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。 B5=[0 0 a b]endfor i=6 j=input(39。j=39。)。k=input(39。k=39。) a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。 B6=[0 0 a b]endfor i=7 j=input(39。j=39。)。k=input(39。k=39。) a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。 B7=[0 0 a b]endfor i=8 j=input(39。j=39。)。k=input(39。k=39。) a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。 B8=[0 0 a b]enddisp(39。誤差方程系數(shù)39。)B=[B1。B2。B3。B4。B5。B6。B7。B8]disp(39。法方程系數(shù)矩陣39。)Nbb=B39。*P*Bdisp(39。法方程逆矩陣39。)A=pinv(Nbb)disp(39。精度估算39。)disp(39。計(jì)算邊長(zhǎng)和方位精度39。)Nbb1=[ ]Nbb2=[ ]r=206265disp(39。P1P4邊權(quán)函數(shù)式的矩陣形式39。)j=input(39。j=39。)。k=input(39。k=39。)a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。c=r*b/S(j,k)。d=r*a/S(j,k)F1=[c a d b]disp(39。P2P3邊權(quán)函數(shù)式的矩陣形式39。)j=input(39。j=39。)。k=input(39。k=39。)a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。c=r*b/S(j,k)。d=r*a/S(j,k)F2=[c a d b]disp(39。P2P4邊權(quán)函數(shù)式的矩陣形式39。)j=input(39。j=39。)。k=input(39。k=39。)a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。c=r*b/S(j,k)。d=r*a/S(j,k)F3=[c a d b]disp(39。計(jì)算協(xié)因數(shù)39。)Q14=F139。*Nbb1*F1Q23=F239。*Nbb2*F2Q24=F339。*Nbb2*F3附錄B當(dāng)進(jìn)行兩次觀測(cè)時(shí)的程序：X=[ ]Y=[ ]disp(39。當(dāng)觀測(cè)兩次時(shí)的權(quán)為39。)P=[ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ]l=[ ]for j=1 for k=2:6 DX(j,k)=X(k)X(j)。 DY(j,k)=Y(k)Y(j)。 S(j,k)=sqrt((X(j)X(k))^2+(Y(j)Y(k))^2) endendfor j=2 for k=3:6 DX(j,k)=X(k)X(j)。 DY(j,k)=Y(k)Y(j)。 S(j,k)=sqrt((X(j)X(k))^2+(Y(j)Y(k))^2) endendfor i=1 j=input(39。j=39。)。k=input(39。k=39。) a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k) B1=[a b a b]endfor i=2 j=input(39。j=39。)。k=input(39。k=39。) a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。 B2=[a b 0 0]endfor i=3 j=input(39。j=39。)。k=input(39。k=39。) a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。 B3=[a b 0 0]endfor i=4 j=input(39。j=39。)。k=input(39。k=39。) a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。 B4=[a b 0 0]endfor i=5 j=input(39。j=39。)。k=input(39。k=39。) a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。 B5=[0 0 a b]endfor i=6 j=input(39。j=39。)。k=input(39。k=39。) a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。 B6=[0 0 a b]endfor i=7 j=input(39。j=39。)。k=input(39。k=39。) a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。 B7=[0 0 a b]endfor i=8 j=input(39。j=39。)。k=input(39。k=39。) a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。 B8=[0 0 a b]enddisp(39。誤差方程系數(shù)39。)B=[B1。B2。B3。B4。B5。B6。B7。B8]disp(39。法方程系數(shù)矩陣39。)Nbb=B39。*P*Bdisp(39。法方程逆矩陣39。)A=pinv(Nbb)disp(39。精度估算39。)disp(39。計(jì)算邊長(zhǎng)和方位精度39。)Nbb1=[ ]Nbb2=[ ]r=206265disp(39。P1P4邊權(quán)函數(shù)式的矩陣形式39。)j=input(39。j=39。)。k=input(39。k=39。)a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。c=r*b/S(j,k)。d=r*a/S(j,k)F1=[c a d b]disp(39。P2P3邊權(quán)函數(shù)式的矩陣形式39。)j=input(39。j=39。)。k=input(39。k=39。)a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。c=r*b/S(j,k)。d=r*a/S(j,k)F2=[c a d b]disp(39。P2P4邊權(quán)函數(shù)式的矩陣形式39。)j=input(39。j=39。)。k=input(39。k=39。)a=DX(j,k)/S(j,k)。b=DY(j,k)/S(j,k)。c=r*b/S(j,k)。d=r*a/S(j,k)F3=[c a d b]disp(39。計(jì)算協(xié)因數(shù)39。)Q14=F139。*Nbb1*F1Q23=F239。*Nbb2*F2Q24=F339。*Nbb2*F340