MATLAB插值函数interp1
来源:互联网 发布:win7 网络打印机脱机 编辑:程序博客网 时间:2024/05/26 07:28
插值法
插值法又称“内插法”,是利用函数f (x)在某区间中已知的若干点的函数值,作出适当的特定函数,在区间的其他点上用这特定函数的值作为函数f (x)的近似值,这种方法称为插值法。如果这特定函数是多项式,就称它为插值多项式。
线性插值法
线性插值法是指使用连接两个已知量的直线来确定在这两个已知量之间的一个未知量的值的方法。
假设我们已知坐标(x0,y0)与(x1,y1),要得到[x0,x1]区间内某一位置x在直线上的值。根据图中所示,我们得到两点式直线方程:
假设方程两边的值为α,那么这个值就是插值系数—从x0到x的距离与从x0到x1距离的比值。由于x值已知,所以可以从公式得到α的值:
同样:
这样,在代数上就可以表示成为:
y = (1 − α)y0 + αy1
或者,
y = y0 + α(y1 − y0)
这样通过α就可以直接得到 y。实际上,即使x不在x0到x1之间并且α也不是介于0到1之间,这个公式也是成立的。在这种情况下,这种方法叫作线性外插—参见 外插值。
已知y求x的过程与以上过程相同,只是x与y要进行交换。
MATLAB中的使用
%{ MATLAB中的插值函数为interp1,其调用格式为: yi= interp1(x,y,xi,'method') 其中x,y为插值点,yi为在被插值点xi处的插值结果;x,y为向量, 'method'表示采用的插值方法,MATLAB提供的插值方法有几种: 'nearest'是最邻近插值, 'linear'线性插值; 'spline'三次样条插值; 'pchip'立方插值.缺省时表示线性插值 注意:所有的插值方法都要求x是单调的,并且xi不能够超过x的范围。%}x = 0:2*pi; y = sin(x); xx = 0:0.5:2*pi; % interp1对sin函数进行分段线性插值,调用interp1的时候,默认的是分段线性插值 y1 = interp1(x,y,xx,'linear'); subplot(2,2,1);plot(x,y,'o',xx,y1,'r') title('分段线性插值') % 临近插值 y2 = interp1(x,y,xx,'nearest'); subplot(2,2,2);plot(x,y,'o',xx,y2,'r'); title('临近插值') %球面线性插值 y3 = interp1(x,y,xx,'spline'); subplot(2,2,3);plot(x,y,'o',xx,y3,'r') title('球面插值') %三次多项式插值法 y4 = interp1(x,y,xx,'pchip'); subplot(2,2,4);plot(x,y,'o',xx,y4,'r'); title('三次多项式插值')
一个例子
%{ MATLAB中的插值函数为interp1,其调用格式为: yi= interp1(x,y,xi,'method') 其中x,y为插值点,yi为在被插值点xi处的插值结果;x,y为向量, 'method'表示采用的插值方法,MATLAB提供的插值方法有几种: 'nearest'是最邻近插值, 'linear'线性插值; 'spline'三次样条插值; 'pchip'立方插值.缺省时表示线性插值 注意:所有的插值方法都要求x是单调的,并且xi不能够超过x的范围。%}%{ 例如:在一 天24小时内,从零点开始每间隔2小时测得的环境温度数据分别为 12,9,9,1,0,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 0
- MATLAB插值函数interp1
- matlab再带的插值函数interp1的四种插值方法
- interp1 函数-------一维数据插值函数
- matlab 插值函数
- MATLAB interp1
- MATLAB interp1
- MATLAB interp1
- Matlab interp1 函数的基本用法简介
- matlab 插值函数理解
- MATLAB 牛顿插值函数
- 一维插值(interp1)和二维插值(interp2)
- [转]一维插值(interp1)和二维插值(interp2)
- 函数插值计算(Matlab程序)
- Matlab基础学习----------------------函数插值
- 【Matlab】图像插值函数interp2理解
- matlab中的interp2插值函数
- 小观matlab插值函数
- 【Matlab】图像插值函数interp2理解
- 219. Contains Duplicate II
- Go语言学习之bufio包(The way to go)
- STL之deque容器的剖析
- JAVA程序设计_图谱
- 枚举操作大全
- MATLAB插值函数interp1
- AngularJS学习总结系列(二)
- docker registry http api v2
- 165. Compare Version Numbers
- SSM框架——详细整合教程(Spring+SpringMVC+MyBatis)
- 处理老大的老大的老大的老大反馈的问题, 发布后居然意外不生效???
- 环境变量的设置
- NYOJ 247 虚拟的城市之旅
- 【Java】Java中的位运算符与、或、非、异或的使用