无聊中，发一个随机分形插值的代码，需要的同志们拿去。原创

[BORLANDSUN]

搞随机分形插值的同志们会用到，需要请收藏，谢谢。
说明：函数用于给一维数据插值，DataIn为待插数据，Dim为待插数据长度，H为Hurst指数，Var为插值方差，Level为插值级数（Level级插值后总的点数为(2^(Level-1))*(Dim-1)+1），DataOut为插值后的数据。

#include "math.h"
#include "Math.hpp"
bool MFIInterpolation(double* DataIn, int Dim, double H, double Var, int Level, double* DataOut)
{
	if(Dim > 1)
	{
// 1. 初始化数据
		int Max = pow(2, Level)*(Dim-1);
		int L = pow(2, Level);								// 初始插值步间隔
		for(int i=0; i<Dim; i++)				// 埴充插值点
		{
			DataOut[i*L] = DataIn[i];
		}

// 2. 插值
		for(int p=0; p<Level; p++)
		{
			int count = pow(2, p)*(Dim-1);              	// 每行插值点数量
			double Scale = pow(1-pow(2, 2*H-2), 0.5);
			for(int i=0; i<=p; i++)
			{
				Scale *= 1/pow(2, H);
			}
			for(int i=0; i<count; i++)
			{
				double D = Scale * RandG(0, Var);
				int temp1 = (i+0.5)*L;
				DataOut[temp1] = (DataOut[i*L] + DataOut[(i+1)*L])/2 + D;
			}
			L /= 2;
		}
		return true;
	}
	else
	{
		return false;
    }
}

使用示例：

double In[3];
In[0] = 0;
In[1] = 80;
In[2] = 150;
Level = 8;
double Out[513];      // 8级插值后总的点数为(2^(Level-1))*(3-1)+1=513
MFIInterpolation(In, 3, 0.6, 100, Level, Out);

[zzbinfo]

热情的顶贴

[paste]

学习，这个是高科技了，用不来的干活。

[samchoy]

俺也来灰常热情的顶一下

[lee0507]

这是干嘛的啊？ 都不懂····呢
解释下 用途 吧 LZ

[xcluii]

数值分析~~~ 谁用bcb啊

[BORLANDSUN]

我用。bcb做数值分析怎么了？