metlab 关于for循环优化的问题

请问一下这段代码运行一下需要12秒可领到要求的时间非常短，在网上查了一下说尽量避免for循环的，那么这段代码需要如何优化才能达到要求呢

  for i=1:length(theta)

           proj = p_back(:,i);

           taxis = (1:size(p_back,1)) - ctrIdx;

           t = x.*costheta(i) + y.*sintheta(i);

           projContrib = interp1(taxis,proj,t(:),'linear');

           img = img + reshape(projContrib,N,N);

       end

可以给一下数据吗，没法测试