Zhenting Xu

now  it's very useful  thank you very much    thanks   

So How to change the code  thank you very much            this initial value should be fixed until to next iteration

My solution is wrong, the correct solution should be u=1 thank you 

%% define the class classdef OptimizationProblemVariables properties u_ini y_ini u y g slack_u slack_y end methods function obj = OptimizationProblemVariables(u_ini,y_ini,u,y,g,slack_u,slack_y) ...

function result = isclose(a, b,rel_tol,abs_tol) if nargin < 3 rel_tol = 1e-9; end if nargin < 4 abs_tol = 0.0; end result = abs(a - b) <= max(rel_tol * max(abs(a), abs(b)), abs_tol); end   func...

%% intialize the parameters in DeePC % DeePC Paramters s=1; % How many steps before we solve again the DeePC problem T_ini = 2;% Size of initial set of data %number of data points used to estimat...

%% Define the new system classdef systems properties sys x0 u y end methods function obj = systems(sys,x0) if nargin < 2 x0 = zeros(1,size(sys.A, 1)); end assert(~isempty(x0) && size(sys.A,1)==l...

%% intialize the parameters in DeePC % DeePC Paramters s=1; % How many steps before we solve again the DeePC problem T_ini = 2;% Size of initial set of data %number of data points used to estimat...

function loss = build_loss(u,y)        P = size(y, 2)    y_ref = ones(size(y)); % Q & R are the weight matrix   loss = norm(y-y_ref, 2)^2; end   function Constraints = build_constraints(u,y) ho...

function obj = build_problem(obj,... lambda_g,lambda_y,lambda_u,lambda_proj) if nargin < 2, lambda_g = 0 ; end if nargin < 3, lambda_y = 0 ; end if nargin < 4, lambda_u = 0 ; end if nargin < 5, l...