传 奇 脚 本 问 题 ， 高 分 悬 赏 。_百度知道小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
高分悬赏一个matlab并行的问题
load 8x8.mat
D=x8S2*27.2;
KBT=0.025;
DPCx2=DPCx*DPCx;
DPCy2=DPCy*DPCy;
a0=17.435242;
b0=30.198726;
gamma=2.7*1e-4;
gamma2=gamma*
x=(-((k-1)/(2*k))):1/k:((k-1)/(2*k));
xi=-0.5/a0:(1/(k*c))/a0:0.5/a0;
=meshgrid(xi,yi);%生成网格采样点
for i=1:1:660;
& & for j=1:1:k
& && &&&for l=1:1:k
& && && &&&Ei(j,l)=D(i,k*j+l-k);& && && && && && && && &
& && &&&end& && && && && && &&&
& & end& &
& & E(:,:,i)=Ei;& &
& & zi=interp2(x,y,Ei,xi,yi,'spline');%插值
& & Z(:,:,i)=
& & =gradient(zi);
& & VX(:,:,i)=WX;
& & VY(:,:,i)=WY;
& & VX(find(isnan(x)==1))=0;%找到NaN的赋值为0
& & VY(find(isnan(x)==1))=0;%找到NaN的赋值为0
& & for a=1:1:k*c+1
& && &&&for b=1:1:k*c+1
& && && && &%b
& && && && &tx=0;
& && && && &ty=0;
& && && && &for m=1:1:k*c+1
& && && && && & for n=1:1:k*c+1
& && && && && && &&&vx=VX(m,n,i);
& && && && && && &&&vy=VY(m,n,i);
& && && && && && &&&vx0=VX(a,b,i);
& && && && && && &&&vy0=VY(a,b,i);
& && && && && && &&&d=gamma/((zi(m,n)-zi(a,b))^2+gamma2);
& && && && && && &&&tx=tx+((1-vx/vx0)*d);& && && && && && && &
& && && && && && &&&ty=ty+((1-vy/vy0)*d);& && && && && && && &
& && && && && & end
& && && && &end
& && && && &Tx(a,b,i)=1/((KBT*DPCx2*2*tx)/(hba*Cx));
& && && && &Ty(a,b,i)=1/((KBT*DPCy2*2*ty)/(hba*Cy));
& && &&&end
& & Taox=Tx(:,:,i);
& & Taoy=Ty(:,:,i);
& & save T
EF=-0.062919;
for i=1:1:660& &
& & fenzixi=0;
& & fenziyi=0;
& & fenmui=0;
& & for j=1:1:(c*k+1)& && &
& && & fenzixji=0;
& && & fenziyji=0;
& && & fenmuji=0;
& && & for l=1:1:(c*k+1)& && && &
& && && & f0(l,j,i)=1/(1+exp(Z(l,j,i)-EF)/KBT);
& && && & fenzixjli=Tx(l,j,i)*VX(l,j,i)*VX(l,j,i)*f0(l,j,i);
& && && & fenziyjli=Ty(l,j,i)*VY(l,j,i)*VY(l,j,i)*f0(l,j,i);
& && && & fenzixji=fenzixji+
& && && & fenziyji=fenziyji+
& && && & fenmuji=fenmuji+f0(l,j,i);
& && && &&&
& && & end
& && & fenzixi=fenzixi+
& && & fenziyi=fenziyi+
& && & fenmui=fenmui+
& & end& &
& & fenzix=fenzix+
& & fenziy=fenziy+
& & fenmu=fenmu+
muix=40*fenzix/
muiy=40*fenziy/
以上是我整个程序的编程，因为数据量是80x80x660，算下来要几天时间，想用实验室的八核计算机并行计算看能不能快一点，求高人指点该怎么改程序和如何设置matlab的参数。
Taox那个是三维数组，i指的就是那个层数。每次就把结果放到一层里面的。
那样的parfor也是过，不过出错显示的信息是里面的Tx无法classified。。
Error: The variable Ei in a parfor cannot be classified.
See Parallel for Loops in MATLAB, &Overview&.
连原文件都给你吧
额，这教研室的网，小木虫的东西都下不动。。。QQ
已经自己解决了。谢谢！
坑爹啊。。。。我还特意让别人下给我了。改进了一晚上。。现在算完大概50分钟，应该可以改进到20分钟以内。
请问楼主现在的计算速度是多少？:tuzi4:诶，有种被骗的赶脚。。。
121*121*660的数据全部算完才50分钟？？？我的改完用parfor算还要大概一天如果八个核算的话。。如果您那个真的是一个小时内的我愿意500金币双手奉上。
。。。八核，一天。。。这 效率太低了
第一次用matlab算东西。新手请见谅。。
没有啊。。你要是一层能算那么快的话我真的很需要啊，现在时间很紧啊。。这个速度能提上来对我很要命的。。为什么你的机器才用了37秒？我这边一层都要几分钟啊。。。你的是什么机器来的。。。
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
登录小木虫
打开微信扫一扫
随时随地聊科研