matlab并行查询,高分悬赏一个matlab并行的问题 - 仿真模拟 - 小木虫 - 学术 科研 互动社区...

CODE:

load 8x8.mat

D=x8S2*27.2;

k=8;

c=10;

KBT=0.025;

hba=1;

DPCx=1;

DPCx2=DPCx*DPCx;

DPCy=1;

DPCy2=DPCy*DPCy;

Cx=1;

Cy=1;

a0=17.435242;

b0=30.198726;

gamma=2.7*1e-4;

gamma2=gamma*gamma;

x=(-((k-1)/(2*k))):1/k(k-1)/(2*k));

y=[((k-1)/(2*k)):-1/k-((k-1)/(2*k)))]';

x=x/a0;

y=y/b0;

xi=-0.5/a01/(k*c))/a0:0.5/a0;

yi=[0.5/b0-1/(k*c))/b0:-0.5/b0]';

[xi,yi]=meshgrid(xi,yi);%生成网格采样点

for i=1:1:660;

i

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)=zi;

[WX,WY]=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

a

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

end

i

Taox=Tx(:,:,i);

Taoy=Ty(:,:,i);

save T

end

EF=-0.062919;

fenzix=0;

fenziy=0;

fenmu=0;

for i=1:1:660

fenzixi=0;

fenziyi=0;

fenmui=0;

for j=1:1c*k+1)

fenzixji=0;

fenziyji=0;

fenmuji=0;

for l=1:1c*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+fenzixjli;

fenziyji=fenziyji+fenziyjli;

fenmuji=fenmuji+f0(l,j,i);

end

fenzixi=fenzixi+fenzixji;

fenziyi=fenziyi+fenziyji;

fenmui=fenmui+fenmuji;

fenmuji;

i

j

end

fenzix=fenzix+fenzixi;

fenziy=fenziy+fenziyi;

fenmu=fenmu+fenmui;

end

muix=40*fenzix/fenmu;

muiy=40*fenziy/fenmu;

save MUI