三維激光掃描儀講解,三維激光掃描儀原理,三維激光掃描儀案例分析

.,講座2,三維激光掃描儀的原理和結(jié)構(gòu),.,復(fù)習(xí),三維激光掃描儀的概念三維激光掃描儀的由來三維激光掃描儀的技術(shù)應(yīng)用和發(fā)展,.,本節(jié)課的內(nèi)容,三維激光掃描儀的原理和結(jié)構(gòu)掌握：坐標(biāo)計(jì)算的原理了解：三維激光掃描儀的結(jié)構(gòu),.,Cyra三維激光掃描系統(tǒng)簡(jiǎn)介,徠卡Cyra三維激光掃描系統(tǒng)組成1、Cyrax2500三維激光掃描儀2、Cyclone3.0軟件。3、筆記本電腦4、三腳架5、電池和充電器6、目標(biāo)牌和電纜,.,二、徠卡Cyra三維激光掃描系統(tǒng)的工作原理,Cyrax2500三維激光掃描儀向目標(biāo)發(fā)射激光脈沖，依次掃描被測(cè)區(qū)域，快速獲得地面景觀的三維坐標(biāo)和反射光強(qiáng)，利用Cyclone3.0軟件進(jìn)行三維建模，生成地面景觀的三維圖象和可量測(cè)點(diǎn)陣數(shù)據(jù)，并可方便地轉(zhuǎn)化為多種輸出格式的圖形產(chǎn)品。,結(jié)構(gòu)和原理,Cyrax三維激光掃描系統(tǒng)由三維激光掃描儀和配套軟件組成。儀器內(nèi)部有一個(gè)激光器，兩個(gè)旋轉(zhuǎn)軸異面且互相垂直的反光鏡。反光鏡由步進(jìn)電機(jī)帶動(dòng)旋轉(zhuǎn)，而激光器發(fā)射的窄束激光脈沖在反光鏡作用下，沿縱向和橫向依次掃過被測(cè)區(qū)域。激光脈沖被物體漫反射后，一部分能量被三維激光掃描儀接收。測(cè)量每個(gè)激光脈沖從發(fā)出到返回儀器所經(jīng)過的時(shí)間，可以計(jì)算出儀器和物體間的距離S。,同時(shí)測(cè)量每個(gè)激光脈沖與儀器固有坐標(biāo)系X軸的夾角，XOY面的夾角，可以由公式1算出被測(cè)物體表面點(diǎn)的三維坐標(biāo)。通常用儀器內(nèi)部坐標(biāo)系統(tǒng)。根據(jù)掃描點(diǎn)的激光反射強(qiáng)度,給反射點(diǎn)匹配顏色。掃描點(diǎn)繪制在屏幕上，組成密集的點(diǎn)云。,.,儀器坐標(biāo)系,.,坐標(biāo)計(jì)算公式,=,.,三、徠卡Cyra三維激光掃描系統(tǒng)的主要特點(diǎn),1、體積小、重量輕，操作簡(jiǎn)單、裝拆便利，具備良好的野外操作性能。2、掃描范圍大、速度快、精度高。3、快速建立三維景觀模型、圖形圖象數(shù)據(jù)一次獲取。,.,Cyrax2500三維激光掃描儀產(chǎn)品參數(shù)1,總體儀器類型高速、高精度激光探測(cè)掃描儀使用界面筆記本電腦掃描儀驅(qū)動(dòng)器伺服馬達(dá)光學(xué)視眼一體化視頻照相機(jī)系統(tǒng)性能單點(diǎn)精度點(diǎn)位6mm/1.5m-50m距離距離1mm角度0.5”形成模型表面的精度2mm,.,Cyrax2500三維激光掃描儀產(chǎn)品參數(shù)2,激光掃描系統(tǒng)激光類型脈沖；專用芯片。顏色綠色安全性2級(jí)（參照CFR1040）對(duì)眼安全，但要避免直射和長(zhǎng)時(shí)間照射眼睛光斑大小6mm/0-50m范圍測(cè)距范圍最大150m推薦1.5m-50m（5%-100%）掃描率每秒1行/每行小于1000個(gè)點(diǎn)每秒2行/每行小于200個(gè)點(diǎn)掃描密度可選性水平方向和垂直方向掃描可以分開選擇可以按一定間隔選取點(diǎn)對(duì)點(diǎn)的掃描方式垂直方向點(diǎn)與點(diǎn)間的最小間隔0.25mm（小于50米）水平方向點(diǎn)與點(diǎn)間的最小間隔0.25mm（小于50米）行掃描每行最多1000個(gè)點(diǎn)列掃描每列最多1000個(gè)點(diǎn),.,Cyrax2500三維激光掃描儀產(chǎn)品參數(shù)3,視野垂直方向最大40角水平方向最大40角三角基座360水平旋轉(zhuǎn)+105/-90垂直方向旋轉(zhuǎn)空間體積19900m3/掃描（精度小于6mm）占據(jù)體積（最大）1592003/掃描（距離小于100m）掃描光雙鏡，隨機(jī)通道由外殼和玻璃罩保護(hù)視頻目標(biāo)480480彩色電力供電AC90-240VAC;50-60HzDC12V功率100W電池密封的鉛酸電池電池容量20C以下供電4小時(shí),.,供掃描與建模的計(jì)算機(jī)平臺(tái),PC系統(tǒng)最低要求推薦處理器200MHz/奔騰500MHz/PRAM64MB256MB硬盤2GB20GB網(wǎng)卡以太網(wǎng)以太網(wǎng)視頻卡SVGA3D圖形加速操作系統(tǒng)WindowsNT4.0WindowsNT4.0顯示器800600；256色1024768；真彩色,.,CYRAX掃描儀的控制與操作,可以建立以下類型的幾何體：平面角射線柱面盒目標(biāo)中心點(diǎn)頂點(diǎn)圓錐球體直線多邊形創(chuàng)建觀察和交互窗口幾何校正根據(jù)激光反射強(qiáng)度獲取彩色掃描圖象用戶定義質(zhì)量檢查方式平面、柱面、球面的區(qū)域擴(kuò)張多次掃描的坐標(biāo)糾正對(duì)目標(biāo)采用平面或球面建模的方式自動(dòng)實(shí)施整體糾正抽取點(diǎn)云圖支持多種坐標(biāo)系統(tǒng)目標(biāo)的插入、復(fù)制和編輯獲取并在計(jì)算機(jī)屏上顯示視頻圖象基本的概念設(shè)計(jì)和二維繪圖工具遙控操作,.,特性和功能1,通過選取線性矩形框來確定目標(biāo)范圍任意選取水平方向和垂直方向的掃描密度完全的“飛行旋轉(zhuǎn)”，縮放和搖動(dòng)，掃描影像的任意旋轉(zhuǎn)，哪怕在掃描的過程中也能做到。給對(duì)象增加顏色和材料屬性影像存儲(chǔ)/恢復(fù)查看掃描儀的位置智能化的3D模型輸出視頻圖象和掃描數(shù)據(jù)圖象同時(shí)顯示3D標(biāo)注可視化客戶許可證掃描圖象和影像的自動(dòng)連接多級(jí)undo/redo操作功能掃描圖自動(dòng)排序診斷信息計(jì)算機(jī)輔助獲取目標(biāo),.,特性和功能2,3D畫線快速將掃描圖轉(zhuǎn)換為格網(wǎng)圖云點(diǎn)圖可永久保存用戶界面：熱鍵、工具條64位雙精度浮點(diǎn)數(shù)據(jù)連續(xù)自動(dòng)存貯客戶/服務(wù)器對(duì)象本底數(shù)據(jù)庫在線幫助項(xiàng)目分層圖軟件使用許可方式靈活多種對(duì)象注釋方式測(cè)量點(diǎn)云和模型中點(diǎn)間斜距數(shù)據(jù)管理建立并管理圖層建立并管理對(duì)象注釋與環(huán)境配套的照明設(shè)施公制或英制測(cè)量單位可選輸出ASC格式的數(shù)據(jù)流，BMP、JPEG格式的圖形文件輸入ASC格式的數(shù)據(jù)流，Riegl，CGP格式的數(shù)據(jù)，通過COE輸入MicroStation格式數(shù)據(jù)。,.,其他三維激光掃描儀的原理和結(jié)構(gòu),Vivid910Optech天寶（法國Mensi）,本節(jié)課結(jié)束,課后要求,.,講座3,三維激光掃描儀的操作方法,.,復(fù)習(xí),三維激光掃描儀的原理和結(jié)構(gòu)坐標(biāo)計(jì)算公式,.,本節(jié)課的內(nèi)容,重要概念：點(diǎn)云、標(biāo)靶、拼合等三維激光掃描儀的使用方法,.,點(diǎn)云,筆記本中的點(diǎn)云圖,.,.,本節(jié)課內(nèi)容,用實(shí)際儀器和筆記本演示：測(cè)量前儀器的準(zhǔn)備和儀器連線方法；掃描軟件的啟動(dòng)；掃描軟件的操作方法；粗掃和精掃以及標(biāo)靶掃描等。,.,Cyclone軟件的使用方法,演示,.,Cyrax掃描儀的控制和操作,計(jì)算機(jī)屏上獲取和顯示視頻圖象遙控測(cè)量用矩形框選擇需要掃描的區(qū)域水平掃描和垂直掃描的密度可以分開設(shè)定完全的“飛行旋轉(zhuǎn)”，縮放和搖動(dòng)，掃描影像的任意旋轉(zhuǎn)，哪怕在掃描的過程中也能作到。自動(dòng)管理掃描順序圖診斷信息計(jì)算機(jī)輔助獲取目標(biāo)自動(dòng)連接掃描影像,.,Cyrax掃描儀的控制和操作,建模、可視化、查詢和瀏覽工具對(duì)云點(diǎn)圖、感光圖、VDB、網(wǎng)點(diǎn)圖、渲染圖及3D模型可以實(shí)行“飛行旋轉(zhuǎn)”、搖晃和縮放以及自由旋轉(zhuǎn)。建立交互式的可視窗口,.,Cyrax掃描儀的控制和操作,測(cè)量和標(biāo)注點(diǎn)云和模型的相關(guān)數(shù)據(jù)：斜距平距高差體積表面積坐標(biāo)參考系統(tǒng)由激光返回信號(hào)的強(qiáng)度生成彩色影像在對(duì)象上添加顏色和材料采用最合適的施工方式動(dòng)態(tài)展現(xiàn)施工各個(gè)階段3D模型用戶可設(shè)置質(zhì)量檢查方式自動(dòng)生成平面和柱面的切面平面、柱面和球面的區(qū)域擴(kuò)展由3D模型派生2D圖,.,Cyrax掃描儀的控制和操作,根據(jù)掃描圖派生2D線畫圖觀察掃描儀所在位置多幅掃描圖的整體拼接采用平面或球面建模自動(dòng)實(shí)施整體糾正對(duì)掃描點(diǎn)自動(dòng)實(shí)施干擾檢查并進(jìn)行可視化的目標(biāo)設(shè)計(jì)根據(jù)格網(wǎng)點(diǎn)生成等高線按照影像數(shù)據(jù)文件格式存儲(chǔ)瞬時(shí)圖根據(jù)點(diǎn)云生成剖面圖、平面圖和斷面圖,.,Cyrax掃描儀的控制和操作,可以創(chuàng)建下面的幾何形狀：平面、柱面、頂點(diǎn)、直線、彎頭、圓錐、圓弧、球面、管狀、錐形管、金屬外型、盒子、立方體、拐角、不規(guī)則多邊形、橢圓、圓、復(fù)線、正方形、矩形、目標(biāo)中心點(diǎn)。,.,儀器直接輸入的格式,ASC格式點(diǎn)位數(shù)據(jù)Riegl通過COE輸入MicroStation格式數(shù)據(jù)CGP,.,儀器直接輸出的格式,AutoCADDXF格式，版本12到2000AutoCADDWG格式，版本12MicroStation數(shù)據(jù)格式（通過MicroStation或COE接入）ASC格式點(diǎn)位數(shù)據(jù)流Cyclone軟件交換格式COEBMPJPEGIGES*PDS*PDMS*AutoPLANT*表示可選,.,間接輸出的格式,PDS，通過DGN（從MicroStation或COE接入）轉(zhuǎn)換PDMS，通過DGN（從MicroStation或COE接入或直接輸入）轉(zhuǎn)換AutoPLANT，通過AutoCADDXF轉(zhuǎn)換,.,本節(jié)結(jié)束,思考題,.,控制標(biāo)靶和目標(biāo)球既是拼合各次掃描點(diǎn)云的基準(zhǔn)點(diǎn)，還是網(wǎng)格坐標(biāo)與儀器內(nèi)部坐標(biāo)轉(zhuǎn)換的媒介，它們的位置要保持穩(wěn)定，不能安裝在變形體上。但掃描儀的視角有限（），無法一次掃描變形體和穩(wěn)固地點(diǎn)的控制標(biāo)靶和目標(biāo)球。因此這是一對(duì)矛盾。解決這個(gè)問題的方法是：在穩(wěn)固地點(diǎn)和變形體上安裝多個(gè)控制標(biāo)靶和目標(biāo)球，增加掃描站數(shù)。儀器由穩(wěn)固地點(diǎn)逐步掃描至變形體，根據(jù)穩(wěn)固地點(diǎn)的控制標(biāo)靶和目標(biāo)球中心坐標(biāo)修正變形體上的目標(biāo)球中心坐標(biāo)。,.,Cyra三維激光掃描系統(tǒng)主要應(yīng)用領(lǐng)域,1、地面景觀形體測(cè)量；2、橋梁改擴(kuò)建工程、橋梁結(jié)構(gòu)測(cè)量；3、大壩和電站工程基礎(chǔ)地形測(cè)量；4、地下工程結(jié)構(gòu)測(cè)量；5、高陡邊坡地形測(cè)量及工程量計(jì)算；6、公路測(cè)量；7、礦山測(cè)量及體積計(jì)算；8、電影特殊效果中的三維景觀制作；9、古文物、建筑測(cè)量和資料保存；10、電廠、化工廠等大型工業(yè)企業(yè)內(nèi)部設(shè)備的測(cè)量；11、管道、線路測(cè)量；12、大比例尺GIS系統(tǒng)和數(shù)據(jù)庫數(shù)據(jù)更新；13、虛擬模型的驗(yàn)證,.,.,.,觀測(cè)實(shí)驗(yàn)1,.,觀測(cè)實(shí)驗(yàn)2,.,掃描點(diǎn)云,.,掃描點(diǎn)云,.,線劃圖,.,.,.,.,.,渲染和紋理,.,.,.,.,講座三維激光掃描儀用于變形監(jiān)測(cè),張國輝遼寧工程技術(shù)大學(xué)測(cè)繪學(xué)院,.,監(jiān)測(cè)對(duì)象和必要性,露天礦開采中的滑坡，井下開采造成的地面沉陷，建（構(gòu)）筑物的變形等，每年都造成大量人員傷亡和財(cái)產(chǎn)損失。,.,已有的變形監(jiān)測(cè)手段,、大地測(cè)量法，建立變形監(jiān)測(cè)網(wǎng)，進(jìn)行邊角網(wǎng)觀測(cè)等、在變形區(qū)域安裝各種傳感器、近景攝影測(cè)量、GPS變形監(jiān)測(cè)、測(cè)量機(jī)器人,.,已有監(jiān)測(cè)手段的優(yōu)缺點(diǎn),、技術(shù)的成熟性、精度和可靠性、成本、效率等,.,新的變形監(jiān)測(cè)設(shè)備和技術(shù),三維激光掃描儀是目前國際上最先進(jìn)的獲取物體三維點(diǎn)陣的設(shè)備。它可以將任何大型的、復(fù)雜的實(shí)體或?qū)嵕暗娜S數(shù)據(jù)完整地采集到計(jì)算機(jī)中，進(jìn)而快速建立目標(biāo)的三維模型并提取線、面、體等各種制圖數(shù)據(jù)，它所采集的三維激光點(diǎn)云數(shù)據(jù)經(jīng)過處理還可應(yīng)用在許多場(chǎng)合。三維激光掃描技術(shù)改變了傳統(tǒng)的單點(diǎn)數(shù)據(jù)采集方式，可以在很短時(shí)間內(nèi)以很小的采樣間隔，獲取實(shí)體表面各點(diǎn)坐標(biāo)，實(shí)現(xiàn)“實(shí)景復(fù)制”,.,三維激光掃描儀應(yīng)用于變形監(jiān)測(cè),優(yōu)勢(shì)：、高效率、數(shù)據(jù)精細(xì)、與近景攝影測(cè)量相比、方便建立模型缺點(diǎn)：、目前價(jià)格高,監(jiān)測(cè)對(duì)象例子,多發(fā)滑坡崩塌事故的某露天礦的高陡邊坡需要變形監(jiān)測(cè)。采用三維激光掃描儀監(jiān)測(cè)變形，需要解決的難題是：如何布設(shè)測(cè)站和控制標(biāo)靶和目標(biāo)球；如何統(tǒng)一各次掃描的坐標(biāo)系統(tǒng)；如何提取變形信息；變形監(jiān)測(cè)誤差的估算；如何顯示變形數(shù)據(jù)等。,.,變形監(jiān)測(cè)方法,測(cè)站地點(diǎn)要選擇在地面堅(jiān)實(shí)、受變形影響小，且離變形體近（m）的地方；布設(shè)專用標(biāo)牌；初始掃描并內(nèi)業(yè)處理；按照觀測(cè)要求，實(shí)施各期掃描并內(nèi)業(yè)處理；提交數(shù)據(jù)和結(jié)論。,.,提取變形信息,由于點(diǎn)云數(shù)據(jù)極其密集，靠視力很難分辨一個(gè)點(diǎn)的變形情況，所以必須對(duì)變形體進(jìn)行處理，以方便提取變形信息。,.,提取方法1,方法一：在變形體表面安裝多個(gè)涂色的球形標(biāo)志，根據(jù)顏色和形狀，可以從掃描點(diǎn)云中將球形標(biāo)志分辨出來，用Cyclone軟件建立球模型并輸出球心坐標(biāo)，通過比較各時(shí)段掃描數(shù)據(jù)中相同球心坐標(biāo)變化來提取變形信息。,.,.,.,提取方法2,方法二：根據(jù)所有點(diǎn)云數(shù)據(jù)建立變形體的數(shù)字地面模型（DTM），讓所有時(shí)段的變形體模型的坐標(biāo)系統(tǒng)一致，直接分析所有模型的變形。對(duì)于滑坡監(jiān)測(cè)，可以采用數(shù)字高程模型（DEM）。,.,DEM,DEM表示的是各坐標(biāo)點(diǎn)處的高程。,.,根據(jù)DEM提取變形信息,由于不同時(shí)間段的DEM并不完全相同，為了比較相同水平坐標(biāo)點(diǎn)的高程變化，需要以初始掃描時(shí)建立的DEM數(shù)據(jù)作為參考，將后面的DEM進(jìn)行內(nèi)插計(jì)算。以該坐標(biāo)相鄰點(diǎn)的高程加權(quán)平均值作為該點(diǎn)高程。最后比較相同點(diǎn)的高程變化來分析變形。,.,變形結(jié)果輸出和顯示,提取出的地形變化可以用Matlab和CAD下自主開發(fā)的軟件來表示。也可以用表格、文本、斷面圖、曲線圖、三維變形曲面圖來表示。,.,精度分析,采用上述方法一來提取變形數(shù)據(jù)，由于模型表面精度為2mm，可以滿足變形監(jiān)測(cè)要求。采用方法二，是將各個(gè)單點(diǎn)數(shù)據(jù)進(jìn)行處理，取其平均值，則精度也要優(yōu)于6mm。參與計(jì)算的點(diǎn)越多精度越好。而用傳統(tǒng)的全站儀來觀測(cè)，其精度在厘米級(jí)。所以兩者沒有明顯的精度區(qū)別。,.,.,.,實(shí)驗(yàn)結(jié)果,實(shí)際量測(cè)地物表面尺寸，與掃描軟件提供的尺寸比較。證明符合良好。,誤差分析,三維激光腳點(diǎn)測(cè)量誤差的影響因素較多，大致可分為三類：儀器誤差、與目標(biāo)物體反射面有關(guān)的誤差、外界環(huán)境條件。儀器誤差是儀器本身性能缺陷造成的測(cè)量誤差，包括激光測(cè)距的誤差；掃描角度測(cè)量的誤差；與目標(biāo)物體反射面有關(guān)的誤差主要包括目標(biāo)物體反射面傾斜的影響和表面粗糙度的影響；外界環(huán)境條件主要包括溫度、氣壓等因素。,.,三維激光掃描儀用于變形監(jiān)測(cè),本節(jié)課的要點(diǎn)：變形監(jiān)測(cè)概念變形監(jiān)測(cè)的方法三維激光掃描儀用于變形監(jiān)測(cè)的實(shí)驗(yàn)對(duì)新技術(shù)的展望,.,變形分析方法,傳統(tǒng)的變形幾何分析主要包括參考點(diǎn)的穩(wěn)定性分析、觀測(cè)值的平差處理和質(zhì)量評(píng)定以及變形模型參數(shù)估計(jì)等內(nèi)容。陳永奇(1988)概括了兩種基本的分析方法，即直接法和位移法。直接法是直接用原始的重復(fù)觀測(cè)值之差計(jì)算應(yīng)變分量或它們的變化率;位移法是用各測(cè)點(diǎn)坐標(biāo)的平差值之差(位移值)計(jì)算應(yīng)變分量。,.,變形數(shù)據(jù)分析方法,對(duì)變形數(shù)據(jù)分析方法研究是極為活躍的，除了傳統(tǒng)的多元回歸分析法以及上述的時(shí)間序列分析法、頻譜分析法和濾波技術(shù)之外，灰色系統(tǒng)理論、神經(jīng)網(wǎng)絡(luò)等非線性時(shí)間序列預(yù)測(cè)方法也得到了一定程度的應(yīng)用。為了彌補(bǔ)單一方法的缺陷，研究多種方法的結(jié)合得到了一定程度的發(fā)展。,.,一般應(yīng)用較多的是單測(cè)點(diǎn)模型，同時(shí)，為顧及觀測(cè)點(diǎn)的整體空問分布特性，多測(cè)點(diǎn)變形監(jiān)控模型也得到了發(fā)展。,.,.,觀測(cè)網(wǎng)選擇和參考點(diǎn)穩(wěn)定性分析,觀測(cè)網(wǎng)一般分為參考網(wǎng)（也叫絕對(duì)網(wǎng)）和相對(duì)網(wǎng)。基準(zhǔn)點(diǎn)是進(jìn)行變形觀測(cè)的起算基準(zhǔn)點(diǎn)。坐標(biāo)系除了提供方位外，更重要的是為系統(tǒng)數(shù)據(jù)處理時(shí)的距離及高差計(jì)算提供基準(zhǔn)。很多三維激光掃描儀沒有對(duì)中裝置。需要由標(biāo)牌來傳遞坐標(biāo)系。,.,對(duì)于參考點(diǎn)穩(wěn)定性檢驗(yàn)，如果一組公共點(diǎn)在兩周期測(cè)量之間內(nèi)部沒有相對(duì)位移，盡管兩期測(cè)量計(jì)算坐標(biāo)時(shí)所依據(jù)的參考系有變化而造成求得的坐標(biāo)有視的變化，但這一組點(diǎn)間的相對(duì)幾何關(guān)系并沒有改變，我們把這樣一組點(diǎn)稱之為穩(wěn)定點(diǎn)組，依據(jù)一組點(diǎn)間的幾何關(guān)系是否保持不變來判別這一組點(diǎn)是否可認(rèn)為它們是相對(duì)穩(wěn)定的。,.,設(shè)第一期觀測(cè)成果中，各基準(zhǔn)點(diǎn)在測(cè)站坐標(biāo)下的坐標(biāo)，第二觀測(cè)周期各點(diǎn)經(jīng)轉(zhuǎn)換后的坐標(biāo)，可得到各點(diǎn)在兩周期觀測(cè)時(shí)刻的坐標(biāo)差。坐標(biāo)差可得各測(cè)量點(diǎn)在兩觀測(cè)周期得到理論位移量,.,變形位移值的計(jì)算,參考點(diǎn)穩(wěn)定性檢測(cè)結(jié)束后，以不同周期觀測(cè)的參考點(diǎn)坐標(biāo)計(jì)算轉(zhuǎn)換參數(shù)，進(jìn)行坐標(biāo)系的轉(zhuǎn)換，從而獲得各個(gè)變形點(diǎn)的統(tǒng)一坐標(biāo)系下的坐標(biāo)。各期觀測(cè)坐標(biāo)求差即可得到變形數(shù)據(jù)。,.,變形觀測(cè)數(shù)據(jù)回歸分析,用Matlab回歸分析：在Matlab中采用多項(xiàng)式曲線擬合方法來實(shí)現(xiàn)兩個(gè)變量之間是這種非線形關(guān)系。函數(shù)Polyfit()根據(jù)一組給定的數(shù)據(jù)在最小平方意義下的擬合多項(xiàng)式的系數(shù)。P=polyfit(x,y,n),x,y是擬合數(shù)據(jù)量，n是返回多項(xiàng)式的階次。,.,代碼,設(shè)置x為時(shí)間/月份，y為累計(jì)位移值，如果用3次多項(xiàng)式來擬合，p=polyfit(x,y,3)中的p為多項(xiàng)式的系數(shù)。先假設(shè)一組數(shù)，計(jì)算x=2時(shí)的值，執(zhí)行polyval(p