用查表法实现对水波的模拟

转载 2007年10月06日 22:51:00
最近在网上看到一个水波摸拟的贴子。觉得算法可以更快一点,因此自己想出这个查表法模拟的程序,大家看看有没有价值。

原理
在一个波长范围内,任意点的位移disp = sin(dist)*force +disp_center;
这里设水波为正弦波,振幅为force,disp_center为振中的位移。
推广到任何点,disp = sin(dist % wave_length)*force + disp_center;
这里wave_length 为波长。可以看出,在离散的数据中,任何点的位移是有限的集合。我们先构造出这个集合,在运算过程中,从集合中选取一个位移值就行了,而不用实时计算每个点的位移。

程序设计如下:取波长为12,任何点的位移由一个12个元素的数组。
int disp[12];


建立一个数组成存贮各点来振中的距离
int dist[600][800];

下面生成一个水波
inline void Build_Disp_Matrix()
{
    for(int i = 0; i < 12; ++i)
    disp[i] = static_cast<int>(sin((double)i*30*DXPI/360)*PRECISION);
}

void Create_Ripple(int center_x, int center_y, int force, int damp, int wave_length)
{

    for(int i = 0 ; i < wave_length; ++i)
    {
        disp[i] = force * disp[i] / PRECISION;
    }

    int dist_x, dist_y;
    for(i = 0; i < SCREEN_HEIGHT; ++i)
    {
        for(int j = 0; j< SCREEN_WIDTH; ++j)
        {
            dist_y = i - center_y;
            dist_x = j - center_x;
            dist[i][j] = static_cast<int>(sqrt(dist_x*dist_x + dist_y*dist_y) + 0.5) % wave_length;
        }
    }

}


渲染水波的代码与原来的代码一样。
void Render_Ripple(LPDIRECTDRAWSURFACE7 lpdds)
{
    DDSURFACEDESC2 ddsd = { sizeof(ddsd) };

    if(disp_index < 0) disp_index = 12;

    --disp_index;

    lpdds->Lock(NULL,&ddsd, DDLOCK_WAIT | DDLOCK_SURFACEMEMORYPTR, NULL);
    short* lpSrcBuffer = static_cast<short*>(pBmpBack->GetBitmapBuffer());

    int pitch = (ddsd.lPitch) >> 1;

    short* lpBuffer = static_cast<short*>(ddsd.lpSurface);

    int left,right,bottom,top;
    int k = 0,l = 0;
    int index_x,index_y;
    int x_off,y_off;
    int pos1,pos2;
    for(int i = 1; i < SCREEN_HEIGHT - 2; ++i)
        for( int j = 1; j < SCREEN_WIDTH - 2; ++j)
        {
            left = right = bottom = top = disp_index;

            left += dist[i][j -1];
            if(left > 12) left -= 12;
            right += dist[i][j+1];
            if(right >12) right -= 12;
            bottom += dist[i+1][j];
            if(bottom > 12) bottom -= 12;
            top += dist[i-1][j];
            if(top >12 ) top -= 12;

            x_off = disp[left] - disp[right];
            y_off = disp[top] - disp[bottom];

            index_x = j + x_off;
            index_y = i + y_off;

            pos1 = i*pitch +j;
            pos2 = index_y * 800 + index_x;
            if(pos2 < 0 || pos2 > 479999) break;

            lpBuffer[pos1] = lpSrcBuffer[pos2];
        }

    pBmpBack->ReleaseBuffer();
    lpdds->Unlock(NULL);
}


不知道这样大家都否看懂。完整的演示程序
 

JavaScript实现的CRC16源码非查表法

  • 2016年01月14日 09:53
  • 646B
  • 下载

CRC32查表法算法实现

CRC参数模型:   Name : "CRC-32"   Width : 32   Poly : 04C11DB7   Init : FFFFFFFF   RefIn : True   R...

QT实现CRC16校验(查表法)

近期产品要做认证,那边需要用CRC校验做串口通信的校验方式,就通过两天的时间改了下做模拟。 比较了CRC校验的几种方法后发现查表法是最快的。 因为字符一共就256个,所以可以将256个字符的CRC...

查表法实现PCM与Alaw、μlaw之间的格式转换 2011

 转载于:http://blog.chinaunix.net/uid-10995602-id-2918727.html 参考地址一:http://yushui213.blogcn.co...

Interger类中进制转换方法的实现原理-通过查表法完成整数的进制转换

前两天没事看Interger类的源码,看到里面封装的整数进制转换的方法的实现挺有意思的。 于是研究了一下,也动手模拟写了一遍还是挺有收获的。 简单的总结就是结合了Java中的位运算以及算法“查表法”...

二次查表法实现真彩色图像YUV/RGB格式转换

作者: 田扬林 许志祥 摘要】 提出了一种基于二次查表实现真彩色图像YUV/RGB格式转换的快速有效算法,并将其应用于H.263编解码中。采用这种二次查表相加法能提高运算速度近4倍而没有任何图像质量...
  • asb2010
  • asb2010
  • 2015年01月23日 14:17
  • 401

查表法8b10b编解码

  • 2014年11月01日 21:37
  • 1.46MB
  • 下载

Ripple游戏水波模拟算法实现

在介绍编程之前,先让我们来回顾一下在高中的物理课上我们所学的关于水波的知识。 水波有如下几个特性:  扩散:当你投一块石头到水中,你会看到一个以石头入水点为圆心所形成的一圈圈的水波,这里,...

vdsp查表法

  • 2012年05月01日 15:50
  • 718B
  • 下载
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:用查表法实现对水波的模拟
举报原因:
原因补充:

(最多只允许输入30个字)