tomography层析成像代码.docxVIP

a=[];a(144,108)=0; p=[];p(108,1)=0; %慢度 t=[];t(144,1)=0; T=[];T(12,12)=0; I=[];I=eye(108); k=[];k(144)=0; %斜率 m=1; %%%%%%%% 设置图形 %%%%%%%%%%%%%%%%%%%%%%%% figure; set(gca,ytick,[0:3:36]); axis([0 45 0 36]); title(单边激发—单边接收：); hold on;grid on;box on; %hold on 图像打开，grid on 网格打开，box on 坐标轴镶边。 %%%%%%%%%%%%%%% 确定矩阵a %%%%%%%%%%%%%%%%%% for S=1:12 for F=1:12 XY=[]; d=[]; k(m)=(F-S)*3/(9*5); x_z=0:45;y_z=k(m)*x_z+(S*3-1.5); plot(x_z,y_z,b); for i=5:5:40