matlab - これを matlab でコーディングするより効率的な方法

Question

だから私はこの論文に取り組んでいます。私は現在、方程式 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)];

私は3つの異なる媒体のためにそれをしなければならないので、jは重要ではありません. しかし、i)「aprime」に基づいて永遠にかかるか、ii)行列に NaN が含まれてしまい、取得できないため、明らかな間違いがあるかどうか疑問に思っていましたeig

現在 aprime はdiff(abs(x-.5),x)

助けてくれてありがとう！