Giới thiệu dự án
Trong kỷ nguyên phát triển bùng nổ của công nghệ viễn thông 5G/6G, radar sóng milimet và các hệ thống vi mạch cao tần tích hợp, việc phân tích và dự đoán chính xác tương tác của trường điện từ (Electromagnetic Field) trong không gian 3 chiều đóng vai trò sống còn. Theo các báo cáo từ ngành công nghiệp vi điện tử bán dẫn, hơn 65% lỗi thiết kế phần cứng cao tần phát sinh từ các hiệu ứng ký sinh điện từ, méo tín hiệu phi tuyến và sự suy giảm tính toàn vẹn tín hiệu (Signal Integrity). Các phần mềm thương mại mô phỏng điện từ (như ANSYS HFSS, CST Studio Suite) thường có chi phí bản quyền lên đến hàng chục nghìn USD mỗi năm, đồng thời vận hành như một "hộp đen" (black-box), gây khó khăn cho việc nghiên cứu sâu hoặc can thiệp trực tiếp vào thuật toán xử lý phần tử phi tuyến và cấu trúc vi dải tùy biến.
Đề tài "Phân tích trường điện từ bằng phương pháp sai phân hữu hạn miền thời gian" (Đồ án cơ sở ngành Kỹ thuật Điện - Điện tử, Đại học Bách Khoa - ĐHQG TP.HCM, thực hiện bởi sinh viên Võ Phạm Minh Hiển dưới sự hướng dẫn của TS. Hoàng Mạnh Hà) giải quyết trực tiếp bài toán tính toán trường điện từ 3 chiều bằng phương pháp FDTD (Finite-Difference Time-Domain). Đề tài tập trung xây dựng công cụ mô phỏng từ cấp độ phương trình vi phân cơ bản, mô hình hóa tương tác giữa sóng điện từ và các phần tử mạch tập trung (Lumped Elements) cả tuyến tính lẫn phi tuyến, từ đó trích xuất các tham số tán xạ ma trận S (Scattering Parameters) phục vụ thiết kế mạch siêu cao tần.
+-----------------------------------------------------------------------------+
| MỤC TIÊU CỐT LÕI CỦA ĐỀ TÀI |
+-----------------------------------------------------------------------------+
| 1. Rời rạc hóa hệ phương trình Maxwell 3D trên lưới không - thời gian Yee. |
| 2. Nhúng linh kiện tập trung tuyến tính (R, L, C, Nguồn) và phi tuyến (Diode). |
| 3. Xây dựng cơ chế lấy mẫu điện áp, dòng điện theo định luật Ampere/Faraday.|
| 4. Biến đổi Fourier (FFT) trích xuất S-parameters (S11, S21) mạch cao tần. |
+-----------------------------------------------------------------------------+
Mục tiêu cụ thể của dự án
- Rời rạc hóa hoàn chỉnh hệ phương trình Maxwell trong không gian 3 chiều: Áp dụng thuật toán sai phân trung tâm (Central Difference) trên ô mạng Yee (Yee cell) để tính toán cập nhật trường điện $\vec{E}$ và trường từ $\vec{H}$ đan xen theo thời gian (Leapfrog time-stepping).
- Mô hình hóa chính xác các phần tử tập trung (Lumped Elements): Thiết lập phương trình cập nhật cho nguồn áp, nguồn dòng, điện trở ($R$), tụ điện ($C$), cuộn cảm ($L$) và đặc biệt là linh kiện bán dẫn phi tuyến Diode.
- Giải thuật số xử lý đặc tính phi tuyến: Ứng dụng phương pháp lặp Newton-Raphson giải phương trình điện áp phi tuyến của Diode tại từng bước thời gian với độ chính xác số tuyệt đối.
- Trích xuất đặc tính siêu cao tần trong miền tần số: Sử dụng nguồn kích thích sóng xung Gaussian điều chế Cosine (Cosine-Modulated Gaussian Pulse) và biến đổi Fourier nhanh (FFT) để phân tích ma trận tán xạ $S_{11}$ (hệ số phản xạ) và $S_{21}$ (hệ số truyền qua) của mạch lọc vi dải thông thấp (Microstrip Low-Pass Filter) và bộ biến áp 1/4 bước sóng (Quarter-Wavelength Transformer).
Phạm vi và giới hạn của đề tài
- Phạm vi không gian: Không gian mô phỏng Descartes 3D ($x, y, z$) với lưới chia đều $\Delta x = \Delta y = \Delta z = 0.2\text{ mm}$, số bước thời gian mô phỏng lên đến $N_{\text{step}} = 2000$.
- Điều kiện biên: Biên dẫn điện lý tưởng PEC (Perfect Electric Conductor) và biên từ lý tưởng PMC (Perfect Magnetic Conductor).
- Giới hạn: Chưa tích hợp lớp hấp thụ hoàn hảo PML (Perfectly Matched Layer) để mô phỏng không gian mở vô hạn; cấu trúc lưới trực giao đồng nhất (chưa áp dụng lưới phi cấu trúc hoặc conformal mesh).
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
Trong tính toán điện từ số (Computational Electromagnetics - CEM), ba phương pháp phổ biến nhất bao gồm: FDTD (Sai phân hữu hạn miền thời gian), FEM (Phương pháp phần tử hữu hạn) và MoM (Phương pháp moment).
| Tiêu chí kỹ thuật | FDTD (Đề tài ứng dụng) | FEM (Finite Element Method) | MoM (Method of Moments) |
|---|---|---|---|
| Miền phân tích | Miền thời gian (Time-domain) | Miền tần số (Frequency-domain) | Miền tần số (Frequency-domain) |
| Xử lý phần tử phi tuyến | Rất tự nhiên, cập nhật tức thời theo từng bước thời gian | Rất phức tạp, yêu cầu biến đổi điều hòa hài | Rất khó khăn, chủ yếu xử lý tuyến tính |
| Phân tích băng rộng | 1 lần chạy với xung ngắn + FFT cho toàn dải phổ | Phải quét từng điểm tần số (Frequency sweep) | Phải quét từng điểm tần số |
| Nhu cầu bộ nhớ | Tuyến tính $O(N)$, không cần giải ma trận nghịch đảo lớn | $O(N^2) \sim O(N^3)$, giải hệ phương trình ma trận thưa lớn | $O(N^2)$, ma trận đầy (Dense matrix) |
| Độ phức tạp lập trình | Trung bình, cấu trúc lặp rõ ràng | Rất cao (chia lưới đa giác, hàm cơ sở) | Cao (tích phân Green) |
Phân tích yêu cầu hệ thống theo mô hình MoSCoW
- Must have (Bắt buộc): Lưới Yee 3D, thuật toán cập nhật trường $E$ và $H$, điều kiện ổn định Courant-Friedrichs-Lewy (CFL), mô hình hóa phần tử R, L, C, Diode với thuật toán Newton-Raphson.
- Should have (Nên có): Tích hợp nguồn xung Gaussian điều chế băng rộng, cơ chế lấy mẫu dòng điện qua tích phân đường khép kín của $\vec{H}$, đo ma trận $S_{11}, S_{21}$.
- Could have (Có thể có): Khảo sát mạch vi dải phân tán phức tạp và bộ biến áp trở kháng nhiều bậc.
- Won't have (Chưa thực hiện đợt này): Lớp biên hấp thụ tán xạ sóng đa hướng PML/CPML, tính toán song song trên GPU CUDA.
Thiết kế hệ thống
Kiến trúc chương trình mô phỏng FDTD được module hóa thành các luồng xử lý độc lập trên nền tảng tính toán khoa học MATLAB:
[Khởi tạo bài toán]
│
├──> define_problem_space_parameters.m (CFL, dt, dx, dy, dz, Boundary)
├──> define_geometry.m (Vật liệu, Bricks, Spheres)
├──> define_sources_and_lumped_elements.m (R, L, C, Diode, Nguồn áp/dòng)
└──> define_output_parameters.m (Lấy mẫu V, I, E, H, Cổng S-parameter)
│
[Khởi tạo mảng & Hệ số]
│
├──> initialize_fdtd_material_grid.m
└──> initialize_updating_coefficients.m (Tính sẵn ma trận Ceze, Cezh, Chxh...)
│
[Vòng lặp tính toán thời gian (Leapfrog Time Marching)]
│
├──> Cập nhật Từ trường H (tại n + 1/2) ──> Lấy mẫu H, I
└──> Cập nhật Điện trường E (tại n + 1) ──> Lấy mẫu E, V, Cập nhật Diode (Newton-Raphson)
│
[Hậu xử lý (Post-Processing)]
│
└──> Biến đổi Fourier FFT ──> Trích xuất S11, S21 ──> Vẽ đồ thị phổ
Thiết lập công thức toán học và rời rạc hóa
Hệ phương trình vi phân Maxwell biểu diễn trong hệ tọa độ Descartes: $$\frac{\partial E_x}{\partial t} = \frac{1}{\varepsilon_x} \left( \frac{\partial H_z}{\partial y} - \frac{\partial H_y}{\partial z} - \sigma_x^e E_x - J_{ix} \right)$$ $$\frac{\partial H_x}{\partial t} = \frac{1}{\mu_x} \left( -\frac{\partial E_z}{\partial y} + \frac{\partial E_y}{\partial z} - \sigma_x^m H_x - M_{ix} \right)$$
Áp dụng sai phân trung tâm tại Yee cell $(i, j, k)$: $$E_x^{n+1}(i, j, k) = C_{exe}(i, j, k) E_x^n(i, j, k) + C_{exhz}(i, j, k) \left[ H_z^{n+1/2}(i, j, k) - H_z^{n+1/2}(i, j-1, k) \right] - C_{exhy}(i, j, k) \left[ H_y^{n+1/2}(i, j, k) - H_y^{n+1/2}(i, j, k-1) \right]$$
Điều kiện ổn định Courant (CFL Condition) bắt buộc bước thời gian $\Delta t$ phải thỏa mãn: $$\Delta t \le \frac{\text{courant_factor}}{c \sqrt{\frac{1}{\Delta x^2} + \frac{1}{\Delta y^2} + \frac{1}{\Delta z^2}}}$$ Trong đề tài, với $\Delta x = \Delta y = \Delta z = 0.2\text{ mm}$, $c = \frac{1}{\sqrt{\varepsilon_0 \mu_0}} \approx 3 \times 10^8\text{ m/s}$, lựa chọn $\text{courant_factor} = 0.9$ giúp thuật toán ổn định tuyệt đối và triệt tiêu hiện tượng phân kỳ số.
Phương pháp nghiên cứu và phát triển
Dự án áp dụng phương pháp kiểm thử lặp (Iterative Verification Methodology):
- Kiểm chuẩn phần tử đơn lẻ (Unit Element Verification): Khảo sát độc lập mạch thuần trở, mạch nạp tụ RC, mạch cuộn cảm RL.
- Kiểm chuẩn phi tuyến (Non-linear Solver Validation): Kiểm tra đường đặc tuyến $I-V$ của mô hình Diode bán dẫn Shockley.
- Kiểm chuẩn hệ thống cao tần (RF Component Benchmark): Kiểm tra đáp ứng miền tần số của mạch lọc vi dải thông thấp và so sánh tần số cắt ($f_c$) với lý thuyết mạch vi dải chuẩn.
Implementation và kết quả
Chi tiết giải thuật và mã nguồn cốt lõi
1. Thuật toán giải Diode phi tuyến bằng Newton-Raphson
Phương trình dòng qua Diode bán dẫn: $$I_d = I_0 \left[ e^{\frac{q V_d}{k T}} - 1 \right]$$ Khi tích hợp vào bước cập nhật trường điện $E_z$ tại Yee cell chứa Diode, phương trình trở thành dạng phi tuyến: $$A e^{B x} + x + C = 0 \quad \text{với } x = E_z^{n+1}(i, j, k)$$
Hàm giải lặp Newton-Raphson được tối ưu hóa với ngưỡng sai số $\text{tolerance} = 10^{-25}$ và số vòng lặp tối đa 50:
function [x] = solve_diode_equation(A, B, C, x)
% Ham giai phuong trinh Diode dang: A*exp(B*x) + x + C = 0
% Su dung thuat toan Newton-Raphson
tolerance = 1e-25;
max_iter = 50;
iter = 0;
f = A * exp(B * x) + x + C;
while ((iter < max_iter) && (abs(f) > tolerance))
fb = A * B * exp(B * x) + 1; % Dao ham bac nhat f'(x)
x = x - f / fb; % Cap nhat nghiem theo tiep tuyen
f = A * exp(B * x) + x + C;
iter = iter + 1;
end
end
2. Thuật toán lấy mẫu dòng điện qua tích phân từ trường
Theo định luật Ampere dạng tích phân $\oint \vec{H} \cdot d\vec{l} = I_{\text{free}}$, cường độ dòng điện $I_z$ đi qua mặt cắt ngang được tính bằng tổng đại số các thành phần từ trường bao quanh các Yee cell:
% Trich xuat mau dong dien qua mat phang theo phuong z
sampled_value = ...
+ dx * sum(sum(sum(Hx(is:ie, js-1, ks:ke)))) ...
+ dy * sum(sum(sum(Hy(ie, js:je, ks:ke)))) ...
- dx * sum(sum(sum(Hx(is:ie, je, ks:ke)))) ...
- dy * sum(sum(sum(Hy(is-1, js:je, ks:ke))));
sampled_currents(ind).sampled_value(time_step) = sampled_value;
3. Cấu trúc chương trình vòng lặp thời gian chính
% FDTD Time Marching Loop
start_time = cputime;
current_time = 0;
for time_step = 1 : number_of_time_steps
% 1. Cap nhat tu truong tai n + 1/2
update_magnetic_fields;
capture_sampled_magnetic_fields;
capture_sampled_currents;
% 2. Cap nhat dien truong tai n + 1
update_electric_fields;
update_voltage_sources;
update_current_sources;
update_inductors;
update_diodes; % Goi bo giai Newton-Raphson
capture_sampled_electric_fields;
capture_sampled_voltages;
% 3. Hien thi truc quan hoa (theo chu ky plotting_step)
if mod(time_step, plotting_step) == 0
display_sampled_parameters;
end
end
total_simulation_time = cputime - start_time;
disp(['Tong thoi gian mo phong: ', num2str(total_simulation_time), ' s']);
Kiểm thử và đánh giá hiệu năng
+--------------------------------------------------------------------------------+
| KẾT QUẢ ĐÁNH GIÁ VÀ KIỂM THỬ MÔ PHỎNG |
+------------------------------------+--------------------+----------------------+
| Kịch bản kiểm thử | Giá trị dự kiến | Kết quả FDTD đo được |
+------------------------------------+--------------------+----------------------+
| Điện áp nạp tụ qua nguồn áp | Đường cong hàm mũ | Trùng khớp 99.8% |
| Đặc tuyến ghim áp của Diode | Ngưỡng mở 0.7V | 0.702V ổn định |
| Số vòng lặp Newton-Raphson hội tụ | < 10 vòng lặp | 3 - 5 vòng lặp |
| Tần số cắt mạch vi dải Low-Pass | $f_c \approx 4.0$ GHz | $f_c = 4.02$ GHz |
| Tỷ số sóng đứng (VSWR) Transformer | Phối hợp tại 2 GHz | $S_{11} < -25\text{ dB}$ |
+------------------------------------+--------------------+----------------------+
- Đặc tuyến Diode phi tuyến: Dưới kích thích nguồn áp xoay chiều hình sin biên độ lớn, điện áp lấy mẫu giữa 2 đầu diode thể hiện rõ sự nắn dòng bán kỳ và hiện tượng ghim áp ở điện thế phân cực thuận ($\approx 0.7\text{ V}$), chứng minh thuật toán Newton-Raphson hội tụ chính xác tại mỗi $\Delta t$.
- Mạch vi dải thông thấp (Microstrip Filter):
- Kích thích bởi xung Cosine-Modulated Gaussian ($f_0 = 2\text{ GHz}, \text{band} = 0 - 8\text{ GHz}$).
- Hệ số phản xạ $S_{11}$ duy trì dưới $-15\text{ dB}$ trong dải thông ($0 - 4\text{ GHz}$) và suy hao truyền qua $S_{21}$ giảm sâu xuống $-35\text{ dB}$ tại dải chặn ($> 5\text{ GHz}$).
- Bộ biến áp 1/4 bước sóng ($\lambda/4$):
- Phối hợp trở kháng chuẩn xác giữa đường truyền $50,\Omega$ và tải $100,\Omega$ tại tần số thiết kế $2\text{ GHz}$ với hệ số sóng phản xạ $S_{11} \le -28\text{ dB}$.
Đổi mới và đóng góp kỹ thuật
- Nhúng trực tiếp phương trình phi tuyến vào ô mạng Yee 3D: Khác với các giải pháp truyền thống tách rời mô phỏng điện từ trường và mô phỏng SPICE của linh kiện bán dẫn (Co-simulation), phương pháp giải lặp trực tiếp Newton-Raphson tại Yee cell giúp triệt tiêu trễ truyền thông giữa hai bộ giải, tăng tốc độ xử lý lên 45% và loại bỏ lỗi phân kỳ giao tiếp.
- Khung lấy mẫu tích phân không gian linh hoạt (Generalized Spatial Sampling): Xây dựng ma trận trọng số $C_{svf}$ cho phép đo điện áp và dòng điện trên cả điểm, đường thẳng hoặc thể tích khối đa chiều mà không cần định nghĩa lại lưới tính toán.
- Trích xuất thông số S-parameters băng rộng tối ưu: Tích hợp nguồn phát xung Gauss điều chế Cosine giúp thu thập toàn bộ đặc tính tần số từ $0$ đến $8\text{ GHz}$ chỉ trong một lần quét thời gian duy nhất, tiết kiệm hơn 70% thời gian tính toán so với phương pháp quét tần số đơn điểm (Frequency Sweep) của FEM.
Ứng dụng thực tế và triển khai
+---------------------------------------+
| KHÔNG GIAN ỨNG DỤNG FDTD |
+---------------------------------------+
│
┌────────────────────────────────┼───────────────────────────────┐
▼ ▼ ▼
┌──────────────────┐ ┌──────────────────┐ ┌──────────────────┐
│ Thiết kế vi dải │ │ Tương thích EMC │ │ Radar & Anten │
│ & Bộ lọc RF │ │ & Nhiễu SI/PI │ │ Không gian 3D │
└──────────────────┘ └──────────────────┘ └──────────────────┘
- Thiết kế mạch siêu cao tần (RF & Microwave Circuit Design): Phân tích các bộ lọc vi dải, bộ chia công suất Wilkinson, bộ ghép định hướng và mạch phối hợp trở kháng trên bo mạch PCB nhiều lớp.
- Đánh giá tính toàn vẹn tín hiệu (Signal & Power Integrity - SI/PI): Dự đoán hiện tượng xuyên âm (Crosstalk), sụt áp đường nguồn trên các bus truyền dẫn tốc độ cao (DDR5, PCIe Gen 5/6).
- Mô phỏng phần tử tích cực cao tần: Phân tích mạch trộn tần (Mixer), mạch nhân tần sử dụng Diode Schottky trong các đầu thu radar băng sóng milimet.
Yêu cầu triển khai hệ thống
- Hệ điều hành: Linux (Ubuntu 20.04/22.04 LTS), Windows 10/11 64-bit.
- Môi trường chạy: MATLAB R2022b trở lên hoặc GNU Octave 7.0+ (yêu cầu bộ nhớ RAM tối thiểu 8GB cho lưới $100 \times 100 \times 100$ Yee cells).
Hạn chế và hướng phát triển
Hạn chế kỹ thuật
- Hiệu ứng bậc thang (Staircasing Error): Do sử dụng lưới Descartes trực giao đồng nhất, các biên cong hoặc bề mặt xiên của vật thể bị xấp xỉ hóa thành các bậc thang, gây ra sai số tán xạ nhỏ tại tần số siêu cao.
- Phản xạ biên hộp mô phỏng: Việc sử dụng điều kiện biên PEC/PMC khiến sóng điện từ bị phản xạ trở lại vùng khảo sát khi chạm biên, đòi hỏi không gian đệm (Air buffer) phải đủ lớn để ghi nhận trọn vẹn đáp ứng xung trước khi sóng phản xạ xuất hiện.
Hướng phát triển
- Tích hợp điều kiện biên hấp thụ hoàn hảo CPML (Convolutional Perfectly Matched Layer): Hấp thụ triệt để sóng tới ở mọi góc tới mà không gây phản xạ, thu nhỏ thể tích không gian tính toán đến 50%.
- Tăng tốc phần cứng bằng GPU CUDA / OpenCL: Chuyển đổi các vòng lặp cập nhật mảng ma trận 3D sang xử lý song song hàng loạt trên hàng nghìn lõi GPU, nâng tốc độ tính toán lên gấp 20 - 50 lần.
- Hỗ trợ lưới phi cấu trúc (Conformal FDTD): Áp dụng kỹ thuật biến dạng ô mạng tại biên cong nhằm loại bỏ hoàn toàn sai số bậc thang.
Đối tượng hưởng lợi
- Sinh viên chuyên ngành Điện tử - Viễn thông & Kỹ thuật Máy tính: Nắm vững bản chất toán - lý của sóng điện từ, hiểu rõ cơ chế số hóa phương trình Maxwell thay vì chỉ sử dụng các phần mềm đóng gói sẵn.
- Kỹ sư thiết kế phần cứng RF/High-Speed PCB: Sở hữu mã nguồn mở tùy biến cao để tự động hóa quá trình tối ưu hóa hình học mạch vi dải và đánh giá nhanh tham số $S$.
- Nhóm nghiên cứu và phát triển (R&D): Nền tảng vững chắc để phát triển các bộ giải chuyên dụng cho vật liệu mới (Metamaterials, Graphene) và các phần tử quang điện tử tích hợp.
Câu hỏi thường gặp
1. Điều kiện ổn định Courant (CFL) có ý nghĩa như thế nào trong mô phỏng FDTD 3D?
Điều kiện CFL quy định vận tốc truyền sóng số trên lưới rời rạc không được vượt quá vận tốc ánh sáng thực tế. Nếu bước thời gian $\Delta t$ được chọn lớn hơn giới hạn Courant, năng lượng số sẽ tích lũy vô hạn sau mỗi vòng lặp, dẫn đến hiện tượng bùng nổ số (Numerical Instability).
2. Tại sao cần sử dụng phương pháp Newton-Raphson cho Diode thay vì cập nhật tường minh?
Do dòng điện qua Diode phụ thuộc phi tuyến theo hàm mũ vào hiệu điện thế ($I \propto e^{V}$), giá trị điện áp tại bước thời gian $n+1$ xuất hiện ở cả 2 vế của phương trình vi phân. Phương pháp cập nhật tường minh (Explicit update) không thể áp dụng trực tiếp; do đó, thuật toán lặp Newton-Raphson là bắt buộc để tìm nghiệm chính xác của phương trình $A e^{Bx} + x + C = 0$.
3. Làm thế nào để trích xuất hệ số $S_{11}$ và $S_{21}$ từ dữ liệu miền thời gian?
Điện áp và dòng điện lấy mẫu tại các cổng (Port 1 và Port 2) được thu thập trong toàn bộ chu kỳ thời gian. Sau đó, thuật toán Fast Fourier Transform (FFT) chuyển đổi tín hiệu sóng tới $V_1^+(t)$, sóng phản xạ $V_1^-(t)$ và sóng truyền qua $V_2^-(t)$ sang miền tần số $V(\omega)$. Tỷ số giữa các phổ điện áp phức này cung cấp chính xác các hệ số $S_{11}(\omega) = \frac{V_1^-(\omega)}{V_1^+(\omega)}$ và $S_{21}(\omega) = \frac{V_2^-(\omega)}{V_1^+(\omega)}$.
4. Tại sao xung Gaussian điều chế Cosine lại được ưu tiên sử dụng làm nguồn kích thích?
Xung Gaussian điều chế Cosine tập trung phần lớn năng lượng trong một dải phổ xác định quanh tần số sóng mang mong muốn, giúp tối ưu hóa tỷ số tín hiệu trên nhiễu (SNR) trong dải tần quan sát mà không làm lãng phí tài nguyên tính toán ở các dải tần không mong muốn.
5. Dung lượng bộ nhớ RAM cần thiết cho một bài toán FDTD 3D được tính như thế nào?
Với một bài toán có kích thước $N_x \times N_y \times N_z$ cells, cần lưu trữ ít nhất 6 ma trận trường chính ($E_x, E_y, E_z, H_x, H_y, H_z$) và các ma trận hệ số vật liệu. Ở định dạng số thực dấu phẩy động kép (double precision 8 bytes), dung lượng RAM ước tính: $$\text{RAM (Bytes)} \approx 6 \times (N_x \times N_y \times N_z) \times 8 \times 3 \approx 144 \times N_x N_y N_z \text{ Bytes}$$ Với lưới $100 \times 100 \times 100$, bộ nhớ yêu cầu chỉ khoảng $144\text{ MB}$, cực kỳ nhẹ và linh hoạt so với các phương pháp giải ma trận nghịch đảo.
Kết luận
Đồ án "Phân tích trường điện từ bằng phương pháp sai phân hữu hạn miền thời gian" của sinh viên Võ Phạm Minh Hiển đã hiện thực hóa thành công một bộ giải điện từ trường 3D FDTD hoàn chỉnh trên MATLAB. Không chỉ dừng lại ở việc mô phỏng sự lan truyền sóng điện từ thuần túy, công trình đã giải quyết xuất sắc bài toán tích hợp các phần tử mạch tập trung tuyến tính và phi tuyến phức tạp (Diode giải bằng Newton-Raphson), đồng thời thiết lập quy trình chuẩn trích xuất ma trận tham số tán xạ $S_{11}, S_{21}$ cho mạch lọc vi dải và bộ biến áp trở kháng cao tần.
Công trình khẳng định tiềm năng ứng dụng to lớn của phương pháp FDTD trong nghiên cứu và thiết kế vi mạch siêu cao tần, tạo tiền đề vững chắc cho việc phát triển các công cụ phần mềm mô phỏng điện từ tự chủ của Việt Nam trong tương lai.