«для» петли я не петли

Question

будет иметь следующий код

Я пытаюсь оптимизировать linprog над данными 8760 набор, но не похоже, чтобы получить петлю будет для каждой строки.
Когда я запускаю его, я получаю размер `A, равный 1x30 (когда это должно быть 8760 на 30).

Кто-нибудь видит, где я ошибался?

Answer 1

Да, каждый раз, когда вы запускаете его вы перезапись с одной строки [PVtech (я, :) WTtech (я, :)]

Попробуйте это: A = [A; PVtech(i,:) WTtech(i,:)]; то есть вертикальная конкатенация

С предраспределения ваш код будет выглядеть следующим образом:

numRows = 8760; 
A = zeros(numRows, 30); 

for i = 1:numRows 

    A(i,:) = [PVtech(i,:) WTtech(i,:)]; 

    b = demand(i); 

    f = [CRF*PVtechcost(i,:).*PVcap(i,:)./PVtech(i,:) CRF*WTtechcost(i,:).*WTcap(i,:)./WTtech(i,:)]; 

    x(i) = linprog(f, A,b,[], [], lb); 

end