GprMa中文说明指导书.doc

GprMa汉字说明指导书 GprMax中GprMax2D使用方法 试验环境： 操作系统：Windows7 软件版本： 参考资料： [1] 一、GprMax2D软件使用 直接运行 ..\\，会出现以下窗口(也能够在命令提醒符窗口输入命令运行) ： 输入文件名 注意：要输入文件全路径；*.in文件只要出现任何语法错误或者路径错误，软件都会自动关闭，不会有任何错误提醒。 出现以下画面(，*.in命令参考前面文章或[1])： 运行完成会发觉..\\、，复制到tools文件夹。 二、数据成像 tools文件夹下有五个m文件：，, , , 。 ，，其余三个以后有空再写。；(参考[1])介电常数；。 gprmax2g函数用于读取GprMax2D软件仿真探地雷达模型生成二进制几何数据。 gprmax2g函数原型： [mesh,header,media] = gprmax2g( filename ) ； media: 存放介质类型，； header: 存放模型几何参数； : 模型名称； : 模型在X轴每次偏移大小(单位:m)； : 模型在Y轴每次偏移大小(单位:m)； : 最大允许时间步长(单位:秒)； : 模型在X轴偏移次数； : 模型在Y轴偏移次数； 例子：如输入文件*.in中定义： #domain: #dx_dy: 那么：=; =; = 1/(c*sqrt(1/^2+1/^2)); (其中c=，为光速，公式参考[1]); =; =; mesh: 存放模型数据，M x N数组，其中M为Y轴方向Yee单元数目，N为X轴方向Yee单元数目；M=，N=； ： filegeo = ; [meshdata,header,media]=gprmax2g(filegeo); figure(1); [MM,NN]=size(meshdata); imagesc((1:NN)*,(1:MM)*,meshdata) axis(equal); xlabel(x(m)); ylabel(y(m)); gprmax函数用于读取GprMax2D与GprMax3D软件仿真探地雷达模型生成二进制波形数据。 gprmax函数原型： [Header, Fields] = gprmax( filename ) ； 1) Header：该变量包含以下组员： : 模型名称； : 迭代次数； : 在X轴每次偏移大小； : 在Y轴每次偏移大小； : 最大允许时间步长； : 仿真次数；等于*.in文件中#analysis:命令第一个参数； 例子：如输入文件*.in中定义： #domain: #dx_dy: #time_window: 12e-9 …… #analysis: 115 b #tx: MyLineSource 12e-9 #rx: #tx_steps: #rx_steps: #end_analysis: …… 那么，=ceil()； =；=； = 1/(c*sqrt(1/^2+1/^2)); =115； =; =; =; =; : : : =； =； =’MyLineSource’； =0；(等于#tx:命令第四个参数) =12e-9；(等于#tx:命令第五个参数) =； =； 2)Fields: 该变量包含以下组员： : 每个波形时间轴。*1； : Z轴方向磁矢量数据。*1* ； : X轴方向电矢量数据。*1* ； : Y轴方向电矢量数据。*1* ； 因为GprMax仿真是基于FDTD算法，以上三者关系为： 详细可参考FDTD算法相关文件。 ： fileout = ; [Header,Fields]=gprmax(fileout); N=1:; %移动次数 Position=*+(N-1)*(*); %天线每次所在位置 Data(:,:)=(:,1,:); %转换数组格式 figure(2); imagesc(Position,*1e9,Data);%画图 colorbar xlabel(Antena Position (m)); ylabel(t(ns)); 试验环境： 操作系统：Windows 7 软件版本：MATLA