putrectangle.m

来自「外国人开发的电磁时域有限差分方法工具包 Electromagnetic Fin」· M 代码 · 共 25 行

function G=PutRectangle(F,x1,y1,x2,y2,Hole)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Electromagnetic Finite-Difference Time-Domain %
% Version 1.20, Release 1                       %
%                                               %
%   (C) Copyright 2005                          %
%   Sharif University of Technology             %
%   School of Electrical Engineering            %
%   All Rights Reserved                         %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

global xCnt yCnt

s=size(F);
sx=s(1);
sy=s(2);
lx = max(1,round(x1));
ux = min(min(xCnt,sx),round(x2));
ly = max(1,round(y1));
uy = min(min(yCnt,sy),round(y2));

F(lx:ux,ly:uy)=Hole;

G=F;