Question

使用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

Answer 1

楼主说得不是很清楚，不知道是要求前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))

追问

你好，谢谢你的解答！就是对M（i，j，g)矩阵的其中i=21,j=30（即M（21，30，g)）元素进行就是求和。我已经把代码附上了，你看下怎么修改？

追答

原来楼主的东西很是很高深的，我以前在做滤波器时，用过椭圆函数
但是很久的事了，当然不说理论，只说代码，我想我前面的回答已经说了
其实，我也不是很有把握
3维矩阵，就是在一根绳子上挂着一串2维矩阵，你的第3维是20吧
每一片是是21*30吧，你是要求1:20的时候，每一片矩阵的所有元素的和，对吗？
最后的结果是一个1*20的向量？稍微改一下：
clear all;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:length(d)
for i=1:length(a)
for j=1:length(b)
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)=sum(sum(M(:,:,g))); %这句改了一下
K12(g)=Mab(g)/((L1*L2)^0.5);
end

追问

也不是很高深啦，是在用完全椭圆积分求解某一耦合系数，虽然你这一改没有报错，但求解的数据是不对的。以我对单个g，如另g=0.01,然后求sum(M(:))的结果不符合的。应该是求和求错了。

要求1:20的时候，每一片矩阵的所有元素的和，是这样的。

追答

我想，你要看看你前面的程序，求得对不对，这是第一
第二，你到底要求解什么，反正每一个2维矩阵都是21*30的
如果你要求这630个元素的和值，就是这个结果
sum是按列求解的，所以需要2个sum
除非你要求解的不是每个630个元素的和值？