Matlab模拟的电磁学时域有限差分法电子书+程序
标签:
•
文件类型: .rar
•
文件大小:
•
下载次数: 1
•
Matlab模拟的电磁学时域有限差分法(电子书+程序)
代码片段和文件信息
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 3184 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Appendix_Afdtd_1d_code.m
文件 836 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Appendix_Ainitialize_plotting_parameters.m
文件 353 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Appendix_Aplot_fields.m
文件 3206 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Appendix_Cpolar_plot_constant_phi.m
文件 2418 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Appendix_Cpolar_plot_constant_theta.m
文件 4529 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_1figure_1_2.m
文件 5155 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_10calculate_domain_size.m
文件 600 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_10define_geometry.m
文件 606 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_10example_10_adefine_geometry.m
文件 1364 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_10example_10_adefine_output_parameters.m
文件 2723 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_10example_10_adefine_problem_space_parameters.m
文件 853 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_10example_10_adefine_sources_and_lumped_elements.m
文件 2640 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_10initialize_thin_wire_updating_coefficients.m
文件 1820 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_10initialize_updating_coefficients.m
文件 4054 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11calculate_and_display_farfields.m
文件 7067 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11calculate_farfields_per_plane.m
文件 506 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11calculate_incident_plane_wave_power.m
文件 661 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11define_sources_and_lumped_elements.m
文件 4680 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11display_transient_parameters.m
文件 250 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_adefine_geometry.m
文件 1352 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_adefine_output_parameters.m
文件 2718 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_adefine_problem_space_parameters.m
文件 661 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_adefine_sources_and_lumped_elements.m
文件 303 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_bdefine_geometry.m
文件 964 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_bdefine_output_parameters.m
文件 2712 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_bdefine_problem_space_parameters.m
文件 663 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_bdefine_sources_and_lumped_elements.m
文件 1826 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_ccapture_sampled_electric_fields.m
文件 292 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_cdefine_geometry.m
文件 1803 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_cdefine_output_parameters.m
............此处省略437个文件信息
% This program demonstrates a one-dimensional FDTD simulation.
% The problem geometry is composed of two PEC plates extending to
% infinity in y and z dimensions parallel to each other with 1 meter
% separation. The space between the PEC plates is filled with air.
% A sheet of current source paralle to the PEC plates is placed
% at the center of the problem space. The current source excites fields
% in the problem space due to a z-directed current density Jz
% which has a Gaussian waveform in time.
% Define initial constants
eps_0 = 8.854187817e-12; % permittivity of free space
mu_0 = 4*pi*1e-7; % permeability of free space
c = 1/sqrt(mu_0*eps_0); % speed of light
% Define problem geometry and parameters
domain_size = 1; % 1D problem space length in meters
dx = 1e-3; % cell size in meters
dt = 3e-12; % duration of time step in seconds
number_of_time_steps = 2000; % number of iterations
nx = round(domain_size/dx); % number of cells in 1D problem space
source_position = 0.5; % position of the current source Jz
% Initialize field and material arrays
Ceze = zeros(nx+11);
Cezhy = zeros(nx+11);
Cezj = zeros(nx+11);
Ez = zeros(nx+11);
Jz = zeros(nx+11);
eps_r_z = ones (nx+11); % free space
sigma_e_z = zeros(nx+11); % free space
Chyh = zeros(nx1);
Chyez = zeros(nx1);
Chym = zeros(nx1);
Hy = zeros(nx1);
My = zeros(nx1);
mu_r_y = ones (nx1); % free space
sigma_m_y = zeros(nx1); % free space
% Calculate FDTD updating coefficients
Ceze = (2 * eps_r_z * eps_0 - dt * sigma_e_z) ...
./(2 * eps_r_z * eps_0 + dt * sigma_e_z);
Cezhy = (2 * dt / dx) ...
./(2 * eps_r_z * eps_0 + dt * sigma_e_z);
Cezj = (-2 * dt) ...
./(2 * eps_r_z * eps_0 + dt * sigma_e_z);
Chyh = (2 * mu_r_y * mu_0 - dt * sigma_m_y) ...
./(2 * mu_r_y * mu_0 + dt * sigma_m_y);
Chyez = (2 * dt / dx) ...
./(2 * mu_r_y * mu_0 + dt * sigma_m_y);
Chym = (-2 * dt) ...
./(2 * mu_r_y * mu_0 + dt * sigma_m_y);
% Define the Gaussian source waveform
time = dt*[0:number_of_time_steps-1].‘;
Jz_waveform = exp(-((time-2e-10)/5e-11).^2);
source_position_index = round(nx*source_position/domain_size)+1;
% Subroutine to initialize plotting
initialize_plotting_parameters;
% FDTD loop
for time_step = 1:number_of_time_steps
% Update Jz for the current time step
Jz(source_position_index) = Jz_waveform(time_step);
% Update magnetic field
Hy(1:nx) = Chyh(1:nx) .* Hy(1:nx) ...
+ Chyez(1:nx) .* (Ez(2:nx+1) - Ez(1:nx)) ...
+ Chym(1:nx) .* My(1:nx);
% Update electric field
Ez(2:nx) = Ceze (2:nx) .* Ez(2:nx) ...
+ Cezhy(2:nx) .* (Hy(2:nx) - Hy(1:nx-1)) ...
+ Cezj(2:nx) .* Jz(2:nx);
Ez(1) = 0; % A
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 3184 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Appendix_Afdtd_1d_code.m
文件 836 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Appendix_Ainitialize_plotting_parameters.m
文件 353 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Appendix_Aplot_fields.m
文件 3206 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Appendix_Cpolar_plot_constant_phi.m
文件 2418 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Appendix_Cpolar_plot_constant_theta.m
文件 4529 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_1figure_1_2.m
文件 5155 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_10calculate_domain_size.m
文件 600 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_10define_geometry.m
文件 606 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_10example_10_adefine_geometry.m
文件 1364 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_10example_10_adefine_output_parameters.m
文件 2723 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_10example_10_adefine_problem_space_parameters.m
文件 853 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_10example_10_adefine_sources_and_lumped_elements.m
文件 2640 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_10initialize_thin_wire_updating_coefficients.m
文件 1820 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_10initialize_updating_coefficients.m
文件 4054 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11calculate_and_display_farfields.m
文件 7067 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11calculate_farfields_per_plane.m
文件 506 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11calculate_incident_plane_wave_power.m
文件 661 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11define_sources_and_lumped_elements.m
文件 4680 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11display_transient_parameters.m
文件 250 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_adefine_geometry.m
文件 1352 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_adefine_output_parameters.m
文件 2718 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_adefine_problem_space_parameters.m
文件 661 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_adefine_sources_and_lumped_elements.m
文件 303 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_bdefine_geometry.m
文件 964 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_bdefine_output_parameters.m
文件 2712 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_bdefine_problem_space_parameters.m
文件 663 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_bdefine_sources_and_lumped_elements.m
文件 1826 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_ccapture_sampled_electric_fields.m
文件 292 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_cdefine_geometry.m
文件 1803 2012-09-07 15:40 MATLAB模拟的电磁学时域有限差分法Chapter_11example_11_cdefine_output_parameters.m
............此处省略437个文件信息
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌抄袭侵权/违法违规的内容, 请发送邮件举报,一经查实,本站将立刻删除。
评论列表(条)