从非零值开始matlab for循环

Question

我在下面的代码中尝试从非零数字开始“number_panels”和“number_turbines”for 循环。

更具体地说，我正在尝试 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 的大小。有谁知道我哪里出错了？

Answer 1

变量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 进行简单索引。

为清楚起见，这里是您提供的代码，并进行了必要的调整，以显示 cmets 中提到的内容的使用。请注意，如果您从5 以外的其他位置开始，则需要修改-4 的number_turbines 索引值。另请注意，您现在需要索引 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