predictions1.m

3个关于预测控制的matlab函数

function [Phi, G]=predictions1(A,B,C,N1,N2,Nu)
Phii=[]; m=[]; Gi = zeros(N2-N1+1, N2-N1+1);    
for i=1:(N2-N1+1)
    Phii=[Phii;C*A^i];        
        for j=1:i %(N2-N1+1)   %Nu        
            if ((i>1) && (j>1))
                Gi(i,j) = Gi(i-1,j-1);
            else
                Gi(i,j)=C*(A^(i-1))*B;
            end
        end
    if (i>100) 
        break
    end
end 
Phi=Phii;
G=Gi(N1:N2,:);

% A = 3
% B = 2
% C = 1
% N1 = 1
% N2 = 5
% Nu = 3
% [Phi, G] = predictionD(A,B,C,N1,N2,Nu)
%