手把手教你做FEKO基础仿真

       FEKO软件是EMSS公司旗下的一款强大的三维全波电磁仿真软件。FEKO 在工业中广泛应用于天线的分析和设计，适用于电台和电视广播、无线系统、蜂窝移动通信系统、遥控无匙开锁系统、轮胎压力监测系统、卫星定位和通信、雷达、RFID 等领域。FEKO 矩量法 (MoM) 求解器广泛用于天线设计，此外，由于这款软件不仅具有模型分解功能（生成和使用等效源），还结合了多层快速多极子算法 (MLFMM) 等全波加速方法，或物理光学 (PO)、射线寻迹几何光学 (RL-GO) 或一致性绕射理论 (UTD) 等渐近方法，因此可以高效地对反射面天线、雷达天线和配有天线罩的天线进行分析。FEKO 还具备适合大型有限阵列的域格林函数法 (DGFM) 等功能，因此，还能准确、高效地对天线阵列进行分析。

       下面通过一个实例对FEKO的基本操作步骤进行介绍，能够得到目标的全方位静态RCS仿真数值，并通过MATLAB实现ffe文件的数据读取。

1. 建立FEKOCAD

打开桌面客户端，点击CADFEKO按钮。

1. 建立模型

点击创建新模型；

点击菜单栏的construct，即可建立对应模型；下面以圆柱体为例；

点击create即可创建。

1. 设置电磁波频率、发射方式与接受方式

左侧栏目为重点设置项目。

右键Frequency，一般设置为单频点或者线性连续频点。分别如下所示

下面以单频点为例进行仿真。点击OK，确认使用单频点。

右键source，设置发射源，右键选择plane wave。

上述设置为垂直极化结果。

注明：Polarisation angle设置为0，则为垂直极化；设置为1，则为水平极化。

右键Requests，设置接收源，右键选择Far fields。设置选项卡Position如下。

设置Advanced选项卡如下。

最后点击Create。

1. 对CAD模型进行网格剖分

点击菜单栏的Mesh。

继续点击左侧Create mesh。酌情选择剖分细腻程度。最后点击mesh。

得到剖分后结果如下，可以看到三角面元网格线。

1. 计算电磁仿真数据

点击Solve/Run

继续点击Feko solver，开始仿真计算。仿真之前会让保存，建立英文文件夹。

保存后，自动开始计算，结果如下。点击OK，完成仿真。

在对应文件夹已经形成相应文件。

1. 通过POSTFEKO进行数据显示

点击菜单栏的POSTFEKO。启动新的界面。

点击新界面的菜单栏中的Import。

点击后选择第一个选项，即ffe文件。在对应窗口中选择刚生成的ffe文件。

彩色线条即为数据显示结果。ffe文件打开如下，可以通过MATLAB读取。

MATLAB读取ffe文件代码如下：

如有代码问题，请加V（Lucky_YQW）交流。

% ---------------------------------------------------------------
%% 目的：从FEKO输出的.ffe文件读取目标RCS、电场数据，生成.txt文件。
% ---------------------------------------------------------------
clear all
clc
close all
fre_Num=1;%%频点个数
azi_Num=361; %%方位角个数
ele_Num=1;%%俯仰角个数
% 查找某个路径下的TXT文件
dir_txt = dir('E:\FEKO数据\*.ffe');
folder = string({dir_txt.folder}');
file = string({dir_txt.name}');
filepath = strcat(folder, '\', file);
TXT_Num = length(filepath);
%将多个数据分别读取并存储
for TAR_order = 1:TXT_Num   %Fixed Position & Velocity
fid1=fopen(filepath(TAR_order),'rt');        %% FEKO7.0计算的RCS_FFE数据
fe = linspace(1e8,12e7,fre_Num);
for i=1:fre_Numfrewind(fid1);%%将文件指针移到文件开头，以便重新读取文件数据n = 18+(i-1)*azi_Num*13;Data1=textscan(fid1,'%f','HeaderLines',n);Data1=Data1{:};RCS_VV(1)=Data1(7);RCS_VH(1)=Data1(8);RCS_total(1)=Data1(9);for j=2:azi_Numfrewind(fid1);%%将文件指针移到文件开头，以便重新读取文件数据m=18+(i-1)*azi_Num*13+(j-1)*13;
data_temp1=textscan(fid1,'%f','HeaderLines',m);data_temp1=data_temp1{:};ele(j)=data_temp1(1);azi(j)=data_temp1(2);EV_r(j)=data_temp1(3);EV_i(j)=data_temp1(4);EH_r(j)=data_temp1(5);EH_i(j)=data_temp1(6);RCS_VV(j)=data_temp1(7);RCS_VH(j)=data_temp1(8);RCS_total(j)=data_temp1(9);endDATA(:,TAR_order)=[RCS_total'];
%     DATA1=[ azi' EV_r' EV_i' EH_r' EH_i' RCS_VV' RCS_VH' RCS_total'];
end
fclose(fid1)
end