plotmodel.m

弹性波动方程正演程序

%这个程序用来绘制一个放在n行m列的数组中的数据,***原始数据顺序***是1列，2列，......m列,文件名为filename
clear;%清除所有工程中变量
a = ones(64,3);
b = a-gray/1.8;

filename='vel_p.dat';%waveform_source_m
n=75;%n行
m=150;%m列

fid=fopen(filename,'r');%打开data文件，把地址付给变量fid，原始数据顺序是1行，2行，......n行
nmdata=fread(fid,[n,m],'float32');%把数据（二进制格式）付给m行n列数组

for i=1:n;
for j=1:m;
if(i<=0);
nmdata(i,j)=0;
else
end
end
end

%axis xy;
imagesc(nmdata);
%b=gray;
%colormap(b);
set(gca,'YDir','reverse','XAxisLocation','top');

%zmax=max(max(nmdata));
%zmin=min(min(nmdata));
%caxis([zmin,zmax]);
ylabel('\fontsize{14} Z/10m');
xlabel('\fontsize{14} X/10m');
text(176,60,['\fontsize{14}','*'])
%axis equal;
axis([0,240,0,100]);