mtseasky
2008-07-25, 23:47
...
a=-5;
jj=5;
N0=1e6;
X=1:0.1:15; %141 ab values
Y=0.3:0.01:0.9; %61 SF2 values
SF2=0.3; %initialize
for i=1:61
ab=1; %initialize
for j=1:141
phan=EUD(enclosed_voxels{jj}.index, phan, ab); % 转换3D矩阵phan, EUD函数中没有循环操作
gEUD=0;
tot=length(enclosed_voxels{jj}.index);
for k=1:tot
gEUD=gEUD+phan(enclosed_voxels{jj}.index(k)).^a;
end
gEUD=(gEUD/tot)^(1/a);
Z(i,j)=exp(-N0*(SF2)^(gEUD/2));
ab=ab+0.1;
end
SF2=SF2+0.01;
end
. . . .
我算上面的程序好慢啊:ft:
请问如何减少一些循环操作
多谢
a=-5;
jj=5;
N0=1e6;
X=1:0.1:15; %141 ab values
Y=0.3:0.01:0.9; %61 SF2 values
SF2=0.3; %initialize
for i=1:61
ab=1; %initialize
for j=1:141
phan=EUD(enclosed_voxels{jj}.index, phan, ab); % 转换3D矩阵phan, EUD函数中没有循环操作
gEUD=0;
tot=length(enclosed_voxels{jj}.index);
for k=1:tot
gEUD=gEUD+phan(enclosed_voxels{jj}.index(k)).^a;
end
gEUD=(gEUD/tot)^(1/a);
Z(i,j)=exp(-N0*(SF2)^(gEUD/2));
ab=ab+0.1;
end
SF2=SF2+0.01;
end
. . . .
我算上面的程序好慢啊:ft:
请问如何减少一些循环操作
多谢