假设有东海平均（120-130E,24-33N）风速某年的时间序列资料为1年，u_t, v_t。一般u_t, v_t的范围分别不超过[-12 12]。

（1）作图时，容易犯的第一个毛病，就是忽视了x和y轴的比例，导致风速矢量的方向看上去不对，譬如：

figure

set(gcf, 'PaperUnits', 'inches');

set(gcf, 'Renderer', 'painters');   % fixmatlab datetick bug in opengl

set(gcf, 'paperposition',[0.4 0.45.0 2.2]);

new_axis = axes('position',[0.10 0.11 0.870.86]);

hwvec = feather(u_t, v_t,'k-');

for i =1:length(hwvec)-1 %the lastline on the axes

 xData = get(hwvec(i),'XData');

 yData = get(hwvec(i),'YData');

 set(hwvec(i),'XData',xData(1:2),'YData',yData(1:2))

end

axis([ 0 365 -12 12])

set(gca,'YTick',[-12:4.0: 12]);

ylabel('Wind (m/s)');

图1 x和y轴比例失当，导致风速矢量看上去方向不合理

正常情况下，东海平均风场应常见偏南风、东南风、西南风、东北风或者西北风，而上图显示基本是正南风和正北风。不是数据有问题，而是数据展示的时候有问题，原因是坐标轴的比例不当。因为feather画图是，自动把x和y的输入都当做一个单位，然后按比例。譬如风速的范围（也就是y轴范围基本是-12到12），而时间范围为1年365天。相对而言，沿x轴方向的风速东分量就看不出来。

（2）解决的方法是加上一个比例，固定x和y轴的比例为1:1。两种方法

Axis equal; %然后

Axis tight

或者

daspect([1 1 1]);

此时的作图结果为：

          图2风速矢量图形的坐标轴比例合理了，但是图形没法看了，因为y轴被明显压扁了

图2的做法比（1）合理，至少风的方向显示符合我们的常识了。譬如250-300天之间的东北风就比图1更合理。然而，图2的结果显然无法放在论文中，因为看不清楚。

（3）解决方法是在放大y轴的数据，比例系数自己选择

图3方法y轴比例，把u, v各放大5倍后的矢量图

ratio        = 5;        % 比例尺

y_range_orig = [-12 12]; %原始比例

y_range      = ratio*y_range_orig;   %新比例

figure

set(gcf, 'PaperUnits', 'inches');

set(gcf, 'Renderer', 'painters');   % fixmatlab datetick bug in opengl

set(gcf, 'paperposition',[0.4 0.45.0 2.2]);

new_axis = axes('position',[0.10 0.11 0.870.86]);

hwvec = feather(ratio*u_t,ratio*v_t,'k-');

set(hwvec,'clip','off'); % 图形在坐标轴之上

daspect([1 1 1]);

for i =1:length(hwvec)-1 %the lastline on the axes

 xData = get(hwvec(i),'XData');

 yData = get(hwvec(i),'YData');

 set(hwvec(i),'XData',xData(1:2),'YData',yData(1:2))

end%for

axis([ 0 365 y_range(1)y_range(2)])

htext    = text(13,10,'(B)');

set(gca,'YTick',[y_range(1):4.0*ratio:y_range(2)],...

'YTicklabel',[y_range_orig(1):4.0:y_range_orig(2)]);

set(gca,'box','on','fontsize',8);

ylabel('Wind (m/s)');

结果合理了。

图3也可以用d:DocumentsCrossShelfFiguregraphicstimeplt.m来画。画图函数timeplt.m，该数据包来自：http://woodshole.er.usgs.gov/operations/sea-mat/timeplt-html/

结果如下：

                             图4. timeplt作的风速分量和矢量图

Happy programming!