请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3710|回复: 1

[程序设计] matlab 画标准差空间场,为什么处理的数据很小,是程序写的不对吗?

[复制链接]

新浪微博达人勋

发表于 2021-11-28 20:41:18 | 显示全部楼层 |阅读模式

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

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x

ii=2020-1901+1;
[m1,n1,l1]=size(pre0);
pre_3=zeros(m1,n1,ii);
pre_4=zeros(m1,n1,ii);
pre_5=zeros(m1,n1,ii);
for k=1:ii
    for i=1:m1
        for j=1:n1
            for im=3:3
                pre_3(i,j,k)=pre_3(i,j,k)+pre0(i,j,(k-1)*12+im);
            end
            for im=4:4
                pre_4(i,j,k)=pre_3(i,j,k)+pre0(i,j,(k-1)*12+im);
            end
            for im=5:5
                pre_5(i,j,k)=pre_5(i,j,k)+pre0(i,j,(k-1)*12+im);
            end
        end
    end
end

year=2020-1901+1;
for i=1:1:year
    pre_3o(:,:,i)=pre_3(:,:,i);
    pre_4o(:,:,i)=pre_4(:,:,i);
    pre_5o(:,:,i)=pre_5(:,:,i);
end

avepre_3=nanmean(pre_3o(:,:,:),3);
avepre_4=nanmean(pre_4o(:,:,:),3);
avepre_5=nanmean(pre_5o(:,:,:),3);
s3=0;s4=0;s5=0;pre_sum0=0;
for i=1:1:year
%     s2_pre3=s3+(pre_3o(:,:,i)-avepre_3)^2;
%     s2_pre4=s4+(pre_4o(:,:,i)-avepre_4)^2;
%     s2_pre5=s4+(pre_5o(:,:,i)-avepre_5)^2;
    s_pre_sum0=pre_sum0+(pre_3o(:,:,i)-avepre_3+pre_4o(:,:,i)-avepre_4+pre_5o(:,:,i)-avepre_5).^2;
end
s_presum=sqrt(s_pre_sum0'/year);


[X,Y]=meshgrid(lon,lat);
contourf(X,Y,s_presum);
colormap(jet);
hcd=colorbar;
set(hcd,'ytick',[0:30:900]);
hold on
load coast   %加载海岸线文件
plot(long,lat,'k-','lineWidth',1);
axis([100 125 15 35]);%调整坐标范围
fnshp_L='D:\Mathlab\R2021a\SHP-shp-文件\省界\省界\bou2_4l.shp';
fnshp_P='D:\Mathlab\R2021a\SHP-shp-文件\省界\省界\bou2_4p.shp';
readL=shaperead(fnshp_L);
readP=shaperead(fnshp_P);
mapshow(readL,'Color','k');
grid on;
hold on;


密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2021-11-29 09:02:00 | 显示全部楼层
不清楚你的变量是什么,计算标准差的话可以直接用内置函数std()的吧
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表