Summing and multiplying matrices of different size

Question

Gbeminiyi Oyedele 2022-11-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1852543-summing-and-multiplying-matrices-of-different-size

评论： Gbeminiyi Oyedele 2022-11-16

采纳的回答： Bjorn Gustavsson

在 MATLAB Online 中打开

So, I have this equation:

The size of C is a 9 by 9 matrice, and X ,Y and sigma parameter is a matrice of size 1 by 9 matrice.

I have this code:

%% Define the function that output classes
% The function will take parameter values, initial condition of the ODE and
% the expected time to run the model.
function [Classes] = model2(para,ICs,maxtime)
opts = odeset('RelTol',1e-4,'AbsTol',1e-3); %tolerance level
tspan = [0:1:maxtime]; %time span
agegroup= 9;
%Shape of the initial condition
All = reshape([ICs.S'; ICs.E'; ICs.I'; ICs.A';ICs.H'; ICs.R'; ICs.D'], 7*agegroup,1);
[t , pop] = ode15s(@(t,y)diff_socialmodel(t,y,para),tspan,All,opts);
%% Define the structure of the output
Classes = struct('S',pop(:,1:1:9),'E',pop(:,10:1:18),'I',pop(:,19:1:27), ...
    'A',pop(:,28:1:36),'H',pop(:,37:1:45),'R',pop(:,46:1:54), ...
    'D',pop(:,55:1:63),'t',t);
%% Solving each the ODE at every time step
function dpop = diff_socialmodel(t,pop,para)
S = pop(1:1:9);
E = pop(10:1:18);
I = pop(19:1:27);
A = pop(28:1:36);
H = pop(37:1:45);
R = pop(46:1:54);
D = pop(55:1:63);
%% Store my output here
dpop = zeros(7*para.agegroup,1);
%% High Income population
dpop(1:1:9) = - (para.beta*para.a*para.l*(para.h.*I+A))*S;
dpop(10:1:18) =  (para.beta*para.a*para.l*(para.h*I+A))*S- para.sigma.*E;
dpop(19:1:27) = para.psi*para.sigma*E + para.epsilon*para.gamma*A - para.eta*I;
dpop(28:1:36) = (1-para.psi)*para.sigma*E - para.gamma*A;
dpop(37:1:45) = para.phi*para.eta*I -para.delta*H;
dpop(46:1:54) = (1-para.epsilon)*para.gamma*A + (1-para.phi)*para.eta*I + para.delta*para.d*H;
dpop(55:1:63) = (1-para.d)*para.delta*H;
end
end

It appears to work but the result is not good. I think I am doing something wrong.

Maybe my code can give a better perspective.

When I change the matrix l, the output yield the same result, which is strange. Each time I change matrix l, I expect a different output.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Davide Masiello 2022-11-16

Please share all of your code.

请先登录，再进行评论。

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

Answer 1

Bjorn Gustavsson 2022-11-16

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1852543-summing-and-multiplying-matrices-of-different-size#answer_1101473

在 MATLAB Online 中打开

Check that you get what you expect from each product and term separately. Then figure out what you should do to fix your discrepancies. It seems likely that the product between the two factors should be an element-wise product:

dpop(1:1:9) = - ((beta*sigma.*S).*(C*(para.h*I+A)));

But that is mainly guess-work.