repeat the code with different value each time

Question

0 个投票

Hi

I have a matrix w= 6x222 . I have 37 window. I want to repeat the code with different w 37 times. First I use the the first 6 columns of w so I use W=w(:,[1:6]); and I have the results then I use the following 6 columns i.e.W=w(:,[7:12]); and so on untile I have the last 6 columns W=w(:,[217:222]);

How can I do this iteration and extract the 37 solutions from each W.

Thanks in advance

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Follow Question

Answer 1

Jan 2023-3-5

编辑：Jan 2023-3-5

在 MATLAB Online 中打开

0 个投票

w = rand(6, 222);
index = 1:6:size(w, 2);
result = zeros(1, numel(index));
for k = 1:numel(index)
    aw = w(:, index(k):index(k)+5);  % A 6x6 submatrix
    result = yourCalculations(aw);
end

I'd avoid the names "w" and "W", because the letters can be confused easily.

A more Matlab'ish way:

ww = reshape(w, 6, 6, []);

Now the calculations might be applicable to the 3D array ww directly without a loop. See e.g. pagemldivide or pagemtimes.

2 个评论
显示无隐藏无

Az.Sa 2023-3-5

编辑：Az.Sa 2023-3-5

在 MATLAB Online 中打开

Thank you very much!

so if my code as follow

clear all
R=readtable('return.1.indexFeb.csv');
R=table2array(R);
w=readmatrix('window.opt.w.05.csv');
W=w(:,[1:6]);
gamma = 2;
Aeq = [1 1 1 1 1 1];
beq = 1;
lb = zeros(6,1);
ub = ones(6,1);
x0=0.1667*ones(1,6);
u = @(x) 1/(1-gamma)*x.^(1-gamma); % CRRA utility 
obj = @(x)-sum(u(x*aw.'*R));
x = fmincon(obj,x0,[],[],Aeq,beq,lb,ub)

so I replace W=w(:,[1:6]) with

R=readtable('return.1.indexFeb.csv');
R=table2array(R);
w=readmatrix('window.opt.w.05.csv');
index = 1:6:size(w, 2);
result = zeros(1, numel(index));
for k = 1:numel(index)
    aw = w(:, index(k):index(k)+5);  % A 6x6 submatrix
    
gamma = 2;
Aeq = [1 1 1 1 1 1];
beq = 1;
lb = zeros(6,1);
ub = ones(6,1);
x0=0.1667*ones(1,6);
%x0=[.25,.25,.25,.25,0,0];
u = @(x) 1/(1-gamma)*x.^(1-gamma); % CRRA utility 
% u = @(x)  1 - exp(-gamma *x);
obj = @(x)-sum(u(x*aw.'*R));
options = optimoptions(@fmincon,'Algorithm','sqp','MaxIterations',1500);
x = fmincon(obj,x0,[],[],Aeq,beq,lb,ub);
end

and this is will calculate my x values corresponding to the w 37 times, correct ?

how can I save the results so I can use it in a different step.

Jan 2023-3-6

在 MATLAB Online 中打开

for k = 1:numel(index)
    ...
    x(k) = ...
end

请先登录，再进行评论。

repeat the code with different value each time

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

2 个评论
显示无隐藏无

更多回答（0 个）

类别

标签

Community Treasure Hunt

repeat the code with different value each time

0 个评论 显示 -2更早的评论 隐藏 -2更早的评论

采纳的回答

2 个评论 显示 无 隐藏 无

更多回答（0 个）

类别

标签

另请参阅

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

2 个评论
显示无隐藏无