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
没有合适的资源?快使用搜索试试~ 我知道了~
程序.rar_RCS_gate6zn_n rcs_散射 时间步_柱 rcs
共64个文件
m:63个
mat:1个
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 135 浏览量
2022-07-15
02:17:56
上传
评论
收藏 348KB RAR 举报
温馨提示
计算一个方柱的散射的波形,并求其RCS。方柱的网格总数是N=500个,时间步为1200步。网格尺寸为lamda/18,时间间隔是ddx/2*c,c为光速。脉冲源位于中心。方柱边长为7*lamda。并与相同尺度的球的散射场进行对比
资源推荐
资源详情
资源评论
收起资源包目录
程序.rar (64个子文件)
程序
create_spheres.m 546B
display_frequency_domain_outputs.m 5KB
calculate_material_component_values.m 5KB
capture_sampled_voltages.m 473B
initialize_display_parameters.m 2KB
calculate_domain_size.m 5KB
display_3D_geometry.m 8KB
initialize_waveforms.m 3KB
define_problem_space_parameters.m 3KB
initialize_capacitor_updating_coefficients.m 2KB
initialize_boundary_conditions.m 1KB
display_problem_space.m 11KB
initialize_sources_and_lumped_elements.m 8KB
capture_sampled_magnetic_fields.m 891B
set_color_axis_scaling.m 1KB
display_animation.m 1KB
plot_e_xy.m 1KB
frequency_to_time_domain.m 578B
post_process_and_display_results.m 145B
define_geometry.m 2KB
initialize_resistor_updating_coefficients.m 2KB
display_sampled_parameters.m 3KB
fdtd_solve.m 736B
object_drawing_functions.m 35KB
time_to_frequency_domain.m 646B
plot_h_xy.m 2KB
update_current_sources.m 645B
update_electric_fields.m 1KB
initialize_animation_parameters.m 6KB
update_voltage_sources.m 645B
initialize_CPML_ABC.m 15KB
update_inductors.m 804B
calculate_frequency_domain_outputs.m 3KB
plot_e_yz.m 1KB
initialize_fdtd_material_grid.m 1KB
display_objects_mesh_in_animation.m 3KB
initialize_diode_updating_coefficients.m 2KB
update_magnetic_fields.m 512B
plot_e_zx.m 1KB
display_transient_parameters.m 4KB
create_PEC_plates.m 1KB
initialize_updating_coefficients.m 2KB
plot_h_yz.m 2KB
run_fdtd_time_marching_loop.m 828B
mesh_data_file.mat 3.85MB
display_material_mesh_gui.m 34KB
initialize_current_source_updating_coefficients.m 2KB
solve_diode_equation.m 402B
create_linear_index_list.m 655B
plot_h_zx.m 2KB
initialize_output_parameters.m 4KB
update_diodes.m 917B
update_electric_field_CPML_ABC.m 3KB
initialize_inductor_updating_coefficients.m 1KB
initialize_fdtd_parameters_and_arrays.m 688B
update_magnetic_field_CPML_ABC.m 3KB
initialize_voltage_source_updating_coefficients.m 2KB
create_bricks.m 670B
recover_a_time_waveform.m 1KB
define_output_parameters.m 5KB
capture_sampled_currents.m 1KB
capture_sampled_electric_fields.m 812B
define_sources_and_lumped_elements.m 2KB
display_material_mesh.m 414B
共 64 条
- 1
资源评论
钱亚锋
- 粉丝: 88
- 资源: 1万+
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功