在matlab 画箭头

最新推荐文章于 2023-07-24 13:12:30 发布
最新推荐文章于 2023-07-24 13:12:30 发布
阅读量960 收藏
点赞数
分类专栏： 图像处理
非常实用的文章
原文地址：在matlab 画箭头 作者：纯情小郎君
完整见链接http://www.mathworks.com/matlabcentral/fx_files/14056/1/content/arrow3_examples.html

ARROW3   EXAMPLES   (R13)
Two-Dimensional Quiver Plots
Example 1. (cf. Joukowski Airfoil Transformation)

t=10*pi/180; u=0.1; v=0.1; r=1.1; tol=8e-2;
s=u+i*v; k=2*r*sin(t); w=exp(i*t);
[x,y]=meshgrid(-2.5:0.1:3.5,-3:0.1:3);
z=x+i*y; z(abs(z-s)<r-tol)=NaN;
f=w*z+exp(-i*t)*r^2./(z-s)+i*k*log(z);
a=0:0.1:2*pi; zc=r*(cos(a)+i*sin(a))+s;
c1=-1.5; c2=2.5; c3=c2-c1;
c=contour(x,y,imag(f),c1:c3/16:c2);
set(gca,'color',0.5*[1 1 1])
hold on, daspect([1 1 1])
fill(real(zc),imag(zc),'y'), colorbar
map=get(gcf,'colormap');
domain=0:1/(size(map,1)-1):1; m=1;
while m<length(c), n=m+c(2,m);
if c(2,m)>9, p=c(:,m+1:n)';
    if p(end,1)>p(1,1), p=flipud(p); end
    ndx=10:10:length(p);
    p2=p(ndx,:); p1=p(ndx-3,:);
    cc=interp1(domain,map,(c(1,m)-c1)/c3);
    set(gca,'ColorOrder',cc)
    arrow3(p1,p2,'0o',0.8)
end, m=n+1;
end, hold off

Example 1. Joukowski Airfoil Transformation
Example 1.

Example 2. (cf. Quiver)

[x,y]=meshgrid(-1:1/21:1);
z=x.*exp(-x.^2-y.^2);
[c,h]=contour(x,y,z);
set(h,'EdgeColor',0.45*[1 1 1])
ndx=1:3:length(x);
x=x(ndx,ndx); y=y(ndx,ndx); z=z(ndx,ndx);
[u,v]=gradient(z,1/7);
p1=[x(:),y(:)]; u=u(:); v=v(:);
m=abs(u+i*v); % gradient magnitude
daspect([1 1 1]), set(gca,'color',0.3*[1 1 1])
hold on, colormap hot, scale=0.4;
arrow3(p1,p1+scale*[u,v],'|',min(1.25*m,0.85))
hold off, h=colorbar;
set(h,'YTickLabel',num2str(str2num(get(h,...
'YTickLabel'))/scale))

Example 2. Surface Gradients
Example 2.

Example 3. (cf. Two-Dimensional Quiver Plots)

[x,y,z]=peaks(-2:0.05:2);
contour(x,y,z,10); h=gca;
ndx=1:4:length(x);
x=x(ndx,ndx); y=y(ndx,ndx); z=z(ndx,ndx);
[u,v]=gradient(z,0.2);
p1=[x(:),y(:)]; u=u(:); v=v(:);
m=abs(u+i*v); % gradient magnitude
hold on, daspect([1 1 1]), scale=0.025;
arrow3(p1,p1+scale*[u,v],'|',0.9*m/max(m))
hold off, title('Dual Colormap')
h1=colorbar; h2=copyobj(h1,gcf);
set(h,'color',0.4*[1 1 1],...
'position',get(h,'position')+[0.05 0 0 0])
% ----------------------------------------------------- Right Colorbar
p=get(h1,'position');
set(h1,'position',[0.87,p(2),0.05,p(4)],...
'YTickLabel',num2str(str2num(get(h1,...
'YTickLabel'))/scale))                               % freeze labels
map=colormap; cdata=reshape(map,size(map,1),1,3);
set(get(h1,'children'),'cdata',cdata)              % freeze colors
set(get(h1,'title'),'string','Gradient')
% ------------------------------------------------------ Left Colorbar
caxis auto, colorbar, colormap autumn
set(h2,'YAxisLocation','left',...
'position',[0.08,p(2),0.05,p(4)])
set(get(h2,'title'),'string','Contour')

Example 3. Dual Colormap
Example 3.
Three-Dimensional Quiver Plots
Example 4. (cf. Quiver3)

[x,y]=meshgrid(-2:0.25:2,-1:0.25:1);
z=x.*exp(-x.^2-y.^2);
surf(x,y,z,'EdgeColor','none')
axis([-2.5 2.5 -1.5 1.5 -1 1]), daspect([1 1 1])
set(gca,'CameraViewAngle',7)
[u,v,w]=surfnorm(x,y,z);
p1=[x(:),y(:),z(:)]; N=[u(:),v(:),w(:)];
hold on
arrow3(p1,p1+0.5*N,'1.5_b')
hold off, colormap spring
light('position',[-1 -1 -1],'style','local')
light('position',[0 1 1]), lighting gouraud

Example 4. Surface Normals
Example 4.

Example 5. (cf. Three-Dimensional Quiver Plots)

vx=2; vy=3; vz=10; a=-32; t=(0:.1:1)';
x=vx*t; y=vy*t; z=vz*t+1/2*a*t.^2;
r=[x,y,z]; v=gradient(r')';
axis([0 3 0 4 -10 2]), pbaspect([2 1 1])
hold on, grid on, view([70 18])
arrow3(r,r+v,'b',0.9)
hold off

Example 5. Velocity Vectors
Example 5.

Cone Plots
Example 6. (cf. Coneplot)

load wind, wind_speed=sqrt(u.^2+v.^2+w.^2);
xmin=min(x(:)); xmax=max(x(:));
ymin=min(y(:)); ymax=max(y(:));
zmin=min(z(:));
hsurfaces=slice(x,y,z,wind_speed,[xmin,xmax],ymax,zmin);
set(hsurfaces,'FaceColor','interp','EdgeColor','none',...
              'AmbientStrength',0.6)
hold on, daspect([2,2,1]), view(30,40), axis tight
xrange=linspace(xmin,xmax,8);
yrange=linspace(ymin,ymax,8);
[cx cy cz]=meshgrid(xrange,yrange,3:4:15);
ui=interp3(x,y,z,u,cx,cy,cz);
vi=interp3(x,y,z,v,cx,cy,cz);
wi=interp3(x,y,z,w,cx,cy,cz);
m=sqrt(ui.^2+vi.^2+wi.^2); m=m(:)/30;
p1=[cx(:),cy(:),cz(:)]; p2=[ui(:),vi(:),wi(:)];
hcones=arrow3(p1,p2,'r',m,3*m,'cone');
set(hcones,'DiffuseStrength',0.8)
hold off, axis off
camproj perspective, camzoom(1.2)
camlight right, lighting phong

Example 6. Wind Speed and Direction
Example 6.

Example 7. (cf. Vector Field Displayed with Cone Plots)

load wind, wind_speed=sqrt(u.^2+v.^2+w.^2);
hiso=patch(isosurface(x,y,z,wind_speed,40),...
          'FaceColor','red','EdgeColor','none');
isonormals(x,y,z,wind_speed,hiso)
hcap=patch(isocaps(x,y,z,wind_speed,40),...
          'FaceColor','interp','EdgeColor','none',...
          'AmbientStrength',0.6);
hold on, colormap hsv, daspect([1,1,1]), view(65,45)
[f verts]=reducepatch(isosurface(x,y,z,wind_speed,30),0.07);
cx=verts(:,1); cy=verts(:,2); cz=verts(:,3);
ui=interp3(x,y,z,u,cx,cy,cz);
vi=interp3(x,y,z,v,cx,cy,cz);
wi=interp3(x,y,z,w,cx,cy,cz);
m=sqrt(ui.^2+vi.^2+wi.^2); m=m(:)/40;
p1=[cx(:),cy(:),cz(:)]; p2=[ui(:),vi(:),wi(:)];
h1=arrow3(p1,p2,'b',m,3*m,'cone');
xrange=linspace(min(x(:)),max(x(:)),10);
yrange=linspace(min(y(:)),max(y(:)),10);
[cx,cy,cz]=meshgrid(xrange,yrange,3:4:15);
ui=interp3(x,y,z,u,cx,cy,cz);
vi=interp3(x,y,z,v,cx,cy,cz);
wi=interp3(x,y,z,w,cx,cy,cz);
m=sqrt(ui.^2 + vi.^2 + wi.^2); m=m(:)/50;
p1=[cx(:),cy(:),cz(:)]; p2=[ui(:),vi(:),wi(:)];
h2=arrow3(p1,p2,'g',m,3*m,'cone');
hold off, axis tight, box on
set(gca,'xtick',[],'ytick',[],'ztick',[])
camproj perspective, camzoom(1.2)
camlight(-45,45), lighting phong

Example 7. Wind Speed and Direction
Example 7.


Feather Plots
Example 8. (cf. Feather)

theta=(-90:10:90)'*pi/180;
n=length(theta); p1=[1:n;zeros(1,n)]';
r=2*ones(n,1); [u,v]=pol2cart(theta,r);
plot([1 n],[0 0],'r')
axis([0 20 -2 2]), daspect([8 2 1])
hold on, grid on
arrow3(p1,p1+[u,v],'r',0.9)
hold off

Example 8. Feather Plot
Example 8.

Example 9. (cf. Plotting Complex Numbers)

t=(0:0.5:10)'; s=0.05+i; Z=exp(-s*t);
n=length(Z); p1=[1:n;zeros(1,n)]';
plot([1 n],[0 0],'b')
axis([0 22 -1 1]), daspect([8 1 1])
hold on, grid on
arrow3(p1,p1+[real(Z),imag(Z)],'b',0.8)
hold off

Example 9. Complex Number Plot
Example 9.

Compass Plots
Example 10. (cf. Compass Plots)

wdir =[45 90 90 45 360 335 360 270 335 270 335 335];
knots=[ 6 6 8 6   5   9   8   8   9 10 14 12];
[x,y]=pol2cart(wdir*pi/180,knots);
polar(0,15), axis(15.9*[-1 1 -1 1])
hold on
arrow3(zeros(length(x),2),[x',y'],'r',1.25)
hold off

Example 10. Wind Speed and Direction
Example 10.

Example 11. (cf. Compass)

Z=eig(randn(20)); m=abs(Z);
R=ceil(max(m)); r=R+0.3;
polar(0,R), axis(r*[-1 1 -1 1])
hold on
arrow3(zeros(length(Z),2),[real(Z),imag(Z)],'b',1.75*m/R)
hold off

Example 11. Complex Number Plot
Example 11.

Reference Frames
Example 12. Coordinate Axes

p=[4 2 1]; axis([-6 6 0 5 -1 1])
pbaspect([2 1.5 1]), view(55,15)
hold on
arrow3(repmat([p(1:2),0],4,1),...
[0 p(2) 0;p(1) 0 0;p;0 0 0],'--o',0,0,0.5)
arrow3(zeros(3),diag([7,5,1]),'o')
arrow3([0 0 0],p,'2.5s',1.5,[],0)
hold off, axis off, camlight left
set(gca,'CameraViewAngle',4)
text(7.1,0,0,'X'), text(0,5.05,0,'Y')
text(0,0,1,'Z','VerticalAlignment','bottom',...
'HorizontalAlignment','center')

Example 12. Coordinate Axes
Example 12.
Example 13. Eigenvectors

n=500; % generate test data
newz=[ 1 1 1]/sqrt(3);
newy=[ 0 -1 1]/sqrt(2);
newx=cross(newy,newz);
R=[newx;newy;newz]; % rotation matrix
U=randn(n,3)*diag([1,2,3])*R;
u=U(:,1)+10; v=U(:,2)+20; w=U(:,3)+30;
% plot eigenvectors with mean origin
plot3(u,v,w,'r.')
axis([0 20 10 29 20 40]), daspect([1 1 1])
view([-70,15]), set(gca,'CameraViewAngle',8)
hold on, grid on
mu=repmat(mean([u,v,w]),3,1);
[V,D]=eig(cov([u,v,w]));
p=9*V'+mu;
arrow3(mu,p,[],1.25,[],0)
p=p+ones(3);
text(p(1,1),p(1,2),p(1,3),'bfU')
text(p(2,1),p(2,2),p(2,3),'bfV')
text(p(3,1),p(3,2),p(3,3),'bfW')
hold off

Example 13. Eigenvectors
Example 13.

Example 14. (cf. Frenet)

z=(0:2/99:2)'; t=2*pi*z; r=2+z;
x=r.*cos(t); y=r.*sin(t);
R=[x,y,z]; dR=gradient(R')';
m=repmat(sqrt(sum(dR.*dR,2)),1,3);
T=dR./m; dT=gradient(T')';
m=repmat(sqrt(sum(dT.*dT,2)),1,3);
N=dT./m; B=cross(T,N);
plot3(x,y,z,'color',0.5*[1 1 1])
axis([-4 4 -4 4 0 2.5])
hold on, grid on, pbaspect([1 1 1])
set(gca,'CameraViewAngle',7)
plot3(0,0,0,'r'), plot3(0,0,0,'color',[0 0.5 0])
plot3(0,0,0,'b')
legend('Curve','Tangent','Normal','Binormal',2)
ndx=1:4:length(x); R=R(ndx,:);
arrow3(R,R+1.4*T(ndx,:),'r',0.9)
arrow3(R,R+1.4*N(ndx,:),'e',0.9)
arrow3(R,R+0.4*B(ndx,:),'b',0.9)
hold off, view(-50,5)
set(gcf,'renderer','zbuffer')

Example 14. Frenet Frames
Example 14.
Named Colors
Example 15. Named Color Table

arrow3('colors',0.4) % 24 named colors with 44 adjustable shades

Example 15. Named Color Table
Example 15.


Color Equivalencies
ColorOrder            Arrow3
Simulink          Arrow3
Color1BlueLightBlueaZure
Color2EvergreenDarkGreenAsparagus
Color3RedOrangekumQuat
Color4Sky blueGrayLight gray
Color5Violet
Color6Pear
Color7Dark gray

胖大星越来越胖
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
用matlab绘制箭头
yyytucj的博客
12-27 1834
3．Annotation，可绘出比较漂亮的箭头，但似乎只能画出二维的箭头，向量长度小于1，并且随对图形放大和平移时剪头不会相应地变化。2. text命令中的\leftarrow和\rightarrow.可用于对图形的标注，但作用仅限于此，不好控制。本例只绘出了一个方向的箭头，若要绘任意向量的剪头，可采用类似的思路,本例仅供参考。4．若要绘制三维剪头且对图形放大和旋转时剪头也随之变化，可自己编程解决。网上的一些代码，可画出漂亮的剪头。用matlab绘制箭头。
Coneplot2:与 coneplot 相同，但使用 quiver5 绘制带有 3d 箭头的箭头-matlab开发
06-01
Coneplot 函数的稍微修改版本，它使用函数 quiver5 以获得更好的 3D 箭头
MATLAB绘制二维和三维箭头
最新发布
刺客
07-24 3215
实践中发现MATLAB没有直接绘制箭头(主要是二维)的函数，所以将自己绘制二维和三维箭头的解决方案和大家分享，如有错误或更高效的方法请大家留言指教。
MATLAB在图中的画一条箭头方向的线
李宏宗的博客
03-21 1819
使用：annotation函数。
matlab绘箭头,matlab绘制箭头
weixin_29200485的博客
03-17 578
一、二维箭头1．调用annotation函数绘制二维箭头annotation函数用来在当前图形窗口建立注释对象(annotation对象)，它的调用格式如下：(1) annotation(annotation_type) % 以指定的对象类型，使用默认属性值建立注释对象。(2) annotation('line',x,y) % 建立从(x(1), y(1))到(x(2), y(2))的线注释对象。...
&nbsp;基于Matlab的辗转相除法
01-30
&nbsp;辗转相除法是整数和多项式理论中求最大公因数和最大公因式的一类重要方法，对于较大的两个整数和次数较高的两个多项式而言，利用辗转相除法手动计算它们的最大公因数和最大公因子运算量非常大，基于减少运算...
MATLAB画箭头
12-26
MATLAB画箭头是MATLAB中的一种常见的图形绘制方式，用于在图像上绘制带箭头的线段。 在MATLAB中，画箭头可以使用自定义函数m_arrow来实现。函数m_arrow的功能是画箭头，需要两个输入参数PointStart和PointEnd，分别...
&nbsp;功率方向元件的MATLAB仿真
01-30
&nbsp;目前，随着舰船综合电力推进系统的快速发展，环形区域配电技术已经成为当前舰船电力推进技术研究的一个重要内容。本文以此为背景，介绍了环形区域配电技术中的相关保护性问题及其重要性，主要分析了其中的方向...
&nbsp;基于MATLAB的飞船测控系统仿真
01-29
&nbsp;为了更好地描述地面监测站对天空中卫星的监测情况，首先建立了椭圆形模型来描述地面监测站的覆盖区域，后来根据飞船轨道的调整逐渐近似为圆形轨道模型。通过对模型的求解，计算出了实现全程监控所需要的测站...
Matlab绘制带箭头坐标轴图形
10-15
Matlab绘制带箭头坐标轴图形,从http://www.matlabsky.com/thread-39948-1-1.html下载
matlab开发-Arrow3version458
08-27
matlab开发-Arrow3version458。从A点到B点画箭头。
matalab 常用图像操作
12-02
matalab 常用图像操作matalab 常用图像操作matalab 常用图像操作matalab 常用图像操作
matlab画箭头函数annotation
xioaming159的博客
12-02 7093
Matlab中箭头的画法annotation画箭头的用法annotation例子 annotation画箭头的用法 Matlab提供了自带的箭头函数annotation，使用格式如下所示。 annotation(‘arrow’,x,y) 上式中x，y的大小要在0~1之间，主要是由于annotation箭头开始和结束的方向并不是对应直角坐标系下的坐标，而是以图像的左下角为（0，0），右上角（1，1）. 因此x，y在非0~1之间都会报错 annotation例子 画一个坐标（0，0）到（2，2）的箭头 an
Simulink 学习（一）
arong1417的博客
06-16 8986
matlab中的simulink仿真学习，是学习笔记，不是经验分享，仅是个人整理学习大佬们的学习资料。
Simulink学习——stateflow快速入门之层次结构、互斥及并行
一大块肥皂的博客
09-21 9195
层级结构 Stateflow中，可以在状态中创建状态，再创建状态，无限套娃。。。如下如图，我在一个状态框里面又画了两个状态框。这里面，外面的框被称为父级，里面的状态被称为子级。其运行顺序如红线所示，里外两个框都有一个大圆点和箭头，这个大黑点与箭头就是stateflow运行的起始位置。 如下图，可以进行不断地套娃： 那么层级结构有何作用？这个很难简单讲清楚，这个我感觉各个使用者的体会可能会有所不同，我这里讲一下我的感受，用Demo举例说明可能更好理解。 如上，一个CD机有多重功..
matlab画一个箭头，箭头两段要落在圆上
weixin_42612405的博客
02-13 517
可以使用 Matlab 的绘图函数 quiver 来绘制箭头，并使用 hold on 和 drawCircle 函数来保证箭头两端落在圆上。 例如： % Define arrow endpointand direction start = [0,0]; direction = [1,1]; endpoint = start + direction; % Plot arrow on top of ...
如何用matlab画出有向图（带箭头）
热门推荐
奔流聚海
09-02 2万+
close all,clear,clc; cm =[       0     10     0     18     0     0     0;      0     0      0     0      8     0     0;      25    0      0     0      0     9     0;      0     0      0     0    ...
利用matlab如何在图形中绘制箭头
iqizheng的专栏
12-12 1万+
一、二维箭头 1．调用annotation函数绘制二维箭头annotation函数用来在当前图形窗口建立注释对象（annotation对象），它的调用格式如下： (1) annotation(annotation_type)  % 以指定的对象类型，使用默认属性值建立注释对象。 (2) annotation('line',x,y)       % 建立从(x(1), y(1))到(x(2),
matlab怎么在曲线上画箭头
05-25
MATLAB中可以使用“quiver”函数在曲线上画箭头。 具体步骤如下： 1. 绘制曲线。 2. 使用“diff”函数计算曲线上每个点的斜率。 3. 使用“quiver”函数在每个点上绘制箭头，箭头的方向与斜率相同。 下面是一个简单的示例代码： ```matlab % 绘制曲线 x = linspace(0, 2*pi, 100); y = sin(x); plot(x, y); % 计算斜率 dx = diff(x); dy = diff(y); slope = dy./dx; % 绘制箭头 hold on; quiver(x(1:end-1), y(1:end-1), dx, dy, 'AutoScale', 'off'); ``` 在这个示例代码中，我们首先绘制了一个正弦曲线。然后，我们使用“diff”函数计算了曲线上每个点的斜率，并将其存储在变量“slope”中。最后，我们使用“quiver”函数在每个点上绘制箭头，箭头的方向与斜率相同。注意，我们在“quiver”函数中设置了“AutoScale”参数为“off”，以确保箭头大小与曲线长度相同。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值