如何求一个三维数组里其中一维的滑动平均

abc111 · 发表于 2016-4-25 21:39:57

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

如何求全球海温360×180×110，求110年海温的滑动平均？

5inklbj · 发表于 2016-4-26 08:46:25

Matlab做循环效率低，能避免就避免。
有时可以用矩阵计算，避免循环，比如这个：

%声明一个SST数组
SST = random('norm',14,10,[360 180 110]);

%% 用矩阵的思想做三点平滑
%三点滑动平均 SST(i) = SST(i-1)*0.25 + SST(i)*0.5 + SST(i+1)*0.25 ;
tmp1 = SST(:,:,1:110-2);
tmp2 = SST(:,:,2:110-1);
tmp3 = SST(:,:,3:110-0);

SST_AfterFilter = tmp1*0.25 + tmp2*0.5 + tmp3*0.25;

oucdaqi · 发表于 2016-4-25 22:13:04

做个格点的循环，smooth函数求滑动平均

weinihou · 发表于 2016-4-26 09:46:20

5inklbj 发表于 2016-4-26 08:46
Matlab做循环效率低，能避免就避免。
有时可以用矩阵计算，避免循环，比如这个：

可以手动点个赞么？反正我是点了~~

abc111 · 发表于 2016-4-26 19:08:16

本帖最后由 abc111 于 2016-4-26 19:13 编辑

5inklbj 发表于 2016-4-26 08:46
Matlab做循环效率低，能避免就避免。
有时可以用矩阵计算，避免循环，比如这个：

你好，我需要做21点滑动平均，这样矩阵是不是写起来会很麻烦？

abc111 · 发表于 2016-4-26 19:11:35

oucdaqi 发表于 2016-4-25 22:13
做个格点的循环，smooth函数求滑动平均

for i=1:110;
cc(:,:,i)=smooth(tte(:,:,i),21);
end

tte是360X180X110的数据
对每个格点的110个数据做21年滑动平均，这样对吗？

5inklbj · 发表于 2016-4-26 20:56:34

abc111 发表于 2016-4-26 19:08
你好，我需要做21点滑动平均，这样矩阵是不是写起来会很麻烦？

还好吧，tmp那部分可以写循环。另外，smooth和我用的平滑方法是不一样的，LZ可以搜搜，看看你究竟想要哪种。

5inklbj · 发表于 2016-4-26 20:56:39

abc111 发表于 2016-4-26 19:08
你好，我需要做21点滑动平均，这样矩阵是不是写起来会很麻烦？

还好吧，tmp那部分可以写循环。另外，smooth和我用的平滑方法是不一样的，LZ可以搜搜，看看你究竟想要哪种。

oucdaqi · 发表于 2016-4-27 00:07:34

对格点做循环呀
smth_sst=zeros(360,180,110);
for i=1:360
   for j=1:180
         smth_sst(i,j,:)=smooth(sst(i,j,:),21);
   end
end
确实循环起来Matlab内存占用很大

丹丹快乐023 · 发表于 2016-11-28 11:33:27

请问楼主做出来了么

		自动登录	找回密码
密码			立即注册

[程序设计] 如何求一个三维数组里其中一维的滑动平均

登录后查看更多精彩内容~

相关帖子