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(
没有合适的资源?快使用搜索试试~ 我知道了~
chapitre4.zip_The Method Method_fdtd MATLAB_researchlgr
共49个文件
m:49个
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 147 浏览量
2022-07-15
21:25:02
上传
评论
收藏 39KB ZIP 举报
温馨提示
The Finite-Difference Time-Domain Method for Electromagnetics with MATLAB Simulations chapter 4
资源推荐
资源详情
资源评论
收起资源包目录
chapitre4.zip (49个子文件)
chapitre4
capture_sampled_voltages.m 487B
initialize_display_parameters.m 2KB
display_3D_geometry.m 9KB
initialize_waveforms.m 519B
initialize_capacitor_updating_coefficients.m 2KB
display_problem_space.m 11KB
initialize_sources_and_lumped_elements.m 9KB
capture_sampled_magnetic_fields.m 915B
set_color_axis_scaling.m 1KB
display_animation.m 1KB
plot_e_xy.m 1KB
post_process_and_display_results.m 70B
initialize_resistor_updating_coefficients.m 2KB
display_sampled_parameters.m 3KB
object_drawing_functions.m 36KB
plot_h_xy.m 2KB
update_current_sources.m 660B
update_electric_fields.m 1KB
initialize_animation_parameters.m 6KB
update_voltage_sources.m 662B
update_inductors.m 827B
plot_e_yz.m 1KB
display_objects_mesh_in_animation.m 3KB
initialize_diode_updating_coefficients.m 2KB
update_magnetic_fields.m 524B
plot_e_zx.m 1KB
display_transient_parameters.m 4KB
initialize_updating_coefficients.m 2KB
plot_h_yz.m 2KB
run_fdtd_time_marching_loop.m 777B
display_material_mesh_gui.m 35KB
initialize_current_source_updating_coefficients.m 2KB
exemple_4_a
define_problem_space_parameters.m 2KB
define_geometry.m 485B
define_output_parameters.m 1KB
define_sources_and_lumped_elements.m 1KB
solve_diode_equation.m 418B
create_linear_index_list.m 677B
plot_h_zx.m 2KB
initialize_output_parameters.m 4KB
update_diodes.m 947B
initialize_inductor_updating_coefficients.m 1KB
initialize_fdtd_parameters_and_arrays.m 711B
initialize_voltage_source_updating_coefficients.m 2KB
define_output_parameters.m 3KB
capture_sampled_currents.m 1KB
capture_sampled_electric_fields.m 835B
define_sources_and_lumped_elements.m 3KB
display_material_mesh.m 425B
共 49 条
- 1
资源评论
我虽横行却不霸道
- 粉丝: 75
- 资源: 1万+
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功