- 积分
- 17
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2021-10-16
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
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;
|
|