MATLAB画风速带有方向的矢量图程序,Matlab向量矢量图

前言

目前在做图像分割时，有相关向量场图的表示，顺便整理一下向量场图的matlab实现方法，整理自：胖大星越来越胖的微博，http://blog.sina.com.cn/s/blog_5e3213f30102v3o8.html，谢谢。

Two-dimensional Quiver Plots 二维向量场图

Joukowski Airfoil Transformation

相关链接：http://www.mathworks.com/matlabcentral/fileexchange/8870-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)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    代码结果图：

Surface Gradients

相关链接：http://www.mathworks.com/help/matlab/ref/quiver.html

代码示例：

[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))    代码结果图：

Dual Colormap

代码示例：

[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