一段插值程序

言深深

闲来无事，写了一小段一维插值程序
其中二维数组x（n,2）有两列，第一列表示横坐标，第二列表示纵坐标；y（n-1,2）存放的是得到新点（蓝点）的坐标，第一列表示横坐标，第二列为纵坐标，n是数据个数；z是输出数组，共有2n-1个需要说明的是，我们插值得到的数值并不能做到完全在理想的曲线上面，只能说是比较近似

subroutine insert_point(x,n,z)
real*8 x(n,2),y(n-1,2),z(2*n-1,2)
y(1,1)=(x(1,1)+x(2,1))/2
y(1,2)=((x(1,2)+x(2,2))/2+(x(2,2)-(x(3,2)-x(2,2))/(x(3,1)-x(2,1))*((x(2,1)-x(1,1))/2.0)))/2
do i=1,n-2
y(i+1,1)=(x(i+1,1)+x(i+2,1))/2
y(i+1,2)=((x(i+1,2)+x(i+2,2))/2+(x(i+1,2)+(x(i+1,2)-x(i,2))/(x(i+1,1)-x(i,1))*((x(i+2,1)-x(i+1,1))/2.0)))/2
enddo
do i=1,n-1
z(2*i,1)=y(i,1)
z(2*i,2)=y(i,2)
z(2*i-1,1)=x(i,1)
z(2*i-1,2)=x(i,2)
enddo
z(2*n-1,1)=x(n,1)
z(2*n-1,2)=x(n,2)
end

言深深

ouhuining 发表于 2011-9-21 22:57

                               
登录/注册后可看大图

我想根据欧洲中心的预报数据，知道任何一个点上的值，这样的程序怎么写呢？
例如我知道A点的坐标，我想知道 ...

可以用双插值方法试试看，具体思想方法如下，但是相关的公式，程序需要自己去琢磨琢磨

其中P，A，m，n表示该点的相关参数（温度降水······）

言深深

zengqing0419 发表于 2011-9-4 18:49

                               
登录/注册后可看大图

不错啊。对于一些零散点插值，你有什么看法吗？

不知道你说的零散点是什么意思哈，可不可以举个例子呢？像随机数一类的吗，还是什么？

言深深

mofangbao 发表于 2011-9-2 16:00

                               
登录/注册后可看大图

你是正规军 绝对的

确实是杂牌的啊，内牛满面啊

mofangbao

你是正规军 绝对的

zengqing0419

不错啊。对于一些零三点插值，你有什么看法吗？

zengqing0419

不错啊。对于一些零散点插值，你有什么看法吗？

ouhuining

我想根据欧洲中心的预报数据，知道任何一个点上的值，这样的程序怎么写呢？
例如我知道A点的坐标，我想知道A点的温度是多少。

                               
登录/注册后可看大图

ouhuining

言深深 发表于 2011-9-22 08:39

                               
登录/注册后可看大图

可以用双插值方法试试看，具体思想方法如下，但是相关的公式，程序需要自己去琢磨琢磨

其中P，A，m，n ...

恩，谢了！我试试~

ouhuining

对公式有点疑问
一般来说，X11=X31，按照公式算法，不论P1，P2，P3，P4的值是什么，只要是在P1P3直线上的点，值都是零。插值就这样？