关于MATLAB计算3维图的向量夹角总是不正确的问题记录

文章目录


问题描述

因为最近在做无人机的一个项目,所以需要画出无人机的轨迹,然后再提取特征值,我这里在计算夹角的时候发现为什么在视觉上明明看的是钝角但是实际计算出来却是锐角的角度。

如下图所示,看起来就是一个钝角,但是计算的结果始终是锐角。

解决方法

后来去请教了一下师姐,才发现是因为向量的方向问题,我想要计算的是上面这一种形式,但是我按照我的写法实际上是下面那一种形式。

改正之前的写法:

matlab 复制代码
vec1 = [x(2)-x(1),y(2)-y(1),z(2)-z(1)];
vec2 = [x(3)-x(2),y(3)-y(2),z(3)-z(2)];

cos_angle =  dot(vec1 ,vec2)/(norm(vec1)*norm(vec2));
angle = acosd(cos_angle);

改正之后的写法:

matlab 复制代码
vec1 = [x(1)-x(2),y(1)-y(2),z(1)-z(2)]; %更改了这里
vec2 = [x(3)-x(2),y(3)-y(2),z(3)-z(2)];

cos_angle =  dot(vec1 ,vec2)/(norm(vec1)*norm(vec2));
angle = acosd(cos_angle);

结果就是钝角了,不得不说,还是师姐更厉害啊。

完整代码

matlab 复制代码
clear;
clc;

x %yourData
y %yourData
z %yourData

vec1 = [x(1)-x(2),y(1)-y(2),z(1)-z(2)];
vec2 = [x(3)-x(2),y(3)-y(2),z(3)-z(2)];

cos_angle =  dot(vec1 ,vec2)/(norm(vec1)*norm(vec2));
angle = acosd(cos_angle);


plot3(x, y, z, '-o', 'LineWidth', 1.5);  % 连续点迹绘制轨迹
hold on;
相关推荐
李长渊哦29 分钟前
深入理解 JavaScript 中的全局对象与 JSON 序列化
开发语言·javascript·json
若水晴空初如梦1 小时前
QT聊天项目DAY06
开发语言·qt
h汉堡2 小时前
C++入门基础
开发语言·c++·学习
HtwHUAT3 小时前
实验四 Java图形界面与事件处理
开发语言·前端·python
鄃鳕3 小时前
QSS【QT】
开发语言·qt
汤姆_5113 小时前
【c语言】深度理解指针4——sizeof和strlen
c语言·开发语言
碎梦归途3 小时前
23种设计模式-结构型模式之外观模式(Java版本)
java·开发语言·jvm·设计模式·intellij-idea·外观模式
muyouking113 小时前
4.Rust+Axum Tower 中间件实战:从集成到自定义
开发语言·中间件·rust
FAREWELL000754 小时前
C#进阶学习(九)委托的介绍
开发语言·学习·c#·委托
我该如何取个名字4 小时前
Mac配置Java的环境变量
java·开发语言·macos