有限差分法电磁场matlab

有限差分法（FDTD）是一种求解电磁场问题的数值方法，可以使用MATLAB进行实现。以下是一个简单的二维TM波的FDTD代码示例：

% 定义常数
c = 3e8; % 光速
dx = 1e-3; % 空间步长
dt = dx/c; % 时间步长
t_end = 1e-8; % 模拟时间
x_num = 100; % x方向网格数
y_num = 100; % y方向网格数

% 初始化场
Ez = zeros(x_num, y_num); % z方向电场
Hy = zeros(x_num, y_num); % y方向磁场

% 模拟主循环
for t = 0:dt:t_end
    % 更新Hy
    for i = 1:x_num-1
        for j = 1:y_num
            Hy(i,j) = Hy(i,j) + (Ez(i,j) - Ez(i+1,j)) * dt/dx;
        end
    end
    
    % 更新Ez
    for i = 2:x_num-1
        for j = 2:y_num-1
            Ez(i,j) = Ez(i,j) + (Hy(i,j) - Hy(i-1,j)) * dt/dx;
        end
    end
    
    % 边界处理（吸收边界）
    Ez(1,:) = 0;
    Ez(x_num,:) = 0;
    Ez(:,1) = 0;
    Ez(:,y_num) = 0;
end

% 可视化结果
imagesc(Ez);