Vectorization of a function with system of ODE's

Question

George Bashkatov 2021-4-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/805196-vectorization-of-a-function-with-system-of-ode-s

回答： Walter Roberson 2021-4-17

I want to modify famp so that it accepts a vector of y values, and reshapes the values to 2 columns, and then computes values for each row, but I don't know how I can do it. Can you help me, please?

function dydz = famp(pa,N0,pe,se,k,y,z)
k=Par.N0*pa/(pa+pe); 
dy1 = y(1).*(pa.*N0+(pa+pe).*(k.*y(1)./(y(1)+y(2)+1))); 
dy2 = y(2).*se.*(k.*y(1)./(y(1)+y(2)+1)); 
 dydz = [dy1; dy2];

What can you advice to read to understand vectorization and optimization in a whole? By the way I know about documentation, but it is not enough.

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

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2021-4-17

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/805196-vectorization-of-a-function-with-system-of-ode-s#answer_677846

在 MATLAB Online 中打开

function dydz = famp(pa,N0,pe,se,k,y,z)
y = reshape(y, [], 2);
k=Par.N0*pa/(pa+pe); 
dy1 = y(:,1).*(pa.*N0+(pa+pe).*(k.*y(:,1)./(y(:,1)+y(:,2)+1))); 
dy2 = y(:,2).*se.*(k.*y(:,1)./(y(:,1)+y(:,2)+1)); 
 dydz = reshape([dy1, dy2],[],1);