Nhập vào chương trình tính tốn trên Fortran 77 theo chiều dài cần mơ phỏng trên EMTP 0m (nhập gần đúng là 0.01 m) , 20m, 40m, 60m, 100m. Các thông số khác được nhập như sau Độ chôn sâu 0.06 m, hằng số điện môi của đất khô (esoil = Ɛr =5), độ dẫn điện của đất ssoil = σe=1.0 mS/cm, bán kính cọc nốt đất 0.02 (m) sai số cho phép tol = 0.001 (1%).
Sau khi có dữ liệu đáp ứng tần số của trở kháng đặc trưng Zc (ω ) và hàm số lan truyền A (ω ) của điện cực nối đất tính tốn bằng chương trình viết trên Fortran 77. Ta sử dụng chương trình viết bằng Matlab để tìm hàm truyền trên miền s của cả Zc (ω ) và A (ω ) và phân tích ra theo cơng thức (3.5) để tìm , , từ đó làm cơ sở dữ liệu cho mô phỏng. Sử dụng các hàm sau trong Matlab [14] :
Hàm FRD : >> frd(response,frequency) hàm thành lập dữ liệu đáp ứng tần số có chứa
đáp ứng (response) và tần số (frequency) để nhập dữ liệu tính tốn từ chương trình viết bằng Fortran ở trên. ( Response chứa dữ liệu dạng số phức Zc (ω ) và A (ω ), frequency chứa dãy tần số từ 1-9 MHZ)
Hàm chgFregUnit : >> chgFreqUnit(sys,newfrequnits) thay đổi đơn vị tần số từ rad/s
sang Hz
Hàm TFEST : >> tfest(data,np,nz) ước tính hàm truyền ( miền s ) từ dữ liệu đáp ứng
tần số frd (data) đã thành lập trước đó với số cực (np) và số zeros (nz) là 25 do tính tốn trong EMTP cho đường dây FD tính với số cực (np) và số zeros (nz)chuẩn là 25.
Hàm POLE :>> P = pole(sys) tìm các poles ( ) trong hàm truyền
Hàm ZERO: >> P = pole(sys) tìm các zero ( ¿ trong hàm truyền
Trong đó sau khi Matlab ước tính được hàm truyền từ dữ liệu đáp ứng tần số mà ta
nhập vào ra dạng (3.25 ) ta sử dụng các phương pháp phân tích của Laplace ngược như sau để phân tích ra thành dạng tổng các phân thức đơn giản.
Dùng hàm POLE ta tìm được . Sau khi tìm được các ta tìm bằng cách chạy vịng lặp trong Matlab sử dụng cơng thức sau.
(3.26)
Sau khi tính tốn ra được , , và đưa các dữ liệu vào file pun trong EMTP. Ta sử mơ hình dường dây phụ thuộc tần số trong EMTP sau đó chọn đường dây. Chọn file dữ liệu đã thực hiện bằng cách vào mục Time – domain model data from Line Data chọn Select
Data File rồi chọn file có đi là .pun chứa dữ liệu đã tính tốn từ Matlab.