Matlab求解插值问题.doc

Matlab求解插值问题.doc

ID:55953538

大小:157.50 KB

页数:5页

时间:2020-06-18

Matlab求解插值问题.doc_第1页
Matlab求解插值问题.doc_第2页
Matlab求解插值问题.doc_第3页
Matlab求解插值问题.doc_第4页
Matlab求解插值问题.doc_第5页
资源描述:

《Matlab求解插值问题.doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、Matlab求解插值问题在应用领域中,由有限个已知数据点,构造一个解析表达式,由此计算数据点之间的函数值,称之为插值。实例:海底探测问题某公司用声纳对海底进行测试,在5×5海里的坐标点上测得海底深度的值,希望通过这些有限的数据了解更多处的海底情况。并绘出较细致的海底曲面图。1、一元插值一元插值是对一元数据点(xi,yi)进行插值。线性插值:由已知数据点连成一条折线,认为相临两个数据点之间的函数值就在这两点之间的连线上。一般来说,数据点数越多,线性插值就越精确。调用格式:yi=interp1(x,y,xi,’linear’)%线性插值zi=inte

2、rp1(x,y,xi,’spline’)%三次样条插值wi=interp1(x,y,xi,’cubic’)%三次多项式插值说明:yi、zi、wi为对应xi的不同类型的插值。x、y为已知数据点。例:已知数据:x0.1.2.3.4.5.6.7.8.91y.3.511.41.61.9.6.4.81.52求当xi=0.25时的yi的值。程序:x=0:.1:1;y=[.3.511.41.61.6.4.81.52];yi0=interp1(x,y,0.025,'linear')xi=0:.02:1;yi=interp1(x,y,xi,'linear');zi

3、=interp1(x,y,xi,'spline');wi=interp1(x,y,xi,'cubic');plot(x,y,'o',xi,yi,'r+',xi,zi,'g*',xi,wi,'k.-')legend('原始点','线性点','三次样条','三次多项式')结果:yi0=0.3500要得到给定的几个点的对应函数值,可用:xi=[0.25000.35000.4500]yi=interp1(x,y,xi,'spline')结果:yi=1.20881.58021.34542、二元插值二元插值与一元插值的基本思想一致,对原始数据点(x,y,z)

4、构造见上面函数求出插值点数据(xi,yi,zi)。单调节点插值函数,即x,y向量是单调的。调用格式1:zi=interp2(x,y,z,xi,yi,’linear’)‘liner’是双线性插值(缺省)调用格式2:zi=interp2(x,y,z,xi,yi,’nearest’)’nearest’是最近邻域插值调用格式3:zi=interp2(x,y,z,xi,yi,’spline’)‘spline’是三次样条插值说明:这里x和y是两个独立的向量,它们必须是单调的。z是矩阵,是由x和y确定的点上的值。z和x,y之间的关系是z(i,:)=f(x,y(

5、i)),z(:,j)=f(x(j),y)即:当x变化时,z的第i行与y的第i个元素相关,当y变化时z的第j列与x的第j个元素相关。如果没有对x,y赋值,则默认x=1:n,y=1:m。n和m分别是矩阵z的行数和列数。例2:已知某处山区地形选点测量坐标数据为:x=00.511.522.533.544.55y=00.511.522.533.544.555.56海拔高度数据为:z=8990878592919693908782929698999591898684828496989592908885848381858081828995969392898686

6、828587989996978885828382858994959392918684888892939495898786838192929697989693958482818485858182808081859093958486819899989796958487808185828384879095868880828184858683828180828788899899979698949287其地貌图为:对数据插值加密形成地貌图。程序:x=0:.5:5;y=0:.5:6;z=[89908785929196939087829296989995918

7、98684828496989592908885848381858081828995969392898686828587989996978885828382858994959392918684888892939495898786838192929697989693958482818485858182808081859093958486819899989796958487808185828384879095868880828184858683828180828788899899979698949287];mesh(x,y,z)%绘原始数据图xi=li

8、nspace(0,5,50);%加密横坐标数据到50个yi=linspace(0,6,80);%加密纵坐标数据到60个[xii,yii]

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。