常用的插值:拉格朗日多项式插值、牛顿插值、分段线性插值、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. 25个python相关的基础概念总结
  2. mybatis mysql查询树形结构_MyBatis collection 集合嵌套查询树形节点
  3. 151. Reverse Words in a String
  4. SonarQube搭建和使用教程
  5. python中lambda的用法
  6. 怎么让打印出来表格的纸下面显示页数_4个使用频率较高的Excel打印技巧,办公室打印文件时常常会用到!...
  7. 【R】数据框data.frame的基本操作【1】
  8. WP8.1学习系列(第二章)——Toast通知
  9. 常见面试算法:Logistic回归、树回归
  10. EL表达式+JSTL,forEach的两种用法
  11. html组件做成圆角,css3制作圆角按钮
  12. Modelica学习笔记1
  13. java上传图片裁剪_上传图像之裁剪
  14. 分享下nas使用技巧,以及最简单的远程开机(用HS-200S做演示)
  15. 玻纤效应对skew的影响(一)
  16. Minecraft 1.12.2模组开发(十八) 自定义附魔
  17. 【转】manifest
  18. computed vue 不 触发_vue computed 无法触发 set
  19. mediasoup中nack的调用机制
  20. UWB的定位算法(简单详细易懂)

热门文章

  1. PyCharm Python迁移项目
  2. 精通SpringBoot---整合RabbitMQ消息队列
  3. MySQL 字符串转in/double类型——CAST/CONVERT函数的用法
  4. Oracle对表空间、用户、用户权限的操作
  5. The 6th Zhejiang Provincial Collegiate Programming Contest-ProblemA:Second-price Auction
  6. [Android]BaseExpandableListAdapter实现可折叠列表
  7. win32api window2con 模块 -系统注册表操作
  8. 转:谈谈Unicode编码,简要解释UCS、UTF、BMP、BOM等名词
  9. asp.net 初步入门使用正则抓取网页信息
  10. Mr.J--Bootstrap使用