弹性波动方程正演程序

				%这个程序用来绘制一个放在n行m列的数组中的数据,***原始数据顺序***按列序,文件名为filename
				clear;%清除所有工程中变量
				filename='snapX.H@';
				%filename1='';
				n=75;%n行
				m=150;%m列
				nshot=1;%总炮数
				nsample=401;%采样点数
				nnn=5;%绘制第nnn炮中!!!!!!!!!
				mmm=150;%第mmm采样点对应的波场!!!!!!!!!!!!!
				kkk=(nnn-1)*nsample+mmm-1;
				fid=fopen(filename,'r');%打开data文件，把地址付给变量fid，原始数据顺序是1列，2列，......n列
				%fid1=fopen(filename1,'r');
				for i=1:kkk
				nmdata=fread(fid,[m,n],'float32');
				%nmdata1=fread(fid1,[m,n],'float32');%把数据（二进制格式）付给m行n列数组（原始数据顺序读取,按列顺序,所以变成m行n列数组）
				end
				nmdata=fread(fid,[n,m],'float32');%把数据（二进制格式）付给m行n列数组
				%nmdata1=fread(fid1,[n,m],'float32');
				%nmdatanew=nmdata1-nmdata;
				imagesc(nmdata);
				%imagesc(nmdatanew);
				colorbar;