顺心随笔(二)

好久没有写博客了。

现在回头看看,真不知道当时是不是只是有一腔热血,去跟着一个离自己比较远的专业去刷题,打比赛。毕竟是自己的选择,虽然我不想以它作为一个职业,但是我仍然很喜欢那种全神贯注的感觉。

路都是自己选的,现在已经脱坑,我希望在自己的路上好好走下去。做每一个决定我都不后悔,毕竟,ACM给了我一段不同的经历,每个人的世界都可以是很丰富的!

现在,大三的自己,被考研和保研包围,想努力出一些比较硬的成果,给自己一些信心,一些去和他们竞争的勇气!

虽然有压力,但内心还是很快乐的,我享受青春的每一秒钟,我享受精力充沛的自己,我享受早起的清晨,我享受每天醒来知道身边有一个花儿一般的她,我享受自己付出可以收获的感觉。

加油,为这一段只有音乐和白开水的时光!为自己的青春年华!毕竟人只能活一次,我的大学时光只有一次,苦也好,乐也好,我愿意全盘接受,来吧!来得更猛烈些吧!!!

根据提供的引用[1],刚体绕顺心转动是指刚体绕过其质心的一条固定轴旋转,且旋转角速度大小保持不变。在Matlab中,可以使用`quiver3`函数来可视化刚体绕顺心转动的现象。具体步骤如下: 1.定义刚体的质心坐标和旋转轴的方向向量,以及旋转角速度大小。 ```matlab % 刚体质心坐标 x0 = 0; y0 = 0; z0 = 0; % 旋转轴方向向量 u = 1; v = 1; w = 1; % 旋转角速度大小 omega = 1; ``` 2.定义旋转轴的起点和终点坐标,并使用`quiver3`函数绘制旋转轴。 ```matlab % 旋转轴起点坐标 x_start = x0 - u/2; y_start = y0 - v/2; z_start = z0 - w/2; % 旋转轴终点坐标 x_end = x0 + u/2; y_end = y0 + v/2; z_end = z0 + w/2; % 绘制旋转轴 quiver3(x_start, y_start, z_start, u, v, w, 'LineWidth', 2, 'MaxHeadSize', 0.5); hold on; ``` 3.定义刚体的顶点坐标,并将其绕顺心旋转一定角度后的坐标计算出来。 ```matlab % 刚体顶点坐标 vertices = [1 1 1; 1 1 -1; 1 -1 1; 1 -1 -1; -1 1 1; -1 1 -1; -1 -1 1; -1 -1 -1]; % 将顶点坐标平移到刚体质心处 vertices(:,1) = vertices(:,1) + x0; vertices(:,2) = vertices(:,2) + y0; vertices(:,3) = vertices(:,3) + z0; % 计算旋转后的顶点坐标 R = vrrotvec2mat([u v w omega]); % 旋转矩阵 vertices_rotated = (R * vertices')'; % 顶点坐标矩阵乘以旋转矩阵 ``` 4.将旋转后的刚体顶点坐标绘制出来。 ```matlab % 绘制旋转后的刚体 patch('Vertices', vertices_rotated, 'Faces', [1 2 4 3; 5 6 8 7; 1 2 6 5; 3 4 8 7; 1 3 7 5; 2 4 8 6], 'FaceColor', 'r', 'FaceAlpha', 0.5); axis equal; ``` 完整的Matlab代码如下: ```matlab % 刚体质心坐标 x0 = 0; y0 = 0; z0 = 0; % 旋转轴方向向量 u = 1; v = 1; w = 1; % 旋转角速度大小 omega = 1; % 旋转轴起点坐标 x_start = x0 - u/2; y_start = y0 - v/2; z_start = z0 - w/2; % 旋转轴终点坐标 x_end = x0 + u/2; y_end = y0 + v/2; z_end = z0 + w/2; % 绘制旋转轴 quiver3(x_start, y_start, z_start, u, v, w, 'LineWidth', 2, 'MaxHeadSize', 0.5); hold on; % 刚体顶点坐标 vertices = [1 1 1; 1 1 -1; 1 -1 1; 1 -1 -1; -1 1 1; -1 1 -1; -1 -1 1; -1 -1 -1]; % 将顶点坐标平移到刚体质心处 vertices(:,1) = vertices(:,1) + x0; vertices(:,2) = vertices(:,2) + y0; vertices(:,3) = vertices(:,3) + z0; % 计算旋转后的顶点坐标 R = vrrotvec2mat([u v w omega]); % 旋转矩阵 vertices_rotated = (R * vertices')'; % 顶点坐标矩阵乘以旋转矩阵 % 绘制旋转后的刚体 patch('Vertices', vertices_rotated, 'Faces', [1 2 4 3; 5 6 8 7; 1 2 6 5; 3 4 8 7; 1 3 7 5; 2 4 8 6], 'FaceColor', 'r', 'FaceAlpha', 0.5); axis equal; ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值