Question

matlab中替代for循环

clear all;
clc;
C=0;T=1;m=1;P=1;
figure;
for tau=-3:0.01:3;
fun=@(t) ((gaussian(t,tau,T,P,m,C)+gaussian(t,0,T,P,m,C)).*(gaussian(t,tau,T,P,m,-C)+gaussian(t,0,T,P,m,-C))).^2;
acf=integral(fun,-10,10);
plot(tau,acf);grid on;
hold on;
end
以上是一个计算脉冲自相关的函数编程，但是计算效率比较低，想要请问一下是不是可以将for循环替换掉，毕竟matlab里面最高效的是利用数组和矩阵，感激不尽！

Answer 1

1. 把画图的语句懂for循环里去掉，改成记录每次循环的结果，最后一起画。因为画图语句是无论如何不可能加速的
   

2. 因为你要算积分，而且被积函数中只有一个参数tau是需要变化的，建议你另写一下函数：

   fun0 = @(tau, t)(...); %写在for外面

   fun = @(t)(fun(tau, t)); % 这个写在for里面

   这样也许快一点

3. 在前两步的基础上进一步去掉for。请参看integral函数的帮助文档（用doc integral），里面提到了同时积多个函数的写法。这一步优化已经可能需要改变gaussian函数的写法了，我没法再说详细。如果你愿意继续优化，可以联系我。