我有一个数组(1,8760)代表每小时负载。每年的负荷是相同的,我想进行20年的计算。所以我尝试使用repmat,toeplitz和trill,以便每次都有一个带有最新负载的矩阵。我的主要目标是使用几何方案计算移动平均值,基本上只计算下面代码中的每个F.问题是我不能使用toeplitz这么大的阵列。还有其他方法吗?
Q=rand(1,8760);
QB=repmat(Q,1,20);
QG=toeplitz(QB);
QC=tril(QG);
F49=mean(QC(50:175200,49:50),2);
F50=mean(QC(54:175200,51:54),2);
F51=mean(QB(62:175200,55:62),2);
F52=mean(QB(78:175200,63:78),2);
F53=mean(QB(110:175200,79:110),2);
F54=mean(QB(174:175200,111:174),2);
F55=mean(QB(302:175200,175:302),2);
F56=mean(QB(558:175200,303:558),2);
F57=mean(QB(1070:175200,559:1070),2);
F58=mean(QB(2094:175200,1071:2094),2);
F59=mean(QB(4142:175200,2095:4142),2);
F60=mean(QB(8238:175200,4143:8238),2);
F61=mean(QB(16430:175200,8239:16430),2);
F62=mean(QB(32814:175200,16431:32814),2);
F63=mean(QB(65582:175200,32815:65582),2);
F64=mean(QB(131118:175200,65583:131118),2);