geoplot.m

非连续波导时间域有限差分法仿真（TD-FDM）

% geometry plot

%figure;

Nx = L/hx+1;
Ny = H/hy+1;
Nix = l/hx-1;
Niy = round((H-2*h)/hy+1);
Nh = h/hy+1;
node = 2*Nx*Ny + Nix*Niy;

MAX_X = 2*Nx+Nix;
MAX_Y = Ny;

X = [xo; L+xo; L+xo; L+xo+l; L+xo+l; 2*L+xo+l];
Y = [yo; yo; yo+h; yo+h; yo; yo];
%line(X,Y);
%hold on;

X = [xo; L+xo; L+xo; L+xo+l; L+xo+l; 2*L+xo+l];
Y = [yo+H; yo+H; yo+H-h; yo+H-h; yo+H; yo+H];
%line(X,Y);
%hold on;

%axis([0 28 -1 6]);

% plot node schema (spacing discretization)
xval = zeros(node,1);
yval = zeros(node,1);
for ii=1:Nx
    for jj=1:Ny
        xval(jj+(ii-1)*Ny) = xo+(ii-1)*hx;
        yval(jj+(ii-1)*Ny) = yo+(jj-1)*hy;
    end
end
for ii=1:Nix
    for jj=1:Niy
        xval(Nx*Ny+jj+(ii-1)*Niy) = xo+hx*Nx+(ii-1)*hx;
        yval(Nx*Ny+jj+(ii-1)*Niy) = yo+hy*(Ny-Niy)/2+(jj-1)*hy;
    end
end
for ii=1:Nx
    for jj=1:Ny
        xval(Nx*Ny+Nix*Niy+jj+(ii-1)*Ny) = xo+Nx*hx+Nix*hx+(ii-1)*hx;
        yval(Nx*Ny+Nix*Niy+jj+(ii-1)*Ny) = yo+(jj-1)*hy;
    end
end
%plot(xval,yval,'.');
%axis([-2e-2 27e-2 -1e-2 4e-2]);