作者fyq (推销员)
看板MATLAB
标题[运算] 矩阵求取固定间格平均问题
时间Mon Jul 29 17:45:32 2013
各位好,不好意思,
这应该是matlab里面很基本的问题,但因为我是新手,
虽然爬过文了却仍然无法理解这样的问题该如何执行。
首先我有一段ADCP(开船时藉由都普勒流速器测得的船下海流资料,共60层)的资料,
我取其中371*60的矩阵
371对应到另外一个矩阵是船所开行的纬度。
我藉由pcolor和contour画出这段航段的流速剖面
figure():clf;
>> gca=pcolor(AnFLatDeg(1110:1480),y,(SerEmmpersec(1110:1480,1:60)/1000)');
%AnFLatDeg是纬度资料,y为深度,SerEmmpersec为u速度
>> axis ij
>> caxis([-0.5 2])
>> set(gca, 'LineStyle','none')
>> hold on
>>
contour(AnFLatDeg(1110:1480),y,(SerEmmpersec(1110:1480,1:60)/1000)',
'k','linewidth',3);
如图
https://www.facebook.com/photo.php?fbid=583108701739297&set=a.189176164465888.61098.100001205170374&type=1&theater
但是,由於这样的contour很不平滑,太多杂讯
我需要将资料每0.04度左右做一个平均
变成一个平滑的曲线
(如果pcolor也可以平滑当然也很棒)
於是我写了一个loop(我真的很不会写loop)
>> for k=1:37
my1(k)=mean(y1(10*(k-1)+1:10*k));
my2(k)=mean(y2(10*(k-1)+1:10*k));
end
y1是SerEmmpersec(1110:1480,1:60)/1000
y2为AnFLatDeg(1110:1480)
这边事先尝试做每十个资料点做一个平均 因为资料长度为371 所以我取37段
但是出来的
并不是37*60这样的矩阵
而是
1*37.....
真是谢谢各位能够帮忙我理解这个问题并且做出成功的每0.04度的平均
麻烦各位多多帮忙!
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.68.57
1F:推 mitsuyo:内插 interp 07/29 22:22