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