tomography层析成像代码

1、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:12for F=1:12XY=;d=;k(m)=(F-S)*3/(9*5);x_z=0:45;y_z=k(

2、m)*x_z+(S*3-1.5); plot(x_z,y_z,b);for i=5:5:40XY(1,i/5)=i;XY(2,i/5)=k(m)*i+(1.5+3*(S -1); end mm=9;for j=(min(S,F)*3):3:(max(S,F)*3 -3)XY(1,mm)=(j-(1.5+3*(S-1)/k(m);XY(2,mm)=j;mm=mm+1;endXY(1,mm)=45;XY(2,mm)=(F-1)*3+1.5;for i=1:mm -1for j=i+1:mmif XY(1,i)=XY(1,j)&XY(2,i)=XY(2,j) XY(:,j)=;mm=mm -1;en

3、dendendfor i=1:mm d(i)=(XY(1,i)A2+(XY(2,i(1.5+(S-1)*3)A2)A(0.5);enddd,nn=sort(d);for i=1:mmif i=1a(m,(S-1)*9+1)=dd(1);elseif k(m)=0if mod(XY(1,nn(i),1)=0a(m,(9*floor(XY(2,nn(i)/3)+floor(XY(1,nn(i)/5)=dd(i) -dd(i-1);elseif mod(XY(2,nn(i),1)=0a(m,(9*(floor(XY(2,nn(i)/3) -1)+floor(XY(1,nn(i)/5)+1)=dd(i

4、) -dd(i-1);endelseif k(m)0if mod(XY(1,nn(i),1)=0a(m,(9*floor(XY(2,nn(i)/3)+floor(XY(1,nn(i)/5)=dd(i) -dd(i-1);elseif mod(XY(2,nn(i),1)=0a(m,(9*(floor(XY(2,nn(i)/3)+floor(XY(1,nn(i)/5)+1)=dd(i) -dd(i-1);endendendm=m+1;endendfor i=1:108p(i,1)=1/4;endp(20,1)=1/3;p(21,1)=1/3;p(70,1)=1/5;p(71,1)=1/5;P1=reshape(p,9,12);P2=rot90(P1);figure;mag=imagesc(P2);t=a*p;for i=1:12for j=1:12T(i,j)=t(i-1)*12+j);endenddisp(单边激发单边接收的走势矩阵：)；disp(矩阵的第i行第j列所表示第i个激发点第j个接收点);y=a;x=0.5*I;S=a*a+x;s1=inv(S