- 积分
- 106
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2022-7-9
- 最后登录
- 1970-1-1
|
发表于 2023-5-9 21:31:15
|
显示全部楼层
大佬写出来了吗?我也是自己写不懂对不对 for i =1:144
for j =1:72
if (~isnan(GPCC_pre(i,j,1)))
% Calculate CDF for model output and observed data
ObservedData=reshape(GPCC_pre(i,j,277:636),360,1);
ModelData1=Historypre{fileindex};
ModelData=reshape(ModelData1(i,j,277:636),360,1);
FurtureData=reshape(GCM_preSSP126(i,j,661:1020),360,1);
F_m_c = cdf("Normal",ModelData);
F_o_c = cdf("Normal",ObservedData);
F_m_p = cdf("Normal",FurtureData);
delta=F_o_c.^-1.*(F_m_p.*FurtureData)-F_m_c.^-1.*(F_m_p.*FurtureData);
GCM_preSSP126new(i,j,:)=GCM_preSSP126(i,j,661:1020)+reshape(delta,1,1,360);
end
end
end
|
|