for i=1:4 p(i) = ps(i) .* (1 - phi(i)) + phi(i) .* pf(i); c(i) = cs(i) .* (1 - phi(i)) + phi(i) .* cf(i); lamdpor(i) = lamds(i) .* (1 - phi(i)) + phi(i) .* lamdaf(i); pcpor(i) = (1-phi(i)).*ps(i).*cs(i)+phi(i) .* pf(i).* cf(i); pcf(i)=pf(i).*cf(i); vt(i) = u(i) .* pcf(i) ./ pcpor(i); % apor(i) = lamdpor(i) ./ c(i) ./ p(i); apor(i) = lamdpor(i) ./ pcpor(i); end向量化

可以将该段代码向量化，使用 MATLAB 中的矩阵运算来实现循环：

p = ps .* (1 - phi) + phi .* pf;
c = cs .* (1 - phi) + phi .* cf;
lamdpor = lamds .* (1 - phi) + phi .* lamdaf;
pcpor = (1 - phi) .* ps .* cs + phi .* pf .* cf;
pcf = pf .* cf;
vt = u .* pcf ./ pcpor;
apor = lamdpor ./ pcpor;

这样不仅可以减少代码的行数，也可以提高代码的运行速度。