版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
主講:41余巧花編輯員:36林燕芳
31任婷
56繆曉丹材料收集員:34黃娟
35胡珊
36郭麗娜
42蔣匯穎數(shù)值分析期中報(bào)告
-----------第四組
插值和擬和插值已知有n+1個(gè)節(jié)點(diǎn)(xj,yj),j=0,1,…,n,其中xj互不相同,節(jié)點(diǎn)(xj,yj)可看成由某個(gè)函數(shù)y=f(x)產(chǎn)生;構(gòu)造一個(gè)相對(duì)簡(jiǎn)單的函數(shù)y=P(x);使P通過(guò)全部節(jié)點(diǎn),即P(xk)=yk,
k=0,1,…,n;用P(x)作為函數(shù)f(x)的近似。yi=interp1(x,y,xi,'method')插值方法被插值點(diǎn)插值節(jié)點(diǎn)xi處的插值結(jié)果‘nearest’
:最鄰近插值‘linear’
:線性插值;‘spline’
:三次樣條插值;‘cubic’
:立方插值。缺省時(shí):線性插值。
注意:所有的插值方法都要求x是單調(diào)的,并且xi不能夠超過(guò)x的范圍。用MATLAB作一維插值計(jì)算機(jī)床加工問(wèn)題用程控銑床加工機(jī)翼斷面的下輪廓線時(shí)每一刀只能沿x方向和y方向走非常小的一步。表3-1給出了下輪廓線上的部分?jǐn)?shù)據(jù)但工藝要求銑床沿x方向每次只能移動(dòng)0.1單位.這時(shí)需求出當(dāng)x坐標(biāo)每改變0.1單位時(shí)的y坐標(biāo)。試完成加工所需的數(shù)據(jù),畫出曲線.例:求解機(jī)床加工問(wèn)題x0=[035791112131415];y0=[01.21.72.02.12.01.81.21.01.6];x=0:0.1:15;y=interp1(x0,y0,x,'spline');plot(x0,y0,'k+',x,y,'r')gridon
要求x0,y0單調(diào);x,y可取為矩陣,或x取行向量,y取為列向量,x,y的值分別不能超出x0,y0的范圍。z=interp2(x0,y0,z0,x,y,’method’)被插值點(diǎn)插值方法插值節(jié)點(diǎn)被插值點(diǎn)的函數(shù)值‘nearest’
最鄰近插值‘linear’
雙線性插值‘cubic’
雙三次插值缺省時(shí),雙線性插值用MATLAB作網(wǎng)格節(jié)點(diǎn)數(shù)據(jù)的插值(二維)
cz
=griddata(x,y,z,cx,cy,‘method’)
要求cx取行向量,cy取為列向量。被插值點(diǎn)插值方法插值節(jié)點(diǎn)被插值點(diǎn)的函數(shù)值‘nearest’
最鄰近插值‘linear’
雙線性插值‘cubic’
雙三次插值'v4‘Matlab提供的插值方法缺省時(shí),雙線性插值用MATLAB作散點(diǎn)數(shù)據(jù)的插值計(jì)算例:航行區(qū)域的警示線某海域上頻繁地有各種噸位的船只經(jīng)過(guò)。為保證船只的航行安全,有關(guān)機(jī)構(gòu)在低潮時(shí)對(duì)水深進(jìn)行了測(cè)量,表3-8是他們提供的測(cè)量數(shù)據(jù):表3-8.水道水深的測(cè)量數(shù)據(jù)x 129.0140.0103.588.0185.5195.0105.5 y 7.5141.523.0147.022.5137.585.5 z 4868688 x 157.5107.577.081.0162.0162.0117.5 y -6.5-81.03.056.5-66.584.0-33.5 z 9988949 航行區(qū)域的警示線其中(x,y)為測(cè)量點(diǎn),z為(x,y)處的水深(英尺),水深z是區(qū)域坐標(biāo)(x,y)的函數(shù)z=z(x,y),船的噸位可以用其吃水深度來(lái)反映,分為4英尺、4.5英尺、5英尺和5.5英尺4檔。航運(yùn)部門要在矩形海域(75,200)×(-50,150)上為不同噸位的航船設(shè)置警示標(biāo)記。請(qǐng)根據(jù)測(cè)量的數(shù)據(jù)描述該海域的地貌,并繪制不同噸位的警示線,供航運(yùn)部門使用。Matlab求解x=[129140103.588185.5195105.5157.5107.57781162162117.5];y=[7.5141.52314722.5137.585.5-6.5-81356.5-66.584-33.5];z=[-4-8-6-8-6-8-8-9-9-8-8-9-4-9];cx=75:0.5:200;cy=-70:0.5:150;cz=griddata(x,y,z,cx,cy','cubic');meshz(cx,cy,cz),xlabel('X'),ylabel('Y'),zlabel('Z')figure(2),contour(cx,cy,cz,[-5-5]);gridon,holdonplot(x,y,'+')xlabel('X'),ylabel('Y')曲線擬合
已知一組(二維)數(shù)據(jù),即平面上n個(gè)點(diǎn)(xi,yi)i=1,…n,
尋求一個(gè)函數(shù)(曲線)y=f(x),
使f(x)
在某種準(zhǔn)則下與所有數(shù)據(jù)點(diǎn)最為接近,即曲線擬合得最好。
+++++++++xyy=f(x)(xi,yi)ii為點(diǎn)(xi,yi)與曲線y=f(x)的距離最常用的方法是線性最小二乘擬和多項(xiàng)式擬合對(duì)給定的數(shù)據(jù)(xj,yj),j=0,1,…,n;選取適當(dāng)階數(shù)的多項(xiàng)式,如二次多項(xiàng)式g(x)=ax^2+bx+c;使g(x)盡可能逼近(擬合)這些數(shù)據(jù),但是不要求經(jīng)過(guò)給定的數(shù)據(jù)(xj,yj);1.多項(xiàng)式f(x)=a1xm+…+amx+am+1擬合指令:a=polyfit(x,y,m)2.多項(xiàng)式在x處的值y的計(jì)算命令:
y=polyval(a,x)輸出擬合多項(xiàng)式系數(shù)a=[a1,…,am,am+1]’
(數(shù)組)輸入同長(zhǎng)度數(shù)組X,Y擬合多項(xiàng)式
次數(shù)多項(xiàng)式擬合指令即要求出二次多項(xiàng)式:中的使得:例對(duì)下面一組數(shù)據(jù)作二次多項(xiàng)式擬合x(chóng)i00.10.20.30.40.50.60.70.80.91yi-0.4471.9783.286.166.167.347.669.589.489.3011.22)計(jì)算結(jié)果:A=[-9.8108,20.1293,-0.0317]解:用多項(xiàng)式
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 人教版八年級(jí)數(shù)學(xué)上冊(cè)課時(shí)優(yōu)化訓(xùn)練-14.2.1平方差公式含答案
- 人教版數(shù)學(xué)九年級(jí)下冊(cè)全冊(cè)教案
- 《狼尾草屬牧草窖貯技術(shù)規(guī)程》編制說(shuō)明
- 2025高考物理步步高同步練習(xí)選修3第四章 原子結(jié)構(gòu)和波粒二象性粒子的波動(dòng)性和量子力學(xué)的建立含答案
- 家園合作培養(yǎng)幼兒良好行為習(xí)慣的研究
- 產(chǎn)教融合實(shí)訓(xùn)基地建設(shè)研究
- 公司負(fù)責(zé)人勞動(dòng)合同簡(jiǎn)易模板
- 共價(jià)鍵模型 隨堂檢測(cè) 下學(xué)期高二化學(xué)魯科版(2019)選擇性必修2
- 年輕教師家訪心得體會(huì)(3篇)
- 新編教師師德師風(fēng)學(xué)習(xí)心得體會(huì)800字
- 壞死性筋膜炎
- 五問(wèn)房地產(chǎn)“以舊換新”
- 新修訂公司法專題講座課件
- 幼兒園廚房安全培訓(xùn)內(nèi)容
- 膽源性肝膿腫的護(hù)理查房
- 中國(guó)加速康復(fù)外科臨床實(shí)踐指南(2021)解讀
- 商場(chǎng)開(kāi)荒保潔方案計(jì)劃
- 計(jì)量監(jiān)管基礎(chǔ)知識(shí)培訓(xùn)課件
- 醫(yī)保政策培訓(xùn)計(jì)劃制度
- 與醫(yī)保有關(guān)的信息系統(tǒng)相關(guān)材料-模板
- 綜合實(shí)踐活動(dòng)(中秋節(jié))
評(píng)論
0/150
提交評(píng)論