Question

我正在寻找摆脱以下循环的方法（在theta上）：

for i=1:1:length(theta)

  V2_ = kV2*cos(theta(i)); 
  X = X0+V2_;
  Y = Y0-V2_*(k1-k2);
  Z = sqrt(X.^2-Z0-4*V2_.*(k.^2*D1+k1));


  pktheta(:,i)=exp(-t/2*V2_).*(cosh(t/2*Z)+...
            Y./((k1+k2)*Z).*sinh(t/2*Z));  


end

其中X0，Y0，Z0和kV2取决于矢量k（相同大小）。 t，D1，k1和k2是数字。由于我必须多次经历这个循环，我怎样才能加速它？感谢

Answer 1

试试这个 -

N = numel(theta);

V2_ = kV2*cos(theta(1:N));
X0 = repmat(X0,[1 N]);
Y0 = repmat(Y0,[1 N]);
Z0 = repmat(Z0,[1 N]);

X = X0 + V2_;
Y = Y0-V2_*(k1-k2);
Z = sqrt(X.^2-Z0-4.*V2_ .* repmat(((1:N).^2)*D1 + k1.*ones(1,N),[size(X0,1) 1]));
pktheta = exp(-t/2*V2_).*(cosh(t/2*Z) + Y./((k1+k2)*Z).*sinh(t/2*Z));

如果有人可以用它发布，肯定BSXFUN必须更快。

在matlab中矢量化

1 个答案: