0

所以我正在写这篇论文。我目前正在尝试对等式 20 之后给出的矩阵进行编码。

我有以下代码：

cx = 1/(1+aprime^2);
cp = 1/period * int(cx*exp(-i*m*K*x),x,0,period);

d_x = aprime/(1+aprime^2);
dp = 1/period * int(d_x*exp(-i*m*K*x),x,0,period);
alpham = (k * sinteta) + m*K;
betam = sqrt((2*pi/lbd)^2 - subs(alpham,m)^2);
beta0 = subs(betam,0);
alpha0 = subs(alpham,0);


for p = -Mmax:Mmax;
    for r = -Mmax:Mmax;

        A((p+Mmax)+1,(r+Mmax)+1,j) = subs(alpham,r)*subs(dp,m,(p-r));

        B(p+Mmax+1,r+Mmax+1,j) = subs(cp,m,p-r);

        C(p+Mmax+1,r+Mmax+1,j) = -subs(alpham,p)*subs(alpham,r)*subs(cp,m,p-r)  +  k^2 * n(j)^2 *kronDel(p,r);

        D(p+Mmax+1,r+Mmax+1,j) = subs(alpham,p)*subs(dp,m,p-r);

    end
end
T(:,:,j) = [A(:,:,j) , B(:,:,j)
    C(:,:,j) , D(:,:,j)];

j 是微不足道的，因为我必须为 3 种不同的介质做它。但我想知道是否有任何明显的错误，因为要么i）它需要永远基于'aprime'或ii）我最终在矩阵中得到NaN，这不允许我采取eig

目前aprime是diff(abs(x-.5),x)

谢谢你的帮助！

4

0 回答 0