常用的插值:拉格朗日多项式插值、牛顿插值、分段线性插值、Hermite 插值和三次样条插值。

1.拉格朗日插值法

function y=lagrange(x0,y0,x);

n=length(x0);m=length(x);for i=1:m

z=x(i);

s=0.0;for k=1:n

p=1.0;for j=1:nif j~=k

p=p*(z-x0(j))/(x0(k)-x0(j));

end

end

s=p*y0(k)+s;

end

y(i)=s;

end

测试:

>> x0=1:5;

>> y0=x0.^2;

>> x=2.5;

>> y=lagrange(x0,y0,x)

y =

6.2500

>> plot(x0,y0)

>> hold on

>> plot(x,y,‘*‘)

2.牛顿插值

function yi=newton(x,y,xi)

n=length(x);

m=length(y);if n~=m

error(‘The lengths of X ang Y must be equal!‘);return;

end

Y=zeros(n);

Y(:,1)=y‘;

for k=1:n-1

for i=1:n-kif abs(x(i+k)-x(i))

error(‘the DATA is error!‘);return;

end

Y(i,k+1)=(Y(i+1,k)-Y(i,k))/(x(i+k)-x(i));

end

end

yi=0;for i=1:n

z=1;for k=1:i-1z=z*(xi-x(k));

end

yi=yi+Y(1,i)*z;

end

测试:

x0=1:5;

y0=x0.^2;

x=2.5;

y=newton(x0,y0,x)

plot(x0,y0)

hold on

plot(x,y,‘*‘)

y =

6.2500

3.分段线性插值

y=interp1(x0,y0,x,‘method‘)

method 指定插值的方法,默认为线性插值。其值可为:

‘nearest‘ 最近项插值

‘linear‘ 线性插值

‘spline‘ 逐段3 次样条插值

‘cubic‘ 保凹凸性3 次插值。

所有的插值方法要求 x0 是单调的。

测试:

x0=1:5;

y0=x0.^2;

x=2.5;

y=interp1(x0,y0,x,‘linear‘)

plot(x0,y0)

hold on

plot(x,y,‘*‘)

y =

6.5000

4. 埃尔米特(Hermite)插值

function y=hermite(x0,y0,y1,x); %y1是导数值

n=length(x0);m=length(x);for k=1:m

yy=0.0;for i=1:n

h=1.0;

a=0.0;for j=1:nif j~=i

h=h*((x(k)-x0(j))/(x0(i)-x0(j)))^2;

a=1/(x0(i)-x0(j))+a;

end

end

yy=yy+h*((x0(i)-x(k))*(2*a*y0(i)-y1(i))+y0(i));

end

y(k)=yy;

end

测试:

>> x0=1:5;

y0=x0.^2;

x=2.5;

y1=[1 1 1 1 1];

y=hermite(x0,y0,y1,x)

plot(x0,y0)

hold on

plot(x,y,‘*‘)

y =

7.4228

5.样条插值

Matlab 中三次样条插值也有现成的函数:

(1)y=interp1(x0,y0,x,‘spline‘);

(2)y=spline(x0,y0,x);

(3)pp=csape(x0,y0,conds);

y=ppval(pp,x)

说明:

csape 的返回值是pp 形式,要求插

值点的函数值,必须调用函数ppval。

pp=csape(x0,y0):使用默认的边界条件,即Lagrange 边界条件。

pp=csape(x0,y0,conds)中的conds 指定插值的边界条件,其值可为:

‘complete‘ 边界为一阶导数,即默认的边界条件

‘not-a-knot‘ 非扭结条件

‘periodic‘ 周期条件

‘second‘ 边界为二阶导数,二阶导数的值[0, 0]。

‘variational‘ 设置边界的二阶导数值为[0,0]。

测试:

x0=1:5;

y0=x0.^2;

x=2.5;

y1=interp1(x0,y0,x,‘spline‘)

y2=spline(x0,y0,x)

pp=csape(x0,y0,‘second‘)

y3=ppval(pp,x)

plot(x0,y0)

hold on

plot(x,y1,‘*‘)

plot(x,y2,‘*‘)

plot(x,y3,‘*‘)

y1 =

6.2500

y2 =

6.2500

pp =

包含以下字段的 struct:

form: ‘pp‘

breaks: [1 2 3 4 5]

coefs: [4×4 double]

pieces: 4

order: 4

dim: 1

y3 =

6.2321

6.二维插值

一维插值:节点为一维变量,插值函数是一元函数(曲线)。

二维插值:节点是二维的,插值函数就是二元函数,即曲面。

6.1插值节点为网格节点

(1)z=interp2(x0,y0,z0,x,y,‘method‘)

(2)pp=csape({x0,y0},z0,conds,valconds);  %三次样条插值

z=fnval(pp,{x,y})

6.2 插值节点为散乱节点

原文:http://www.cnblogs.com/youngsea/p/7498743.html

matlab实现查值,[数学建模(六)]使用MATLAB实现插值相关推荐

  1. 一个实例用matlab分别画出拟合和插值,数学建模培训之四--拟合与插值专题(hzd).doc...

    您所在位置:网站首页 > 海量文档 &nbsp>&nbsp高等教育&nbsp>&nbsp微积分 数学建模培训之四--拟合与插值专题(hzd).doc6 ...

  2. 【数学建模】基于matlab武汉地铁2号线路线地图动态模拟【含Matlab源码 1092期】

    一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[数学建模]基于matlab武汉地铁2号线路线地图动态模拟[含Matlab源码 1092期] 点击上面蓝色字体,直接付费下载,即可. 获取代 ...

  3. 【数学建模】基于matlab船舶三自由度MMG模型【含Matlab源码 1925期】

    ⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[数学建模]基于matlab船舶三自由度MMG模型[含Matlab源码 1925期] 点击上面蓝色字体,直接付费下载,即可. 获取代码方式 ...

  4. 数学建模专栏 | 开篇:如何备战数学建模竞赛之 MATLAB 编程

    作 者 简 介 卓金武,MathWorks中国高级工程师,教育业务经理,在数据分析.数据挖掘.机器学习.数学建模.量化投资和优化等科学计算方面有多年工作经验,现主要负责MATLAB校园版业务.曾2次获 ...

  5. 【数学建模】基于matlab动态水波仿真【含Matlab源码 2056期】

    ⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[数学建模]基于matlab动态水波仿真[含Matlab源码 2056期] 点击上面蓝色字体,直接付费下载,即可. 获取代码方式2: 付费 ...

  6. 【数学建模】基于matlab zernike泽尼克多项式仿真【含Matlab源码 1953期】

    一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[数学建模]基于matlab zernike泽尼克多项式仿真[含Matlab源码 1953期] 点击上面蓝色字体,直接付费下载,即可. 获取 ...

  7. 【数学建模】基于matlab三维海浪模型仿真【含Matlab源码 1159期】

    一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[数学建模]基于matlab三维海浪模型仿真[含Matlab源码 1159期] 点击上面蓝色字体,直接付费下载,即可. 获取代码方式2: 付 ...

  8. 数学建模——一维、二维插值模型详解Python代码

    数学建模--一维.二维插值模型详解Python代码 一.一维插值 # -*-coding:utf-8 -*- import numpy as np from scipy import interpol ...

  9. 一篇文章带你认识数学建模中的二维插值

    本篇文章主要是认识数学建模中的二维插值 二维插值的具体计算可参考:一篇文章带你搞定二维插值的 MATLAB 计算 文章目录 一.引言 二.网格节点的概念 三.散乱节点的概念 四.二维插值的常见方法 1 ...

最新文章

  1. 描述一下Spring框架的作用和优点?
  2. java面试题(开发框架)
  3. 通过checkbox选择以逗号拼接删除字符串
  4. Axis 实现的SOAP附件的传输
  5. java 时间l_java中获取系统时间并格式化输出
  6. unity图片拖不进去_Unity UGUI实现简单拖拽图片功能
  7. python装饰器的使用
  8. DevExpress v18.1新版亮点——WinForms篇(六)
  9. 网奇iwms插件之“我浏览过的文章”
  10. 网页服务器目录属性,IIS 6.0 F1:网站属性 - 目录安全性选项卡
  11. jquery初级视频教程
  12. 广播系统可以跨服务器控制吗,数字广播系统技术参数和要求.doc
  13. Scrapy爬取二手房信息
  14. Quartus 在工具栏设置功能快捷方式
  15. mysql网络投票系统设计_学生在线投票表决系统设计与开发(JSP,MySQL)
  16. HarmonyOS无法添加outlook日历(报错该账户名称已被使用)
  17. (js队列,堆栈) (FIFO,LIFO)
  18. 骑行318、 2016.7.15
  19. 干货分享 | FMEA何时做?谁来做?
  20. C程序设计(谭浩强)--学习笔记

热门文章

  1. 项目Beta冲刺Day3
  2. 关于sqoop增量导入oracle的数据到HDFS的出现的错误总结
  3. 使用DOM操纵样式表
  4. Centos6.x服务器配置jdk+tomcat+mysql环境
  5. 【缅怀妈妈系列诗歌】之八:妈妈,我不会忘记
  6. java基础—计算: 2012-3-17到2012-4-6中间有多少天?
  7. 单页面与多页面的优缺点
  8. Django分页器和自定义分页器
  9. python练习题-day2
  10. ORA-12519: TNS:no appropriate service handler found 解决方法