请帮忙看看如何提高这段循环积分代码的效率
因为前面要做符号运算，得到表达式，然后后面要用一个嵌套的二次循环积分计算，因此其代码效率是相当低的，用profview分析，发现主要是inline占用了大量的时间。由于inline/eval的低效，在循环中调用inline或者eval将占用了绝大部分时间，但这里积分却又不得不调用，下式中的fun=inline(vectorize(fai))可以直接写成fun=vectorize(fai)，但实际上在积分时依然会用到eval。
由于int是求精确解析解的，效率自然更加低效。
这个代码运行10多个小时也没完，还不一定出结果，最后还可能出错。有没有什么办法解决这个问题？有没有更好的积分办法？
syms x; num=0; N=200; LL=sym(zeros(1,N));
for i=0:N-1 %这里用maple生成一个拉盖尔函数数组供下面调用
LL(i+1)=maple(['L(',num2str(i),',1,x)']);
end
for j=0:N-1, for i=j:N-1
num=num+1; fai=exp(-x)*x^m*LL(i+1)*LL(j+1);
fun=inline(vectorize(fai)); %这里的inline可以不写，实际quadl时会自然调用
N1=quadl(fun,0,l0);
% 下面还有一些处理代码，略去
end; end
敬请各位指点。
上面还有个初始条件： m=1;
[ 本帖最后由 ChaChing 于 2010-5-29 00:41 编辑 ]
页:
[1]