Trích mô hình mạch tương đương từ dữ liệu S-parameter đo được: vector fitting

File Touchstone đo từ máy phân tích mạng là dữ liệu trực tiếp nhất về một linh kiện thụ động — nhưng nó chỉ là một bảng số rời rạc theo tần số. Bộ mô phỏng mạch miền thời gian như SPICE không chạy thẳng được trên bảng đó một cách ổn định: nó cần một mô hình liên tục, có cực và thặng dư, có thể viết thành mạch. Vì vậy bước trích mô hình mạch tương đương từ dữ liệu S-parameter đo được là cầu nối bắt buộc giữa phòng đo và phòng thiết kế.

Phương pháp được dùng rộng rãi nhất cho bước này là vector fitting. Bài này giải thích thuật toán làm gì, vì sao một mô hình khớp tốt vẫn có thể làm mô phỏng phát nổ, quy trình kiểm chứng mô hình trước khi giao cho người thiết kế, và vai trò của máy phân tích mạng trong chất lượng mô hình cuối cùng.

Vì sao phải trích mô hình mạch tương đương thay vì dùng thẳng bảng số

Một file S-parameter mô tả đáp ứng tại các tần số đã đo, không hơn. Khi đưa vào mô phỏng quá độ, trình mô phỏng phải nội suy giữa các điểm, ngoại suy ra ngoài dải đo (đặc biệt về phía DC), rồi chuyển sang miền thời gian. Mỗi bước đều có thể sinh ra hành vi phi vật lý: đáp ứng xuất hiện trước kích thích, năng lượng sinh ra từ một linh kiện thụ động, hoặc dao động không tắt.

Một mô hình hữu tỉ — tổng các phân thức cực/thặng dư — giải quyết gốc rễ vấn đề: nó xác định tại mọi tần số, có thể chuyển trực tiếp sang dạng không gian trạng thái, và từ đó tổng hợp thành mạch RLC có nguồn phụ thuộc mà bất kỳ trình mô phỏng SPICE nào cũng chạy được. Antonini (IEEE Trans. EMC, 45(3), 2003) trình bày bốn cách tổng hợp mạch tương đương như vậy từ đáp ứng miền tần số.

Trước khi fitting, dữ liệu đầu vào phải sạch. Các phép kiểm tính thụ động, tính thuận nghịch và tính nhân quả của file Touchstone — được IEEE 370-2020 đưa thành thước đo chất lượng — đã được phân tích trong bài kiểm tra chất lượng file Touchstone: tính thụ động, nhân quả và ngoại suy DC. Fitting một file hỏng chỉ cho ra một mô hình hỏng được làm trơn đẹp hơn.

Vector fitting hoạt động thế nào

Vector fitting được Gustavsen và Semlyen công bố năm 1999 (IEEE Trans. Power Delivery, 14(3), tr. 1052–1061). Mô hình có dạng tổng các số hạng c_n/(s − a_n) cộng hằng số d và số hạng s·h; các cực a_n và thặng dư c_n hoặc là số thực, hoặc đi thành cặp liên hợp phức.

Cái khó của bài toán là nó phi tuyến theo các cực. Ý tưởng của vector fitting là biến nó thành hai bước tuyến tính, cả hai đều với cực đã biết:

1. Chọn cực khởi đầu. Bài báo gốc khuyến nghị các cặp cực liên hợp phức có phần ảo phân bố tuyến tính trên dải tần quan tâm, phần thực bằng một phần trăm phần ảo. Cực khởi đầu thực có thể làm bài toán tuyến tính bị suy biến.
2. Định vị lại cực. Nhân hàm cần khớp với một hàm trọng số σ(s) chưa biết, biểu diễn σ(s) cũng bằng các cực khởi đầu, rồi giải bình phương tối thiểu. Các không điểm của σ(s) chính là tập cực mới, tốt hơn. Cực nào rơi sang nửa phải mặt phẳng phức thì đảo dấu phần thực.
3. Lặp lại bước định vị cho tới khi cực hội tụ, rồi giải một bài toán tuyến tính cuối cùng để lấy thặng dư.

Hai cải tiến quan trọng sau đó: relaxed vector fitting (Gustavsen, 2006) định nghĩa lại hàm trọng số để cải thiện hội tụ khi dữ liệu có nhiễu — đúng tình huống của dữ liệu đo; và fast vector fitting (Deschrijver và cộng sự, IEEE MWCL, 2008) giảm mạnh chi phí tính cho hệ nhiều cổng. Với ma trận S nhiều cổng, mọi phần tử dùng chung một tập cực và chỉ khác nhau ở thặng dư.

Ổn định, nhân quả và thụ động: ba điều kiện để mô hình chạy được

Một mô hình có sai số khớp nhỏ chưa chắc dùng được. Có ba điều kiện vật lý phải kiểm riêng:

• Ổn định và nhân quả. Khi mọi cực có phần thực không dương, mô hình vector fitting đồng thời ổn định và nhân quả; điều kiện này gần như bắt buộc cho mô phỏng miền thời gian (Triverio, chương “Vector Fitting”, arXiv 1908.08977). Bước đảo dấu cực không ổn định trong thuật toán gốc bảo đảm điều này.
• Thụ động. Với tham số tán xạ, mô hình thụ động khi giá trị suy biến lớn nhất của ma trận S không vượt quá 1 tại mọi tần số thực. Đây là điểm hay bị bỏ qua: áp dụng vector fitting lên đáp ứng của một hệ thụ động vẫn có thể cho ra mô hình không thụ động do sai số xấp xỉ và sai số số học.
• Ổn định khi ghép vào mạch. Gustavsen và Semlyen (IEEE Trans. Power Systems, 16(1), 2001) chỉ ra các ví dụ mô hình chỉ có cực ổn định nhưng vi phạm tiêu chí thụ động vẫn làm mô phỏng mất ổn định khi ghép với phần mạch còn lại. Cực ổn định là điều kiện cần, không đủ.

Vì vậy sau fitting luôn có bước kiểm và cưỡng bức tính thụ động. Hai họ phương pháp phổ biến là kiểm tra trị riêng thuần ảo của ma trận Hamilton để xác định dải tần vi phạm (Grivet-Talocia, IEEE TCAS-I, 2004), và nhiễu loạn thặng dư để kéo mô hình về vùng thụ động với thay đổi đáp ứng nhỏ nhất. Một số triển khai chỉ nhiễu loạn thặng dư, không đụng hằng số, để giữ nguyên điểm DC.

Quy trình kiểm chứng mô hình trước khi giao cho người thiết kế

Kiểm chứng mô hình là bước phân biệt một mô hình “khớp đẹp trên đồ thị” với một mô hình dùng được. Quy trình gợi ý:

1. Kiểm dữ liệu vào. Chạy các phép kiểm thụ động, thuận nghịch, nhân quả trên file đo. Kiểm dải tần đủ thấp: máy phân tích mạng thường không đo tới DC, nên điểm DC phải được bổ sung sau; tần số bắt đầu càng thấp thì ngoại suy DC càng đáng tin.
2. Chọn bậc mô hình. Đáp ứng trơn kiểu thông thấp thường chỉ cần 1–3 cực thực; dữ liệu có cộng hưởng cần số cặp cực phức tương ứng. Quá ít cực thì khớp kém; quá nhiều cực thì mô hình học cả nhiễu đo và có thể sinh cộng hưởng giả ở ngoài dải khớp.
3. Đánh giá sai số khớp bằng sai số trung bình bình phương (RMS) trên toàn dải và trên từng phần tử của ma trận S, không chỉ nhìn đường S21.
4. Kiểm và cưỡng bức thụ động, rồi tính lại sai số sau cưỡng bức — bước cưỡng bức làm thay đổi đáp ứng, và thay đổi đó phải nhỏ.
5. Kiểm miền thời gian. Mô phỏng đáp ứng xung và đáp ứng bậc của mô hình, so với biến đổi ngược của dữ liệu đo. Tìm dấu hiệu đáp ứng trước kích thích hoặc dao động không tắt.
6. Kiểm mạch tổng hợp trong chính trình mô phỏng đích. Mạch tổng hợp có thể chứa điện trở âm — bình thường về mặt toán học, nhưng là lý do phải chạy thử trong đúng công cụ sẽ dùng.

Không có quy tắc phổ quát về mật độ điểm tần số tối thiểu. Điều chắc chắn là quá ít điểm sẽ làm hệ phương trình thiếu ràng buộc và fitting thất bại, còn các cộng hưởng hẹp phải có đủ điểm đo phủ lên đỉnh. Hãy chọn IF bandwidth và số điểm quét trên máy phân tích mạng với mục đích fitting ngay từ đầu, thay vì dùng lại một file đo cho mục đích khác.

Dữ liệu S-parameter tốt bắt đầu từ máy phân tích mạng

Chất lượng mô hình bị chặn trên bởi chất lượng dữ liệu S-parameter. Ba yếu tố phía phần cứng ảnh hưởng trực tiếp: dải tần xuống thấp (quyết định ngoại suy DC), dải tần lên cao (quyết định mô hình có đúng ở sườn tín hiệu nhanh không), và dải động (quyết định nhiễu trong các phần tử có suy hao lớn như S21 của đường truyền dài hay S31 xuyên âm).

Keysight N5245B PNA-X có dải từ 900 Hz hoặc 10 MHz tới 50 GHz, 2 hoặc 4 cổng — phù hợp đặc tính hoá đường mạch và đầu nối tốc độ cao, nơi cần cả dải thấp lẫn dải cao trong cùng một phép đo nhiều cổng.

Keysight E5080B ENA dải 9 kHz tới 20 GHz hoặc 100 kHz tới 53 GHz, 2 hoặc 4 cổng, dải động 140 dB — lựa chọn cân bằng cho linh kiện thụ động và bộ lọc, với điểm bắt đầu 9 kHz rất có lợi cho ngoại suy DC. Cho trạm đo gọn hoặc dây chuyền, Keysight P5024A thuộc dòng Streamline USB, dải 9 kHz tới 20 GHz, điều khiển bằng máy tính ngoài, dải động tối thiểu 140 dB từ 50 MHz tới 3 GHz.

Ở phía phần mềm, PathWave ADS của Keysight có công cụ Broadband SPICE Model Generator: chuyển dữ liệu tham số mạng đo hoặc mô phỏng của linh kiện thụ động thành mô hình mạch tương đương dạng phân thức cực/thặng dư, có tuỳ chọn thuật toán relaxation, cưỡng bức thụ động và xuất netlist cho HSPICE, Spectre, SPICE3. Ghi chú phát hành ADS 2026 Update 1.0 nêu công cụ này được cải thiện độ chính xác khớp ở DC và tần số thấp. Techmaster cung cấp Keysight W4321EP, một gói phần mềm thiết kế PathWave; thành phần gói được cấu hình theo quy trình mô phỏng, nên hãy nêu rõ nhu cầu trích mô hình khi yêu cầu báo giá.

Nếu mục tiêu là mô hình hành vi cho bộ đệm vào/ra số thay vì mô hình mạng thụ động, cách tiếp cận khác hẳn — xem bài trích mô hình IBIS từ dữ liệu đo.

Kết luận

Trích mô hình mạch tương đương bằng vector fitting là quy trình ba tầng: dữ liệu đo sạch, thuật toán khớp hội tụ, và kiểm chứng vật lý — ổn định, nhân quả, thụ động — trước khi mô hình rời phòng đo. Tầng thứ ba hay bị bỏ qua nhất, vì một mô hình khớp đẹp trên đồ thị vẫn có thể làm mô phỏng miền thời gian mất ổn định. Và toàn bộ chuỗi phụ thuộc vào tầng thứ nhất: sai số của máy phân tích mạng sẽ được fitting trung thành biến thành một phần của mô hình.

Techmaster Electronics là đại lý uỷ quyền của Keysight tại Việt Nam, vận hành phòng hiệu chuẩn theo ISO/IEC 17025 cho máy phân tích mạng và bộ chuẩn hiệu chuẩn cơ khí/điện tử — nền tảng để dữ liệu S-parameter có truy xuất nguồn gốc trước khi đưa vào trích mô hình. Để được tư vấn cấu hình hoặc đặt lịch hiệu chuẩn, vui lòng liên hệ hotline 0908 173 345.

Biên soạn bởi Khanh Nguyen — Marketing Manager, Techmaster Electronics JSC.