GAMES101-Assignment8

一、总览

1.1 连接绳子的约束

在rope.cpp 中, 实现Rope 类的构造函数。这个构造函数应该可以创建一个
新的绳子(Rope) 对象，该对象从start 开始，end 结束，包含num_nodes 个节点。也就是如下图所示：
在这里插入图片描述

每个结点都有质量，称为质点；质点之间的线段是一个弹簧。通过创建一系列的质点和弹簧，你就可以创建一个像弹簧一样运动的物体。
pinned_nodes 设置结点的索引。这些索引对应结点的固定属性(pinned attribute)应该设置为真（他们是静止的）。对于每一个结点，你应该构造一个Mass对象，并在Mass对象的构造函数里设置质量和固定属性。（请仔细阅读代码，确定传递给构造函数的参数）。你应该在连续的两个结点之间创建一个弹簧，设置弹簧两端的结点索引和弹簧系数k，请检查构造函数的签名以确定传入的参数。
运行./ropesim。你应该可以看到屏幕上画出绳子，但它不发生运动。

1.2 显式/半隐式欧拉法

胡克定律表示弹簧连接的两个质点之间的力和他们之间的距离成比例。也就是：
在这里插入图片描述

在Rope::simulateEuler 中, 首先实现胡克定律。遍历所有的弹簧，对弹簧两端的质点施加正确的弹簧力。保证力的方向是正确的！对每个质点，累加所有的弹簧力。
一旦计算出所有的弹簧力，对每个质点应用物理定律：
在这里插入图片描述

运行./ropesim。仿真应该就开始运行了，但是只有3 个结点，看起来不够多。在application.cpp 文件的最上方，你应该可以看到欧拉绳子Verlet绳子的定义。改变两个绳子结点个数（默认为3个），比如16或者更多。
运行./ropesim -s 32 来设置仿真中每帧不同的仿真步数。尝试设置较小的
值和较大的值（默认值为64）。

1.3 显式Verlet

Verlet 是另一种精确求解所有约束的方法。这种方法的优点是只处理仿真中
顶点的位置并且保证四阶精度。和欧拉法不同，Verlet 积分按如下的方式来更新下一步位置：
在这里插入图片描述

除此之外，我们可以仿真弹簧系数无限大的弹簧。不用再考虑弹簧力，而是用解约束的方法来更新质点位置：只要简单的移动每个质点的位置使得弹簧的长度保持原长。修正向量应该和两个质点之间的位移成比例，方向为一个质点指向另一质点。每个质点应该移动位移的一半。
只要对每个弹簧执行这样的操作，我们就可以得到稳定的仿真。为了使运动更加平滑，每一帧可能需要更多的仿真次数。

1.4 阻尼

向显示Verlet方法积分的胡克定律中加入阻尼。现实中的弹簧不会永远跳动，因为动能会因摩擦而减小。阻尼系数设置为0.00005, 加入阻尼之后质点位置更新如下：
在这里插入图片描述

应该修改的函数是:

  • rope.cpp 中的Rope::rope(…)
  • rope.cpp 中的void Rope::simulateEuler(…)
  • rope.cpp 中的void Rope::simulateVerlet(…)

二、安装依赖

$ sudo apt install libglu1-mesa-dev freeglut3-dev mesa-common-dev
$ sudo apt install xorg-dev

三、参考答案

3.1 rope.cpp中的Rope::rope(…)

实现1.1 连接绳子的约束

    Rope::Rope(Vector2D start, Vector2D end, int num_nodes, float node_mass, float k, vector<int> pinned_nodes){// TODO (Part 1): Create a rope starting at `start`, ending at `end`, and containing `num_nodes` nodes.for(int i = 0; i < num_nodes; i++){// 创建质点Vector2D pos = start + (end - start) * (double(i) / double(num_nodes));masses.push_back(new Mass(pos, node_mass, false));}for(int i = 0; i < num_nodes - 1; i++){// 创建弹簧springs.push_back(new Spring(masses[i], masses[i + 1], k));}
//        Comment-in this part when you implement the constructorfor (auto &i : pinned_nodes) {masses[i]->pinned = true;}}

3.2 rope.cpp 中的void Rope::simulateEuler(…)

实现1.2 显式/半隐式欧拉法

  • 胡克定律
    在这里插入图片描述
    for (auto &s : springs){// TODO (Part 2): Use Hooke's law to calculate the force on a nodeauto len = (s->m1->position - s->m2->position).norm();s->m1->forces += -s->k * (s->m1->position - s->m2->position) / len * (len - s->rest_length);s->m2->forces += -s->k * (s->m2->position - s->m1->position) / len * (len - s->rest_length);}
  • 显式/隐式欧拉法
    在这里插入图片描述
    if (!m->pinned){// TODO (Part 2): Add the force due to gravity, then compute the new velocity and position// 显式欧拉法：下一个时刻位置用当前速度计算auto a = m->forces / m->mass + gravity;m->position += m->velocity * delta_t;//根据当前速度更新位置m->velocity += a * delta_t;//更新下一时刻速度// 隐式欧拉法：下一个时刻位置用下一时刻速度计算// auto a = m->forces / m->mass + gravity;// m->velocity += a * delta_t;//更新速度// m->position += m->velocity * delta_t;//根据下一时刻速度更新位置 // TODO (Part 2): Add global damping}
  • 完整代码
    void Rope::simulateEuler(float delta_t, Vector2D gravity){for (auto &s : springs){// TODO (Part 2): Use Hooke's law to calculate the force on a nodeauto len = (s->m1->position - s->m2->position).norm();s->m1->forces += -s->k * (s->m1->position - s->m2->position) / len * (len - s->rest_length);s->m2->forces += -s->k * (s->m2->position - s->m1->position) / len * (len - s->rest_length);}for (auto &m : masses){if (!m->pinned){// TODO (Part 2): Add the force due to gravity, then compute the new velocity and position// 显式欧拉法：下一个时刻位置用当前速度计算auto a = m->forces / m->mass + gravity;m->position += m->velocity * delta_t;//根据当前速度更新位置m->velocity += a * delta_t;//更新下一时刻速度// 隐式欧拉法：下一个时刻位置用下一时刻速度计算// auto a = m->forces / m->mass + gravity;// m->velocity += a * delta_t;//更新速度// m->position += m->velocity * delta_t;//根据下一时刻速度更新位置 // TODO (Part 2): Add global damping}// Reset all forces on each massm->forces = Vector2D(0, 0);}}

3.3 rope.cpp 中的void Rope::simulateVerlet(…)

实现1.3 显式Verlet

  • 显式Verlet
    在这里插入图片描述
    void Rope::simulateVerlet(float delta_t, Vector2D gravity){for (auto &s : springs){// TODO (Part 3): Simulate one timestep of the rope using explicit Verlet （solving constraints)auto len = (s->m1->position - s->m2->position).norm();s->m1->forces += -s->k * (s->m1->position - s->m2->position) / len * (len - s->rest_length);s->m2->forces += -s->k * (s->m2->position - s->m1->position) / len * (len - s->rest_length);}for (auto &m : masses){if (!m->pinned){Vector2D temp_position = m->position;// TODO (Part 3.1): Set the new position of the rope massauto a = m->forces / m->mass + gravity;m->position = temp_position + (temp_position - m->last_position) + a * delta_t * delta_t;m->last_position = temp_position;// TODO (Part 4): Add global Verlet damping}m->forces = Vector2D(0, 0);}}

3.4 增加阻尼

  • rope.cpp 中的void Rope::simulateVerlet(…)
    向显示Verlet方法积分的胡克定律中加入阻尼。现实中的弹簧不会永远跳动，因为动能会因摩擦而减小。阻尼系数设置为0.00005, 加入阻尼之后质点位置更新如下：
    在这里插入图片描述
    void Rope::simulateVerlet(float delta_t, Vector2D gravity){for (auto &s : springs){// TODO (Part 3): Simulate one timestep of the rope using explicit Verlet （solving constraints)auto len = (s->m1->position - s->m2->position).norm();s->m1->forces += -s->k * (s->m1->position - s->m2->position) / len * (len - s->rest_length);s->m2->forces += -s->k * (s->m2->position - s->m1->position) / len * (len - s->rest_length);}for (auto &m : masses){if (!m->pinned){Vector2D temp_position = m->position;// TODO (Part 3.1): Set the new position of the rope massauto a = m->forces / m->mass + gravity;//m->position = temp_position + (temp_position - m->last_position) + a * delta_t * delta_t;//m->last_position = temp_position;// TODO (Part 4): Add global Verlet dampingdouble damping_factor = 0.00005;m->position = temp_position + (1 - damping_factor) * (temp_position - m->last_position) + a * delta_t * delta_t;m->last_position = temp_position;}m->forces = Vector2D(0, 0);}}
  • rope.cpp 中的void Rope::simulateEuler(…)
    向欧拉方法中加入阻尼，直接使用−kdv作为阻尼，而不是相对速度
    void Rope::simulateEuler(float delta_t, Vector2D gravity){for (auto &s : springs){// TODO (Part 2): Use Hooke's law to calculate the force on a nodeauto len = (s->m1->position - s->m2->position).norm();s->m1->forces += -s->k * (s->m1->position - s->m2->position) / len * (len - s->rest_length);s->m2->forces += -s->k * (s->m2->position - s->m1->position) / len * (len - s->rest_length);}for (auto &m : masses){if (!m->pinned){// TODO (Part 2): Add the force due to gravity, then compute the new velocity and position// 显式欧拉法：下一个时刻位置用当前速度计算// auto a = m->forces / m->mass + gravity;// m->position += m->velocity * delta_t;//根据当前速度更新位置// m->velocity += a * delta_t;//更新下一时刻速度// 隐式欧拉法：下一个时刻位置用下一时刻速度计算// auto a = m->forces / m->mass + gravity;// m->velocity += a * delta_t;//更新速度// m->position += m->velocity * delta_t;//根据下一时刻速度更新位置 // TODO (Part 2): Add global dampingfloat kd = 0.005;auto a = m->forces / m->mass + gravity - kd * m->velocity / m->mass;m->velocity += a * delta_t; m->position += m->velocity * delta_t; //隐式欧拉法}// Reset all forces on each massm->forces = Vector2D(0, 0);}}

四、编译

$ mkdir build
$ cd build
$ cmake ..
$ make

通过查看rope.cpp中的Rope::rope(…)的引用，可以定位到application.cpp文件，

// Create two ropes 
// 通过修改构造函数的第三个参数指定质点个数(默认三个质点，这里改成了10个)
ropeEuler = new Rope(Vector2D(0, 200), Vector2D(-400, 200), 10, config.mass,config.ks, {0});
ropeVerlet = new Rope(Vector2D(0, 200), Vector2D(-400, 200), 10, config.mass,config.ks, {0});

附件

作业7压缩包

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.mzph.cn/news/642529.shtml

如若内容造成侵权/违法违规/事实不符，请联系多彩编程网进行投诉反馈email:809451989@qq.com，一经查实，立即删除！

相关文章

洛谷P3392 涂国旗-----暴力美学的演绎

洛谷P3392 涂国旗-----暴力美学的演绎

# 涂国旗 ## 题目描述 某国法律规定&#xff0c;只要一个由 $N \times M$ 个小方块组成的旗帜符合如下规则&#xff0c;就是合法的国旗。&#xff08;毛熊&#xff1a;阿嚏——&#xff09; - 从最上方若干行&#xff08;至少一行&#xff09;的格子全部是白色的&#xff1b…
阅读更多...
【论文阅读】Automated Runtime-Aware Scheduling for Multi-Tenant DNN Inference on GPU

【论文阅读】Automated Runtime-Aware Scheduling for Multi-Tenant DNN Inference on GPU

该论文发布在 ICCAD’21 会议。该会议是EDA领域的顶级会议。 基本信息 AuthorHardwareProblemPerspectiveAlgorithm/StrategyImprovment/AchievementFuxun YuGPUResource under-utilization ContentionSW SchedulingOperator-level schedulingML-based scheduling auto-searc…
阅读更多...
CGAL 网格连通聚类

CGAL 网格连通聚类

文章目录 一、简介二、实现代码三、实现效果参考资料一、简介 这里的思路其实与点云的欧式聚类非常类似，区别在于点云的欧式聚类是通过搜索半径对点云进行聚类，至于基于连通性网格面片聚类则是通过面片的邻近关系对面片进行聚类，大致的过程与欧式聚类是相同的： 首先，需要指…
阅读更多...
Rustdesk 中VP8 / VP9 / AV1 是什么？

Rustdesk 中VP8 / VP9 / AV1 是什么？

环境&#xff1a; Rustdesk1.1.9 VP8 / VP9 / AV1 问题描述&#xff1a; VP8 / VP9 / AV1 是什么&#xff1f; 解决方案&#xff1a; 1.VP8、VP9和AV1是视频编解码器&#xff0c;用于压缩和解压缩视频数据。它们是由Google和Alliance for Open Media&#xff08;AOM&#…
阅读更多...
011 变量

011 变量

变量的定义 注意事项 作用域 类变量&#xff1a;定义在类中的变量&#xff0c;类的整个生命周期内可用实例变量&#xff1a;定义在类中的变量&#xff0c;类的实例可以使用局部变量&#xff1a;定义在方法中的变量&#xff0c;仅在方法内可以使用 class Dog{// 类变量static …
阅读更多...
Oracle 19c RAC集群管理 ---------关键参数以及常用命令

Oracle 19c RAC集群管理 ---------关键参数以及常用命令

Oracle 19c RAC集群管理 ---------关键参数 Oracle 19C RAC 参数最佳实践 --开启强制归档 ALTER DATABASE FORCE LOGGING; --设置 30分钟 强制归档 ALTER SYSTEM SET ARCHIVE_LAG_TARGET1800 SCOPEBOTH SID*; --设置期望undo保持时间3h ALTER SYSTEM SET UNDO_RETENTION21600…
阅读更多...
系统学习Python——警告信息的控制模块warnings：常用函数-[warnings.simplefilter]

系统学习Python——警告信息的控制模块warnings：常用函数-[warnings.simplefilter]

分类目录&#xff1a;《系统学习Python》总目录 函数 warnings.simplefilter(action, categoryWarning, lineno0, appendFalse)在警告过滤器种类列表中插入一条简单数据项。函数参数的含义与filterwarnings()相同&#xff0c;但不需要正则表达式&#xff0c;因为插入的过滤器总…
阅读更多...
机器学习 | 深入理解并掌握核心概念

机器学习 | 深入理解并掌握核心概念

在如今数字化时代的浪潮下&#xff0c;机器学习已经成为人工智能领域的璀璨明星。它像一面魔镜&#xff0c;赋予计算机系统学习和改进的能力&#xff0c;让机器能够从海量数据中提取规律、预测未来&#xff0c;甚至做出智能决策。本 专栏 将带您踏上机器学习的奇妙之旅&#xf…
阅读更多...
题解：CF1920E. Counting Binary Strings

题解：CF1920E. Counting Binary Strings

题解&#xff1a;CF1920E. Counting Binary Strings 题意简述 题目链接&#xff1a;Problem - E - Codeforces。 洛谷翻译&#xff1a;Counting Binary Strings - 洛谷。 思路解析 假设我们有一个01串str&#xff08;设里面有z个“1”&#xff09;&#xff0c;我们要求它里…
阅读更多...
Ubuntu (Linux) 下创建软链接（即符号链接，相当于windows下的快捷方式）方法

Ubuntu (Linux) 下创建软链接（即符号链接，相当于windows下的快捷方式）方法

Ubuntu (Linux) 下创建软链接&#xff08;即符号链接&#xff0c;相当于windows下的快捷方式&#xff09;方法 使用创建软链接的命令 #命令格式如下。注意&#xff1a;请使用绝对路径&#xff0c;否则链接可能失效 ln -s <源文件或目录的绝对路径> <符号链接文件&am…
阅读更多...
ImportError: cannot import name ‘complete_to_chordal_graph‘解决

ImportError: cannot import name ‘complete_to_chordal_graph‘解决

文章目录 ImportError: cannot import name complete_to_chordal_graph解决杂话问题解决注 ImportError: cannot import name complete_to_chordal_graph’解决 杂话 大伙儿应该都用过NetworkX这个库函数吧&#xff0c;是一个很方便的图相关的库&#xff0c;但是如果你的版本…
阅读更多...
IDEA更改页面不重启

IDEA更改页面不重启

IDEA更改页面不重启 idea若依 修改包名 idea 1、修改IDEA设置 File -> Settings -> Build Execution Deployment -> Build Project automatically 勾选2、勾选Running Ctrl Shift Alt / 然后选择 Registry&#xff0c;勾上 Compiler.autoMake.allow.when.app.runn…
阅读更多...
算法训练营第五十五天|392.判断子序列 115.不同的子序列

算法训练营第五十五天|392.判断子序列 115.不同的子序列

目录 Leetcode392.判断子序列Leetcode115.不同的子序列 Leetcode392.判断子序列 文章链接&#xff1a;代码随想录 题目链接&#xff1a;392.判断子序列 思路&#xff1a;正常动规记录&#xff0c;重点在于后面采用判断 class Solution { public:bool isSubsequence(string s, …
阅读更多...
【MySQL进阶】视图_存储过程_存储函数_触发器

【MySQL进阶】视图_存储过程_存储函数_触发器

文章目录 视图基本介绍视图操作视图创建视图查询视图修改视图删除 存储过程基本介绍基本操作存储语法变量IF语句参数传递CASEWHILEREPEATLOOP游标 存储函数触发器基本介绍基本操作 总结 视图 基本介绍 视图概念&#xff1a;视图是一种虚拟存在的数据表&#xff0c;这个虚拟的表…
阅读更多...
算法训练营第五十六天|583. 两个字符串的删除操作 72. 编辑距离

算法训练营第五十六天|583. 两个字符串的删除操作 72. 编辑距离

目录 Leetcode583. 两个字符串的删除操作Leetcode72. 编辑距离 Leetcode583. 两个字符串的删除操作 文章链接&#xff1a;代码随想录 题目链接&#xff1a;583. 两个字符串的删除操作 思路&#xff1a;直接记录需要改&#xff08;增或删&#xff09;几个&#xff0c;也就是求不…
阅读更多...
seafile+onlyoffice集成部署

seafile+onlyoffice集成部署

准备docker-compose version: "3" services:db:image: mariadb:10.11container_name: seafile-mysqlrestart: alwaysenvironment:- MYSQL_ROOT_PASSWORDsea123456 # Requested, set the roots password of MySQL service.- MYSQL_LOG_CONSOLEtruevolumes:- ./data/…
阅读更多...
透明拼接屏在汽车领域的应用

透明拼接屏在汽车领域的应用

随着科技的进步&#xff0c;透明拼接屏作为一种新型的显示技术&#xff0c;在汽车领域的应用越来越广泛。尼伽小编将围绕透明拼接屏在汽车本身、4S店、展会、工厂等方面的应用进行深入探讨&#xff0c;并展望未来的设计方向。 一、透明拼接屏在汽车本身的应用 车窗显示&#x…
阅读更多...
openssl3.2/test/certs - 011 - primary client-EKU root: croot-cert

openssl3.2/test/certs - 011 - primary client-EKU root: croot-cert

文章目录 openssl3.2/test/certs - 011 - primary client-EKU root: croot-cert概述笔记END openssl3.2/test/certs - 011 - primary client-EKU root: croot-cert 概述 openssl3.2 - 官方demo学习 - test - certs 笔记 \file my_openssl_linux_log_doc_011.txt \note open…
阅读更多...
Python入门（一）

Python入门（一）

anaconda安装 官网&#xff1a;https://www.anaconda.com下载 jupyter lab 简介&#xff1a; 包含了Jupyter Notebook所有功能。 JupyterLab作为一种基于web的集成开发环境&#xff0c;你可以使用它编写notebook&#xff0c;操作终端&#xff0c;编辑markdown文本&#xf…
阅读更多...
Linux中文件属性的获取(stat、chmod、Istat、fstat函数的使用)

Linux中文件属性的获取(stat、chmod、Istat、fstat函数的使用)

修改文件权限 函数如下&#xff1a; chmod/fchmod函数用来修改文件的访问权限: #include <sys/stat.h> int chmod(const char *path, mode_t mode); int fchmod(int fd, mode_t mode); 成功时返回0&#xff1b;出错时返回EOF 注意&#xff1a;在vmware和windows共享的文…
阅读更多...
最新文章

Copyright @ 2022~2023