一、概念与应用

1、概念

插值法又称“内插法”,是利用函数f (x)在某区间中已知的若干点的函数值,作出适当的特定函数,在区间的其他点上用这特定函数的值作为函数f (x)的近似值,这种方法称为插值法。如果这特定函数是多项式,就称它为插值多项式(百度百科)

插值法作用——有效预测未知点

已知点近似计算未知点的近似计算方法,即构造一个多项式函数,使其通过所有已知点,然后用求得的函数预测位置点

讲解matlab中插值函数的

2、应用

    常常需要根据已知的函数点进行数据、模型的处理和分析,而通常情况下现有的数据是极少的,不足以支撑分析的进行,这时就需要使用一些数学的方法,“模拟产生”一些新的但又比较靠谱的值来满足需求。

二、matlab插值函数

MATLAB中的插值函数为interp1,

语法:  yi= interp1(x,y,xi,'method')

其中x,y为插值点,yi为在被插值点xi处的插值结果;x,y为向量, 'method'表示采用的插值方法,MATLAB提供的插值方法有几种: 'method'是最邻近插值, 'linear'线性插值; 'spline'三次样条插值; 'cubic'三次插值.缺省时表示线性插值

【注意:所有的插值方法都要求x是单调的,并且xi不能够超过x的范围,对于非单调x,采用x,y分别构造关于t的参数方程】

【注意:所有的插值方法都要求x是单调的,并且xi不能够超过x的范围,对于非单调x,采用x,y分别构造关于t的参数方程】

【注意:所有的插值方法都要求x是单调的,并且xi不能够超过x的范围,对于非单调x,采用x,y分别构造关于t的参数方程】

(重要事情说3遍)

例如:在一 天24小时内,从零点开始每间隔2小时测得的环境温度数据分别为

12,9,9,10,18 ,24,28,27,25,20,18,15,13,

推测中午12点(即13点)时的温度.

x=0:2:24;
       y=[12   9   9   10   18  24   28   27   25   20  18  15  13];

a=13;
      y1=interp1(x,y,a,'spline')

结果为:  27.8725

若要得到一天24小时的温度曲线,则:

xi=0:1/3600:24;

yi=interp1(x,y,xi, 'spline');

plot(x,y,'o' ,xi,yi)

命令1 interp1
功能 一维数据插值(表格查找)。该命令对数据点之间计算内插值。它找出一元函数f(x)在中间点的数值。其中函数f(x)由所给数据决定。
x:原始数据点
Y:原始数据点
xi:插值点
Yi:插值点
格式
(1)yi = interp1(x,Y,xi) 
返回插值向量yi,每一元素对应于参量xi,同时由向量x 与Y 的内插值决定。参量x 指定数据Y 的点。
若Y 为一矩阵,则按Y 的每列计算。yi 是阶数为length(xi)*size(Y,2)的输出矩阵。
(2)yi = interp1(Y,xi) 
假定x=1:N,其中N 为向量Y 的长度,或者为矩阵Y 的行数。
(3)yi = interp1(x,Y,xi,method) 
用指定的算法计算插值:
’nearest’:最近邻点插值,直接完成计算;
’linear’:线性插值(缺省方式),直接完成计算;
’spline’:三次样条函数插值。对于该方法,命令interp1 调用函数spline、ppval、mkpp、umkpp。这些命令生成一系列用于分段多项式操作的函数。命令spline 用它们执行三次样条函数插值;
’pchip’:分段三次Hermite 插值。对于该方法,命令interp1 调用函数pchip,用于对向量x 与y 执行分段三次内插值。该方法保留单调性与数据的外形;
’cubic’:与’pchip’操作相同;
’v5cubic’:在MATLAB 5.0 中的三次插值。

-------------------------------------------详解-----------------------------------------

一、pchip全称(piecewise cubic hermite interpolation polynomial)分段三次艾米特插值多项式

infomation : matlab 指令 help pchip

语法:

yi = pchip(x,y,xi)——返回与 xi中的查询点对应的插值 p 的向量。p 的值由 x 和 y 的保形分段三次插值确定。

输入参数:

x,y为采样点,维度相同,xi为内插点,是 pchip 计算出的插值函数值 p 的 x 坐标

返回:yi为xi对应的输出点

pp = pchip(x,y)——返回一个分段多项式结构体以用于 ppval 和样条实用工具 unmkpp

输入参数:

x,y为采样点,维度相同

返回:pp结构体

以下是help pchip信息

example

x = -3:3;
y = [-1 -1 -1 0 1 1 1];
t = -3:.01:3;
p = pchip(x,y,t);
s = spline(x,y,t);
plot(x,y,'o',t,p,'-',t,s,'-.')
legend('data','pchip','spline','Location','SouthEast')

pchip与spline的区别:

  • 样条插值曲线产生光滑的结果,如二阶导连续
  • 样条插值曲线产生更精确的结果,如果数据包含光滑函数的值
  • pchip没有超调量,如果数据不光滑,更少的震荡
  • pchip建立起来更少代价
  • 两种方法评估代价相等

源解释:

曲线坐标系与直角坐标系转换(一)——基础:matlab插值函数简介相关推荐

  1. OpenCV的图像直角坐标系转极坐标系的函数warpPolar()详解,并附自己写的实现直角坐标系转极坐标系的MATLAB代码

    [草稿]详解OpenCV实现图像直角坐标系转极坐标系的函数warpPolar(),并附自己写的实现直角坐标系转极坐标系的MATLAB代码 有些时候我们需要把图像或矩阵从直角坐标系(笛卡尔坐标系)转换到 ...

  2. 东北大学数据科学基础(MATLAB)-笔记

    (本地用typora写的md笔记,导入图片得建一个picGo图床,嫌麻烦,所以如有需要带完整图片的笔记,参见: 课件,参见:东北大学数据科学基础MATLAB学习PPT.zip) 文章目录 一. 零散知 ...

  3. 【图像转换】基于matlab灰度图像转换彩色图像【含Matlab 1233期】

    一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[图像转换]基于matlab灰度图像转换彩色图像[含Matlab 1233期] (https://download.csdn.net/dow ...

  4. 【自学】零基础MATLAB入门笔记

    [自学]零基础MATLAB入门笔记 (一) 第一章 MATLAB的新建工程和基本操作与处理 文章目录 [自学]零基础MATLAB入门笔记 (一) 前言 一.MATLAB简介 二.软件下载 三.新建工程 ...

  5. 视频教程-三十八课时零基础matlab精通优化算法-Matlab

    三十八课时零基础matlab精通优化算法 图像和算法等领域有多年研究和项目经验:指导发表科技核心期刊经验丰富:多次指导数学建模爱好者参赛. 宋星星 ¥100.00 立即订阅 扫码下载「CSDN程序员学 ...

  6. MATLAB基础-matlab的搜索路径

    MATLAB基础-matlab的搜索路径 matlab的搜索路径为:变量->内部函数->程序文件(当前文件目录优先于搜索路径文件目录) 变量优先于内部函数:当变量名与内部函数重复时,优先调 ...

  7. 《对冲基金建模与分析基于MATLAB》简介及PDF下载

    转 <对冲基金建模与分析--基于MATLAB>简介及PDF下载 内容简介 本书是关于用MATLAB对对冲基金进行建模和分析的入门读物.在对对冲基金的基本概念.分类.相关工具和指标系统介绍的 ...

  8. MATLAB R2018a简介

    MATLAB R2018a简介 MathWorks及其产品概述 MATLAB是矩阵实验室(Matrix Laboratory)的简称,是一种用于算法开发.数据可视化.数据分析及数值计算的高级技术计算语 ...

  9. TypeScript基础入门 - 函数 - 简介

    2019独角兽企业重金招聘Python工程师标准>>> 转载 TypeScript基础入门 - 函数 - 简介 项目实践仓库 https://github.com/durban89/ ...

  10. wkhtmltopdf:wkhtmltopdf(将html转换成pdf的利器)简介、安装、使用方法详细攻略

    wkhtmltopdf:wkhtmltopdf(将html转换成pdf的利器)简介.安装.使用方法详细攻略 目录 wkhtmltopdf简介 wkhtmltopdf安装 wkhtmltopdf使用方法 ...

最新文章

  1. ZooKeeper 基本介绍
  2. mysql主备数据库配置_MySQL双主互备配置
  3. 泛函编程(19)-泛函库设计-Parallelism In Action
  4. 逻辑回归原理(python代码实现)
  5. html ios 定位吗,iOS对html简单解析
  6. 'adb' 不是内部或外部命令,也不是可运行的程序或批处理文件
  7. SAP Spartacus避免不能向下兼容的breaking changes
  8. Parhaps you are running on a JRE rather than a JDK?
  9. Matlab协方差矩阵的计算原理
  10. Python 从入门到进阶
  11. python3.8使用aiml总结
  12. 训练loss不下降原因总结
  13. 360 极速模式 css,请问 chrome 和 360 极速模式 渲染页面不一样
  14. Word2010页眉添加两条横线
  15. 教你炒股票25:吻,MACD、背弛、中枢
  16. HDLC(High-Level Data Link Control)(对比以太网)
  17. OTA三种差分算法安装使用。bsdiff,hdiffpatch,xdelta3
  18. 推荐一款强大的在线编译器
  19. 流行的ORM框架简介
  20. 游戏开发中常用的数学知识---矩阵(一)

热门文章

  1. 国产高清卫星影像时代来了,打造中国版“谷歌地球”!
  2. java线程安全的计数器_Java线程安全计数器的简单代码示例
  3. Builder中使用Access数据库
  4. (安徽省地图)行政区划图高清矢量cdr|pdf(详细版2021年)
  5. 空间计量经济学 matlab,空间计量经济学基于MATLAB的应用分析
  6. fastreport java 集成_报表生成器配置FastReport.Net环境
  7. ActivityManager Displayed 源码位置
  8. 网站丨平淡的生活里增添一点幸福感
  9. mcafee完全彻底卸载方法
  10. Odin WindowEditor使用体会