看不到颤抖的情节

24日视图(30天)
기홍박
기홍박 2023年5月3日在8:06
编辑: 鸡向量 2023年5月3日在28
我正在试图让一个代码路径规划的势场。
我定义势能(u_att u_rep)然后力(fxa fya, fxr,马其顿)通过使用“梯度”。
尽我所知,力矢量(fxa fya, fxr,马其顿)正确地计算,
但是当我阴谋力矢量与“颤抖”,我什么也看不见。
我如何解决这个问题?
代码下面附势场
谢谢你!
关闭所有;清晰;clc;
图,;
轴(100 0 100[0]),轴广场,盒子;
包含(“x”),ylabel (“y”);
obslist = [60 60 15;
20 40 7;
70年30 5];
x_start = 10 [10];
x_goal = (90 - 90);
情节(x_start (1) x_start (2),“处方”,“线宽”2);
情节(x_goal (1) x_goal (2),“波”,“线宽”2);
ang =(0:π/ 50:2 *π)';
xpu = cos (ang);您什么=罪(ang);
i = 1:尺寸(obslist, 1)
xp = obslist(我,1)+ obslist(我,3)* xpu;
yp = obslist(我,2)+ obslist(我,3)*您什么;
填充(xp, yp,“y”);
情节(xp, yp,“k”);
结束
%以上代码从韩科院2023年春季ME652讲座
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
片= 1;
(X, Y) = meshgrid(0:片:100,0:片:100);
Ka = 5;
u_att = 0.5 * Ka *√(X-x_goal (1))。^ 2 + (Y-x_goal (2)) ^ 2);
[fxa, fya] =梯度(u_att,切片,切片);
基米-雷克南= 500;
u_rep = 0(大小(X));
i = 1:尺寸(obslist, 1)
u_rep = u_rep + 0.5。* Kr * 1. /√(X-obslist(我,1))。^ 2 + (Y-obslist(我,2))^ 2)- 1. /√((obslist(我,3)+ 2)))^ 2;
结束
(fxr,马其顿)=梯度(u_rep,切片,切片);
fX = - fxr fxa;
财政年度= -马其顿fya;
箭袋(X, Y,外汇、财政年度“k”)
%轮廓(X, Y, u_att + u_rep, 25)

接受的答案

鸡向量
鸡向量 2023年5月3日在剩
编辑:鸡向量 2023年5月3日在28
因为他们有太多的箭头和重叠,导致自动缩放不会显示任何东西。
设置规模 “关闭” :
箭袋(X, Y,外汇、财政年度0)
和调整的值 外汇 财政年度 相应的行动。
在评论中你可以看到一个例子。
2的评论
기홍박
기홍박 2023年5月3日在九11
谢谢你的帮助!

登录置评。

答案(1)

Chunru
Chunru 2023年5月3日8:23
图,;
轴(100 0 100[0]),轴广场,盒子;
包含(“x”),ylabel (“y”);
obslist = [60 60 15;
20 40 7;
70年30 5];
x_start = 10 [10];
x_goal = (90 - 90);
情节(x_start (1) x_start (2),“处方”,“线宽”2);
情节(x_goal (1) x_goal (2),“波”,“线宽”2);
ang =(0:π/ 50:2 *π)';
xpu = cos (ang);您什么=罪(ang);
i = 1:尺寸(obslist, 1)
xp = obslist(我,1)+ obslist(我,3)* xpu;
yp = obslist(我,2)+ obslist(我,3)*您什么;
填充(xp, yp,“y”);
情节(xp, yp,“k”);
结束
%以上代码从韩科院2023年春季ME652讲座
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
片= 1;
(X, Y) = meshgrid(0:片:100,0:片:100);
Ka = 5;
u_att = 0.5 * Ka *√(X-x_goal (1))。^ 2 + (Y-x_goal (2)) ^ 2);
[fxa, fya] =梯度(u_att,切片,切片);
基米-雷克南= 500;
u_rep = 0(大小(X));
i = 1:尺寸(obslist, 1)
u_rep = u_rep + 0.5。* Kr * 1. /√(X-obslist(我,1))。^ 2 + (Y-obslist(我,2))^ 2)- 1. /√((obslist(我,3)+ 2)))^ 2;
结束
(fxr,马其顿)=梯度(u_rep,切片,切片);
fX = - fxr fxa;
财政年度= -马其顿fya;
%颤抖(X, Y,外汇、财政年度“k”)
箭袋(X, Y,外汇、财政年度0)
%检查你的代码,外汇和财政年度有正值,自动定量
%的工作
max (fX(:))、马克斯(财政年度(:))
ans =正
ans =正
1评论
기홍박
기홍박 2023年5月3日在9:12
谢谢你的帮助!

登录置评。

类别

找到更多的在编程帮助中心文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!

翻译的