弹性波动方程正演程序

				%这个程序用来绘制一个放在n行m列的数组中的数据,***原始数据顺序***是1列，2列，......m列,文件名为filename
				
				clear;%清除所有工程中变量
				M=2250
				filename='wavelet.H@'
				fid=fopen(filename,'r');%打开data文件，把地址付给变量fid，原始数据顺序是1行，2行，......n行
				nmdata=fread(fid,'float32')
				plot(nmdata);
				colorbar;