[求助]求助各位高手MATLAB运行速度的问题

求助各位高手MATLAB运行速度的问题 Vt \g9-[  
7p\&D?  
我编写的一个程序运行的网格有80*180，运行的时间步数为1000，按道理我的程序运行完成应该不用多久的，可是却用了270秒钟！ g"Hl 30o  
>);M\,1\I  
同样的都是fdtd的程序我运行别人的一个程序网格500*200，时间步数也有4000。可是这个程序仅仅用了360秒！为什么差距会那么大？我的 .. -dF (_ %C  
A@?-
"=h}  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

Matlab循环的效率很低

对，尽量不要用循环，能用数组的用数组。

for p=Imin:Imax                           %hx hy 的fdtd迭代 ~Al3Dv9x  
     for j=Jmin:Jmax-1 5A5t
 
         T1=ez(p,j+1)-ez(p,j); "+`u ]  
         hx(p,j)=cay(ob(p,j),ob(p,j+1))*hx(p,j)-cby(ob(p,j),ob(p,j+1))/377*T1; Btr>ek  
     end lfd-!(tXD  
end Jy"\_Vvl  

<tsexsw  
T1=ez(Imin:Imax,Jmin+1:Jmax)-ez(Imin:Imax,Jmin:Jmax-1); i|,}y`C#  
vr?u=_%Z  
hx(Imin:Imax,Jmin:Jmax-1)=cay(ob(Imin:Imax,Jmin:Jmax-1),ob(Imin:Imax,Jmin+1:Jmax)).*hx(Imin:Imax,Jmin:Jmax-1)-... jLO$[c`;  
    -cby(ob(Imin:Imax,Jmin:Jmax-1),ob(Imin:Imax,Jmin+1:Jmax))/377.*T1; q(R|3l^6T  
pV\>?  
像这种就是矩阵运算吗？可是这样改的时候size（cay(ob(Imin:Imax,Jmin:Jmax-1),ob(Imin:Imax,Jmin+1:Jmax))）就会超过10000*10000 /[/{m]  
MATLAB就会提示out of memory了，这该怎么办呢？

for p=Imin:Imax                           %hx hy 的fdtd迭代 'Ox "YE  
     for j=Jmin:Jmax-1 ZFH-srs{  
         T1=ez(p,j+1)-ez(p,j); >o`+j$j  
         hx(p,j)=cay(ob(p,j),ob(p,j+1))*hx(p,j)-cby(ob(p,j),ob(p,j+1))/377*T1; agruS'c g  
     end #4"eQ*.*"  
end eLgq )  
Uugq.'>  
T1=ez(Imin:Imax,Jmin+1:Jmax)-ez(Imin:Imax,Jmin:Jmax-1); 31#jLWY'0  
UQYHR+  
hx(Imin:Imax,Jmin:Jmax-1)=cay(ob(Imin:Imax,Jmin:Jmax-1),ob(Imin:Imax,Jmin+1:Jmax)).*hx(Imin:Imax,Jmin:Jmax-1)-... 1gt 7My  
    -cby(ob(Imin:Imax,Jmin:Jmax-1),ob(Imin:Imax,Jmin+1:Jmax))/377.*T1; X&rsWk  
\UM&|yk:  
像这种就是矩阵运算吗？可是这样改的时候size（cay(ob(Imin:Imax,Jmin:Jmax-1),ob(Imin:Imax,Jmin+1:Jmax))）就会超过10000*10000 R-,L"Vv  
MATLAB就会提示out of memory了，这该怎么办呢？ |)0Ta9~  

for p=Imin:Imax                           %hx hy 的fdtd迭代 7i02M~*uS  
     for j=Jmin:Jmax-1 Y,,Z47% E  
         T1=ez(p,j+1)-ez(p,j); 8WpNlB+:{  
         hx(p,j)=cay(ob(p,j),ob(p,j+1))*hx(p,j)-cby(ob(p,j),ob(p,j+1))/377*T1; ''!j:49  
     end hBaG*J{  
end OUO'w6m!  
lg ,%  
T1=ez(Imin:Imax,Jmin+1:Jmax)-ez(Imin:Imax,Jmin:Jmax-1); "s_S!;w@  
vgg)f~  
hx(Imin:Imax,Jmin:Jmax-1)=cay(ob(Imin:Imax,Jmin:Jmax-1),ob(Imin:Imax,Jmin+1:Jmax)).*hx(Imin:Imax,Jmin:Jmax-1)-... Vj8-[ww!  
    -cby(ob(Imin:Imax,Jmin:Jmax-1),ob(Imin:Imax,Jmin+1:Jmax))/377.*T1; Vu4LC&q  
U#w0E G  
像这种就是矩阵运算吗？可是这样改的时候size（cay(ob(Imin:Imax,Jmin:Jmax-1),ob(Imin:Imax,Jmin+1:Jmax))）就会超过10000*10000 Buq(L6P9r  
MATLAB就会提示out of memory了，这该怎么办呢？

这个，我还没遇到这种情况。是不是你的计算机的配置不够好啊，内存多大？matlab应该不会限制数组的大小把

不会啊我查看的memory命令明明还有内存可用的，我的内存是2G的足够了才是！ T?jN/}qg  

而且很奇怪的是我定义的cay数组才是2*2的可是size（cay(ob(Imin:Imax,Jmin:Jmax-1),ob(Imin:Imax,Jmin+1:Jmax))）却是变得那么大了呢？