![MATLAB插值實驗報告--數(shù)學(xué)實驗(共11頁)_第1頁](http://file3.renrendoc.com/fileroot_temp3/2022-2/22/2907a179-5a8f-4f50-84fb-94de9b25971d/2907a179-5a8f-4f50-84fb-94de9b25971d1.gif)
![MATLAB插值實驗報告--數(shù)學(xué)實驗(共11頁)_第2頁](http://file3.renrendoc.com/fileroot_temp3/2022-2/22/2907a179-5a8f-4f50-84fb-94de9b25971d/2907a179-5a8f-4f50-84fb-94de9b25971d2.gif)
![MATLAB插值實驗報告--數(shù)學(xué)實驗(共11頁)_第3頁](http://file3.renrendoc.com/fileroot_temp3/2022-2/22/2907a179-5a8f-4f50-84fb-94de9b25971d/2907a179-5a8f-4f50-84fb-94de9b25971d3.gif)
![MATLAB插值實驗報告--數(shù)學(xué)實驗(共11頁)_第4頁](http://file3.renrendoc.com/fileroot_temp3/2022-2/22/2907a179-5a8f-4f50-84fb-94de9b25971d/2907a179-5a8f-4f50-84fb-94de9b25971d4.gif)
![MATLAB插值實驗報告--數(shù)學(xué)實驗(共11頁)_第5頁](http://file3.renrendoc.com/fileroot_temp3/2022-2/22/2907a179-5a8f-4f50-84fb-94de9b25971d/2907a179-5a8f-4f50-84fb-94de9b25971d5.gif)
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)
文檔簡介
1、精選優(yōu)質(zhì)文檔-傾情為你奉上新鄉(xiāng)學(xué)院數(shù)學(xué)與信息科學(xué)系實驗報告實驗項目名稱 插值實驗 所屬課程名稱 數(shù)學(xué)實驗 實 驗 類 型 綜合性實驗 實 驗 日 期 班 級 學(xué) 號 姓 名 成 績 一、實驗概述:【實驗?zāi)康摹空莆沼肕atlab插值的方法,了解拉格朗日插值、線性插值、樣條插值的基本思想,了解三種網(wǎng)格節(jié)點數(shù)據(jù)的插值方法的基本思想,了解掌握用Matlab計算一維差值和二維插值的方法?!緦嶒炘怼?.一位握又性插值、樣條插值的基本思想,了解三種網(wǎng)格節(jié)點數(shù)據(jù)的差222222222222222222222222222222222222222222222222222222222222222222222222
2、222222222222222222222222 拉格朗日(Lagrange)插值。已知函數(shù)f(x)在n+1個點x0,x1,xn處的函數(shù)值為 y0,y1,yn 求一n次多項式函數(shù)Pn(x),使其滿足: Pn(xi)=yi,i=0,1,n.解決此問題的拉格朗日插值多項式公式如下其中Li(x) 為n次多項式:稱為拉格朗日插值基函數(shù)【實驗環(huán)境 】Matlab 7.0.1 Windows 8 專業(yè)版 二、 實驗內(nèi)容:【實驗方案】1.利用interp2命令,以及畫圖命令。2.會用各種插值命令。3.熟悉并掌握拉格朗日插值?!緦嶒炦^程】(實驗步驟、記錄、數(shù)據(jù)、分析)1.山區(qū)地貌:在某山區(qū)測得一些地點的高程如
3、下表:(平面區(qū)域1200x 4000,1200y 3600),試作出該山區(qū)的地貌圖和等高線圖,并對幾種插值方法進(jìn)行比較 36003200280024002000160012001480 1500 1550 1510 1430 1300 1200 9801500 1550 1600 1550 1600 1600 1600 15501500 1200 1100 1550 1600 1550 1380 10701500 1200 1100 1350 1450 1200 1150 10101390 1500 1500 1400 900 1100 1060 9501320 1450 1420 1400
4、1300 700 900 8501130 1250 1280 1230 1040 900 500 700y/x1200 1600 2000 2400 2800 3200 3600 4000解:原始圖程序:x=1200:400:4000;y=1200:400:3600;z=1480 1500 1550 1510 1430 1300 1200 980;1500 1550 1600 1550 1600 1600 1600 1550;1500 1200 1100 1550 1600 1550 1380 1070;1500 1200 1100 1350 1450 1200 1150 1010;1390
5、1500 1500 1400 900 1100 1060 950;1320 1450 1420 1400 1300 700 900 850;1130 1250 1280 1230 1040 900 500 700;mesh(x,y,z);title('原始圖'),xlabel('x'),ylabel('y'),zlabel('z') 最鄰近插值圖:xi,yi=meshgrid(1200:10:4000,1200:10:3600);zi=interp2(x,y,z,xi,yi,'nearest');mesh(xi,y
6、i,zi)title('最鄰近插值圖'),xlabel('x'),ylabel('y'),zlabel('z') 線性插值:xi,yi=meshgrid(1200:10:4000,1200:10:3600);zi=interp2(x,y,z,xi,yi,'linear');mesh(xi,yi,zi)title('線性插值'),xlabel('x'),ylabel('y'),zlabel('z') 三次樣條插值:xi,yi=meshgrid(1200:
7、10:4000,1200:10:3600);zi=interp2(x,y,z,xi,yi,'spline');mesh(xi,yi,zi)title('三次樣條插值'),xlabel('x'),ylabel('y'),zlabel('z') 立方差值:xi,yi=meshgrid(1200:10:4000,1200:10:3600);zi=interp2(x,y,z,xi,yi,'cubic');mesh(xi,yi,zi)title('立方插值'),xlabel('x
8、9;),ylabel('y'),zlabel('z') 等高線圖:x=1200:400:4000;y=1200:400:3600;z=1480 1500 1550 1510 1430 1300 1200 980;1500 1550 1600 1550 1600 1600 1600 1550;1500 1200 1100 1550 1600 1550 1380 1070;1500 1200 1100 1350 1450 1200 1150 1010;1390 1500 1500 1400 900 1100 1060 950;1320 1450 1420 1400
9、1300 700 900 850;1130 1250 1280 1230 1040 900 500 700;meshc(x,y,z);title('等高線'),xlabel('x'),ylabel('y'),zlabel('z') 2. 針對這兩個函數(shù),分別用用Lagrange插值選取11個基點計算插值(ych),并把插值函數(shù)和原函數(shù)都畫在一個坐標(biāo)系內(nèi)(注意是否有龍格震蕩現(xiàn)象)。解:M文件: function y=lag(x0,y0,x)for i=1:length(x) s=0; for k=1:length(x0) l=1;
10、for j=1:length(x0) if j=k l=l*(x(i)-x0(j)/(x0(k)-x0(j); end end s=l*y0(k)+s; end y(i)=s;end命令:x0=-6:1.2:6;y01=1./(1+x0.2);y02=1./(1+25*x0.2);x=-6:0.01:6;y=lag(x0,y01,x);z=lag(x0,y02,x);y1=1./(1+x.2);y2=1./(1+25*x.2);subplot(1,2,1)plot(x,y),title('g(x)=1/(1+x2)')hold onplot(x,y1,'r')subplot(1,2,2)plot(x,z),title(&
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 2025年蛋撻皮合作協(xié)議書
- 2025年無機(jī)械動力飛機(jī)合作協(xié)議書
- 2025年九年級下學(xué)期語文教學(xué)工作總結(jié)標(biāo)準(zhǔn)范文(二篇)
- 2025年中山市店鋪出租合同(2篇)
- 2025年中小學(xué)走讀生安全責(zé)任協(xié)議模板(三篇)
- 2025年二年級教師心得體會例文(6篇)
- 2013-2022年北京市中考真題物理試題匯編:磁現(xiàn)象章節(jié)綜合
- 2025年個人客戶信息保密協(xié)議范文(2篇)
- 倉儲裝修終止協(xié)議樣本
- 文化產(chǎn)業(yè)基地裝修合同
- HYT 235-2018 海洋環(huán)境放射性核素監(jiān)測技術(shù)規(guī)程
- 中國香蔥行業(yè)市場現(xiàn)狀分析及競爭格局與投資發(fā)展研究報告2024-2034版
- 消化系統(tǒng)常見疾病康復(fù)
- 婦科惡性腫瘤免疫治療中國專家共識(2023)解讀
- 2024年浪潮入職測評題和答案
- 小班數(shù)學(xué)《整理牛奶柜》課件
- 皮膚感染的護(hù)理診斷與護(hù)理措施
- 中考語文真題雙向細(xì)目表
- 2024年江蘇省對口單招英語試卷及答案
- 藥品集采培訓(xùn)課件
- 高中物理考試成績分析報告
評論
0/150
提交評論