使用matlab 对矩阵M(i,j,g)的前面两维进行求和，结果表示为Md(g),用Md(g)=sum(M(:,:,g))不对，如何修改？

clc;
format long
L1=3.151*10^(-6);
L2=2.660*10^(-6);
u0=4*pi*10^(-7);
a1=2.525:0.55:13.525;
b1=0.5075:0.018:1.0295;
a=0.001.*a1 ;
b=0.001.*b1 ;
d=0.001:0.001:0.02;
for g=1:1:20
for i=1:1:21
for j=1:1:30
k(i,j,g)=(4*a(i)*b(j)/((a(i)+b(j))^2+d(g)^2))^0.5;
[K(i,j,g),E(i,j,g)]=ellipke(k(i,j,g));
M(i,j,g)=u0*((a(i)*b(j))^0.5)*((2/k(i,j,g)-k(i,j,g))*K(i,j,g)-2/k(i,j,g)*E(i,j,g));
end
end
Mab(g)=1.1*sum(M(:,:,g)) ;/这步错了
K12(g)=Mab(g)/((L1*L2)^0.5);
end

楼主说得不是很清楚，不知道是要求前2维什么形式的和呢？
是“每一片”矩阵的元素的和？还是？希望说清楚些，
不知是不是你说的意思：
clear all;clc;
a=[1,2;3,4];
b=ones(2,2);
c=magic(2);
d=cat(3,a,b,c);
Md=sum(sum(d))

同问。。。