matlab

matlab

matlab 索引超出矩阵维度

现在在学习matlab复现公式，出现了一些问题，想请大神帮忙看一看。

function time = producedelay(Sortmatix1,aerfa_1_1,F,f)
global B W W1 N0 d M OMEGA task omega Sortmatix0 power1 power2

%-------------------统计每个用户占用了几个信道--------------------------%
%--------------对应子载波的分配存在了Sortmatix1中---------------------%
%--------------对应子载波的功率存在了Sortmatix2中---------------------%
d1=Sortmatix0.Sortmatix1;
Sortmatix2=d10;
bb=zeros(1,12);
for m=1:8
midle_a=[d1(1,m) d1(2,m) d1(3,m) d1(4,m) d1(1,m+8) d1(2,m+8) d1(3,m+8) d1(4,m+8) d1(1,m+16) d1(2,m+16) d1(3,m+16) d1(4,m+16)];
bb=waterfillingup(midle_a,power1);
for n=1:4
nn=mod(n,4);
if (nn0) nn=4;
end
Sortmatix2(nn,m)=bb(1,n);
end
for n=5:8
nn=mod(n,4);if (nn0) nn=4; end
Sortmatix2(nn,m+8)=bb(1,n);
end
for n=9:12
nn=mod(n,4);if (nn0) nn=4; end
Sortmatix2(nn,m+16)=bb(1,n);
end
end%上行功率分配
%--------------------------
bb=zeros(1,32);
for m=1:3 %%%下行m是基站数，上行不是
mm=(m-1)*8;
midle_a=[d1(5,mm+1) d1(5,mm+2) d1(5,mm+3) d1(5,mm+4) d1(5,mm+5) d1(5,mm+6) d1(5,mm+7) d1(5,mm+8) d1(6,mm+1) d1(6,mm+2) d1(6,mm+3) d1(6,mm+4) d1(6,mm+5) d1(6,mm+6) d1(6,mm+7) d1(6,mm+8) d1(7,mm+1) d1(7,mm+2) d1(7,mm+3) d1(7,mm+4) d1(7,mm+5) d1(7,mm+6