object_drawing_functions.m

this a 3d fdtd code using matlab and it can be used as a reference for other structures. Dear sir,

function object_drawing_functions(obj)

if strcmp(obj.obj_type, 'voltage_source')
    draw_a_voltage_source(obj);
end

if strcmp(obj.obj_type, 'current_source')
    draw_a_current_source(obj);
end

if strcmp(obj.obj_type, 'resistor')
    draw_a_resistor(obj);
end

if strcmp(obj.obj_type, 'capacitor')
    draw_a_capacitor(obj);
end

if strcmp(obj.obj_type, 'inductor')
    draw_an_inductor(obj);
end

if strcmp(obj.obj_type, 'diode')
    draw_a_diode(obj);
end

if strcmp(obj.obj_type, 'sampled_voltage')
    draw_a_sampled_voltage(obj);
end

if strcmp(obj.obj_type, 'sampled_current')
    draw_a_sampled_current(obj);
end

if strcmp(obj.obj_type, 'thin_wire')
    draw_a_thin_wire(obj);
end

if strcmp(obj.obj_type, 'sampled_electric_field')
    draw_a_sampled_electric_field(obj);
end

if strcmp(obj.obj_type, 'sampled_magnetic_field')
    draw_a_sampled_magnetic_field(obj);
end

if strcmp(obj.obj_type, 'axis_outside_domain')
    draw_axis_outside_domain(obj);
end

if strcmp(obj.obj_type, 'axis_at_origin')
    draw_axis_at_origin(obj);
end

%===================================================================
function draw_a_current_source(obj)

ix = obj.min_x; iy = obj.min_y; iz = obj.min_z; 
ax = obj.max_x; ay = obj.max_y; az = obj.max_z;

switch (obj.direction)
case ('xn')
    vx = [ax ix]; vy = 0.5*[(iy+ay) (iy+ay)]; vz = 0.5*[(iz+az) (iz+az)];
    len = ax - ix;
    ir = len/20;
    [v1, f1] = a_line([ix ix+0.2*len], vy, vz, ir);
    [v2, f2] = a_line([ix+0.8*len ax], vy, vz, ir); 
    if ((az-iz)<(ay-iy)) dn = 'z'; else dn = 'y'; end
    cx = (vx(1)+vx(2))/2;
    cy = (vy(1)+vy(2))/2;
    cz = (vz(1)+vz(2))/2;
    [v3, f3] = a_circle(cx, cy, cz, len*0.3, dn, ir);    
    [v4, f4] = a_vector([cx+len/10 cx-len/10], [cy cy], [cz cz], len/16, len/8, 0.5);
case ('xp')
    vx = [ix ax]; vy = 0.5*[(iy+ay) (iy+ay)]; vz = 0.5*[(iz+az) (iz+az)];
    len = ax - ix;
    ir = len/20;
    [v1, f1] = a_line([ix ix+0.2*len], vy, vz, ir);
    [v2, f2] = a_line([ix+0.8*len ax], vy, vz, ir); 
    if ((az-iz)<(ay-iy)) dn = 'z'; else dn = 'y'; end
    cx = (vx(1)+vx(2))/2;
    cy = (vy(1)+vy(2))/2;
    cz = (vz(1)+vz(2))/2;
    [v3, f3] = a_circle(cx, cy, cz, len*0.3, dn, ir);    
    [v4, f4] = a_vector([cx-len/10 cx+len/10], [cy cy], [cz cz], len/16, len/8, 0.5);
case ('yn')
    vx = 0.5*[(ix+ax) (ix+ax)]; vy = [ay iy]; vz = 0.5*[(iz+az) (iz+az)];
    len = ay - iy;
    ir = len/20;
    [v1, f1] = a_line(vx, [iy iy+0.2*len], vz, ir);
    [v2, f2] = a_line(vx, [iy+0.8*len ay], vz, ir); 
    if ((ax-ix)<(az-iz)) dn = 'x'; else dn = 'z'; end
    cx = (vx(1)+vx(2))/2;
    cy = (vy(1)+vy(2))/2;
    cz = (vz(1)+vz(2))/2;
    [v3, f3] = a_circle(cx, cy, cz, len*0.3, dn, ir);    
    [v4, f4] = a_vector([cx cx],[cy+len/10 cy-len/10], [cz cz], len/16, len/8, 0.5);
case ('yp')
    vx = 0.5*[(ix+ax) (ix+ax)]; vy = [iy ay]; vz = 0.5*[(iz+az) (iz+az)];
    len = ay - iy;
    ir = len/20;
    [v1, f1] = a_line(vx, [iy iy+0.2*len], vz, ir);
    [v2, f2] = a_line(vx, [iy+0.8*len ay], vz, ir); 
    if ((ax-ix)<(az-iz)) dn = 'x'; else dn = 'z'; end
    cx = (vx(1)+vx(2))/2;
    cy = (vy(1)+vy(2))/2;
    cz = (vz(1)+vz(2))/2;
    [v3, f3] = a_circle(cx, cy, cz, len*0.3, dn, ir);    
    [v4, f4] = a_vector([cx cx],[cy-len/10 cy+len/10], [cz cz], len/16, len/8, 0.5);
case ('zn')
    vx = 0.5*[(ix+ax) (ix+ax)]; vy = 0.5*[(iy+ay) (iy+ay)]; vz = [az iz];
    len = az - iz;
    ir = len/20;
    [v1, f1] = a_line(vx, vy, [iz iz+0.2*len], ir);
    [v2, f2] = a_line(vx, vy, [iz+0.8*len az], ir); 
    if ((ax-ix)<(ay-iy)) dn = 'x'; else dn = 'y'; end
    cx = (vx(1)+vx(2))/2;
    cy = (vy(1)+vy(2))/2;
    cz = (vz(1)+vz(2))/2;
    [v3, f3] = a_circle(cx, cy, cz, len*0.3, dn, ir);    
    [v4, f4] = a_vector([cx cx], [cy cy], [cz+len/10 cz-len/10], len/16, len/8, 0.5);
case ('zp')
    vx = 0.5*[(ix+ax) (ix+ax)]; vy = 0.5*[(iy+ay) (iy+ay)]; vz = [iz az];
    len = az - iz;
    ir = len/20;
    [v1, f1] = a_line(vx, vy, [iz iz+0.2*len], ir);
    [v2, f2] = a_line(vx, vy, [iz+0.8*len az], ir); 
    if ((ax-ix)<(ay-iy)) dn = 'x'; else dn = 'y'; end
    cx = (vx(1)+vx(2))/2;
    cy = (vy(1)+vy(2))/2;
    cz = (vz(1)+vz(2))/2;
    [v3, f3] = a_circle(cx, cy, cz, len*0.3, dn, ir);    
    [v4, f4] = a_vector([cx cx], [cy cy], [cz-len/10 cz+len/10], len/16, len/8, 0.5);
end
v = [v1; v2; v3; v4];
f = f1;
f = [f; max(max(f)) + f2];
f = [f; max(max(f)) + f3];
f = [f; max(max(f)) + f4];
patch('Vertices',v,'Faces',f,'facecolor','m', ...
    'facealpha',1,'edgecolor','none');

% draw boundaries
v = [0 0 0; 1 0 0; 1 1 0; 0 1 0; 0 0 1; 1 0 1; 1 1 1; 0 1 1];
vx = v(:,1) * (ax - ix) + ix;
vy = v(:,2) * (ay - iy) + iy;
vz = v(:,3) * (az - iz) + iz;
f = [1 2 6 5; 2 3 7 6; 3 4 8 7; 4 1 5 8; 1 2 3 4; 5 6 7 8];
patch('Vertices',[vx vy vz],'Faces',f,'facecolor','none','LineWidth',1, ...
    'facealpha',0.3,'edgecolor','m','linestyle','-');

% show label
text(ix,iy,iz,obj.obj_id,'VerticalAlignment','bottom', ...
    'BackgroundColor','w','edgecolor','m','fontsize',12);

%===================================================================
function draw_a_voltage_source(obj)

ix = obj.min_x; iy = obj.min_y; iz = obj.min_z; 
ax = obj.max_x; ay = obj.max_y; az = obj.max_z;

switch (obj.direction)
case ('xn')
    vx = [ax ix]; vy = 0.5*[(iy+ay) (iy+ay)]; vz = 0.5*[(iz+az) (iz+az)];
    len = ax - ix;
    ir = len/30;
    [v1, f1] = a_line([ix ix+0.2*len], vy, vz, ir);
    [v2, f2] = a_line([ix+0.8*len ax], vy, vz, ir); 
    if ((az-iz)<(ay-iy)) dn = 'z'; else dn = 'y'; end
    cx = (vx(1)+vx(2))/2;
    cy = (vy(1)+vy(2))/2;
    cz = (vz(1)+vz(2))/2;
    [v3, f3] = a_circle(cx, cy, cz, len*0.3, dn, ir);    
    [v4, f4] = a_plus(cx-len/10, cy, cz, len/5, dn, ir);
    [v5, f5] = a_minus(cx+len/10, cy, cz, len/5, dn, ir, 'x');
case ('xp')
    vx = [ix ax]; vy = 0.5*[(iy+ay) (iy+ay)]; vz = 0.5*[(iz+az) (iz+az)];
    len = ax - ix;
    ir = len/30;
    [v1, f1] = a_line([ix ix+0.2*len], vy, vz, ir);
    [v2, f2] = a_line([ix+0.8*len ax], vy, vz, ir); 
    if ((az-iz)<(ay-iy)) dn = 'z'; else dn = 'y'; end
    cx = (vx(1)+vx(2))/2;
    cy = (vy(1)+vy(2))/2;
    cz = (vz(1)+vz(2))/2;
    [v3, f3] = a_circle(cx, cy, cz, len*0.3, dn, ir);    
    [v4, f4] = a_plus(cx+len/10, cy, cz, len/5, dn, ir);
    [v5, f5] = a_minus(cx-len/10, cy, cz, len/5, dn, ir, 'x');
case ('yn')
    vx = 0.5*[(ix+ax) (ix+ax)]; vy = [ay iy]; vz = 0.5*[(iz+az) (iz+az)];
    len = ay - iy;
    ir = len/30;
    [v1, f1] = a_line(vx, [iy iy+0.2*len], vz, ir);
    [v2, f2] = a_line(vx, [iy+0.8*len ay], vz, ir); 
    if ((ax-ix)<(az-iz)) dn = 'x'; else dn = 'z'; end
    cx = (vx(1)+vx(2))/2;
    cy = (vy(1)+vy(2))/2;
    cz = (vz(1)+vz(2))/2;
    [v3, f3] = a_circle(cx, cy, cz, len*0.3, dn, ir);    
    [v4, f4] = a_plus(cx, cy-len/10, cz, len/5, dn, ir);
    [v5, f5] = a_minus(cx, cy+len/10, cz, len/5, dn, ir, 'y');
case ('yp')
    vx = 0.5*[(ix+ax) (ix+ax)]; vy = [iy ay]; vz = 0.5*[(iz+az) (iz+az)];
    len = ay - iy;
    ir = len/30;
    [v1, f1] = a_line(vx, [iy iy+0.2*len], vz, ir);
    [v2, f2] = a_line(vx, [iy+0.8*len ay], vz, ir); 
    if ((ax-ix)<(az-iz)) dn = 'x'; else dn = 'z'; end
    cx = (vx(1)+vx(2))/2;
    cy = (vy(1)+vy(2))/2;
    cz = (vz(1)+vz(2))/2;
    [v3, f3] = a_circle(cx, cy, cz, len*0.3, dn, ir);    
    [v4, f4] = a_plus(cx, cy+len/10, cz, len/5, dn, ir);
    [v5, f5] = a_minus(cx, cy-len/10, cz, len/5, dn, ir, 'y');    
case ('zn')
    vx = 0.5*[(ix+ax) (ix+ax)]; vy = 0.5*[(iy+ay) (iy+ay)]; vz = [az iz];
    len = az - iz;
    ir = len/30;
    [v1, f1] = a_line(vx, vy, [iz iz+0.2*len], ir);
    [v2, f2] = a_line(vx, vy, [iz+0.8*len az], ir); 
    if ((ax-ix)<(ay-iy)) dn = 'x'; else dn = 'y'; end
    cx = (vx(1)+vx(2))/2;
    cy = (vy(1)+vy(2))/2;
    cz = (vz(1)+vz(2))/2;
    [v3, f3] = a_circle(cx, cy, cz, len*0.3, dn, ir);    
    [v4, f4] = a_plus(cx, cy, cz-len/10, len/5, dn, ir);
    [v5, f5] = a_minus(cx, cy, cz+len/10, len/5, dn, ir, 'z');
case ('zp')
    vx = 0.5*[(ix+ax) (ix+ax)]; vy = 0.5*[(iy+ay) (iy+ay)]; vz = [iz az];
    len = az - iz;
    ir = len/30;
    [v1, f1] = a_line(vx, vy, [iz iz+0.2*len], ir);
    [v2, f2] = a_line(vx, vy, [iz+0.8*len az], ir); 
    if ((ax-ix)<(ay-iy)) dn = 'x'; else dn = 'y'; end
    cx = (vx(1)+vx(2))/2;
    cy = (vy(1)+vy(2))/2;
    cz = (vz(1)+vz(2))/2;
    [v3, f3] = a_circle(cx, cy, cz, len*0.3, dn, ir);    
    [v4, f4] = a_plus(cx, cy, cz+len/10, len/5, dn, ir);
    [v5, f5] = a_minus(cx, cy, cz-len/10, len/5, dn, ir, 'z');
end
v = [v1; v2; v3; v4; v5];
f = f1;
f = [f; max(max(f)) + f2];
f = [f; max(max(f)) + f3];
f = [f; max(max(f)) + f4];
f = [f; max(max(f)) + f5];
patch('Vertices',v,'Faces',f,'facecolor','g', ...
    'facealpha',1,'edgecolor','none');

% draw boundaries
v = [0 0 0; 1 0 0; 1 1 0; 0 1 0; 0 0 1; 1 0 1; 1 1 1; 0 1 1];
vx = v(:,1) * (ax - ix) + ix;
vy = v(:,2) * (ay - iy) + iy;
vz = v(:,3) * (az - iz) + iz;
f = [1 2 6 5; 2 3 7 6; 3 4 8 7; 4 1 5 8; 1 2 3 4; 5 6 7 8];
patch('Vertices',[vx vy vz],'Faces',f,'facecolor','none','LineWidth',1, ...
    'facealpha',0.3,'edgecolor','g','linestyle','-');

% show label
text(ix,iy,iz,obj.obj_id,'VerticalAlignment','bottom', ...
    'BackgroundColor','w','edgecolor','g','fontsize',12);

%===================================================================
function draw_a_resistor(obj)

ix = obj.min_x; iy = obj.min_y; iz = obj.min_z; 
ax = obj.max_x; ay = obj.max_y; az = obj.max_z;

t = 0.1 * ones(1,11);
[X,Y,Z] = cylinder(t,18);
X(5,:) = X(5,:) + 0.2;
X(7,:) = X(7,:) - 0.2;

switch (obj.direction)
case ('x')
    len = ax - ix;
    vx = ix; vy = 0.5*(iy+ay); vz = 0.5*(iz+az);
    X = X*len; Y = Y*len; Z = Z*len;    
    if ((az-iz)<(ay-iy)) 
        [f,v] = surf2patch(Z+vx,X+vy,Y+vz,Z);       
    else
        [f,v] = surf2patch(Z+vx,Y+vy,X+vz,Z);       
    end
case ('y')
    len = ay - iy;
    vx = 0.5*(ix+ax); vy = iy; vz = 0.5*(iz+az);
    X = X*len; Y = Y*len; Z = Z*len;    
    if ((ax-ix)<(az-iz))
        [f,v] = surf2patch(Y+vx,Z+vy,X+vz,Z);       
    else
        [f,v] = surf2patch(X+vx,Z+vy,Y+vz,Z);       
    end
case ('z')
    len = az - iz;
    vy = 0.5*(iy+ay); vx = 0.5*(ix+ax); vz = iz;
    X = X*len; Y = Y*len; Z = Z*len;    
    if ((ax-ix)<(ay-iy)) 
        [f,v] = surf2patch(Y+vx,X+vy,Z+vz,Z);       
    else
        [f,v] = surf2patch(X+vx,Y+vy,Z+vz,Z);       
    end
end
patch('Vertices',v,'Faces',f,'facecolor','k', ...
    'facealpha',1,'edgecolor','none');

% draw boundaries
v = [0 0 0; 1 0 0; 1 1 0; 0 1 0; 0 0 1; 1 0 1; 1 1 1; 0 1 1];
vx = v(:,1) * (ax - ix) + ix;
vy = v(:,2) * (ay - iy) + iy;
vz = v(:,3) * (az - iz) + iz;
f = [1 2 6 5; 2 3 7 6; 3 4 8 7; 4 1 5 8; 1 2 3 4; 5 6 7 8];
patch('Vertices',[vx vy vz],'Faces',f,'facecolor','none','LineWidth',1, ...
    'facealpha',0.3,'edgecolor','k','linestyle','-');

% show label
text(ix,iy,iz,obj.obj_id,'VerticalAlignment','bottom', ...
    'BackgroundColor','w','edgecolor','k','fontsize',12);

%===================================================================
function draw_a_capacitor(obj)

ix = obj.min_x; iy = obj.min_y; iz = obj.min_z; 
ax = obj.max_x; ay = obj.max_y; az = obj.max_z;

t = 0.1 * ones(1,11);
[X,Y,Z] = cylinder(t,18);
X(5,:) = X(5,:) + 0.2;
X(7,:) = X(7,:) - 0.2;

switch (obj.direction)
case ('x')
    len = ax - ix;
    vx = ix; vy = 0.5*(iy+ay); vz = 0.5*(iz+az);
    ir = len/20;
    [v1, f1] = a_line([ix ix+0.4*len], [vy vy], [vz vz], ir);
    [v2, f2] = a_line([ix+0.6*len ax], [vy vy], [vz vz], ir);
    [v3, f3] = a_line([ix+0.4*len ix+0.42*len], [vy vy], [vz vz], 4*ir);
    [v4, f4] = a_line([ix+0.58*len ix+0.6*len], [vy vy], [vz vz], 4*ir);
case ('y')
    len = ay - iy;