如何用MATLAB 编写interpl插值函数

yi = interp1(x,y,xi,method)

已知样本点坐标x,y，求xi处的函数值yi，插值方法是method。

function [输出变量] = 函数名称(输入变量）

% 注释

% 函数体

如下所示，是编写的一个求1到n之和的求和函数 eg_sum，按照上述格式，编写代码如下并保存文件，注意文件命名时不能以数字开头：

function [s]  = eg_sum( n )

% calculate the sum of 1..n

s = 0

for i = 1:n

s = s + i

end

回到控制台Command Window（主界面）,注意要使用之前保存的函数文件名来调用函数，如下所示：

>>clear

>>a = eg_sum(10)

运行结果： a = 55

扩展资料：

对于函数f(x)，常常不仅知道它在一些点的函数值，而且还知道它在这些点的导数值。这时的插值函数P(x)，自然不仅要求在这些点等于f(x)的函数值，而且要求P(x)的导数在这些点也等于f(x)的导数值。这就是埃尔米特插值问题，也称带导数的插值问题。

从几何上看，这种插值要寻求的多项式曲线不仅要通过平面上的已知点组，而且在这些点（或者其中一部分）与原曲线“密切”，即它们有相同的斜率。可见埃尔米特插值多项式比起一般多项式插值有较高的光滑逼近要求。

参考资料来源：百度百科-插值

a=[3.946 05.193 505.762 1006.311 1506.795 2007.257 2507.704 3008.170 350...

8.603 4009.073 4509.511 50010.001 55010.435 60010.932 65011.400 700...

11.938 75012.428 80012.995 85013.560 90014.246 95015.063 1000]

x=a(:,2)

y=a(:,1)

x2=1:60:1200

y2=interp1(x,y,x2, 'pchip')

plot(x,y,'bo',x2,y2,'r*')

legend('千分尺数据','分段插值数据','Location', 'Best')

grid on

xlabel('x 容量值(L)'),ylabel('y 电流值（mA）')

---