我在下面的代码中尝试从非零数字开始循环“number_panels”和“number_turbines”。
更具体地说,我正在尝试 3000 到 4000 个“number_panels”,面板间隔为 500 和从 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
当我这样做时,我得到一个 size(hourly_deficit) = 4001,9,24,2 而我期望并试图达到 3,4,24,2 大小。有谁知道我哪里出错了?
最佳答案
变量 number_panels
的值从 3000 开始,而不是 0 或 1。因此,当您使用该变量作为索引值对矩阵进行索引时,Matlab 认为您需要第 3001 个索引,因此给您一个包含 3000 个零的矩阵,第 3001 个零设置为您要求的值。
如果您跟随循环进入下一个循环,number_panels
的值将变为 3500。根据您的代码,您现在的索引为 3501。这意味着从 3002 到 3500 的所有位置都将用零填充,而 3501 将设置为您给它的任何值。
相同的逻辑适用于 number_turbines
,唯一的区别是您将按 1 进行索引,而不是像使用 number_panels
时那样按 500 进行索引。
如果您想恢复到您期望的大小矩阵,您需要修改调用索引值的方式。这可以通过多种方式完成。你可以在 for 循环中有一个计数器,或者你可以使用模数数学。当您使用的步长不是 1
时,模数数学运算效果不佳。当您的索引值是起始索引的倍数时,它也不起作用。
您必须弄清楚在那个领域什么最适合您。尤其是因为您想使用的步长不是 1。但是对于从 5 到 8 的 number_turbines
,您可以使用 number_turbines - 5 + 1
进行简单索引或更简洁的 number_turbines-4
。
为清楚起见,这里是您提供的代码,其中进行了必要的调整以显示评论中提到的内容的用法。请注意,如果您从 5
开始,则需要修改 number_turbines
索引值的 -4
。另请注意,您现在需要索引 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
关于matlab - 从非零值开始 matlab for 循环,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/11615140/