求问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 楼主
我想要这样的图 @做有野心的人
2017年04月26日 02点04分 4
这个像是mesh函数得到的图,不太懂,帮顶
2017年04月26日 07点04分
回复 做有野心的人 你也不懂,那懂的也不多了,哎
2017年04月26日 13点04分
回复 做有野心的人 :你好,请问gprmax2d可以弄成环形天线的么,还是必须是线天线
2017年05月02日 03点05分
回复 做有野心的人 :你好,那个雷克子波什么波形,自定义的波形有什么编写格式么
2017年05月02日 03点05分
1