求问gprmax3d的用法 - gprmax吧

level 5

大学叫兽W 楼主

#medium: 6.0 0.0 0.0 0.001 1.0 0.0 concrete
--------------------------------------
#domain: 0.6 1.3 0.65
#dx_dy_dz: 0.01 0.01 0.01
#time_window: 12e-9
--------------------------------------
#box: 0 0 0 0.6 1.3 0.5 concrete
#box: 0.2 0.55 0.1 0.4 0.75 0.3 free_space
--------------------------------------
#hertzian_dipole: 1.0 600e6 ricker MyDipole
#analysis: 21 second.out b
#tx: x 0.3 0.125 0.55 MyDipole 0.0 12e-9
#rx: 0.3 0.375 0.55
#tx_steps: 0.0 0.04 0.0
#rx_steps: 0.0 0.04 0.0
#end_analysis:
--------------------------------------
#geometry_file: second.geo
#messages: y
这是我的in文件，为什么出来的图形不对

2017年04月24日 11点04分 1

level 5

大学叫兽W 楼主

这是画出来的图形为什么会成这个样子

2017年04月24日 11点04分 2

做有野心的人

你用了surf函数绘图，就是这个样子啊。你想要什么样子的图 [啊]

2017年04月25日 08点04分

level 5

大学叫兽W 楼主

clear all;clc;
name = '..\windows\second.geo';
[mesh,ID,header,media] = gprmax3g(name);
figure
p1=patch(isosurface(mesh,1));
set(p1,'facecolor',[1,0,0],'edgecolor','none');
p2=patch(isosurface(mesh,2));
set(p2,'facecolor',[0,1,0],'edgecolor','none');
[Header, Fields] = gprmax('..\windows\second.out');
radargram(:,:) = Fields.ez(:,1,:);
t=1e9*Fields.t;
Y = Header.dy*Header.rx:Header.dy*Header.RxStepY:Header.dy*Header.rx+(Header.dy*Header.RxStepY*(Header.NSteps-1));
figure;
imagesc(Y,t,radargram);
title('Ez');
ylabel('Time (ns)');
xlabel('Scan length (m)');
figure;
[x,y]=meshgrid(Y,t);
surf(x,y,radargram,'EdgeColor','none')
view([1,1,0.5]);
title('Ez' );
ylabel('Time (m)');
xlabel('Scan length (m)');
这是我的matlab程序

2017年04月24日 11点04分 3

level 5

大学叫兽W 楼主