【正文】 4 N16 N17 N19 N20 N22 N23 N24 N26 N27 N29 N30 N32 N35 N36 最大殘差 續(xù)表412最小殘差 內(nèi)符合精度 v/m外符合精度 12節(jié)點(diǎn)殘差走勢(shì)圖(1) 在二次曲面擬合中，用9個(gè)基點(diǎn)進(jìn)行擬合時(shí)，內(nèi)符合精度和外符合精度都是最小。所以在進(jìn)行多項(xiàng)式曲面擬合時(shí)選擇9個(gè)基點(diǎn)時(shí)其擬合效果最好(2) 在多面函數(shù)擬合中，可以看出隨著節(jié)點(diǎn)的增加，內(nèi)符合精度在變小，同一節(jié)點(diǎn)中隨著平滑因子變大，內(nèi)符合精度和外符合精度都在減小。當(dāng)平滑因子選擇適當(dāng)時(shí)，GPS 水準(zhǔn)聯(lián)測(cè)點(diǎn)越多、顯著點(diǎn)越多時(shí)，擬合精度就越高。(3) 在所取的數(shù)據(jù)中，二次曲面擬合的最大殘差和最小殘差都比多面函數(shù)的較小。二次曲面的外符合精度明顯比多面函數(shù)的外符合精度低，所以，在本數(shù)據(jù)模型中二次曲面的擬合效果比較好。(1) 本章首先介紹了GPS水準(zhǔn)擬合的精度評(píng)定，然后給出了用MATLAB計(jì)算GPS水準(zhǔn)擬合的流程圖。(2) 第三節(jié)和第四節(jié)分別進(jìn)行兩個(gè)案例，一是狹長(zhǎng)線性區(qū)域，二是丘陵區(qū)域。對(duì)狹長(zhǎng)線性區(qū)域分別進(jìn)行了多項(xiàng)式曲線擬合和多面函數(shù)擬合。對(duì)丘陵區(qū)域分別進(jìn)行了多項(xiàng)式曲線擬合，多項(xiàng)式曲面擬合和多面函數(shù)擬合。(3) 對(duì)各個(gè)模型的數(shù)據(jù)進(jìn)行了具體的分析，得出了一些比較合理的結(jié)論。25南京工業(yè)大學(xué)本科生畢業(yè)設(shè)計(jì)(論文)第五章 總結(jié)與展望本文主要從理論和實(shí)際應(yīng)用兩方面研究GPS高程擬合方法，經(jīng)過對(duì)比分析之后，總結(jié)出各種擬合模型的優(yōu)劣，為在工程中取得更好的擬合效果提供參考。高程異常的高精度獲取是構(gòu)建區(qū)域似大地水準(zhǔn)面，實(shí)現(xiàn)GPS高程向正常高（水準(zhǔn)高）轉(zhuǎn)換的關(guān)鍵。本文在總結(jié)歸納常規(guī)的高程異常擬合方法的基礎(chǔ)上，著重討論和分析了多項(xiàng)式曲線擬合，多項(xiàng)式曲面擬合和多面函數(shù)擬合。本文通過分析比較，并編制帶有可改變擬合方法中關(guān)鍵參數(shù)的計(jì)算機(jī)程序，將其作為研究的工具，通過改變這些參數(shù)進(jìn)行反復(fù)、大量的試算后得出一些有益的經(jīng)驗(yàn)和結(jié)論：(1) 對(duì)于多項(xiàng)式擬合，二次曲面就能得到較為滿意的結(jié)果，增大多項(xiàng)式曲面的次數(shù)則意味著需要擬合數(shù)據(jù)中的已知點(diǎn)個(gè)數(shù)就更多，反而降低了多項(xiàng)式擬合方法的適應(yīng)性，建議采用二次曲面擬合比較好。(2) 多面函數(shù)擬合中平滑因子是決定擬合精度高低的關(guān)鍵參數(shù)，根據(jù)實(shí)際的工程實(shí)例的不同，平滑因子值的大小不定，很難找到其中的規(guī)律，相當(dāng)復(fù)雜。(3) MATLAB作為具有強(qiáng)大數(shù)據(jù)處理功能的工具是可以進(jìn)行GPS水準(zhǔn)擬合，并且精度良好，穩(wěn)定性可靠。在傳統(tǒng)的以數(shù)學(xué)模型為基礎(chǔ)的擬合方式中，由于數(shù)學(xué)模型的影響，必然存在模型誤差，如何對(duì)數(shù)學(xué)模型誤差的大小進(jìn)行定量的分析，并對(duì)擬合結(jié)果進(jìn)行模型誤差改正，將會(huì)使擬合精度得到很大的提高，擴(kuò)展 GPS高程擬合的應(yīng)用前景。(1) 在這幾種擬合方法中由于精力和時(shí)間問題，沒有細(xì)致研究擬合點(diǎn)的點(diǎn)位分布對(duì)GPS高程擬合的精度的影響。(2) 在面積更大、地形條件更復(fù)雜的測(cè)區(qū)，多面函數(shù)擬合方法擬合GPS高程的效果有待進(jìn)一步的應(yīng)用分析。(3) 本文只是在已有觀測(cè)數(shù)據(jù)的基礎(chǔ)上做的試驗(yàn)研究，即進(jìn)行的只是GPS高程擬合模型的選擇與應(yīng)用分析，并未對(duì)影響GPS高程擬合結(jié)果的原始數(shù)據(jù)的質(zhì)量控制進(jìn)行探討。應(yīng)從GPS大地高測(cè)量、正常高水準(zhǔn)測(cè)量、擬合模型選擇等多方面誤差源作進(jìn)一步的探討與研究，以確保GPS高程擬合的質(zhì)量。29致謝本文是在導(dǎo)師李明峰教授和助教的精心指導(dǎo)和悉心關(guān)懷下完成的，在我的學(xué)業(yè)和論文的研究工作中無不傾注著導(dǎo)師辛勤的汗水和心血。導(dǎo)師嚴(yán)謹(jǐn)?shù)闹螌W(xué)態(tài)度、淵博的知識(shí)、無私的奉獻(xiàn)精神使我深受啟迪。從尊敬的導(dǎo)師身上，我不僅學(xué)到了扎實(shí)、寬廣的專業(yè)知識(shí)，也學(xué)會(huì)了做人的道理。在此，我要向?qū)熤乱宰钪孕牡母兄x和最崇高的敬意！同時(shí)，我也要感謝在日常學(xué)習(xí)和生活中給我莫大幫助的王永明和李春進(jìn)助教，他們嚴(yán)謹(jǐn)?shù)膶W(xué)態(tài)度和人生理念，深深影響著我，督促我不斷進(jìn)步，并將使我終生受益。他聯(lián)系相關(guān)部門為我尋找論文所需數(shù)據(jù)，解決我論文后期最主要的問題，感謝他給予我的幫助和支持！感謝包國(guó)鋒、朱桂林兩位同學(xué)的關(guān)懷和幫助，感謝梁偉亮與我同甘共苦、共同奮斗。衷心感謝我的父母家人以及女朋友的父母家人! 在他們的全力支持和幫助下，才使我順利完成學(xué)業(yè)，在此向他們表示衷心的感謝和祝福！最后，衷心地感謝在百忙之中評(píng)閱論文和參加答辯的各位專家、學(xué)者！南京工業(yè)大學(xué)本科生畢業(yè)設(shè)計(jì)(論文)參考文獻(xiàn)[1] 盧輝,[J].測(cè)繪科學(xué),2009,34(2)[2] 姚東,王愛民,馮峰,[M].北京:人民郵電出版社,2000.[3] 許波,[M].北京:清華大學(xué)出版社,2000.[4] 石博強(qiáng),[M].上海:中國(guó)鐵道出版社,2005.[5] 謝劭峰,[J].測(cè)繪與空間地理信息,2005,28(6)[6] 徐紹銓,張華海,楊志強(qiáng),[M].武漢:武漢大學(xué)出版社,2005[7] 曹先革,劉金鵬,[J].地理空間信息,2009,7(2)[8] 蔡群,[J].科技資訊,2006,32[9] 的GPS水準(zhǔn)擬合方法的應(yīng)用與比較[J].測(cè)繪與空間地理信息,2010,33(6)[10] 陳本富,王貴武,沈慧, 的數(shù)據(jù)處理方法在GPS高程擬合中的英語[J].昆明理工大學(xué)學(xué)報(bào),2009,34(5)[11] 吳偉盛,[J].科技資訊,2008,28[12] 劉萬選,[J].測(cè)繪與空間地理信息,2009,32(3) [13] 侯東陽,[J].大地測(cè)量與地球動(dòng)力學(xué),2010,30(6) [14] 李秀海,[J]. 測(cè)繪與空間地理信息,2010,33(1) [15] 胡伍生,華錫生,[J]. 測(cè)繪學(xué)報(bào),2002,31(2)[16] 陳建寶,[J]. 昆明大學(xué)學(xué)報(bào),2006,17(4)[17] [J]. 測(cè)繪學(xué)報(bào),2003,(7)：1114.[18] [J].人民黃河, 2010(9)[19] [J].礦山測(cè)量,2009(3)[20] 黃聲享，郭英起， 在測(cè)量工程中的應(yīng)用[M].北京:測(cè)繪出版社，200731南京工業(yè)大學(xué)本科生畢業(yè)設(shè)計(jì)(論文)附錄多項(xiàng)式曲線擬合MATLAB程序x=[ ] %已知點(diǎn)z=[ ] %已知點(diǎn)高程異常B=xlsread() %待求點(diǎn)和高程異常x1=B(:,2) %待求點(diǎn)z1=B(:,1) %待求點(diǎn)擬合前的高程異常p=polyfit(x,z,1)z2=polyval(p,x1)%擬合后的高程異常v=z1z2多項(xiàng)式曲面擬合MATLAB程序A=xlsread()。 %已知點(diǎn)的x，y和高程異常E=xlsread()。 %待求點(diǎn)的X,Y和擬合前高程異常a=A(:,1)。x0=0。for i=1:8 x0=x0+a(i,1)。endx0=x0/8 %求已知點(diǎn)X的平均值b=A(:,2)。y0=0。for i=1:8 y0=y0+b(i,1)。endy0=y0/8 %求已知點(diǎn)Y的平均值p=zeros(8,8)。for n=1:16 B(n,1)=E(n,1)x0。 B(n,2)=E(n,2)y0。 Z(n)=E(n,3)。 M(n,1)=1。 M(n,2)=B(n,1)。 M(n,3)=B(n,2)。 M(n,4)=B(n,1)^2。 M(n,5)=B(n,1)*B(n,2)。 M(n,6)=B(n,2)^2。 d(n)=sqrt(B(n,1)^2+B(n,2)^2)。 P(n,n)=1/(d(n))^2。endX=inv(M39。*P*M)*M39。*P*Z39。 %求系數(shù)V=M*Xv=E(:,3)V多面函數(shù)擬合MATLAB程序function q=q(x1,y1,x2,y2)k=1q=sqrt((x1x2)^2+(y1y2)^2+k^2)。A=xlsread() %已知數(shù)據(jù)B=xlsread() %待求數(shù)據(jù)x=B(:,1)y=B(:,2)z=B(:,3)X=A(:,1)Y=A(:,2)Z=A(:,3)for m=1:16 for n=1:8 W(m,n)=[q(x(m,1),y(m,1),X(n,1),Y(n,1))] endend for m=1:16 for n=1:8 F(m,n)=[q(X(m,1),Y(m,1),X(n,1),Y(n,1))]%系數(shù)矩陣 endend H=W*inv(F39。*F)*F39。*Zv=zH33南京工業(yè)大學(xué)本科生畢業(yè)設(shè)計(jì)(論文)35