plotseis.m

弹性波动方程正演程序

%这个程序用来绘制一个放在n行m列的数组中的数据,***原始数据顺序***是1列，2列，......m列,文件名为filename

clear;%清除所有工程中变量
filename='seismicZ.H@';%waveform_record_cal_1 waveform_record_obs
n=401;%n行
m=46;%m列
noshot=1;

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