电磁仿真CST场分布MATLAB后处理代码

% CST场监视器导出数据后处理代码 % CST数据导出(以采样点个数导出,不要以步长导出)： % (1) 选择Arrows箭头模式导出多个Z切面数据; % (2) 选择Contour热图模式,定义Fields on Plane并输入Z切面坐标导出单个Z切面数据. clc;clear; filename=('C:\Users\PC\Desktop\num1.25.txt'); Mx=11;Ny=32;%阵列单元尺寸 nx=220; ny=640; %导出数据的采样点数 a=importdata(filename);D=a.data; % X=D(1:nx,1);Y=D(1:nx:length(D)/nz,2); %以距离长度为坐标 X=(11/nx:11/nx:Mx);Y=(32/ny:32/ny:Ny); %以单元间隔为坐标 %=====================导入数据中包含多个Z轴坐标点时选用===================== % CutZ=find(D(:,3)==1); %查找所选定的Z轴坐标截面在数据D中行位置 % re_Ex=D(min(CutZ):max(CutZ),4); ExRe=reshape(re_Ex,nx,ny)'; %Z轴截面上xy平面的电场实部 % re_Ey=D(min(CutZ):max(CutZ),7); ExIm=reshape(re_Ey,nx,ny)'; %Z轴截面上xy平面的电场虚部 %=====================导入数据中只含特定个Z轴坐标时选用===================== re_Ex=D(1:length(D),4); ExRe=reshape(re_Ex,nx,ny)'; %Z轴截面上xy平面的电场实部 re_Ey=D(1:length(D),7); ExIm=reshape(re_Ey,nx,ny)'; %Z轴截面上xy平面的电场虚部 %========================================================================= clear ('a','D'); %清除占内存的巨大初始导入数据 Amplitude=sqrt(ExRe.^2+ExIm.^2); Phase=atand(ExIm./ExRe); subplot(1,3,1); surf(X,Y,ExRe);title('Real'); set(gca,'xtick',[0:1:11],'ytick',[0:4:32],'xlim',[0,11],'ylim',[0,32]);%设置坐标间隔&范围 shading interp; %去掉黑色网格 取光滑 view(0,90); %view(方位角,俯仰角) subplot(1,3,2); surf(X,Y,Amplitude);title('Amplitude'); set(gca,'xtick',[0:1:11],'ytick',[0:4:32],'xlim',[0,11],'ylim',[0,32]);%设置坐标间隔&范围 shading interp; %去掉黑色网格 取光滑 view(0,90); %view(方位角,俯仰角) subplot(1,3,3); surf(X,Y,Phase);title('Phase'); set(gca,'xtick',[0:1:11],'ytick',[0:4:32],'xlim',[0,11],'ylim',[0,32]);%设置坐标间隔&范围 shading interp; %去掉黑色网格 取光滑 view(0,90); %view(方位角,俯仰角)