matlab循环积分矢量化问题

findpass

1. 
   
2. function H3jmn=H3jmn(T0,lamb,Gamma,pjun,N)
   
3. H3jmn=zeros(N,N,N);
   
4. for jj=1:N
   
5.     for n=1:N
   
6.         for m=1:N
   
7.             fH3=@(x)(Cnx(x,m,T0,lamb,Gamma,pjun).*dBnxdx(x,n,T0,lamb,Gamma,pjun)+...
   
8.                 Bnx(x,n,T0,lamb,Gamma,pjun).*dCnxdx(x,m,T0,lamb,Gamma,pjun)).*...
   
9.                 WB(x,T0,lamb).*Bnx(x,jj,T0,lamb,Gamma,pjun);
   
10.             H3jmn(m,n,jj)=double(quadl(fH3,0,1));
    
11.         end
    
12.     end
    
13. end
    

复制代码

就是T0,lamb,Gamma,pjun,N是参数，然后计算H3jmn是个三维数组，就是一定要（m,n,j)位置上分别对应值，怎么才能把这个for的循环去掉呢？:L

messenger

参考此贴，http://forum.simwe.com/thread-792777-1-1.html

findpass

谢谢messenger