side1_aa.m

来自「用fdtd方法 模拟集成光学离子交换中扩散离子浓度的分布 进而求得折射率的分布。」· M 代码 · 共 9 行

function sd1p=side1_point(mat,i,j,dt,d,db,alf,h)
sd1p=dt.*(...
    d(i,j)./(h^2.*(1-alf(i,j).*mat(i,j))).*(mat(i+1,j)+mat(i,j+1)+mat(i,j-1)-3*mat(i,j)+alf(i,j)./(1-alf(i,j).*mat(i,j)).*((mat(i+1,j)-mat(i,j)).^2+(mat(i,j+1)-mat(i,j)).^2))...
    +...
    (d(i+1,j)-d(i,j))./h.*(1-mat(i,j))./(1-alf(i,j).*mat(i,j)).*((mat(i+1,j)-mat(i,j))./h.*(1+alf(i,j)./(1-alf(i,j).*mat(i,j))))...
    +...
    (db(i+1,j)-db(i,j))./h.*mat(i,j).*(1-alf(i,j))./(1-alf(i,j).*mat(i,j)).*((mat(i+1,j)-mat(i,j))./h.*(1-alf(i,j).*(1-mat(i,j))./(1-alf(i,j).*mat(i,j))))...
    );