弹性波动方程正演程序

				%这个程序用来绘制一个放在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				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]);