c语言实现matlab函数,C语言实现matlab的interp2()函数

项目要用到matlab中的Vq = interp2(X,Y,V,Xq,Yq)函数，即把一个已知经纬度和对应值的矩阵，插值变换到一个给定经纬度网格中，也就是对给定网格填值，需要用到插值，这里使用双线性内插法。

*(这只是一个初步完成代码，仅供参考)

以下是对应C代码和测试程序：

//************************************

// 函数名称: inter_linear()

// 函数说明：计算两点之间某一给定点的值，(x0,y0)->(x1,y1),已知x(x0

// 返 回 值: double

// 参 数: double x0,y0,x1,y1/*in*/ 两个点的坐标(x0,y0)(x1,y1)

// 作 者：WSS

// 作成日期：2019/09/05

//************************************

double inter_linear(double x0, double y0, double x1, double y1, double x)

{

double a0, a1, y;

a0 = (x - x1) / (x0 - x1);

a1 = (x - x0) / (x1 - x0);

y = a0*y0 + a1*y1;

return y;

}

//************************************

// 函数名称: interp2d()

// 函数说明：二维插值，同matlab的interp2()功能

// 返 回 值: double

// 参 数: x,y分别为长度为m和n的向量(一维数组)，z为矩阵(对应的二维数组(m,n))

// a,b分别为长度为asize和bsize的向量(一维数组)，out_result为矩阵(对应的二维数组(asize,bsize))

// 作 者：WSS

// 作成日期：2019/09/05

//************************************

int interp2d(double *x, double *y, double *z, int m, int n, double *a, double *b, int asize,int bsize,double *out_result)

{

//a,b是待插值的矩阵行和列，遍历a，取出行，遍历b取出列，将插值计算后的值存储在out_result中

double pointa = ;

double pointb = ;

double w1, w2,w;

const int nu_val = -;

int tempi = ;

int tempj =;

for (int i = ; i < asize; ++i)

{

for (int j = ; j <