从非零值开始matlab for循环

Question

我有下面的代码,我试图从非零数字开始循环"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尺寸.有谁知道我哪里出错了？

Answer 1

变量的值number_panels从3000而不是0或1开始.因此,当您使用该变量索引矩阵作为索引值时,Matlab认为您需要第3001个索引,从而为您提供一个3000个零的矩阵,并且设置了3001个你问的问题.

如果您按照循环进入下一个循环,则值number_panels变为3500.您现在根据代码索引为3501.这意味着从3002到3500的所有位置都将填充零,3501将被设置为您给它的任何值.

同样的逻辑适用于number_turbines 唯一的区别是你将像1一样索引1而不是500 number_panels.

如果要返回大小矩阵,您需要修改调用索引值的方式.这可以通过多种方式完成.您可以在for循环中使用计数器,也可以使用模数学.当您使用非步长时,模数数学不能很好地工作1.当你到达索引值是起始索引的倍数时,它也不起作用.

你必须在那个领域找出最适合你的方法.特别是因为您希望使用不是1的步长.但是对于number_turbines从5到8 的步长,您可以使用简单索引number_turbines - 5 + 1或更简洁number_turbines-4.

为清楚起见,以下是您提供的代码以及必要的调整,以显示评论中提到的内容的使用.请注意,您将需要修改-4的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