[OpenGL] 水面波动场景模拟 - 基于Gerstnder波实现

2024-03-03 22:50

本文主要是介绍[OpenGL] 水面波动场景模拟 - 基于Gerstnder波实现,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

【更新】我的新博客:www.ryuzhihao.cc,当然这个csdn博客也会更新

              本文在新博客中的链接:点击打开链接

1. 时间:2017/3/19、大三下学期

2. 参考文献:

[1] John Hany的博客:http://johnhany.net/2014/02/water-rendering-with-gerstner-wave/

[2] 论文:基于Gerstner波模型的海洋卷浪实时仿真(作者:黄威 同济大学)  


3. 首先给出我实现的效果:

             

            

            

 


4. 关于Gerstner波的基本原理:

        经典的Gerstner模型从动力学的角度描述了海浪各质点的运动.Gerstner波是最早的用于计算机图形学海浪仿真的方法,1986年被Fournier首次引入计算机图形图像领域.

4.1 单个Gerstner波的数学模型:

                     
      上式为只考虑时刻t单个波的情况。即:原质点为(x0,y0),随着时间t,该点做波幅 r波数k 角频率w 周期运动


4.2 多个Gerstner叠加波的数学模型

线性海浪理论认为,海浪是由多个不同波幅和不同角频率的波线性迭加而成.因此, 我们采用由上式得出的多个波的三维离散形式不明白原理的,可以直接拿去用:-D:
                          
       式中 ,定义xoz平面为水面静止时的平面,x轴正向朝右,z轴正向朝外,y轴为垂直方向,正向朝上,海面上的每一个点(x,y,z)都绕其静止时的位置(x0,y0,z0)圆周运动.aij是组成波的波幅;ki是波数; ω i 是角频率,并在xoz水平面上沿x轴成θj角方向传播;n为频率的划分个数;m为方向的划分个数。


4.3 最后一步:关于Gerstner参数的选择:

        令k=(ω×ω)/g (其中g=9.8为重力加速度)


5. 相关代码:

存储方式:使用顶点数组存储,并用三角网格绘制出来。

我设定的一组比较好用的参数:

    n = 3;  // 频率划分个数m = 2; // 方向的划分个数double _thetas[m] = {0.38,1.42};  // 传播方向double _amplitudes[n][m] ={  //波幅0.2,0.2,0.3,0.50,0.2,0.6,};double _omegas[n] = {3.27,3.31,3.42};     //  角频率double _ks[n] = {1.091,1.118,1.1935};  // 波数

定义一个方便计算的三维点类

class Point3
{
public:double x,y,z;Point3():x(0),y(0),z(0){}Point3(double x,double y,double z):x(x),y(y),z(z){}
};

获取(x,y,z)点在time时刻的Gerstner波对应的值,剩下的工作只需要,调用Gerstner()函数对顶点数组中每个顶点都定时计算对应的值然后绘制出来即可:

Point3 Water::Gerstner(double x, double y, double z)
{for(int i=0; i<n; i++){for(int j=0; j<m; j++){x-= cos(thetas[j])*amplitudes[i][j]*sin(waveNums[i]*(x*cos(thetas[j])+z*sin(thetas[j]))-omegas[i]*time);y+= amplitudes[i][j]*cos(waveNums[i]*(x*cos(thetas[j])+z*sin(thetas[j]))-omegas[i]*time);z-= sin(thetas[j])*amplitudes[i][j]*sin(waveNums[i]*(x*cos(thetas[j])+z*sin(thetas[j]))-omegas[i]*time);}}return Point3(x,y,z);
}


这篇关于[OpenGL] 水面波动场景模拟 - 基于Gerstnder波实现的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/771178

相关文章

Python实现终端清屏的几种方式详解

《Python实现终端清屏的几种方式详解》在使用Python进行终端交互式编程时,我们经常需要清空当前终端屏幕的内容,本文为大家整理了几种常见的实现方法,有需要的小伙伴可以参考下... 目录方法一:使用 `os` 模块调用系统命令方法二:使用 `subprocess` 模块执行命令方法三:打印多个换行符模拟

SpringBoot+EasyPOI轻松实现Excel和Word导出PDF

《SpringBoot+EasyPOI轻松实现Excel和Word导出PDF》在企业级开发中,将Excel和Word文档导出为PDF是常见需求,本文将结合​​EasyPOI和​​Aspose系列工具实... 目录一、环境准备与依赖配置1.1 方案选型1.2 依赖配置(商业库方案)二、Excel 导出 PDF

Python实现MQTT通信的示例代码

《Python实现MQTT通信的示例代码》本文主要介绍了Python实现MQTT通信的示例代码,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一... 目录1. 安装paho-mqtt库‌2. 搭建MQTT代理服务器(Broker)‌‌3. pytho

使用zip4j实现Java中的ZIP文件加密压缩的操作方法

《使用zip4j实现Java中的ZIP文件加密压缩的操作方法》本文介绍如何通过Maven集成zip4j1.3.2库创建带密码保护的ZIP文件,涵盖依赖配置、代码示例及加密原理,确保数据安全性,感兴趣的... 目录1. zip4j库介绍和版本1.1 zip4j库概述1.2 zip4j的版本演变1.3 zip4

python生成随机唯一id的几种实现方法

《python生成随机唯一id的几种实现方法》在Python中生成随机唯一ID有多种方法,根据不同的需求场景可以选择最适合的方案,文中通过示例代码介绍的非常详细,需要的朋友们下面随着小编来一起学习学习... 目录方法 1:使用 UUID 模块(推荐)方法 2:使用 Secrets 模块(安全敏感场景)方法

Spring StateMachine实现状态机使用示例详解

《SpringStateMachine实现状态机使用示例详解》本文介绍SpringStateMachine实现状态机的步骤,包括依赖导入、枚举定义、状态转移规则配置、上下文管理及服务调用示例,重点解... 目录什么是状态机使用示例什么是状态机状态机是计算机科学中的​​核心建模工具​​,用于描述对象在其生命

Spring Boot 结合 WxJava 实现文章上传微信公众号草稿箱与群发

《SpringBoot结合WxJava实现文章上传微信公众号草稿箱与群发》本文将详细介绍如何使用SpringBoot框架结合WxJava开发工具包,实现文章上传到微信公众号草稿箱以及群发功能,... 目录一、项目环境准备1.1 开发环境1.2 微信公众号准备二、Spring Boot 项目搭建2.1 创建

IntelliJ IDEA2025创建SpringBoot项目的实现步骤

《IntelliJIDEA2025创建SpringBoot项目的实现步骤》本文主要介绍了IntelliJIDEA2025创建SpringBoot项目的实现步骤,文中通过示例代码介绍的非常详细,对大家... 目录一、创建 Spring Boot 项目1. 新建项目2. 基础配置3. 选择依赖4. 生成项目5.

nginx -t、nginx -s stop 和 nginx -s reload 命令的详细解析(结合应用场景)

《nginx-t、nginx-sstop和nginx-sreload命令的详细解析(结合应用场景)》本文解析Nginx的-t、-sstop、-sreload命令,分别用于配置语法检... 以下是关于 nginx -t、nginx -s stop 和 nginx -s reload 命令的详细解析,结合实际应

Linux下删除乱码文件和目录的实现方式

《Linux下删除乱码文件和目录的实现方式》:本文主要介绍Linux下删除乱码文件和目录的实现方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录linux下删除乱码文件和目录方法1方法2总结Linux下删除乱码文件和目录方法1使用ls -i命令找到文件或目录