📄 pec.m
字号:
% 该函数确定导体边界条件
function [Exc,Ezc]=PEC(Nx,Ny,Nz,dx,dy,dz,W,S,TT,NH)
Exc(1:Nx,1:Ny,1:Nz)=1; Ezc(1:Nx,1:Ny,1:Nz)=1;
Nxsta=round((Nx-(2*W+S)/dx)/2); Nxend=round(Nx/2+(2*W+S)/(2*dx));
Exc(Nxsta:Nxsta+W/dx,NH:NH+TT/dy,1:Nz)=0;
Exc(Nxend-W/dx:Nxend,NH:NH+TT/dy,1:Nz)=0;
Ezc(Nxsta:Nxsta+W/dx,NH:NH+TT/dy,1:Nz)=0;
Ezc(Nxend-W/dx:Nxend,NH:NH+TT/dy,1:Nz)=0;
Exc(:,1,:)=0; %接地板
Ezc(:,1,:)=0;
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -