有没有学过微波电磁波方面会c++编程的大佬，这边有个fdtd相关请教一下？

Question

1关于电磁波入射的fdtd问题，想用这个公式改一下周期边界条件我应该怎么去改

void CFDTD::BoudaryExchangeE(){
    register int i,j,k;
    //PEC
    for(k=0;k<Nz+1;k++){
        for(j=1;j<=Ny/2+1;j++){
            Ey[k][j][1] = 0.0;
            Ey[k][j][Nx/2+1] = 0.0;
            Ez[k][j][1] = 0.0;
            Ez[k][j][Nx/2+1] = 0.0;
        }
    }
    //PEC
    for(k=0;k<Nzi+1;k++){
        for(j=1;j<=Ny/2+1;j++){
            Ey_inc[k][j][1] = 0.0;
            Ey_inc[k][j][Nx/2+1] = 0.0;
            Ez_inc[k][j][1] = 0.0;
            Ez_inc[k][j][Nx/2+1] = 0.0;
        }
    }

}

用图中的公式展开成sin cos形式改写

Answer 1

https://blog.csdn.net/qq_43946615/article/month/2020/05