params.Q_x=1000*eye(params.nstates); % Weighting matrix Q for state cost params.R_u=0.1*eye(params.nctrl); % Weighting matrix R for energy cost params.umin=-2*ones(params.nctrl,1); % minimum value of ...