各位好，不好意思，
這應該是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
推
07/29 22:22, , 1F
07/29 22:22, 1F