張志偉 張劍輝

(南京橋新建設(shè)有限公司，江蘇南京 211103)

1 概述

隨著科學(xué)技術(shù)的迅猛發(fā)展和我國現(xiàn)代化進(jìn)程的不斷加快，城市各類高層、超高層建筑物和重要建筑物日益增多。隨著建筑物施工的進(jìn)展，荷載不斷增加，在地基基礎(chǔ)和上部結(jié)構(gòu)的共同作用下，建筑物將發(fā)生沉降。為了建筑物的施工使用安全，建筑物施工過程中，應(yīng)加強(qiáng)沉降觀測及時(shí)掌握建筑工程沉降程度及沉降趨勢，處理沉降觀測結(jié)果，并采取合理的方法科學(xué)預(yù)測建筑物的沉降變形，為勘察、設(shè)計(jì)、施工等部門提供詳盡的一手資料，指導(dǎo)合理的施工工序，保證建筑物的正常使用壽命和建筑物的安全性。在沉降觀測數(shù)據(jù)處理中，建立沉降量預(yù)測模型、繪制沉降曲線是建筑物沉降觀測的重要內(nèi)容。Matlab有一個(gè)功能強(qiáng)大的曲線擬合工具箱Curve Fitting Tool，使用方便。在Matlab完成繪制沉降曲線后，數(shù)據(jù)分析人員即能根據(jù)實(shí)際沉降曲線和Matlab工具箱提供的多種類型的線性、非線性曲線擬合形式，實(shí)現(xiàn)沉降量預(yù)測模型的建立。本文根據(jù)南京江寧區(qū)某小區(qū)其中一棟建筑物的沉降觀測數(shù)據(jù)為例，對Matlab在建筑物沉降觀測數(shù)據(jù)處理中的應(yīng)用進(jìn)行了研究。

2 M atlab Curve Fitting Tool工具箱簡介

Curve Fitting Tool工具箱提供用于曲線擬合數(shù)據(jù)的圖形工具和函數(shù)。使用該工具箱可進(jìn)行數(shù)據(jù)分析、比較候選模型和刪除偏值等。本文涉及的主要功能有:用于曲線擬合的圖形工具。

使用自定義方程求解線性和非線性模型。

2.1 曲線擬合類型

Curve Fitting Tool工具箱提供的擬合類型有:Custom Equations:用戶自定義的函數(shù)類型;Exponential:指數(shù)函數(shù)，有2種類型:a*exp(b*x)，a*exp(b*x)+c*exp(d*x);Fourier:傅立葉函數(shù)，有7種類型，基礎(chǔ)型是a0+a1*cos(x*w)+b1*sin(x*w);Gaussian:高斯函數(shù)，有8種類型，基礎(chǔ)型是a1*exp(－((xb1)/c1)^2);Interpolant:插值函數(shù)，有 4種類型，linear，nearest neighbor，cubic spline，shape-preserving;Polynomial:多項(xiàng)式函數(shù)，有9 種類型，linear，quadratic，cubic，4-9th degree;Power:冪函數(shù)，有2 種類型:a*x^b，a*x^b+c;Rational:有理數(shù)函數(shù)，分子、分母共有的類型是 linear，quadratic，cubic，4-5th degree，此外，分子還包括 constant型;Smoothing Spline:平滑函數(shù);Sum of Sin Functions:正弦曲線函數(shù)，有8種類型，基礎(chǔ)型是a1*sin(b1*x+c1);Weibull:只有1種，a*b*x^(b－1)*exp(－a*x^b)

2.2 應(yīng)用

1)點(diǎn)擊Matlab左下角的Start→Toolboxes→Curve Fitting;2)點(diǎn)擊Curve Fitting Tool進(jìn)入曲線擬合工具箱界面“Curve Fitting Tool”;3)在Data Sets頁面里的X Data選項(xiàng)中選取x向量，Y Data選項(xiàng)中選取y向量，此時(shí)Curve Fitting Tool窗口中顯示出這一數(shù)據(jù)組的散點(diǎn)分布圖;4)選擇好所需的擬合曲線類型及其子類型，并進(jìn)行相關(guān)設(shè)置。——如果是非自定義的類型，根據(jù)實(shí)際需要點(diǎn)擊“Fit options”按鈕，設(shè)置擬合算法、修改待估計(jì)參數(shù)的上下限等參數(shù);——如果選 Custom Equations，點(diǎn)擊“New”按鈕，彈出自定義函數(shù)等式窗口，有“Linear Equations線性等式”和“General Equations構(gòu)造等式”兩種標(biāo)簽。然后點(diǎn)擊OK。即可顯示出擬合結(jié)果(見圖1)，非常方便。

3 應(yīng)用實(shí)例

南京江寧區(qū)某小區(qū)共有26棟高層住宅樓，其中第26棟建筑物的沉降點(diǎn)布置見圖2，其中3，7，11，14，17，21沉降觀測數(shù)據(jù)見表1。

圖1 Curve Fitting Tool擬合結(jié)果圖

圖2 26棟沉降點(diǎn)布置圖

表1 沉降觀測數(shù)據(jù)表mm

為比較擬合模型的穩(wěn)定性，選用 3，7，11，14，17，21 6 個(gè)點(diǎn)進(jìn)行擬合。根據(jù)“Table of Fits”框中顯示的SSE(誤差平方和):擬合數(shù)據(jù)和原始數(shù)據(jù)對應(yīng)點(diǎn)的誤差的平方和。SSE越接近0，說明模型選擇和擬合更好，數(shù)據(jù)預(yù)測也越成功;RMSE(均方根):回歸系統(tǒng)的擬合標(biāo)準(zhǔn)差;R-square(相關(guān)系數(shù)):其正常取值范圍為0～1，越接近1，表明方程的變量對y的解釋能力越強(qiáng)，這個(gè)模型對數(shù)據(jù)的擬合也越好。選用擬合函數(shù)為二項(xiàng)指數(shù)擬合，其擬合公式為y=a·ebx+c·edx。根據(jù)表1的數(shù)據(jù)，得擬合結(jié)果見表2。其中3號點(diǎn)的沉降量散點(diǎn)和擬和曲線圖見圖1。

表2 擬合結(jié)果表

從表2可以看出，這6個(gè)點(diǎn)的擬合函數(shù)系數(shù)相近，說明所選模型較穩(wěn)定。

4 結(jié)語

本文利用Matlab Curve Fitting Tool工具箱對建筑物沉降觀測數(shù)據(jù)處理進(jìn)行分析處理，建立預(yù)測模型。由于選擇模型是根據(jù)SSE，RMSE和R-square的數(shù)值比較而得，因此所建立的預(yù)測模型更具科學(xué)性，較符合實(shí)際情況。能更好的為指導(dǎo)施工、驗(yàn)證設(shè)計(jì)計(jì)算參數(shù)、檢驗(yàn)施工質(zhì)量、建成后安全運(yùn)營以及工程質(zhì)量事故判斷與處理提供必要的數(shù)據(jù)支撐。

［1］陳 健.MATLAB在變形監(jiān)測數(shù)據(jù)處理中的應(yīng)用［J］.城市勘測，2009(2):130-133.

［2］基于MATLAB的高層建筑沉降變形監(jiān)測數(shù)據(jù)處理［J］.測繪與空間地理信息，2010，32(5):3-5.

［3］石博強(qiáng)，滕貴法，李海鵬，等.Matlab數(shù)學(xué)計(jì)算范例教程［M］.北京:中國鐵道出版社，2004.

［4］張智星.MATLAB程序設(shè)計(jì)與應(yīng)用［M］.北京:清華大學(xué)出版社，2002:4.

［5］王 劍，張書畢，史先領(lǐng)，等.MATLAB工具箱在測繪數(shù)據(jù)處理中的應(yīng)用［J］.四川測繪，2006，29(1):17-20.

［6］毛濤濤，王正林，王 玲，等.精通MATLABGUI設(shè)計(jì)［M］.北京:電子工業(yè)出版社，2008.