我有下面的代码,我试图从非零数字开始循环“number_panels”和“number_turbines”。
更具体地说,我正在尝试使用500个面板间隔的3000到4000“number_panels” 具有一个涡轮机间隔的5到8“number_turbines” 即
number_of_days = 2;
for number_panels = 3000:500:4000 % range of PV panel units examined
for number_turbines = 5:8 % range of wind turbine units examined
for h=1:24 %# hours
for d = 1:number_of_days %# which day
n = h + 24*(d-1);
% hourly_deficit_1(...,..., h, d)= Demand(n)-(PV_supply(n)... %
hourly_deficit(number_panels + 1, number_turbines + 1, h,d) = hourly_annual_demand(n) - (hourly_annual_PV(n)*number_panels) - (hourly_annual_WT(n)*number_turbines);% hourly power deficit (RES supply with demand)
if hourly_deficit(number_panels + 1, number_turbines + 1, h,d)< 0 % zero out negative hourly deficit values (this is power surplus from RES)
hourly_deficit(number_panels + 1, number_turbines + 1, h,d) = 0;
end
当我这样做时,我得到一个尺寸(hourly_deficit)= 4001,9,24,2,而我期待并试图达到3,4,24,2尺寸。有谁知道我哪里出错了?
答案 0 :(得分:5)
变量number_panels
的值从3000而不是0或1开始。因此,当您使用该变量索引矩阵作为索引值时,Matlab认为您需要第3001个索引并因此为您提供一个矩阵是3000零,3001是你所要求的。
如果你按照循环进入下一个循环,number_panels
的值就变为3500.你现在根据代码索引为3501。这意味着从3002到3500的所有位置都将填充零,3501将被设置为您给它的任何值。
同样的逻辑同样适用于number_turbines
唯一的区别是,您将使用number_panels
进行索引,而不是像1
那样使用500。
如果你想回到大小矩阵,你需要修改你调用索引值的方式。这可以通过多种方式完成。您可以在for循环中使用计数器,也可以使用模数学。当您使用的步长不是number_turbines
时,模数数学效果不佳。当你到达索引值是起始索引的倍数时,它也不起作用。
你必须找出在这个领域最适合你的东西。特别是因为你想使用不是1的步长。但对于从5到8的number_turbines - 5 + 1
,你可以使用number_turbines-4
或更简洁的-4
进行简单索引。 / p>
为清楚起见,以下是您提供的代码以及必要的调整,以显示评论中提到的内容的使用。请注意,如果您从number_turbines
以外的其他地方开始,则需要修改5
索引值的number_panels
。另请注意,您现在需要对number_of_days = 2;
number_panels = 3000:500:4000;
for idx_number_panels = 1:length(number_panels) % range of PV panel units examined
for number_turbines = 5:8 % range of wind turbine units examined
for h=1:24 %# hours
for d = 1:number_of_days %# which day
n = h + 24*(d-1);
% hourly_deficit_1(...,..., h, d)= Demand(n)-(PV_supply(n)... %
hourly_deficit(idx_number_panels , number_turbines -4, h,d) = hourly_annual_demand(n) - (hourly_annual_PV(n)*number_panels(idx_number_panels)) - (hourly_annual_WT(n)*number_turbines);% hourly power deficit (RES supply with demand)
if hourly_deficit(idx_number_panels, number_turbines -4, h,d)< 0 % zero out negative hourly deficit values (this is power surplus from RES)
hourly_deficit(idx_number_panels, number_turbines -4, h,d) = 0;
end
end
end
end
end
向量进行索引,因为它不是循环值。
希望这有帮助!
{{1}}